Taglibro de H

誤差を含んだXでの回帰をもうちょっとかんがえてみる。今回は2000回計算を繰り返して傾きの値の分布を推定する。


実際のXの値で回帰したときと、誤差を含んだXの値で直線回帰したときの傾きの値。


実際のXの値で回帰したときと、誤差を含んだXの値をつかって下のモデルでMCMC推定したときの傾きの値。

なんだかほとんどおなじになっている。

続きを読む

預金金利引き上げ(強制)で景気刺激とかいっている人がまあ大臣ということで、まあなんですね。

# ほかの人々ですけど、代替医療への傾倒とかEMへの傾倒なんかも、まあなんですね。

カワセミ


コサギ


エナガ


コガモ

カワセミは南浅川で2羽、城山川で1羽を確認。

続きを読む

コサギ

風が冷たかった。

誤差のあるxとポアソン分布するyとの関係をかんがえる。

こういうデータ。

set.seed(1234)
n <- 100
sd.x <- 0.1
beta <- 0
beta.x <- 0.5
x0 <- rnorm(n, 0, 1)
x <- x0 + rnorm(n, 0, sd.x)
log.y <- beta + beta.x * x0
y <- rpois(n, exp(log.y))

glmではこうなる。

> summary(glm(y ~ x, family = poisson))

Call:
glm(formula = y ~ x, family = poisson)

Deviance Residuals: 
       Min          1Q      Median          3Q         Max  
-1.8005747  -1.1532144   0.0004878   0.3378818   2.6383498  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  0.09312    0.09769   0.953    0.340    
x            0.41419    0.08584   4.825  1.4e-06 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 119.886  on 99  degrees of freedom
Residual deviance:  97.729  on 98  degrees of freedom
AIC: 261.55

Number of Fisher Scoring iterations: 5

続きを読む

アップル - ダウンロードから。

WidgetTerm
Dashboardウィジェットのターミナル。ちょっとシェルを使いたいときとかに便利そう。

Leopard ライフサイエンス辞書
定番のライフサイエンス辞書の2010年版。

今日の高尾山。陵南大橋から。

Model IIの回帰をMCMCでやってみる。

続きを読む

O'ReillyのNutshellシリーズ(日本語版だとクイックレファレンスシリーズ)のR本。Amazonに注文していたものがとどいた。前回、プログラミング方面の本のようだなどと書いてしまったが、プログラミングだけではなく、グラフィックから各種統計処理からBioconductorまで幅ひろく網羅されている。巻末のR Referenceが108ページもあったりする。





表紙の鳥は オウギワシ Harpia harpyjaであるとの説明が版元ウェブページにあり。

今年も集計してみる → 12月と2009年に見た鳥

続きを読む