全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
23429 12
2014-08-16
R自身包含var.test和bartlett.test方差齐性检验,分别对应两组和多组数据。在car包中有levene.test(也是SPSS的默认方差齐性检验方法)。对同一组数据,var.test和bartlett.test检验p值接近,而与levene.test相差甚远。数据是用的Robert A. Muenchen的“R for SAS and SPSS users(2nd edit.)”的mydata100
mydata100.xls
大小:(35 KB)

 马上下载


这个问题我直接发邮件问Robert A. Muenchen,他回复说暂时也想不起来问题出在哪里,但如果有时间会帮忙查查资料找找原因。

请教论坛的朋友们,有知道问题出在哪里的吗?谢谢!


结果如下所示:
> var.test(formula = posttest~gender,data = mydata100)

        F test to compare two variances

data:  posttest by gender
F = 0.7279, num df = 46, denom df = 51, p-value = 0.2763
alternative hypothesis: true ratio of variances is not equal to 1
95 percent confidence interval:
0.4132401 1.2929517
sample estimates:
ratio of variances
         0.7278693

> bartlett.test(formula = posttest~gender,data = mydata100)

        Bartlett test of homogeneity of variances

data:  posttest by gender
Bartlett's K-squared = 1.196, df = 1, p-value = 0.2741

> library(car)
> leveneTest(y = mydata100$posttest,group = mydata100$gender)
Levene's Test for Homogeneity of Variance (center = median)
      Df F value Pr(>F)
group  1  0.4308 0.5131
      97


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2014-8-16 17:04:59
不奇怪啊,算法的原理不一样啊,前两者是对原始数据的方差进行检验的,leveneTest是对方差模型的残差进行组间齐性检验的,当然结果不一样啦,呵呵。一般认为是要求残差的方差齐,所以一般的统计软件都做的是leveneTest。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-8-17 12:22:06
yangming98 发表于 2014-8-16 17:04
不奇怪啊,算法的原理不一样啊,前两者是对原始数据的方差进行检验的,leveneTest是对方差模型的残差进行组 ...
涨知识了,非常感谢!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-9-5 10:56:53
你好,我安装你的脚本做levenetest,总是出现   Error: unexpected symbol in "levenetest(y=OTUANOVAOUT richness"    这是怎么回事呢,求解答
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-9-5 11:05:17
火飞猫 发表于 2017-9-5 10:56
你好,我安装你的脚本做levenetest,总是出现   Error: unexpected symbol in "levenetest(y=OTUANOVAOUT r ...
我的数据是OUT  richness是因变量,Treatment是自变量,储存的文件名是OTUANOVA
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2018-3-4 11:45:55
请教各位:奇怪的是我用r里的levenetest做方差齐性分析得到的结果和SPSS分析的结果都不一样!我的两组数据x1<-c(1,1,1); X2<-C(1,0.9,0.9)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群