如何测量R中的曲线拟合?

时间:2014-04-14 20:53:13

标签: r time-series curve-fitting forecasting

我基于拟拟数族进行曲线拟合,现在我想估计曲线拟合的性能以及预测误差。我该怎么做? 我把我的代码放在这里和情节

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(抱歉,我无法在此处显示图片)

我很感激任何建议和程序

0 个答案:

没有答案