我有若干城市的时间序列数据,使用lme4这个包分析影响因素时出现几个问题
1、glmer做时间序列数据时能否像gee模型一样指定工作相关矩阵‘ar1’之类的,如果可以该加在哪?
2、模型的结果在输出summary的时候随机效应只显示了随机截距的方差,没显示残差的方差,在哪里输出?
3、使用MuMIn包的r.squaredGLMM输出R2时,输出了三组R2,三组之间有何区别?
ggg1 <- glmer(case(1|id),data = sj,offset = log(SUM),family="poisson")urban2+(1|id),data = sj,offset = log(SUM),family="poisson")
ggg2 <- glmer(case
ggg3 <- glmer(caseurban2+ccc+(1|id),data = sj,offset = log(SUM),family="poisson")urban2*ccc+(1|id),data = sj,offset = log(SUM),family="poisson")
ggg4 <- glmer(case
summary(ggg1)
anova(ggg2,ggg1)
r.squaredGLMM(ggg4,ggg1)
Random effects:
Groups Name Variance Std.Dev.
id (Intercept) 0.7456 0.8635
Number of obs: 390, groups: id, 30
++++++++++++++++++++++++++++++++++
R2m R2c
delta 0 0.19335449
lognormal 0 0.34531503
trigamma 0 0.06466862