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
ってこと？具体的なコード（再現性のあるサンプルコード）を示してくれないと、
なんで出来ないのかこちらには分からない。
674:１３２人目の素数さん
10/07/04 01:35:42
nnet使ってるけど、これって出力層のユニット数は変更できないの？
675:名無しさん＠そうだ選挙に行こう
10/07/11 17:05:40
fit <- lm(Y ~ X1 + X2 + X3, df)
anova(fit)
で出てくるF値は、t値の二乗になっていないですよね。
このF値は統計学的にどういう意味のあるF値なんですか。
676:名無しさん＠そうだ選挙に行こう
10/07/11 20:31:09
XPSP3でソースコードからshlibをenableにしてRをインストールする方法を知ってる人いませんか？
Rtoolsを利用して通常のtar→makeまではできたのですが、UNIXでの./configure --enable-R-shlib のようなコマンドがわからず困っています…
677:名無しさん＠そうだ選挙に行こう
10/07/11 20:39:55
ＳＡＳ＞＞ＳＡＳＳ＞Ｒ＞えくせろ
678:名無しさん＠そうだ選挙に行こう
10/07/11 21:41:45
SPSSだろ
679:１３２人目の素数さん
10/07/13 10:02:28
>>677
少なくともあなたがコンピュータに無縁なのはよく分かった
680:１３２人目の素数さん
10/07/14 21:58:23
こんばんは、質問させてください。
factanalを使って1因子の因子分析をしているのですが、初期解を求めようとすると
「以下にエラー solve.default(cv) :
システムは数値的に特異です：条件数の逆数 = 2.13584e-17 」
というエラーが出てしまいます。
prcompを使うと正しく計算されます。
しかし「前回課題で使用したデータ」を「factanalを使って1因子分析」するのが条件なので、
関数やデータをいじるわけにもいきません。
どうしたら解決できるのでしょうか。
681:１３２人目の素数さん
10/07/14 22:06:13
psppというのもあるよ。
682:１３２人目の素数さん
10/07/14 22:18:12
>>681
Rに比べると怪しげに見える。
よく知らんけど、Rはgccと同じくらいのレベルで標準を確立したりしていないか？
683:１３２人目の素数さん
10/07/14 23:13:21
>>680
それってなぜエラーが出るか考えさせる問題ってことはないよね？
684:１３２人目の素数さん
10/07/15 11:33:10
>>681
まだプロジェクトが続いていたのか。
昔試したときには、GUI化される前のSPSS Baseのコピーって
かんじだったけど。公式ページを見たら、psppもGUI化されたみたいだな。
685:１３２人目の素数さん
10/07/17 20:20:27
which関数の条件文ってAND条件使えないの？
例えば、
x<-1:10
which(x<5&&x>2)
ってしたら
integer(0)
ってでてくる
686:１３２人目の素数さん
10/07/17 20:25:44
&は1個でいい
687:１３２人目の素数さん
10/07/20 10:16:24
>>685
&と&&は意味が違うので注意って入門書に書いていなかった？
688:１３２人目の素数さん
10/08/01 12:35:35
>>684
>GUI化
結構いい線いってると思うけど、本物と比べると出来ないことが多すぎて
代わりにはなりませんな。
689:１３２人目の素数さん
10/08/01 13:09:12
batファイルからperlを呼び出して前化工して、
次にRscriptを呼び出して結果を出力して、
後始末して、ということをしてるけど、
batファイルというのがどうも不格好なので、
rファイルだけで全部やらせられない?
690:１３２人目の素数さん
10/08/04 23:37:19
データフレームの特定の列をfactorにしたい場合って、こう書くしかないのかな？
data$colulmn <- as.factor( data$column );
691:１３２人目の素数さん
10/08/05 22:22:59
age
692:１３２人目の素数さん
10/08/11 15:48:15
sage