« [R]例2.3 直流発電量の予測と予測区間(「線形回帰分析」(朝倉書店)pp.72-75) | トップページ | [R]平均勤続年数と所定内賃金(「線形回帰分析」(朝倉書店)pp.116-118) »

2023年11月 5日 (日)

[R]例2.4 年齢と最高血圧(「線形回帰分析」(朝倉書店)pp.75-82)その1

5つのモデルについて、lm関数を使用して回帰パラメーターを推定したところまで。

> dtf <- read.csv("table2_5.csv", header = TRUE)
> xxi <- as.double(dtf$AGE)
> yyi <- as.double(dtf$BP)
> # (AGE, BP)
> r <- lm(yyi ~ 1 + xxi)
> print(r)
Call:
lm(formula = yyi ~ 1 + xxi)
Coefficients:
(Intercept) xxi
112.3167 0.4451
> # (AGE, log(BP))
> r <- lm(log(yyi) ~ 1 + xxi)
> print(r)
Call:
lm(formula = log(yyi) ~ 1 + xxi)
Coefficients:
(Intercept) xxi
4.732231 0.003365
> # (log(AGE), log(BP))
> r <- lm(log(yyi) ~ 1 + log(xxi))
> print(r)
Call:
lm(formula = log(yyi) ~ 1 + log(xxi))
Coefficients:
(Intercept) log(xxi)
4.411 0.127
> # (AGE^2 / 10000, log(BP))
> r <- lm(log(yyi) ~ 1 + I(xxi ^ 2 / 10000))
> print(r)
Call:
lm(formula = log(yyi) ~ 1 + I(xxi^2/10000))
Coefficients:
(Intercept) I(xxi^2/10000)
4.7942 0.3863
> # (AGE^1.825 / 10000, log(BP))
> r <- lm(log(yyi) ~ 1 + I(xxi ^ 1.825 / 10000))
> print(r)
Call:
lm(formula = log(yyi) ~ 1 + I(xxi^1.825/10000))
Coefficients:
(Intercept) I(xxi^1.825/10000)
4.7883 0.8272

« [R]例2.3 直流発電量の予測と予測区間(「線形回帰分析」(朝倉書店)pp.72-75) | トップページ | [R]平均勤続年数と所定内賃金(「線形回帰分析」(朝倉書店)pp.116-118) »

R(本の計算を再現)」カテゴリの記事

コメント

コメントを書く

コメントは記事投稿者が公開するまで表示されません。

(ウェブ上には掲載しません)

« [R]例2.3 直流発電量の予測と予測区間(「線形回帰分析」(朝倉書店)pp.72-75) | トップページ | [R]平均勤続年数と所定内賃金(「線形回帰分析」(朝倉書店)pp.116-118) »