search for: nlmecontrol

...llo all. I'm doing a simulation study where I will be making use of the 'nlme' package. I want to loosen up the convergence criteria so that I increase the likelihood of convergence (potentially at the cost of obtaining slightly less than ideal results). The parameters in the function nlmeControl() control the convergence criteria. These default values can be modified to make convergence critera more or less stringent. In an effort to get more converging results, I first increased the iterations for 'maxIter', 'pnlsMaxIter', 'msMaxIter', and 'niterEM' by a fa...

All, I'm trying to create an augPred plot in the nlme library, similar to the plot on p.43 of Pinheiro & Bates (Mixed Effects Models in S and S-Plus) for their Pixel data. My data structure is the same as the example but I still get the error msg below. > comp.adj.UKV <- groupedData(adj.UKV ~ Time | Patient_no/Lisinopril, data = comp.adj.UKV.frm, order.groups = F) >

...A fixed DESCRIPTORS ONLY.csv",sep=",",header=TRUE) #Load in my data set ge$study = as.factor(ge$study) #Force study to be a factor variable gelist = nlsList(RBA~SSBiLin(Molecular.Volume,nu,a,b,c,x0)|study,data=ge) #This is where I get the error contr = nlmeControl(pnlsTol=0.3) ge0.nlme=nlme(gelist,control=contr,verbose=F) ##########Errors####################### > gelist = nlsList(RBA~SSBiLin(Molecular.Volume,nu,a,b,c,x0)|study,data=ge) Error in qr(attr(rhs, "gradient")) : NA/NaN/Inf in foreign function call (arg 1) In addition: Warning messag...

...light), data=light, fixed = a.max + a.qe + lcp ~ trt, random = a.max + a.qe + lcp ~ 1 | bench/line, start = list(fixed=c(17.4305, 0.182444, 0.00928341, -0.00057221, 44.8384, 8.67678)), method="ML", verbose=T, control = nlmeControl(maxIter=250, msMaxIter=200, pnlsMaxIter=20, gradHess=TRUE, returnObject=TRUE, niterEM=100)) What I'd like to do is to compare reduced models with something like random = list(a.max ~ 1 | bench, a.qe + lcp ~ 1 | bench/line) i.e., to drop the nested term on...

...~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 Veterinary Science C01 JL Shute Building Camden University of Sydney

....a2<https://maps.google.com/?q=1,A.a2&entry=gmail&source=g>=0,A.a3=0,B=1,B.a2=0,B.a3=0,C=0,C.a2=0,C.a3=0,D=1,D.a2<https://maps.google.com/?q=1,D.a2&entry=gmail&source=g>=0,D.a3=0) try({withTimeout(nlme(model=y~fpl.B.range(x,exp(1),A,B,C,D), control=nlmeControl(maxIter=50,pnlsMaxIter=7,msMaxIter=50,niterEM=25), data=dat, na.action=na.omit, fixed=myFormula,random=list(f.block=pdSymm(A+B+C+D~1)), start=INIT),timeout=10)}) [[alternative HTML version deleted]]

...t; > Values > > getData.lme(nlme) Extract lme Object Data > > nlme(nlme) Nonlinear Mixed-Effects > > Models > > nlme.nlsList(nlme) NLME fit from nlsList > > Object > > nlmeControl(nlme) Control Values for nlme Fit > > nlmeObject(nlme) Fitted nlme Object > > nlmeStruct(nlme) Nonlinear Mixed-Effects > > Structure > > predict.nlme(nlme) P...

Hi, I was very much hoping someone could help me with the following. I am trying to convert some SAS NLMIXED code to NLME in R (v.2.1), but I get an error message. Does anyone have any suggestions? I think my error is with the random effect "u" which seems to be parametrized differently in the SAS code. In case it's helpful, what I am essentially trying to do is estimate parameters

...2=0,C.a3=0,D=1,D.a2=0,D.a3=0) for (i in 1:100) { print(paste("Iteration ",i,"...this will stall soon")) set.seed(i) dat$y <- dat$y+rnorm(nrow(dat), mean = 0, sd = 0.1) try({withTimeout(nlme(model=y~fpl.B.range(x,exp(1),A,B,C,D), control=nlmeControl(maxIter=50,pnlsMaxIter=7,msMaxIter=50,niterEM=25), data=dat, na.action=na.omit, fixed=myFormula,random=list(f.block=pdSymm(A+B+C+D~1)), start=INIT),timeout=10)}) } [[alternative HTML version deleted]]