Poisson模型——R软件

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我想在新数据data2中计算Poisson模型的AIC值，如何算呢？下面是我计算到新数据残差的代码：

#Poisson模型
#估计模型参数的数据
dat=read.csv("data1.csv",header=T)
#检验模型的新数据（预测）
da=read.csv("data2.csv",header=T)

poiss <- glm(y ~ x1+x2+x3+x4, data = dat, family = poisson)
summary(poiss)
logLik(poiss)
AIC(poiss)
prepoi<-predict(poiss, newdata=da)
residpoi=da$y-prepoi

谢谢了！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：poisson Iss r软件 son Data 模型 软件 poisson

AIC公式
[urlhttp://www.pinggu.org/bbs/viewth ... p;page=1#pid4216379][/url]

好的
呵呵不错谢谢

您好！我不知道如何在上面的代码上添加这公式，您帮我添进去吧！谢谢！

x <- c(rep(0,4),rep(1,3),rep(0,2),
rep(1,3),rep(0,2),rep(1,3),0,1,0)
n <- c(4,3,4,6,10,14,12,2,3,4,5,10,4,
5,9,5,10,4,7,4)
nx <- data.frame(n,x)
nx.glm <- glm(n ~ x, family="poisson", data=nx)
summary(nx.glm)
Call:
glm(formula = n ~ x, family = "poisson", data = nx)

Deviance Residuals:
     Min        1Q    Median        3Q       Max  
-1.75393  -0.49739   0.05042   0.46525   1.68649  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)   
(Intercept)   1.3610     0.1601   8.499  < 2e-16 ***
x             0.7908     0.1931   4.096  4.2e-05 ***
---
Signif. codes:  0 ｀***＇ 0.001 ｀**＇ 0.01 ｀*＇ 0.05 ｀.＇ 0.1 ｀ ＇ 1

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 32.616  on 19  degrees of freedom
Residual deviance: 14.502  on 18  degrees of freedom
AIC: 89.951

Number of Fisher Scoring iterations: 4

##自行计算AIC
AIC = -2*sum(dpois(nx.glm$y, nx.glm$fitted.values, log=TRUE)*nx.glm$prior.weights)+2*nx.glm$rank
AIC   # 89.95149