2009/10/13

Rで単回帰分析

# 回帰分析
test.lm <- lm(y ~ x, data = data.frame(testdata))

# 結果の表示
# 回帰式:y = Ax + B のAとBを表示したり、AとBのF統計量のP値を表示したり
summary(test.lm)

# AとBの取り出し
A = test.lm$coefficients[2]
B = test.lm$coefficients[1]

# F統計量のP値の取り出し
summary(test.lm)$coefficient[, "Pr(>|t|)"]

# AのP値
summary(test.lm)$coefficient[2, "Pr(>|t|)"]
# BのP値
summary(test.lm)$coefficient[1, "Pr(>|t|)"]

# 散布図上に回帰直線を描画
abline(test.lm)

0 件のコメント:

コメントを投稿