統計解析フリーソフト　Ｒ　【第３章】

統計解析フリーソフト　Ｒ　【第３章】at MATH

統計解析フリーソフト　Ｒ　【第３章】 - 暇つぶし2ch■コピペモード
□スレを通常表示
□オプションモード
□このｽﾚｯﾄﾞのURL
■項目テキスト

672:１３２人目の素数さん
 10/06/23 14:34:40 
質問です。ご回答お願いします。 
 
10行2列のマトリックスAと20行2列のマトリックスBを使って 
回帰分析を行いたいんですが、エラーが出ます。 
 
具体的には、A,Bの一列目を従属変数(濃度)、二列目を独立変数(スペクトル)にしてます。 
関数lmを使ってAで作った検量線に、関数predictを使ってBの二列目を代入したい 
(スペクトルの値から濃度を予測したい)と考えているのですが、 
「データ数が違う」という旨の警告が出てできません。 
 
分かりにくいかもしれませんが、宜しくお願いします。

673:１３２人目の素数さん
 10/06/23 15:33:28 
>>672 
Rjpwikiの方が親切に回答してもらえると思うぞ。 
> A <- data.frame(matrix(1:20,nrow=10)) 
> head(A) 
  X1 X2 
1  1 11 
2  2 12 
3  3 13 
4  4 14 
5  5 15 
6  6 16 
> B <- data.frame(matrix(runif(40),nrow=20)) 
> head(B) 
         X1        X2 
1 0.6353921 0.2679418 
2 0.2917369 0.3124313 
3 0.2857031 0.2795522 
4 0.1951349 0.4510532 
5 0.4426567 0.4894139 
6 0.5981174 0.3142484 
> predict(lm(X1~X2,A),B) 
        1         2         3         4         5         6         7         8  
-9.732058 -9.687569 -9.720448 -9.548947 -9.510586 -9.685752 -9.718648 -9.481105  
        9        10        11        12        13        14        15        16  
-9.477939 -9.628372 -9.697301 -9.711742 -9.717192 -9.791366 -9.910725 -9.650475  
       17        18        19        20  
-9.903119 -9.638157 -9.547609 -9.051490  
ってこと？具体的なコード（再現性のあるサンプルコード）を示してくれないと、 
なんで出来ないのかこちらには分からない。

次ページ