« 石狩湾新港 2009.0125 | トップページ | カラスのねぐら »

glm() 自動Akaike weight プログラム

生態屋さんで統計やっている人が一度は見たことあるであろう有名サイト『飯島の雑記帳

ここで、glm()のakaike weightを求めるプログラムが配布されている。

しかし、自分はこのプログラムを試してみたことがあるのだが、エラーを吐き出してうまく動かなかった。特に使わなかったので、長らく放置していたのだが、久しぶりにRネタの記事を書くためにちょっと書きなおしてみる。

といっても使うのはTaglibro de Hにある説明変数シラミつぶし関数である。

「glm_aw.r」をダウンロード

↑をクリックしてファイルダウンロードしてください。

使い方
glm.aw(y, indeps, dis, offset=NULL,data=NULL)
y=目的変数
indeps=説明変数 c("x1","x2")のように入れていく
dis=分布形
offset=offset項(なくてもよい)

set.seed(13)
n <- 40
x1 <- runif(n, 0, 1)
x2 <- x1 + rnorm(n, 0, 0.2)
x3 <- x1 + rnorm(n, 0, 0.2)
x4 <- x1 + rnorm(n, 0, 0.2)
x5 <- x1 + rnorm(n, 0, 0.2)
z  <- x1 + rnorm(n, 0, 1) #offset用
y  <- x1 * 10 + rnorm(n, 0, 1)
dat<- data.frame(y, x1, x2, x3, x4, x5, z)

result<- glm.aw(y, c("x1","x2","x3","x4","x5"),gaussian,z, dat)
head(result)

   as.character.tm.  vec.aic     dAIC akaike.weight
2            x1 + 1 144.7587 0.000000    0.23329566
4       x1 + x2 + 1 145.9729 1.214172    0.12713133
10      x1 + x4 + 1 146.2318 1.473145    0.11169080
6       x1 + x3 + 1 146.7564 1.997742    0.08592162
18      x1 + x5 + 1 146.7587 1.999999    0.08582473
12 x1 + x2 + x4 + 1 147.4386 2.679901    0.06109047

2008.1.28 追記
引数dataを追加。

|

« 石狩湾新港 2009.0125 | トップページ | カラスのねぐら »

統計/R」カテゴリの記事

コメント

この記事へのコメントは終了しました。

トラックバック


この記事へのトラックバック一覧です: glm() 自動Akaike weight プログラム:

« 石狩湾新港 2009.0125 | トップページ | カラスのねぐら »