RパッケージadehabitatHSと関数を使用して悪いテスト結果が得られますcompana()
。
コンパナテストでは、生息地の利用可能性と動物による生息地の使用状況を比較します。このテストは組成分析を実行し、ウィルクのラムダとp値を生成します。
コンパナテストの規定された形式と内容の要件を満たす経験的データ(16匹の動物、8つの生息地カテゴリ)でテストを実行していますが、テストは次のように返されますlambda = NaN, p-value = 1.0
。
NA
NaNの標準的な原因についてデータを調べましたが、値の欠落( s)、負の比率、行の合計が100%を超えるなどの問題は見られません。さらに、完全なデータセット(16行8列)を使用するとエラーが発生しますが、サブセット(1列が欠落)を使用するとエラーは発生しません。7つ以下しか含まれていない限り、どの列が除外されているかは問題ではないようです。
私のデータの詳細は以下の通りです。どんな助けや提案も大歓迎です。
プラットフォーム:OSX 10.7.4(Lion)Rバージョン:R 2.15.1 GUI 1.52 Leopardビルド64ビット(6188)
Rコード:
使用された生息地の割合(マトリックス):
used.m
Atr Cre DW Dist Gra Mar Mes Opn
1 2.76 19.31 4.14 17.24 0.00 28.28 0.00 0.00
2 15.12 1.16 3.49 5.81 2.33 32.56 25.58 3.49
3 57.35 13.24 1.47 0.00 6.62 0.00 8.09 12.50
4 41.18 0.98 2.94 0.00 3.92 0.00 29.41 18.63
5 9.52 17.01 23.13 0.00 0.00 37.41 3.40 0.00
6 2.54 67.01 0.00 0.51 0.51 6.09 3.05 0.00
7 10.45 15.67 5.22 5.97 0.00 0.00 33.58 10.45
8 6.34 0.00 0.00 4.23 0.00 0.00 0.00 21.83
9 0.66 18.54 1.32 39.74 0.00 8.61 2.65 0.66
10 20.75 0.00 0.00 0.00 0.00 0.00 0.00 49.06
11 20.14 18.06 12.50 8.33 0.00 5.56 22.22 9.72
12 0.00 0.00 0.00 58.04 0.00 0.00 13.99 0.00
13 44.23 4.81 0.00 0.00 0.00 0.00 28.85 19.23
14 13.04 2.61 1.74 2.61 0.00 19.13 6.09 25.22
15 35.96 3.51 0.00 2.63 0.00 1.75 25.44 0.88
16 0.00 0.00 0.00 35.14 0.00 7.21 0.00 8.11
利用可能な生息地の割合(マトリックス):
avail.m
Atr Cre DW Dist Gra Mar Mes Opn
1 1.11 39.18 15.50 0.29 0.00 30.18 4.15 3.43
2 18.96 1.02 0.87 7.40 4.36 37.80 6.17 10.56
3 65.75 3.85 0.57 0.00 12.36 0.00 9.97 3.43
4 20.52 5.44 2.01 0.00 0.08 0.00 10.99 44.22
5 11.16 43.55 10.34 0.00 0.00 18.11 9.60 3.17
6 25.24 9.88 2.77 2.44 4.95 12.90 11.44 2.19
7 5.47 21.13 6.67 0.57 0.00 0.10 30.73 9.59
8 6.23 0.00 0.00 4.85 0.00 0.00 0.00 5.72
9 0.85 1.97 0.09 33.52 0.00 5.58 1.73 0.84
10 18.97 0.00 0.00 0.00 0.00 0.00 0.00 5.59
11 0.16 44.41 19.16 3.53 0.00 1.28 21.65 4.13
12 6.67 3.89 0.31 13.63 0.00 0.00 20.56 0.00
13 21.40 2.38 2.66 0.00 0.00 0.00 21.23 46.96
14 42.62 9.36 1.49 3.89 0.43 3.61 20.79 5.80
15 33.32 12.90 4.17 4.09 0.15 3.96 12.32 4.94
16 5.65 0.00 0.00 1.36 0.00 0.59 0.00 9.71
adehabitatHSパッケージをロードします
library(adehabitatHS)
完全なデータセット(すべての行、すべての8列)でcompanaテストを実行します
fulldata.comp <- compana(used.m, avail.m)
チェック結果
fulldata.comp$test
Lambda P
NaN 1
部分データ(すべての行、列1〜7)でcompanaを実行します
partdata.comp <- compana(used.m[, 1:7], avail.m[, 1:7])
チェック結果
partdata.comp$test
Lambda P
0.3319731 0.6280000