同一数据集中的多个Brown-Forsythe检验和paircomp

Question

我的实验有36个比较，3个隔室，每个隔室21个数据点。数据集具有2268个Obs。总共。第一行看起来像这样：

Comparison  Value   Compartment 
A   -0.715126087    1  
A   -0.481391603    2 
A   0.374693449 1  
A   null    2  
A   0.857450232 1  
A   null    2  
A   -1.608637992    3
A   -1.670859336    3  
A   -1.731618976    3
...    
AJ  -5.84850106 3

我运行ANOVA 每次比较并测试了假设，但在某些比较中有违反情况。我想运行非参数的Brown-Forsythe F *测试，然后进行成对比较。我的尝试：

if(!require(onewaytests)){install.packages("onewaytests")} # for Brown-Forsythe Test
if(!require(tidyverse)){install.packages("tidyverse")}
if(!require(broom)){install.packages("broom")}

fm3 <- my_data %>% 
  group_by(Comparison) %>% 
  do(multitst = paircomp(bf.test(Value ~ Compartment, data = .)))
fm3 %>% tidy(multitst)

问题：未创建fm3，因为它在遇到第一个非统计显着性比较后便立即停止运行并出现错误。这是我得到的错误：

paircomp.owt（bf.test（Value〜Compartment， ：由于存在差异，因此无法进行成对比较 没有统计学意义（alpha = 0.05）。

希望：在这样的表中有一个输出：

Comparison  term  comparison  estimate  conf.low  conf.high  adj.p.value
1   A   Compartment 1-2 3.360531    2.9495551   3.7715075   2.023126e-11
2   A   Compartment 3-2 3.537098    3.1293925   3.9448041   2.023126e-11
3   A   Compartment 3-1 0.176567    -0.2382312  0.5913652   5.779274e-01
...
108 AJ  Compartment 3-1 4.50731363  2.85507469  6.1595525749    1.954819e-07

有人可以帮忙吗？预先感谢！

Answer 1

您可以包含ifelse（）语句以在比较不重要时填写NA，在比较不重要时返回结果。 我不熟悉Brown-Forsythe测试，但是在这种情况下，我最喜欢的是随机测试。我猜这个决定应该基于使用这种特殊的非参数测试会丢失多少信息。