我基于拟拟数族进行曲线拟合,现在我想估计曲线拟合的性能以及预测误差。我该怎么做? 我把我的代码放在这里和情节
h=23
daneS = cumsum(unitsSold[,22])
dopasowanie = function(daneS, h)
{
n = length(daneS)
standaryzacja = 1/daneS[h]
uczenie = 1:h
testowanie = (h+1):n
daneS = daneS*standaryzacja
model = glm(daneS[uczenie]~uczenie, family=quasibinomial)
prognoza = predict(model, newdata=data.frame(uczenie=testowanie))
return(list(parametry=model$coefficients, prognoza=prognoza/standaryzacja))
}
szeregNieskum=units_sold[,22]
szeregSkumulowane=cumsum(szeregNieskum)
ab=dopasowanie(szeregSkumulowane, 22)$parametry
dopasowane = (1/(1+exp(-ab[1]-ab[2]*1:30)))*szeregSkumulowane[22]
### Plot based on amount of sales
plot(szeregNieskum, type="l", xlim=c(1,30),ylim=c(-100,400))
lines(c(diff(dopasowane)), col="red")
### Plot based on sales dynamics
plot(diff(szeregNieskum), type="l", xlim=c(1,30),ylim=c(-100,400))
lines(c(diff(diff((dopasowane)))), col="red")
这是情节http://imageshack.com/a/img834/2521/g6sq.png(抱歉,我无法在此处显示图片)
我很感激任何建议和程序