Displaying 2 results from an estimated 2 matches for "mlkyld".
Did you mean:
medyld
2004 Aug 19
0
nlme R vs S plus
...ations. Is this because the nlme function
is slightly different and if so, is there a way to solve this problem?
The code I am using is below:
W3<-deriv(~A*ti^exp(logB)*exp(-exp(logC)*ti),c("A","logB","logC"),
function(ti,A,logB,logC){})
Lact.nlme<-{nlme(model=MLKYLD~W3(DIM,A,logB,logC),
fixed=A+logB+logC~1,
random=A+logB+logC~1|ID,
data=LacData.x0, start=c(13.41143,log(0.152792),log(0.002494)), control=nlmeControl(msMaxIter=200),
verbose=T)
}
Any help or advice would be greatly appreciated,
Evelyn
--
Evelyn Hall
PhD Student
Faculty of Vet...
2005 Apr 11
1
extracting correlations from nlme
Hi,
I would like to know how (if) I can extract some of the information from
the summary of my nlme.
at present, I get a summary looking something like this:
> summary(fit.nlme)
Nonlinear mixed-effects model fit by maximum likelihood
Model: MLKYLD ~ W4(DIM, logA, B, C)
Data: ADHIS.x0
AIC BIC logLik
265314 265401.6 -132647
Random effects:
Formula: list(logA ~ 1 , B ~ 1 , C ~ 1 )
Level: ID
Structure: General positive-definite
StdDev Corr
logA 0.480476962 logA B
B 0.140839992 -0.855
C 0....