R中RCB设计的方差齐性检验

R中RCB设计的方差齐性检验,r,statistics,R,Statistics,我有下面的问题,我想对一个完全随机的块设计应用方差齐性检验。 我试过左旋,但有两个因素不能接受 库(车) 图书馆(tibble) tb=tibble(“治疗”=因子(x=rep(1:3,c(4,4,4)), 标签=c(“第1部分”、“第2部分”、“第3部分”), “bloque”=系数(x=c(代表(1:4,3)), 标签=c(“东北”、“西北”、“东南”、“西南”), “值”=c(250350219375, 400,525,390,580, 275,340,200,310)) 在这里,我应用

我有下面的问题,我想对一个完全随机的块设计应用方差齐性检验。 我试过左旋,但有两个因素不能接受

库(车)
图书馆(tibble)
tb=tibble(“治疗”=因子(x=rep(1:3,c(4,4,4)),
标签=c(“第1部分”、“第2部分”、“第3部分”),
“bloque”=系数(x=c(代表(1:4,3)),
标签=c(“东北”、“西北”、“东南”、“西南”),
“值”=c(250350219375,
400,525,390,580,
275,340,200,310))
在这里,我应用了levene测试

levenest(值~treat*bloque,tb)#无显示结果

对于RCBD的方差同质性是否还有其他统计检验?

每个组组合只有1个观察值:

table(tb$treat,tb$bloque)
         
          NE NW SE SW
  disen 1  1  1  1  1
  disen 2  1  1  1  1
  disen 3  1  1  1  1
不能使用n=1计算方差(即,最佳拟合是单个值)

如果至少有一个组的数据n>1,则它可以工作:

mtcars$cyl = factor(mtcars$cyl)
mtcars$gear = factor(mtcars$gear)

with(mtcars,table(gear,cyl))
    cyl
gear  4  6  8
   3  1  2 12
   4  8  4  0
   5  2  1  2

leveneTest(mpg ~ gear*cyl,mtcars)

Levene's Test for Homogeneity of Variance (center = median)
      Df F value  Pr(>F)  
group  7  2.1926 0.07172 .
      24                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

每个组组合只有1个观察:

table(tb$treat,tb$bloque)
         
          NE NW SE SW
  disen 1  1  1  1  1
  disen 2  1  1  1  1
  disen 3  1  1  1  1
不能使用n=1计算方差(即,最佳拟合是单个值)

如果至少有一个组的数据n>1,则它可以工作:

mtcars$cyl = factor(mtcars$cyl)
mtcars$gear = factor(mtcars$gear)

with(mtcars,table(gear,cyl))
    cyl
gear  4  6  8
   3  1  2 12
   4  8  4  0
   5  2  1  2

leveneTest(mpg ~ gear*cyl,mtcars)

Levene's Test for Homogeneity of Variance (center = median)
      Df F value  Pr(>F)  
group  7  2.1926 0.07172 .
      24                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

完美,坦克你完美,坦克你