特に、ブートストラップオプションと、それが の扱いにくい出力ではなくリストを返すという事実が好きなので、 lawstatlevene.test
パッケージの関数を使用して等分散性をテストしたいと思います。私のデータセットからNAをデフォルトで除外する必要があることは明らかです。以下、元のデータを提供します。car::leveneTest
lawstat::levene.test
testset.logcount<-c(6.86923171973098, 6.83122969386706, 7.30102999566398,7.54282542695918,6.40823996531185, 6.52891670027766, 6.61278385671974, 6.71933128698373, 6.96567197122011, 6.34242268082221, 6.60205999132796, 6.69897000433602, 6.6232492903979, 6.54157924394658, 6.43136376415899, 6.91381385238372,6.44715803134222, 6.30102999566398, 6.10037054511756, 6.7481880270062,NA, 4.89762709129044,5.26951294421792, 5.12385164096709, 5.11394335230684, 4.43136376415899, 5.73957234445009, 5.83250891270624, 5.3451776165427, 5.77887447200274, 5.38524868240322, 5.75127910398334, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA)
testset.treat<-structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("CTL","TRM"), class = "factor")
実行するlawstat::levene.test(y=testset.logcount,group=testset.treat)
と、次のエラー メッセージが表示されます: Error in contrasts<-
( *tmp*
, value = contr.funs[1 + isOF[nn]]) : 対比は、2 つ以上の水準を持つ因子にのみ適用できます
私によると、testset.treat には明らかに 2 つのレベルがあります。
また、leveneTest(y=testset.logcount,group=testset.treat)
またはfligner.test(x=testset.logcount,g=testset.treat)
両方を使用すると、エラーなしで実行されます。
lawstat::levene.testでこの特定のエラーが発生した理由を見つけることができませんでした。誰かが私を助けてくれることを願っています。
x86_64-w64-mingw32/x64 プラットフォーム (Windows 7、64 ビット) で R 3.0.0 を実行しています。