我无法轻易找到答案的快速编码问题-搜索了R帮助并在线上找到了相当多的地方。
我正在跟踪a post regarding prediction vs. confidence intervals。
有一行代码如下所示:
pd <- transform(pd, Fitted = ilink(fit), Upper = ilink(fit + (2 * se.fit)),
Lower = ilink(fit - (2 * se.fit)))
我对fit
和参数Fitted
所做的事情有些迷惑。我认为它们已用于预测模型输出,但是我看不到帖子中显示的代码中的任何地方定义了fit
。
答案 0 :(得分:4)
Fitted
是要在pd
数据框中创建的新变量的名称(Upper
和Lower
也是如此)fit
是pd
数据框中已经存在的变量(se.fit
也是如此)不用看博客文章,我很确定pd
是从predict(., se.fit=TRUE)
生成的:?predict.lm
是
如果“ se.fit”为“ TRUE”,则包含以下组件的列表为 返回:
fit:如上所述的向量或矩阵
se.fit:预测均值的标准误
residual.scale:残留标准偏差
df:残差的自由度
如果使用的是GLM,则应查看?predict.glm()
,它会给出相似(但不相同)的答案。如它所说
type:所需预测的类型。默认的 [即type =“ link”] 线性预测因子;另一种是“响应” 响应变量的大小。
如果您将type="link"
用于具有logit链接的二项式模型,则会得到对数奇数(logit)范围的响应,但不限于(0,1)。