1回答
我有一个包含mod <- lmer(value ~ time + study + (1|participant))形式的因子的模型，其中我对研究系数并不真正感兴趣，但由于各种原因，我希望将其保持不变。但我想根据mod进行预测，而不受研究系数的影响。我已经对学习因子进行了总和编码，所以我真的只需要一种方法来删除学习系数。 从系数手动计算预测似乎有点麻烦。我想一种替代方法是将研究系数</e
浏览 18提问于2020-11-22得票数 0
2回答
谁能告诉我为什么从具有随机斜率的lmer模型中提取的斜率系数与从适合同一数据集的lmList模型中提取的斜率系数之间存在偏差？
谢谢..。
浏览 5提问于2014-03-15得票数 1
我对Julia语言(在Windows7机器上使用版本0.6.0 (2017-06-19 13:05 UTC)官方发行版x86_64-w64-mingw32 )非常陌生。我有R背景，发现对于非常大的数据集，混合模型的实现速度很慢(n>2,000,000，p> 100)。因此，我寻找替代方案，当涉及到评估时间时，Julia似乎闪电般的速度。我想在这里提出的问题是关于dmbate的MixedModels.jl包。尽管与lme4相
浏览 2提问于2017-08-02得票数 1
我对R还是新手，我正试图在数据集中的90列上循环一个混合模型。我的数据集看起来像下面的数据集，但是有90个预测器，而不是7个，我需要在连续模型中将其评估为固定效果。然后，我需要存储模型输出(系数和P值)，以最终构造一个图，总结每个预测器的大小效应。我知道lme4混合模型对P值估计的讨论。+ (1|id) + (1|experiment), data = mydata)
我的
浏览 4提问于2022-08-27得票数 1
回答已采纳
我正在进行面板回归，并决定改用lm()，因为plm()对于测试数据(以及linearmodels中的linearmodels)没有一个很好的predict()函数，对于我来说，作为计量经济学新手，lme4语法是不直观的。我想使用lm并预测未知的数据。我的lm()方程看起来像这样，作者是一个固定的效果很明显它
浏览 3提问于2022-03-11得票数 -1
回答已采纳
1回答
我一直在研究国际足联的一些数据，我想看看运球和控球之间的关系。我在下面创建了一个线性回归，然而，当我使用.score方法获得R2时，我得到的输出为-1.9345。这显然是错误的，因为R2不能大于1。有人能用外行的术语解释一下我哪里错了吗？
此外，我还注意到，当我执行metrics.r2_score(x_test,y_test))时，我得到了正确的R2。y_trai
浏览 0提问于2021-05-20得票数 0
1回答
我正试图预测在我的二项式数据上运行的glmer模型的值(x轴中的天数)。总活着和全部死亡都是计数数据。这是我的模型，以及下面相应的步骤。，我们已经对过分散进行了解释。my.dredge.models<-get.models(dredge.models)
然后，我们进行了一个模型平均值，以平均系数的最佳拟合模型。silly<-model.avg(my.d
浏览 0提问于2014-01-29得票数 6
正如标题中所指出的，我正试图生成一个预测光栅，描述使用的相对概率。为了创建一个可重复的示例，我使用了来自MaungaWhau包的maxlike数据。library(maxlike)data(MaungaWhau)
我们可以制作两个光栅层，一个光栅堆栈，以及一个SpatialPoints对象。我的第一个问题是如果确认。raster包vignette声明“Raster对象中的名称应该与模型所期望的</
浏览 4提问于2016-03-09得票数 1
回答已采纳
1回答
当拟合混合效果模型(或任何其他模型)时，使用新数据或新模型参数预测“反事实”拟合值通常是有用的。在我的情况下，我想用新的随机效果来预测拟合值。lme4包允许用户插入新参数，以替换merMod类的beta、theta和sigma插槽。这将允许一个人在反事实固定效应系数下预测拟合值。我如何才能对随机效果估计做同样的事情？我的第一个想法是直接修改merMod对象
浏览 0提问于2016-03-12得票数 3
1回答
我已经建立了一个线性回归模型，我得到了训练集和验证集的两个不同的分数，这意味着什么？这是我得到的训练数据得分: 0.8802644749581061，验证数据得分: 0.7124061878004602reg = LinearRegression().fit(X_train, y_train)这是从下面这行代码中得到
浏览 0提问于2019-10-22得票数 1
有一个R函数(通过nlme包) ，它包含一个具有公共模型的lm对象列表，它根据分组因子对数据进行分区。这有很多模型通常具有的有用方法，比如预测。glm是否有类似的功能？我确实看到了，但与lmlist不同的是，它似乎更像是一个花哨的容器，没有任何方法如预测。
也就是说，glm的代码是否相当于以下代码？
浏览 0提问于2016-05-19得票数 1
回答已采纳
2回答
线性回归不返回所有系数
、、、
我正在对所有的预测因子进行线性回归(我有384个预测因子)，但是从总结中只能得到373个系数。我想知道为什么R不返回所有系数，以及如何得到所有384个系数？
浏览 0提问于2018-04-22得票数 0
回答已采纳
1回答
模型的r-反测试
、、
我必须评估我的模型的预测能力。我想做什么？我回归的结果34565自由度残差标准误差: 271.6倍R-平方: 0.9552，校正R-平方: 0.9545 F
浏览 3提问于2019-10-14得票数 2
回答已采纳
1回答
R中的相关与回归
、、
我想在R中使用快速食品dataset fromopenintro`包执行以下任务。( a)为Sonic、Subway和Taco Bell中所有项目的卡路里、total_fat、糖和钙之间的关系创建一个相关矩阵，忽略na.omit()中缺少的值。建立一个回归模型，根据卡路里、钠和蛋白质来预测一家餐馆是McDonalds还是Subway。
( c)对饱和脂肪、纤维和糖的热量进行回归预测。根据标准化回归系数，确定最优<em
浏览 1提问于2022-02-25得票数 0
回答已采纳
1回答
根据文档，它展示了如何导入python中不标准的R包。幸运的是，他们做了我需要的例子，那就是lme4from rpy2.robjects import FloatVectorimport rpy2.rinterface as rinterface
base = import
浏览 3提问于2018-01-18得票数 0
1回答
考虑lme4包中的睡眠研究数据，如下所示。该研究包含18名受试者，他们在不同的日期重复测量反应(反应就是反应)。library("lme4") Reaction Days Subject2 258.7047 1fit1 = lmer(Reaction ~ Days + (1 | Subject), data = sleepstudy) 我们可以使用"ranef(fi
浏览 38提问于2021-11-19得票数 0
回答已采纳
我使用leaps库中的regsubsets来执行最佳子集选择。我需要将它生成的系数与我在模拟数据时指定的“真实”系数进行比较(通过比较，即它们之间的差值的平方，以及总和的平方根)，对于每个预测器。由于regsubsets生成了16个不同的模型，因此我使用一个循环自动完成此操作。它可以工作，但当我从具有x个预测器的</em
浏览 12提问于2021-02-13得票数 0
回答已采纳
1回答
利用rpy2调整随机效应
、、、、
我使用lmer模型()来模拟不同国家不同产品的价格弹性。在使用历史数据对模型进行培训之后，有时对其中一些产品的弹性是正的(根据定义，它应该是负的，或者仅仅是由于业务限制)。所以我需要手动调整一些系数，只有那些没有意义的系数。/country)
model = lmerTest.lmer(model_str, data = df)
在这个问题中，也解决了同样的问题，但在
浏览 14提问于2022-11-23得票数 0
1回答
从lme4模型中获取预测器的类
、、、、
在对lme4模型进行拟合之后，我想知道如何从terms(fit)[[3]]获得预测器的类。注意：所有东西都必须从模型中提取出来。library(lme4)
h <- read.csv('https://raw.githubusercontent.com/hkil/m/master/h.cs
浏览 3提问于2020-05-21得票数 1
回答已采纳
我在model.avg()包中使用MuMIn函数平均了三个GLMM模型。我现在想用这个平均模型来预测。但是，当我调用predict() R时，我似乎认为我想要predict.merMod()，因此会抛出一个警告消息，因为我请求了标准错误，这些错误不是由predict.merMod()提供的。相反，我想调用predict.averaging()，这是对类averaging模型的适当调用，它将为模型预测提供标准错误。如果我试图通过指定predict.averaging()或
浏览 8提问于2020-05-12得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
