我有4点火车,我用" nls"适合列车组,然后预测包含2个点的测试集的响应。然而,"预测"命令总是返回列车集的值。
The code is attached below:
##Following is train set
HD = c(714,715,716.6,717.6)
p_l = c(0.5,0.1,0.05826374, 0.005982334)
##Fitting with nls
raw_data = data.frame(HD,p_l)
exp_fit = nls(p_l~exp(a+b*HD),data = raw_data,trace = T,start = list(a = 0,b
= 0))
##Following is test set
HD_test = c(718.2,719.17)
p_l_test = predict(exp_fit,newdata = HD_test)
答案 0 :(得分:1)
您需要将数据框或命名列表作为newdata
传递。
使用与原始预测变量相同的名称(HD
)。否则,newdata
被视为缺失,并返回原始训练数据中的拟合值。
来自nls docs:
<强> newdata 强>
一个命名列表或数据框,用于查找要预测的变量。如果缺少newdata,则返回原始数据点的拟合值。
HD_test = data.frame(HD = c(718.2,719.17)) # wrap in data frame, name "HD"
p_l_test = predict(exp_fit, newdata = HD_test)
p_l_test
[1] 0.0009320202 0.0002184518