-
1 votesanswersviews
R中面板数据分析的大问题
我在R中进行面板数据分析存在一个主要问题 . 我正在尝试为我的分析做一个合并的OLS,固定效果和随机效应模型 . 我使用的数据集如下所示(这只是一个例子,我的数据集包含超过2900个数据): Y X1 X2 X3 X4 X5 X6 ID Year 0,03 0 0 0,05 0,01 0,00 -0,03 ... -
1 votesanswersviews
由于负方差,plm包的随机效应估计误差
我正在尝试估计一个面板数据模型,该模型具有随机效应Swamy和Ahora对旅游目的地市场份额数据的转换 . 请在https://github.com/Joseperles/Statistical-questions中找到附加的数据集 该模型根据以下代码估算: random<-plm(ld_Share~ld_Gdp+ld_Gdppc+ld_Gkf+ld_Cpi+ld_Fdi+ld_Exrate... -
31 votesanswersviews
在R中,使用qqmath或dotplot绘制来自lmer(lme4包)的随机效果:如何使它看起来很花哨?
qqmath函数使用lmer软件包的输出产生很好的随机效应图 . 也就是说,qqmath非常适合绘制层次模型中的截距及其在点估计周围的误差 . 下面使用名为Dyestuff的lme4包中的内置数据,给出了lmer和qqmath函数的一个示例 . 代码将使用ggmath函数生成分层模型和一个漂亮的图 . library("lme4") data(package = "l... -
1 votesanswersviews
GAM与mgcv包随机效应的方差 - 协方差矩阵
随机效应与lme4包的随机效应和方差 - 协方差矩阵提取如下: library(lme4) fm1 <- lmer(Reaction ~ Days + (1|Subject), sleepstudy) fm1.rr <- ranef(fm1,condVar=TRUE) fm1.pv <- attr(rr[[1]],"postVar") 我想知道如何用mgcv... -
1 votesanswersviews
测试R中固定效应和随机效应之间相互作用的重要性 . 正确的语法?
我正在使用R中的lmer进行混合模型,我对某些编码有点困惑 . 我测量了35个近交系(基因型)在果蝇中的雄性和雌性 Health 状况超过三个区块 . 我的响应变量是'适应度',n = 10个人/性别/行/块测试 . 性是固定的,块是随机的,并且嵌套在块内是随机的 . 我主要对性和线之间的相互作用感兴趣 . 因此我的模型看起来像 m1<-lmer(FitnessCured~Sex+(1|Bl... -
0 votesanswersviews
随机和固定效果之间的3种相互作用
我有一组看起来像这样的数据: rep stage line temp surv 1 L 149 18 0.6 2 L 149 18 0.7 3 L 149 18 0.25 1 A 149 18 1 2 A 149 18 ... -
26 votesanswersviews
如何在混合效果模型中获得系数及其置信区间?
在 lm 和 glm 模型中,我使用函数 coef 和 confint 来实现目标: m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous coef(m) confint(m) 现在我使用lme4包中的 lmer 函数为模型添加了随机效果 - 使用了混合效果模型 . 但是,函数 coef 和 confint...