アールメカブ

• バックアップ一覧
• 差分 を表示
• ソース を表示
• バックアップ を表示
• R_Adj.r.squared へ行く。
  • 1 (2007-10-08 (月) 13:12:01)

• 追加された行はこの色です。
• 削除された行はこの色です。

[[日録]]
[[Baayen>R_Baayen]] から．
調整済み &mathml(R^2); が大きな値でも，残差からはモデルが適切ではないと判断される例．
#ref(baayen248.png,center,nowrap,nolink,残差からモデルは適切ではない)

# p.249
 # p.249
 alice.g <- alice.growth@data$data
 par(mfrow = c(2,2), bg = "white")
 plot(log(alice.g$Tokens), log(alice.g$Types)) #

 alice.g.lm <- lm( log(alice.g$Types) ~ 
   log(alice.g$Tokens))
     log(alice.g$Tokens))
 abline(alice.g.lm , col = "darkgrey")

 # r <- round(summary(alice.g.lm)$adj.r.squared,dig = 3)
 mtext(expression(paste(R^2, "= 0.998")))
 
 r <- round(summary(alice.g.lm)$adj.r.squared,dig = 3)
 mtext(substitute(R^2 == r.swap, list(r.swap = r))) 
 plot(log(alice.g$Tokens), resid(alice.g.lm))
 abline(h= 0)
 mtext("Residuals")

 z <- zipf.fnc(alice, plot = T)
 plot(log(z$rank), log(z$frequency), type = "S")
 z.lm <- lm(log(z$frequency) ~ log(z$rank))
 abline(z.lm, col= "darkgrey")
 r <- round(summary(z.lm)$adj.r.squared,dig = 3)
 mtext(substitute(R^2 == r.swap, list(r.swap = r))) 
 plot( log(z$rank), resid(z.lm))
 abline(h= 0)
 mtext("Residuals")