@font-face {
font-family: "ＭＳ 明朝";
}@font-face {
font-family: "Verdana";
}@font-face {
font-family: "Verdana";
}@font-face {
font-family: "Monaco";
}p.MsoNormal, li.MsoNormal, div.MsoNormal { margin: 0cm 0cm 0.0001pt;
font-size: 10pt; font-family: Verdana; }.MsoChpDefault { font-family: Cambria;
}div.WordSection1 { page: WordSection1; }
Greetings,
I am trying
to add a cubic spline to an NLME biological growth model. After looking in the
R-help archive, I’ve tried the functions ns() and bs() from the splines-package
and rcs() from the Design-package:
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
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
