Displaying 1 result from an estimated 1 matches for "ans51g".
Did you mean:
ans1
2004 Jan 14
2
Generalized least squares using "gnls" function
...rged parameter estimates. I also
get the same answers with "gnls" without specifying the "weights"
argument.
However, when I attempt to use the "gnls" function and try to estimate
the variance function, as a power function, I get the following error
message:
> ans51g <- gnls(log(b51) ~ p0 + p1/(1 + exp(-(log(dose)-p2)/p3))^p4,
start=list(p0=3,p1=1,p2=4,p3=2,p4=1.5),control=gnlsControl(tol=1.e-
07),weights=varPower())
Error in eval(expr, envir, enclos) : Object "." not found
>
What am I doing wrong here and how can I do a GLS analysis with a...