自身の前後で地価の説明変数の影響がどう変化したかなという興味の元統計を勉強し始めました、最近。それでRにてlm()で重回帰分析してsummary()をしようとしています。
その途中でけっつまづいたので助けてください!
aが被説明変数でb1,b2が説明変数として
yy<-read.csv("tes.csv")
a<-yy[,1]
b1<-yy[,2]
b2<-yy[,3]
b3<-yy[,4]
y<-data.frame(a,b1,b2,b3)
y.lm<-lm(a~.,data=y)
『summary(y.lm)最終的にはこれをやりたいのだがlm()のとこで』
以下にエラー lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
因子保管モードを変更するのは不正です
追加情報: 警告メッセージ:
In model.response(mf, "numeric") :
因子応答を持つ type="numeric" の使用は無視されるでしょう
とかえされるんですがなにがいけないんですかね ?
a <- c(148, 160, 159, 153, 151, 140, 158, 137, 149, 160)
b1 <- c(41, 49, 45, 43, 42, 29, 49, 31, 47, 47)
b2 <- c(41, 159, 153, 151, 140, 31, 47, 47, 47, 47)
b3 <- c(41, 159, 153, 151, 140, 31, 47, 47, 47, 47)
y.lm<-lm(a~.,data=y)
summary(y.lm)
これはいけるんですがね。ちなみにcsvのなかは後者のような文字列があるのですが...
助けてください!
それ以外にもRで壁にあたった時自分で対処できるようになれそうなサイトorぐぐり方も指南していただけると今後が幸いです。
No.1ベストアンサー
- 回答日時:
読み込んでいる yyのデータが数字として認識されてないのでは?
実際 書かれている後半の例は確かに動きます b2,b3 が相関が1なので b3はNAになってますが、これは、多重共形性というか 行列式が0でそれでわるからですよね。これは今関係ないです。
たぶん
tes.csv に問題があって 数として認識されてなくて、factor か 文字として認識されているんだと思いますが。
それと間違いではないですが yy は すでにdata.frame なので わざわざ yを作る必要はないかと
Rのサイトですが
cse.naro.affrc.go.jp/takezawa/r-tips/r/24.html
が便利です(今日はメンテナンスでみれませんが)
> a <- c(148, 160, 159, 153, 151, 140, 158, 137, 149, 160)
> b1 <- c(41, 49, 45, 43, 42, 29, 49, 31, 47, 47)
> b2 <- c(41, 159, 153, 151, 140, 31, 47, 47, 47, 47)
> b3 <- c(41, 159, 153, 151, 140, 31, 47, 47, 47, 47)
> y<-data.frame(a,b1,b2,b3)
> y
a b1 b2 b3
1 148 41 41 41
2 160 49 159 159
3 159 45 153 153
4 153 43 151 151
5 151 42 140 140
6 140 29 31 31
7 158 49 47 47
8 137 31 47 47
9 149 47 47 47
10 160 47 47 47
> y.lm<-lm(a~.,data=y)
> y.lm
Call:
lm(formula = a ~ ., data = y)
Coefficients:
(Intercept) b1 b2 b3
108.03255 0.96957 0.02844 NA
> summary(y.lm)
Call:
lm(formula = a ~ ., data = y)
Residuals:
Min 1Q Median 3Q Max
-5.9392 -1.5571 -0.5075 2.5066 5.0608
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 108.03255 7.33188 14.735 1.59e-06 ***
b1 0.96957 0.18281 5.304 0.00112 **
b2 0.02844 0.02309 1.232 0.25777
b3 NA NA NA NA
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 3.6 on 7 degrees of freedom
Multiple R-squared: 0.8504, Adjusted R-squared: 0.8077
F-statistic: 19.9 on 2 and 7 DF, p-value: 0.001294
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
イベントスタディ分析
-
4つのサイコロを同時に振って...
-
2つのサイコロの目の合計が偶数...
-
一般常識を教えてください。1割...
-
五分を6回連続で外すのはなん...
-
重複順列
-
数学の質問です。 1〜6までの番...
-
数学の問題で4C0の答えを教えて...
-
統計学の問題なんですが・・・
-
2つのサイコロを投げた時、目...
-
サイコロ5個を同時に投げて任意...
-
確率で「試行の独立」「事象の...
-
大中小3個のサイコロを同時に投...
-
1個のサイコロを3回続けて投げ...
-
3個のサイコロを同時に1回振る...
-
75%を3回連続で引かない確率
-
男子4人と女子4人が輪の形にな...
-
3つのサイコロを振って出た目の...
-
数学の確率の超有名問題で、サ...
-
異なる4つのサイコロを同時に投...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
近似直線の傾き同士の相関
-
クロスカントリー分析・クロス...
-
イベントスタディ分析
-
ローリング回帰とはなんですか?
-
<計量経済学>定数項のない回...
-
統計ソフトRの重回帰分析で困っ...
-
一次同時(次?)について質問...
-
回帰分析の結果
-
回帰分析の時に対数をとる意味は?
-
重回帰分析におけるP値とデー...
-
多変量解析の手法
-
重回帰のダミー変数について
-
最小2乗法→回帰分析→?
-
相関係数の特徴
-
回帰分析の際のデータの指数化
-
重回帰分析をして仮説検定を行...
-
回帰分析:理論値と実績値の差...
-
エクセルの回帰分析について
-
4つのサイコロを同時に振って...
-
2つのサイコロの目の合計が偶数...
おすすめ情報