Displaying 1 result from an estimated 1 matches for "logfv".
Did you mean:
logfp
2002 Aug 29
0
fitting nlme model inside a function
I'm not an experienced user of the 'nlme' package but I ran into a curious
problem recently that I couldn't find on the mailing list anywhere.
I have simple function which calls 'groupedData' and 'nlme':
g <- function(dataset) {
unig <- groupedData(logfv ~ Vs+v1+v2+v3+PHA| EQ, data=dataset)
modelfv <- nlme(logfv~a + c*log(Vs/a1)+ b1*( v1 +
v2*(Vs-300)^2/(300-120)^2 )*log(PHA/0.1)
+ b2*( v3 + v2 + (-1*v2)*(Vs-300)^2 /(300-120)^2)*log(PHA/0.1),
data=unig, fixed=a1+c+b1+b2~1, random=a~1|EQ ,...