[B! glm] syou6162のブックマーク

syou6162 id:syou6162

glmに関するsyou6162のブックマーク (2)

2005-08-27
顕微鏡でみた小さな生物たちの博物館。見ていて飽きない。クンショウモとかいいよね。割礼がエイズ感染率を下げるかも？という報告を受けて、その是非を語る Slate の記事。書いたのは自分の子供2人にも生まれてすぐに割礼を受けさせた女性記者で、「セックス時の快感はどうなのよ？」という（親として心配されるにはなんとも複雑な）疑問について、大人になってから割礼を受けた読者からの意見を求めている。 Rでロジスティック回帰を行なうときには一般化線形モデルの関数 glm() を使う。ロジスティックモデルの場合、リンク関数はロジットなので glm() の引数に正しくは family=binomial(link="log it") を使うのだけど、分布にbinomial（二項分布）を指定した場合のデフォルトはロジットなので (link="log it") は省略していい。またもや Kutner et al. "
syou6162 2008/04/25
glm

r
リンク
Brownleeの排気損失率データ
出処:回帰分析の基礎、早川毅著、朝倉出版 (経済企画庁調査局、消費者動向調査による) >年度<-c(1966:1984) >普及率<-c(0.003,0.016,0.054,0.139,0.263,423,0.611,0.758,0.859,0.903,0.937,0.954,0.978,0.978,　0.982,0.985,0.989,0.988,0.992) > tv<-glm(普及率~年度,family=binomial) 関数glmによる当てはめ値はfittedで返すことができる。図2にカラーテレビの普及率の実測値と関数glmを用いたロジスティック回帰モデルの予測値の折れ線プロットを示す。 >plot(年度,普及率,type="l") >lines(年度,fitted(tv),lty=2,col="red",lwd=2) >legend(1975,0.5,c("実測
syou6162 2008/04/13
r

glm
リンク
1