gpt4 book ai didi

R stargazer,lme4和lmerTest不兼容

转载 作者:行者123 更新时间:2023-12-04 12:00:24 33 4
gpt4 key购买 nike

R新手在这里。我在使用lmerTeststargazer时遇到问题。我在这里按照教程进行操作,以使stargazer与R中的lme4一起使用。

http://svmiller.com/blog/2015/02/quasi-automating-the-inclusion-of-random-effects-in-rs-stargazer-package/

我没有任何问题地运行了这个示例。

library(lme4)
library(stargazer)
data(cake)
summary(M1 <- lmer(angle ~ temp + (1 | replicate) + (1|recipe:replicate), cake, REML= FALSE))
summary(M2 <- lmer(angle ~ factor(temperature) + (1 | replicate) + (1|recipe:replicate), cake, REML= FALSE))
anova(M1,M2)
stargazer(M1, M2, style="ajps", title="An Illustrative Model Using Cake Data", dep.var.labels.include = FALSE,
covariate.labels=c( "Temperature (Continuous)", "Temperature (Factor $<$ 185)", "Temperature (Factor $<$ 195)", "Temperature (Factor $<$ 205)", "Temperature (Factor $<$ 215)", "Temperature (Factor $<$ 225)")
)

在这种情况下,如果我包括 lmerTest包,则 stargazer不再起作用。
library(lme4)
library(lmerTest)
library(stargazer)
data(cake)
summary(M1 <- lmer(angle ~ temp + (1 | replicate) + (1|recipe:replicate), cake, REML= FALSE))
summary(M2 <- lmer(angle ~ factor(temperature) + (1 | replicate) + (1|recipe:replicate), cake, REML= FALSE))
anova(M1,M2)
stargazer(M1, M2, style="ajps", title="An Illustrative Model Using Cake Data", dep.var.labels.include = FALSE,
covariate.labels=c( "Temperature (Continuous)", "Temperature (Factor $<$ 185)", "Temperature (Factor $<$ 195)", "Temperature (Factor $<$ 205)", "Temperature (Factor $<$ 215)", "Temperature (Factor $<$ 225)")
)

Error in objects[[i]]$zelig.call :
$ operator not defined for this S4 class

我确实很想使用 stargazer,但是我的实验要求使用 merModLmerTest对象,但 stargazer不支持该对象。有人知道解决方法吗?将 merModLmerTest对象转换为兼容的 lmerMod对象有多困难?

最佳答案

这是一个简单的解决方法。将lmerTest的输出(位于merModLmerTest类中)转换为lmerMod类。这将与观星者兼容。

class(model) <- "lmerMod"

关于R stargazer,lme4和lmerTest不兼容,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/31319030/

33 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com