Displaying 1 result from an estimated 1 matches for "lddata".
Did you mean:
olddata
2009 Jun 03
0
fitting polynomial, for integration.
...I plot this the generated line is not
representative of the data. The polynomial intercept is correct, however
the gradient is well below the actual dataset, and does not follow the
same contour.
A code snippet and output are provided below. Any help would be much
appreciated.
Thanks,
Blaise.
lddata.mdl <- lm(lddata[,2]~poly(lddata[,1],9,raw=TRUE))
lddata.ply <- polynom(lddata.mdl$coefficients);
lddata.int <- integral(lddata.ply);
plot(lddata[,1],lddata[,2],pch=".")
lines(predict(lddata.ply,timeseq),col=4,lwd=2,lty=1)
lines(predict(lddata.mdl,newdata=list(timeseq=x)),col=4,...