我最近在学习使用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'
请各位高手帮帮忙，看看错在哪里？谢谢！！