Displaying 3 results from an estimated 3 matches for "glns".
Did you mean:
glms
2012 Mar 21
1
nlme error on dimensions in multiplication
...mpleNmleError.RData")
load(file=con) #ds1, fm1, ctrl
close(con)
ds1 <- subset(ds1, select=c("ColNoM","minT","resp"))
head(ds1)
plot( resp ~ minT, col=rainbow(length(unique(ColNoM)))[ColNoM], ds1)
lines( fitted.values(fm1) ~ minT, ds1 ) # starting values from glns fit
ds2 <- ds1
# ds2 <- unique(ds1) # does not resolve the error
tmp.fit1 <- nlme(
resp ~ exp(beta0l) + beta1 * exp(-exp(beta2l)*minT)
,fixed=list(beta0l+beta1+beta2l~1)
,random= beta0l ~ 1 | ColNoM
, weights =varPower(fixed=0.4, form=~resp)
, start=list(fixed = coef(fm1))
#, meth...
2002 Sep 10
9
R 1.6.0 beta
...llcheck can take a long time, especially in the boot and nlme
packages).
The current version in known to pass the checks for our Linux/x86
systems but has two known problems on Solaris/gcc: An optimizer bug in
gcc 3.0.3 causing "cluster" to fail in mona.f and a convergence issue
with the glns example in nlme.
--
O__ ---- Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk) FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-...
2003 Jul 24
0
nls.control in gnls
...I'm using gnls instead of nls. I want to
get rid of the start= option, but I keep getting the following error
message:
Error in nls(formula = y ~ spot.shape.fct(xcord, ycord, background,
spotintensity, :
number of iterations exceeded maximum of 50
>From the help page and the glns code I can see that the initial start
values are improved by estimating the parameters in a model without the
correlation structures that are possible with gnls. This is quite
reasonable and will quite likely improve the convergence succes and time
once a correlation is introduced.
Obviously t...