Displaying 1 result from an estimated 1 matches for "gradnm".
Did you mean:
gradim
2011 Mar 18
0
adding a cubic spline to a model
...kage:
?
M7.W<-nlme(model= RW ~ (R*(RPrevYear/K)^m)*ns(Year,4)*exp(m*((1-RPrevYear^2/K^2)/2)), data=data, fixed=R+K+m~1, random= R+K~1|IDPlot, start= c(R=3, K=90, m=0.1), correlation=corARMA(p=2, q=0, form= ~Year))
?
All of them
give the same error-message:
?
Error in
X[, fmap[[nm]]] <- gradnm :
? number of items to replace is not a multiple
of replacement length
?
Am I using
the ns() function wrongly or should I use another function?
?
Thanks for
any hints,
?
Ellen,
Msc
student, KULeuven, Belgium