使用broom::augment函数获得响应比例下的预测值
我正在使用R拟合glm模型，可以使用predict.glm(object=fm1, type="response")得到响应比例尺下的预测值，其中fm1是拟合的模型。我想知道如何使用broom包中的augment函数得到响应比例尺下的预测值。以下是我的最小工作示例。Dilution <-...
在使用caret的train()函数训练randomForest对象并使用公式进行预测时出现错误。
在一台64位Linux机器上，使用 R 3.2.0 版本、caret 6.0-41 版本和 randomForest 4.6-10 版本。当尝试对使用公式训练的randomForest对象应用caret包中的train()函数进行预测时，使用predict()方法会返回错误。但是，如果使用ra...
"predict.glm(, type="terms")" 实际上是做什么的？（涉及IT技术）
我对R中的predict.glm函数的工作方式感到困惑。根据帮助文档， “terms”选项返回一个矩阵，给出模型公式中每个项在线性预测量上的拟合值。 因此，如果我的模型形式为f(y) = X * beta，则命令 predict(model, X, type='terms') 预计会...
在R中的回归（逻辑回归）：寻找特定y值（结果）的x值（预测变量）
我已经拟合了一个逻辑回归模型，该模型从mtcars数据集中的mpg预测二进制结果vs。下面显示了绘图结果。如何确定任何特定vs值的mpg值？例如，我想找出当vs的概率为0.50时mpg的值是多少。感谢任何人提供的帮助！model <- glm(vs ~ mpg, data = mtcar...
glmer logit - 用概率刻度测量的交互作用效应（使用`predict`复制“effects”）
我正在使用lme4包运行glmer logit模型。我对各种二元和三元交互作用效应及其解释感兴趣。为简化问题，我只关心固定效应系数。 我设法编写了一段代码来计算并绘制这些效应的logit比例，但我在将它们转换为预测概率比例时遇到了困难。最终，我希望复制“effects”包的输出结果。 该示...
Python ARIMA模型中的外生变量如何在样本外进行预测
我正在尝试使用Python中的statsmodels ARIMA包预测时间序列，并包括一个外生变量，但无法确定在预测步骤中插入外生变量的正确方法。请参见此处的文档。 import numpy as np from scipy import stats import pandas as pd ...
使用`lm`的线性模型：如何获取预测值之和的预测方差
我正在对多个预测变量的线性模型进行预测值求和，就像下面的例子一样，并希望计算这个总和的组合方差、标准误差和可能的置信区间。lm.tree <- lm(Volume ~ poly(Girth,2), data = trees) 假设我有一组名为Girths的集合：newdat <- ...
如何处理面板数据回归中的NA值？
我正在尝试在包含NA的数据上预测拟合值，这是基于plm生成的模型。以下是一些示例代码： require(plm) test.data <- data.frame(id=c(1,1,2,2,3), time=c(1,2,1,2,1), y=c(1,3,5,10,8), x=c(1,...
如何从线性模型（lm）中预测x值
我有这样的数据集：x <- c(0, 40, 80, 120, 160, 200) y <- c(6.52, 5.10, 4.43, 3.99, 3.75, 3.60) 我使用lm()计算了一个线性模型：model <- lm(y ~ x) 我想知道如果我有新的y值，例如yn...
在R中决定glm逻辑回归模型的阈值
我有一些带有预测变量和二元目标的数据。例如：df <- data.frame(a=sort(sample(1:100,30)), b= sort(sample(1:100,30)), target=c(rep(0,11),rep(1,4),rep(0,...