「統計学関連なんでもあり」の過去ログ--- 048

No.23131　回帰曲線の描き方　　【コロン】　2021/09/06(Mon) 14:22

青木先生，お世話になっております。

今回の質問は以前お尋ねしたものと同じです。ワード検索での私の検索の仕方がまずいのかなかなかヒットせず，過去ログを一つ一つ丁寧に見ていったつもりですが，それでも見つけることが出来ませんでしたので，失礼を承知で同じ質問をさせてください。

以下のデータに対して，回帰曲線を描きたいのです。出典は以下のものです。2次，3次をどのように描けばいいのでしょうか。

小塩真司.2007.『実践形式で学ぶSPSSとAmosによる心理・調査データ解析』東京図書（pp.32-37）

motivation	score
6.13	5.59
2.30	2.78
4.66	6.31
7.28	4.49
3.30	7.52
4.12	7.34
2.49	2.77
6.35	6.68
8.23	2.36
2.41	4.56
2.55	3.19
7.26	2.38
5.60	8.74
5.40	6.84
5.27	6.93
2.51	4.66
6.05	6.06
3.04	5.21
6.66	3.63
5.74	5.83
7.67	3.39
6.59	6.35
6.69	4.09
7.64	4.67
4.06	6.77
7.19	5.30
4.90	7.15
4.87	7.00
8.05	2.69
3.68	6.87
2.88	5.49
5.11	5.47
7.35	2.77
3.21	5.18
3.88	7.11
2.91	3.09
6.03	8.18
5.18	7.95
4.50	4.79
3.26	7.26

No.23132　Re: 回帰曲線の描き方　　【aoki】　2021/09/06(Mon) 23:30

データがデータフレーム df に付値されているとします

以下でよろしいかと思います

a = lm(score ~ motivation + I(motivation^2), data=df)
x = seq(min(df$motivation), max(df$motivation), length=500)
y2 = predict(a, newdata=list(motivation=x))
b = lm(score ~ motivation + I(motivation^2) + I(motivation^3), data=df)
y3 = predict(b, newdata=list(motivation=x))
plot(score ~ motivation, data=df, pch=19, col="red")
lines(y2 ~ x, data=df, col="blue")
lines(y3 ~ x, data=df, col="brown")
legend("topright", c("order 2", "order 3"), col=c("blue", "brown"), lty=1)

No.23134　Re: 回帰曲線の描き方　　【コロン】　2021/09/07(Tue) 12:50

青木先生

お忙しい中，ありがとうございました。また同じ質問をしてしまい申し訳ございませんでした。