自身の前後で地価の説明変数の影響がどう変化したかなという興味の元統計を勉強し始めました、最近。それで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で質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# gcc のコンパイラオプションについて 2 2022/12/21 17:11
- 経済学 IS-LM曲線とAD-AS曲線の相関関係を教えていただけませんか? 5 2022/11/27 12:23
- 国産車 レクサスが日本でもLMを発売するようですが、価格予想として1600万程度と噂が? しかも素人の位置づ 4 2023/05/07 21:50
- 統計学 加重最小二乗法=①「変数を自然対数変換」=②「誤差項の分散の逆数を重み付け」? 8 2022/11/26 11:15
- 経済学 経済学について 1 2023/05/14 21:13
- その他(趣味・アウトドア・車) 【時計】ベルト幅について(現行:カルティエ タンクLM) 2 2022/04/03 20:29
- 経済学 is-lm-bp分析と国際収支の件で質問 1 2023/05/09 22:17
- 経済学 大学のマクロ経済の質問です。 2 2023/01/30 12:56
- Visual Basic(VBA) 【VBA】特定の文字で改行(次の行)に行きたい。 3 2022/04/11 17:20
- その他(プログラミング・Web制作) データ解析ソフトRでのファイル入力read.csvがエラーになります 7 2022/03/27 22:11
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
イヤーダミーについて
-
BINGが間違えた、とっても簡単...
-
重複順列
-
4つのサイコロを同時に振って...
-
確率の問題 数学と実生活と
-
統計学の問題なんですが・・・
-
一般常識を教えてください。1割...
-
調査結果は、標準誤差率が何%...
-
3σについて教えてください(基...
-
6人が円形のテーブルを囲んで座...
-
P(A|B)などの読み方
-
統計学的に信頼できるサンプル...
-
5人の人々を3人と2人のグループ...
-
サイコロの目が10回目に「1」が...
-
00~99、AA~ZZの組み合わせっ...
-
a.b.c.d.eの5個から3個を選んで...
-
数学の質問です。 一枚の硬貨を...
-
サイコロを6回振って1回も1が...
-
効果量のマイナス表示に関して
-
サイコロを三つ同時に振るとき ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
近似直線の傾き同士の相関
-
クロスカントリー分析・クロス...
-
イベントスタディ分析
-
計量経済の問題
-
重回帰分析におけるP値とデー...
-
GDP統計の「寄与度」、「寄与率」
-
回帰分析の時に対数をとる意味は?
-
因子分析における因子得点の解...
-
相関係数の特徴
-
回帰分析の際のデータの指数化
-
回帰分析の目的変数と説明変数...
-
夏に気温が上がるとビールの売...
-
統計ソフトRの重回帰分析で困っ...
-
ローリング回帰とはなんですか?
-
回帰分析:理論値と実績値の差...
-
モンテカルロシュミレーション...
-
4つのサイコロを同時に振って...
-
重複順列
-
統計学の問題なんですが・・・
-
2つのサイコロの目の合計が偶数...
おすすめ情報