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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我想问下，我用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)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：方差齐 Bartlett Bartle 基础数据 seed 百度

有人吗？求助！在线等！

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##

机械唯物主义123 发表于 2016-12-24 01:07
x

终于有人回复了！谢谢你，可是我运行你的代码运行不出来。图在下面

运行不出来

已经运行出来了啊？

少打了几个#号，不要把#号间的代码复制进去