lme4::glmer与Stata的melogit命令比较
最近我一直在尝试对相对较大的数据集拟合很多随机效应模型，比如说大约有50,000人（或更多）在最多25个时间点上进行观察。由于样本量很大，我们需要调整很多预测变量 - 也许是50个左右的固定效应。我使用R中的lme4 :: glmer 将模型拟合到二元结果，每个主题都有随机截距。我不能具体说明...
R stargazer、lme4和lmerTest不兼容性问题
我是R的新手。我在使用lmerTest和stargazer时遇到了问题。我正在按照此处的教程操作，以使stargazer与R中的lme4一起工作。 http://svmiller.com/blog/2015/02/quasi-automating-the-inclusion-of-rando...
如何在Python的statsmodels线性混合效应模型中拥有多个组？
我正在尝试使用Python的statsmodels线性混合效应模型来拟合具有两个随机截距（例如两个组）的模型。我无法弄清楚如何初始化模型，以便我可以这样做。 以下是示例。我有类似以下数据的数据（取自here）： subject gender scenario attitude ...
在Linux上安装nloptr - 致命错误：没有找到nlopt.h文件或目录。
在Linux（Ubuntu 14.04.4 LTS）上安装nloptr包时，遇到了困难，希望得到帮助。我查看了许多问题，但没有找到解决办法。 由于nloptr安装出现非零退出状态，我无法在R（版本3.3.1 / Rstudio版本0.99.902）中安装lme4软件包。当我尝试时： ...
如何在lmer()中使用update()更新随机部分？
我想使用update()函数更新模型的随机部分，特别是添加一个随机效应。大多数例子（help("update")，help("update.formula")，lme4:用 R 进行混合效应建模）都关注于模型的固定部分。在下面的示例中，如何使用update()从fm0到fm1？ librar...
nlmer 长期数据
我一直在使用“nlme”软件包中的R正畸数据集。只需使用install.packages("nlme");library(nlme);head(Orthodont)查看。该数据集由27个儿童随时间测量下垂体和翼腭裂之间的距离组成。 使用lme4软件包，我可以拟合一个非线性混合效应模型，使用逻...
如何处理线性混合模型（lme4）中的奇异适配？
我正在进行一项关于鸟类的研究，使用变量“巢穴”作为随机变量运行多个线性混合模型。但是，在某些模型中，我遇到了所谓的“奇异拟合”：我的巢穴随机变量的方差和标准误差均为0.00。 一些背景信息：我正在研究野生鸟类在嘈杂环境中生活对一些氧化应激参数的影响。为此，我们为每个巢穴的雏鸟取了一次血样进行...
使用glmer解决警告信息的方法
与许多其他人一样，我在运行使用lme4软件包的glmer函数的模型时遇到了困难。 以下是我的模型： model = glmer(depvar ~ variety*cover+amplitude+time+ (1|pp) + (1|stim), data = datafile, fam...
将一个命名的模型列表传递给anova.merMod
我希望能够将一个命名的模型列表（merMod对象）传递给anova()函数，并在输出中保留模型名称。这在使用mclapply()并行运行一批缓慢的模型（如glmer）时特别有用。我想到的最好方法是对模型列表进行去命名处理，然后使用do.call，但这并不理想，因为我的模型可能会被命名为“mod...
在nlme和lme4中拟合相同的模型
数据来自这里。 library(nlme) dat0 <- read.table("aids.dat2",head=T) dat1 <- dat0[dat0$day<=90, ] # use only first 90-day data dat2 <- dat1[...