[问答] 求助－－使用R 做方差分析并做两两比较 [推广有奖]

0关注
0粉丝

学前班

60%

还不是VIP/贵宾

威望: 0 级
论坛币: 0 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 33 点
帖子: 2
精华: 0
在线时间: 1 小时
注册时间: 2010-6-11
最后登录: 2017-12-13

楼主

zhui126 发表于 2011-4-23 22:35:15 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我最近在学习使用R编程，做一个方差分析程序，程序如下：
library(foreign, pos=4)
Dataset <- read.spss("D:/stat.data/german2010512.sav",
  use.value.labels=TRUE, max.value.labels=Inf, to.data.frame=TRUE)
library(multcomp, pos=4)
library(abind, pos=4)
AnovaModel.1 <- aov(days.mgraine.5 ~ group, data=Dataset)
summary(AnovaModel.1)
numSummary(Dataset$days.mgraine.5 , groups=Dataset$group,
  statistics=c("mean", "sd"))
pairwise.t.test(days.mgraine.5, group, p.adj="bonferroni", paired=T)
具体内容是从一个SPSS数据当中读取数据，然后按病人的头痛天数（days.mgraine.5）和组别（group）进行分析，程序运行结果如下：
> library(foreign, pos=4)
> Dataset <- read.spss("D:/stat.data/german2010512.sav",
+ use.value.labels=TRUE, max.value.labels=Inf, to.data.frame=TRUE)
read.spss("D:/stat.data/german2010512.sav", use.value.labels = TRUE, 里有警告:
  D:/stat.data/german2010512.sav: Unrecognized record type 7, subtype 18 encountered in system file
> library(multcomp, pos=4)
> library(abind, pos=4)
> AnovaModel.1 <- aov(days.mgraine.5 ~ group, data=Dataset)
> summary(AnovaModel.1)
         Df Sum Sq Mean Sq F value  Pr(>F)
group       3 76.7 25.5776  2.6472 0.04858 *
Residuals 435 4203.0  9.6621
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
37 observations deleted due to missingness
> numSummary(Dataset$days.mgraine.5 , groups=Dataset$group,
+ statistics=c("mean", "sd"))
         mean    sd n NA
group A 2.314815 3.911756 108 13
group B 2.027273 1.988310 110  9
group C 2.495495 3.011716 111  7
group D 3.163636 3.229813 110  8
> pairwise.t.test(days.mgraine.5, group, p.adj="bonferroni", paired=T)
运行到最后一句，也就是pairwise.t.test(days.mgraine.5, group, p.adj="bonferroni", paired=T)时，出现错误，提示误于factor(g) : 找不到对象'group'
请各位高手帮帮忙，看看错在哪里？谢谢！！