Displaying 8 results from an estimated 8 matches for "pnlstol".
2003 Sep 16
2
gnls( ) question
...unanswered question:
What exactly does "Warning message: Step halving factor
reduced below minimum in NLS step in: gnls(model = y ~ 5 + ...)"
mean? I have tried to address this by specifying "control =
list(maxIter = 1000, pnlsMaxIter = 200, msMaxIter = 1000,
tolerance = 1e-06, pnlsTol = 1e-04, msTol = 1e-07, minScale =
1e-10, returnObject = TRUE)" in my model calls, but this
does not entirely eliminate the problem (I am running gnls( )
24 separate times on separate data sets).
Much thanks in advance,
david paul
#Constructing Indicator Variables
indicator <- past...
2004 Jul 23
1
nlme parameters in nlmeControl
...ean to increase the default value? Will increasing this value lead to more convergence because the criteria is more liberal?
Is the tolerance criteria met when the change in the log likelihood is reduced below the
specified value, or is it a value other than the log likelihood? What about for the
pnlsTol? That is, what is the quantity whose change must be less than pnlsTol to continue?
Is there anything wrong with reducing 'minAbsParApVar' to zero (so that a parameter can be fixed across all individuals and thus not have a variance)?
Is anything, other than the help file, written about...
2012 Feb 05
1
Covariate model in nlme
...eterized in terms of CL, V1, V2 and Q.
basemodel<-nlme(Conc ~TwoCompModel(CL,Q,V1,V2,Time,ID),
data = data2, fixed=list(CL+Q+V1+V2~1),
random = pdDiag(CL+V1+V2~1),
start=c(CL=log(20),Q=log(252),V1=log(24.9),V2=log(120)),
control=list(returnObject=TRUE,msVerbose=TRUE,
msMaxIter=20,pnlsMaxIter=20,pnlsTol=1),
verbose=TRUE).
Now, I want to include covariates (age and weight) in the base model. The
relationships between CL and covariate is linear and given as the following
equation.
CL=TVCL+a*wt-b*age, where TVCL is mean clearance, a and b are coefficients
of weight and age, respectively.
How do I...
2007 Apr 20
1
nlme trouble
...RIPTORS
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 message:
NaNs...
2002 May 10
0
PNLS step error in nlme
Hi R-help;
In using nlme, I often encounter the following error message:
Error in nlme.formula(model = response ~ ......, data = ...
:
Step halving factor reduced below minimum in PNLS step
Resetting the pnlsTol default value does not solve the problem. Any
suggestions or workarouds for solving the problem?
TIA,
Richard
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", &quo...
2005 Nov 02
1
NLME
Dear All,
Using:-
R : Copyright 2005, The R Foundation for Statistical Computing
Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0
and
Jose Pinheiro, Douglas Bates, Saikat DebRoy and Deepayan Sarkar (2005). nlme: Linear and
nonlinear mixed effects models. R package version 3.1-65.
on a WINDOWS 2000 machine
I am trying to run the scripts from the Mixed Models book and am running into some
2005 Jul 07
2
Problems with nlme: Quinidine example
This concerns the "Clinical Study of Quinidine" example on page 380
of the book "Mixed-Effects Models in S and S-PLUS" by Pinheiro and Bates (2000).
I have tried to reproduce the example, but get an error:
> library(nlme)
> fm1Quin.nlme <- nlme(conc ~ quinModel(Subject, time, conc, dose, interval, lV,
lKa, lCl),
+ data=Quinidine,
+
2005 Dec 12
2
convergence error (lme) which depends on the version of nlme (?)
Dear list members,
the following hlm was constructed:
hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I)
the grouped data object is located at and can be downloaded:
www.anicca-vijja.de/lg/hlm_example.Rdata
The following works:
library(nlme)
summary( fitlme <- lme(hlm) )
with output:
...
AIC BIC logLik
425.3768 465.6087 -197.6884
Random effects: