在logistic回归模型中使用sjplot和子集函数时遇到问题
问题解决了:我将R从4.2.0更新为4.2.1,并且该问题是自主解决的。任何遇到类似问题的人都可以首先尝试一下。 我正在尝试使用“子集”功能根据模型的标…
将随机效应纳入随机森林回归的时间序列数据
我有三个住宅的时间序列空气污染数据(例如PM2.5,CO2,TEMP,室外PM2.5),以及居民以二进制格式记录的活动日记(例如,在活动中进行1时,1时进行了1…
计算边缘效应:为什么ggeffect和ggemeans给出不同的答案?
示例 library(glmmTMB) library(ggeffects) ## Zero-inflated negative binomial model (m <- glmmTMB(count ~ spp + mined + (1|site), ziformula=~s…
对暴露的影响进行建模,该曝光对仅测量一次的结果进行了多次测量。是否可以?
我有一再的测量数据对173名患者,其中暴露于3次,并且一次连续结局。数据集还包括每个测量值的患者年龄,这就是数据的分类。 # Load packages library…
如何从模型结构中删除数据的某些元素
以下是我的数据示例: > dput(mydata1) structure(list(subject = c("E1", "E1", "E1", "E1", "E1", "E1", "E1", "E1", "E1", "E1", "E1", "E1"), blo…
ggeffects :: ggpredict()not nlme :: lme()返回人口级别的预测间隔
我正在尝试从 ggeffects获取总体级别预测间隔(PI):ggpredict()使用 type =“ re” 从 nlme:lme()模型。 GGPREDICT并未返回LME()模型的预期…
R:NLME软件包中具有自相关错误的混合效应模型:如何检查模型ARMA假设?
我正在使用包装nlme: model <- lme(y ~ x, random = ~ 1 | group, data = data, correlation = corARMA(form = ~ x | group, p=1, q=1) 比较AIC值的L…
研究的分层模型 - 我可以在r中使用混合效果吗?
这是一项研究,其中3只小鼠用作对照,而3只是突变体。 &gt; 2000测量是每只鼠标的。我想拟合一个模型如下: lmer(测量值〜鼠标 +(鼠标|处理),dat…
二项式重复数据中的glmmadaptive
我正在使用R中的GLMMADAPTIVE进行重复的二项式数据 。 。 所以我运行了这个模型。我的模型是:固定效果(时间*治疗+持续时间),随机效应仅具有随机…
use.u = true bootmer函数
我有一个关于多级模型的随机效果(BLUP)的置换置信区间的问题。 我目前正在使用 bootmer ,并且有一个参数 use.u = true ,它允许人们将蓝图视为固定…
lme4 :: lmer-使用随机斜率时(任何(bgrad&lt; 0))中的错误
我正在尝试使用 lmer 函数从 lme4 软件包中首次在R中运行线性混合效应模型,而且我一直在遇到我不喜欢的错误t知道如何解释。 现在,我有错误: if中的…
如何在R GGPLOT2中与相互作用拟合混合效应模型回归?
我与之合作的数据集具有互动术语。我想将模型与x轴中的交互项和y轴中的 y.var 一起拟合。我尝试关注这个例子关于如何在ggplot2中复制它太远(即绘图样…
尝试scale_linetype_manual以添加传奇,但不成功
遵循 sell the the the the the the the the the (name =“新传奇标题”)但是传奇标题仍然没有更改,我现在只想根据行类型生成Legend的行绘图,但没…