理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

散布図に直線を描く

f:id:yoshida931:20190524000831p:plain:w400

x <- rnorm(20)  #乱数を使用してるので毎回変わります
y <- c(1:20)
plot(x,y)

xs <- seq(min(x), max(x), length=1000)

# xの影響を含むモデル
fit <- glm(y~x,family = poisson)
lines(xs, exp(fit$coef[1] + xs*fit$coef[2]))

# 切片のみのモデル
fit.null <- glm(y~1,family = poisson)
lines(xs, exp(fit.null$coef[1] + xs*0))