Google 到底想說什麼?
December 9, 2010
一門英烈
December 9, 2010
鴉片有益論
孔子和平獎
袓孫皆俊傑
磕頭相輝映
孔子和平獎
袓孫皆俊傑
磕頭相輝映
BBC 新聞部表示
November 25, 2010
R 打敗 Python?
November 14, 2010
賓拉登原來這般仙風道骨啊
November 9, 2010
GLM, 沒你想的那麼好
April 11, 2010
作個很簡單的R
實驗:
beta1 <- rnorm(4) C <- cbind(rexp(30, 3.2), abs(rnorm(30, 3.2, 1.4)), rgamma(30, 3.2, 1.4)) lambda <- abs(cbind(rep(1, 30), C)%*%beta1) + rnorm(30) y <- rep(NA, 30) for (j in 1:30) y[j] <- rpois(1, lambda[j]) glm(y~C, family="poisson")
你可以看到結果有多離譜。
這代表其實GLM
不是那麼可靠的方法。 只要在mean的部份,有那一點點小小的變動,結果就會整個跑掉,這對想要套套公式生paper的學者們,(尤其是那些非physics, mathematics, engineering的學者們), 這真是一個壞消息:每天用來套的公式原來是那麼的不可靠
可惜東榮公司並沒有投資統計方法論的研究,我只好自己來了,特聘我自己開設巫毒統計學專欄,以饗那些時間多到殺不完的待業,失業或結業人士。