2011-07-27 15 views
0

5つのサイトの密度データの正常性の多変量テストを実行しようとしています。 mvnormtestパッケージのmshapiro.test。各サイトは列で、密度は以下のとおりです。一番上の行をヘッダー(サイト名)として、5列目と5行目です。ここで私は自分のデータをロードする方法である:mshapiro.testを使用しているときにエラーが発生しました。「U []は列数(サンプルサイズ)が3〜5000の行列ではありません」

datafilename="/Users/megsiesiple/Documents/Lisa/lisadensities.csv" 
data.nc5=read.csv(datafilename,header=T) 
attach(data.nc5)` 

データは次のようになります。私はmshapiro.test(呼び出すとき

 B07  B08  B09  B10  M 
1 72571.43 17714.29 3142.86 22571.43 8000.00 
2 44571.43 46857.14 49142.86 16857.14 7142.86 
3 54571.43 44000.00 26571.43 6571.43 17714.29 
4 57714.29 38857.14 32571.43 2000.00 5428.57 

)data.nc5のために、私はこのメッセージが表示されます。Error in mshapiro.test(data.nc5) : U[] is not a matrix with number of columns (sample size) between 3 and 5000

mshapiro.test()を使用してShapiro-Wilkテストを実行するには、データが3から5000の間の数の数値行列でなければならないことがわかります。ただし、.csvを数値のみの行列(つまり、私はサイト名を省略すると)、私はまだエラーが発生します。行列を別に設定する必要はありますか?誰にもこの問題がありましたか? ありがとう!

+2

あなたはそれが不要だやっていることに応じて、私の時間帯に後半っぽいではありません

mat <- data.matrix(data.nc5) mshapiro.test(mat) 

をお試しください'attach(data.nc5)'( 'attach'を使うと混乱するのは簡単です) –

答えて

2

read.csv()でRに数値行列を表示する場合、表示するコードと同様のコードを使用すると、データフレームとして読み込まれ、行列ではありません。

(あなたが再現可能な例を与えないように試験し、それが今;-)

+0

Yeehaw!出来た!次回はその例を示します。画像やその他のものとして貼り付けるかどうかはわかりませんでした。ありがとう、ギャビン。 – Margaret

+0

@マーガレット:dput(data.nc5)から出力された5 x 5オブジェクトの場合、Gavinにチェックマークを付けることを忘れないでください。彼はそれが必要ではありませんが、それは良いマナーです。 –

+1

@マーガレットあなたが答えに満足しているなら、それを受け入れてください(回答の横にある大きな目盛りを確認してください)。これは、提供されたAnswersの1つが問題を解決したのと同じまたは同様の問題を抱えている可能性があることを他人に伝えます。これは、StackOverflowサイトの価値をAnswers/Advice –

関連する問題