No.1ベストアンサー
- 回答日時:
企業でSQCを推進する立場の者です。
こんな感じでしょうか。
これに信頼区間とか予測区間を加えようとすると、
絶対にlmでやった方が楽ですけど・・・。
rm(list=ls()) # メモリー上の残骸を消す
# 適当に点列を作ってプロットする
# 仮に、β=1.5,y切片=3 とする。
x <- rnorm(0:100,5,2.5)
y <- 3 + 1.5 * x + rnorm(0:100,0,0.5)
xy <- data.frame(cbind(x,y))
plot(xy,pch=16,xlim=c(-5,15),ylim=c(-5,25),las=1)
abline(h=0)
abline(v=0)
xy.c <- scale(xy,,scale=F) # 中心化
# 偏差平方和を求める
# XTXで求めるという、定石を使う
SS <- t(as.matrix(xy.c)) %*% (as.matrix(xy.c))
Sxx <- diag(SS)[1] # XTXの左上部分
Sxy <- SS[lower.tri(SS)] # XTXの左下部分
# 回帰のパラメータを求める
beta1 <- Sxy / Sxx
beta0 <- mean(y) - beta1 * mean(x)
# 回帰線を追加する
abline(beta0,beta1)
# 念のため検算
names(xy) <- c("X","Y")
lm(Y~.,data=xy)
beta0
beta1
この回答へのお礼
お礼日時:2017/06/19 18:37
大変分かりやすい解説ありがとうございました。まだ勉強中の身でして、lm関数で細かく見ると何が起きているのかを理解する為、やってみたのですが上手くいかず困っておりました。本当にありがとうございます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
x二乗-3xy+y二乗 この因数分解...
-
絶対値の二乗の思考過程 |x-y|^2
-
acrobat8(standard)で図形を書...
-
(3x-y)(2x+5y)の展開の仕方を簡...
-
0≦z≦xy、x^2+y^2≦a^2 、x≧0 、y...
-
x+y=5 xy=-3 のとき、x二乗-3xy...
-
数学の質問です!! 「2つの素...
-
数学の問題で!
-
なんどもすいません
-
とても急いでいます!
-
xy=0ならばx=0またはy=0 を証明...
-
素イデアルの冪と準素イデアル
-
数学、因数分解
-
高1 数II x+y+z=−1、xy+yz+zx+...
-
x^3+y^3+z^3
-
因数分解
-
合成関数の問題です。sinxy=x+y...
-
命題ー必要十分条件
-
平方根の計算
-
チャート例題185について 画像...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
絶対値の二乗の思考過程 |x-y|^2
-
とても急いでいます!
-
x^3+y^3+z^3
-
x²+xy-4x-y+3 を因数分解して...
-
x二乗-3xy+y二乗 この因数分解...
-
acrobat8(standard)で図形を書...
-
x²-xy-y+xの因数分解って、どの...
-
x+y=5 xy=-3 のとき、x二乗-3xy...
-
数Ⅰ「xとyについて降べきの順に...
-
xy=0ならばx=0またはy=0 を証明...
-
因数分解
-
数式で項のアルファベットの順...
-
高1 数II x+y+z=−1、xy+yz+zx+...
-
【代数学】可換群の証明
-
画像の問題の解説のところで。 ...
-
Q(x+y, x^2+y^2)の存在する範囲...
-
高校一年の数Ⅰに苦戦しています...
-
xy−x−y+1これの因数分解の仕方...
-
ラグランジュの未定乗数を二つ...
-
数学・ベクトル場の線積分・面...
おすすめ情報