0

特に、ブートストラップオプションと、それが の扱いにくい出力ではなくリストを返すという事実が好きなので、 lawstatlevene.testパッケージの関数を使用して等分散性をテストしたいと思います。私のデータセットからNAをデフォルトで除外する必要があることは明らかです。以下、元のデータを提供します。car::leveneTestlawstat::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 を実行しています。

4

3 に答える 3

0

levene.test() をいじってみたところ、問題は欠損値にあると思います。

test <- cbind(testset.logcount, testset.treat)
test <- test[complete.cases(test),] #removing Nas
levene.test(test[,1], test[,2])


        modified robust Brown-Forsythe Levene-type test based on the absolute
        deviations from the median

data:  test[, 1]
Test Statistic = 0.9072, p-value = 0.3487

これは車のレベン テスト (df = 29) と一致するため、欠落している行を自動的に削除する必要があります。

> leveneTest(y=testset.logcount,group=testset.treat)
Levene's Test for Homogeneity of Variance (center = median)
      Df F value Pr(>F)
group  1  0.9072 0.3487
      29     
于 2013-08-12T15:16:10.073 に答える
0

私も同様のプロセスを経ました。目的は、'y' が数値ベクトルで、'g' がデータの係数である場合に、2 番目の要素が消去されないようにすることです。

yg = na.omit(data.frame(y,g))
y = yg[,1]
g = yg[,2]
于 2014-11-13T09:09:58.670 に答える