標本共分散の計算<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />次の行列があります。[35 36 40、38 37 40、48 41 42、50 46 54]ここで、サンプル共分散を計算しようとしています。 手作業で計算する場合、サンプル共分散行列は次のようになります。[40.6875 23.5 28.0、23.5 15.5 21.5、28.0 21.5 34.0]ただし、これをJMPで再現することはできません。 共分散を計算するためのオプションが見つかりましたが、同じ結果が返されません。
誰かが私を正しい方向に向けるのを手伝ってくれる?私はこのソフトウェアに非常に慣れていないし、大学生としてソフトウェアを使ったことは一度もありません。
ありがとうございます。
アンナ
サンプル共分散の計算について下記のJSLコードを試してください(スクリプトウィンドウを開き、貼り付けて実行します)。
sample_Q = 共分散 ( [ 35 36 40 、 38 37 40 、 48 41 42 、 50 46 54 ] ) ;
pop_Q = 共分散 ( [ 35 36 40 、 38 37 40 、 48 41 42 、 50 46 54 ] ) * 3 / 4 ;
見せる ( sample_Q 、 pop_Q ) ;
結果が出るはずです
sample_Q =
[ 54、25 31、33 37、33、
31,33 20,66 28,66、
37,33 28,66 45,33]。
pop_Q =
[ 40,6875 23,5 28、
23,5 15,5 21,5、
28 21,5 34]。
後者はあなたの計算に等しい。両者の差は3/4、すなわち(N-1)/ 1です。あなたが使った式は、実際には(不偏の)標本共分散ではなく、母集団の平均が分かっていれば正しい選択となる母集団共分散のためのものだと思います。Nの代わりにN-1を分母として使用するようにしてください。そうすれば、JMPの公式と同じ共分散推定値になるはずです。