【求助】方差齐性的检验出现报错；每组至少需要有两个观察量

7952

收藏 2016-12-23

我想问下，我用bartlett.test()来做方差齐性的检验，我作出来的A也是数据框类型了，可是为什么就是报错，用老师的又可以
，然后我对比了一下，差别也就只是我的是用随机产生的数做成的数据框，老师的是自带的基础数据框，可是我的题目就是要求用随机产生的数据来做。也查了下百度，还是解决不了。因此来向各位求助！麻烦大家了！

我的代码如下：
set.seed(1)
x<-rnorm(25,mean=100,sd=5)
y<-rnorm(25,mean=100,sd=5)
a<-c(x,y)
A <- data.frame(x,y)
View(A)
bartlett.test(x~y,data=A)#检验方差齐性#

老师的代码如下：
attach(cholesterol)
library(car)
qqPlot(lm(response~trt), main="Q~Q plot",labels=F) #检验正态性#
bartlett.test(response~trt,data=cholesterol) #检验方差齐性#
detach(cholesterol)

附件列表

捕获2.PNG

原图尺寸 8.88 KB

扫码加我拉你入群

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

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

全部回复

Licko123456

2016-12-23 20:43:06

有人吗？求助！在线等！

扫码加我拉你入群

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

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

机械唯物主义123

2016-12-24 01:07:15

x<-rnorm(25,mean=100,sd=5)
y<-rnorm(25,mean=100,sd=5)
a<-c(x,y)
##方法1
library(lmtest)
bptest(y~x)
# studentized Breusch-Pagan test

data:  y ~ x
BP = 3.4403, df = 1, p-value = 0.06362

p>0.05,其次性较好#
##方法2
factorx=rep(c("x","y"),25)
bartlett.test(a~factorx)##

      Bartlett test of homogeneity of variances

data:  a by factorx
Bartlett's K-squared = 0.57304, df = 1, p-value = 0.4491##