我正在尝试使用R
函数预测predict()
中的值,方法是将变量传递到模型中。
我收到以下错误:
Error in eval(predvars, data, env) :
numeric 'envir' arg not of length one
这是我的data frame
,名称为df:
df <- read.table(text = '
Quarter Coupon Total
1 "Dec 06" 25027.072 132450574
2 "Dec 07" 76386.820 194154767
3 "Dec 08" 79622.147 221571135
4 "Dec 09" 74114.416 205880072
5 "Dec 10" 70993.058 188666980
6 "Jun 06" 12048.162 139137919
7 "Jun 07" 46889.369 165276325
8 "Jun 08" 84732.537 207074374
9 "Jun 09" 83240.084 221945162
10 "Jun 10" 81970.143 236954249
11 "Mar 06" 3451.248 116811392
12 "Mar 07" 34201.197 155190418
13 "Mar 08" 73232.900 212492488
14 "Mar 09" 70644.948 203663201
15 "Mar 10" 72314.945 203427892
16 "Mar 11" 88708.663 214061240
17 "Sep 06" 15027.252 121285335
18 "Sep 07" 60228.793 195428991
19 "Sep 08" 85507.062 257651399
20 "Sep 09" 77763.365 215048147
21 "Sep 10" 62259.691 168862119', header=TRUE)
str(df)
'data.frame': 21 obs. of 3 variables:
$ Quarter : Factor w/ 24 levels "Dec 06","Dec 07",..: 1 2 3 4 5 7 8 9 10 11 ...
$ Coupon: num 25027 76387 79622 74114 70993 ...
$ Total: num 132450574 194154767 221571135 205880072 188666980 ...
代码:
model <- lm(df$Total ~ df$Coupon)
> model
Call:
lm(formula = df$Total ~ df$Coupon)
Coefficients:
(Intercept) df$Coupon
107286259 1349
现在,当我运行predict
时,我收到了上面显示的错误。
> predict(model, df$Total, interval="confidence")
Error in eval(predvars, data, env) :
numeric 'envir' arg not of length one
知道我哪里错了吗?
由于
答案 0 :(得分:40)
这里有几个问题:
newdata
的{{1}}参数需要预测变量变量。因此,您应该为predict()
传递值,而不是Coupon
,这是模型中的响应变量。
预测变量需要作为数据框中的命名列传入,以便
Total
知道它所代表的数字代表什么。 (当您考虑更复杂的模型,具有多个预测变量时,对此的需求变得清晰。)
要使其正常工作,您的原始调用应通过predict()
参数传递df
,而不是直接在公式中使用它。 (这样,data
中列的名称将能够匹配公式的RHS上的名称。
随着这些变化的合并,这将有效:
newdata