本笔记中原始数据及代码均来源于李东风先生的R语言教程,在此对李东风先生的无私分享表示感谢。
预测区间
设X取x0,Y的预测值为y0^=a^+b^x0。 置信水平为1-α的预测区间为
其含义是P(y0∈预测区间)=1-α。 λ是标准正态分布双侧α分位数qnorm(1 - alpha/2)。 用predict(lmres)得到yi^,i=1,2,...,n的值, lmres表示回归结果变量。
用predict(lmres, interval="prediction") 同时得到预测的置信区间, 需要的话加入level=选项设定置信度。 回归的置信区间有两种, 上述的区间是“预测区间”(CLI), 使得P(y0∈预测区间)=1-α; 另一种区间称为均值的置信区间(CLM), 使得P(Ey0∈均值置信区间)=1-α。
控制
如果需要把Y的值控制在[yi,yu]范围内,问如何控制X的范围,