私は、私のデータフレームでshapiro-wilkes正規性仮説検定のp値のテーブルを取得することに取り組んできました。これは、カンマ区切りの CSVとしてのデータ フレーム (「mdf1」という名前)です。
R での Shapiro-Wilkes テストには、3 を超えるサンプル サイズが必要です。データ フレーム (「変数」と「サイト」という 2 つの関連因子を含む) をサブセット化するために、次のコードを使用しました。
Z <- as.data.frame(data.table(mdf1)[, list(freq=.N, value=value), by=list(Site,variable)][freq > 3])
これにより、3 より大きい n の「サイト」*「変数」の組み合わせに属するすべての値を含むデータ フレーム「Z」が生成されました。次に、Z をddply
関数に渡して、shapiro-wilkes p のテーブルを取得しようとしました。 -値:
norm2 <- ddply(Z, .(Site, variable), summarize, n=length(value), sw=shapiro.test(value)[2])
このコマンドの結果は次のとおりです。
Error in shapiro.test(val) : all 'x' values are identical
それはどうしてですか?何かご意見は?