Displaying 20 results from an estimated 3000 matches similar to: "step halving factor reduced below minimum"
2003 Dec 17
1
step halving factor reduced below minimum
Dear Splus and R users:
I have a problem in fitting a NLME model, the error message is:
"step halving factor reduced below minimum in PNLS step"
What does it mean and how to fix the problem, could anyone help me
about it?
Any suggestion/help would be greatly appreciated.
Mei
2006 Jul 18
2
Using corStruct in nlme
I am having trouble fitting correlation structures within nlme. I would like to
fit corCAR1, corGaus and corExp correlation structures to my data. I either
get the error "step halving reduced below minimum in pnls step" or
alternatively R crashes.
My dataset is similar to the CO2 example in the nlme package. The one major
difference is that in my case the 'conc' steps are
2009 Jan 04
1
POSIXct and chron issues with tz
Dear All-
I am trying to merge two data files - they have different date formats
and different times zones. I need to match up the date/time of the
datasets and then invoke a conditional statement, such as: if dataC$mph
is >= 12 then keep dataM$co23 for the corresponding time/date stamp.
snippets of data files:
*dataC.txt*
LST in mph Deg DegF DegF2 % volts Deg
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
2005 Jul 17
1
how to solve the step halving factor problems in gnls and nls
Hi R-users,
Could you give me some advice in
solving the problem of such error message from gnls and nls?
## begin error message
"Problem in gnls(y1 ~ glogit4(b, c, m, t, x), data.frame(x..: Step halving
factor reduced below minimum in NLS step "
##and
"Problem in nls(y ~ 1/(1 + exp((xmid - x)/scal)), data = x..: step factor
reduced below minimum "?
Thank you in
2003 Aug 14
1
gnls - Step halving....
Hi all,
I'm working with a dataset from 10 treatments, each
treatment with 30 subjects, each subject measured 5
times. The plot of the dataset suggests that a
3-parameter logistic could be a reasonable function to
describe the data. When I try to fit the model using
gnls I got the message 'Step halving factor reduced
below minimum in NLS step'. I´m using as the initial
values of the
2008 Aug 08
2
aggregate
Dear All-
I have a dataset that is comprised of the following:
doy yr mon day hr hgt1 hgt2 hgt3 co21 co22 co23 sig1 sig2 sig3 dif flag
244.02083 2005 09 01 00 2.6 9.5 17.8 375.665 373.737 373.227 3.698 1.107
0.963 -0.509 PRE
244.0625 2005 09 01 01 2.6 9.5 17.8 393.66 384.773 379.466 15.336 11.033
5.76 -5.307 PRE
244.10417 2005 09 01 02 2.6 9.5 17.8 411.162 397.866 387.755 6.835 5.61
6.728
2008 Jan 31
1
Error handling in nlme call
In some trial simulation work I need to create batch files that will
repeatedly generate pseudoreplicate datasets and then create non-
linear mixed effects models using nlme. Inevitably these models
sometimes fail to converge but I need the batch file to simply move on
to another simulation rather than abort. I am using the try() function
as in
model<-try((nlme(...)))
which handles
2007 Jan 22
0
Branch 'interpreter' - 3 commits - libswfdec/swfdec_bits.c libswfdec/swfdec_swf_decoder.c test/Makefile.am test/swfedit.c test/swfedit_file.c test/swfedit_file.h test/swfedit_tag.c test/swfedit_tag.h test/swfedit_token.c test/swfedit_token.h
libswfdec/swfdec_bits.c | 27 +-
libswfdec/swfdec_swf_decoder.c | 46 ----
test/Makefile.am | 19 +
test/swfedit.c | 81 +++++++
test/swfedit_file.c | 209 ++++++++++++++++++
test/swfedit_file.h | 59 +++++
test/swfedit_tag.c | 61 +++++
test/swfedit_tag.h | 56 ++++
test/swfedit_token.c |
2008 Oct 15
2
Network meta-analysis, varConstPower in nlme
Dear Thomas Lumley, and R-help list members,
I have read your article "Network meta-analysis for indirect treatment
comparisons" (Statist Med, 2002) with great interest. I found it very
helpful that you included the R code to replicate your analysis;
however, I have had a problem replicating your example and wondered if
you are able to give me a hint. When I use the code from the
2006 Feb 17
0
trouble with extraction/interpretation of variance struct ure para meters from a model built using gnls and varConstPower
Works perfectly. Thank you.
-Hugh Rand
-----Original Message-----
From: Spencer Graves [mailto:spencer.graves at pdf.com]
Sent: Sunday, January 15, 2006 6:41 PM
To: Rand, Hugh
Cc: 'r-help at lists.R-project.org'
Subject: Re: [R] trouble with extraction/interpretation of variance
structure para meters from a model built using gnls and varConstPower
How about this:
>
2004 Oct 18
3
manual recreation of varConstPower using new fixed effects variables in nlme
Hello, I am trying to design new variance structures
by using fixed effects variables in combination with
the VarPower function. That is, I would like to
create and evaluate my own variance function in the
data frame and then incorporate it into the model
using varPower, with value=.5.
As a start, I am trying to recreate the function of
VarConstPower by introducing two new variables in the
2009 Jun 22
1
The gradient of a multivariate normal density with respect to its parameters
Does anybody know of a function that implements the derivative (gradient) of
the multivariate normal density with respect to the *parameters*?
It?s easy enough to implement myself, but I?d like to avoid reinventing the
wheel (with some bugs) if possible. Here?s a simple example of the result
I?d like, using numerical differentiation:
library(mvtnorm)
library(numDeriv)
f=function(pars, xx, yy)
2010 Oct 15
2
How to extract parameter estimates of variance function from lme fit
Dear R-Users,
I have a question concerning extraction of parameter estimates of
variance function from lme fit.
To fit my simulated data, we use varConstPower ( constant plus power
variance function).
fm<-lme(UPDRS~time,data=data.simula,random=~time,method="ML",weights=varConstPower(fixed=list(power=1)))
I extract the results of this function by using the following codes:
2008 Feb 24
0
problem with ML estimation
dear list,
as a part my problem. I have to estimate some parameters using ML
estimation. The form of the likelihood function
is not straight forward and I had to use a for loop to define the function.
I used "optim" to maximise the result but
was not sure of the programme.
To validate my results, I tried to write a function to obtain the MLE of a
bivariate normal in the same manner.
On
2012 Jan 03
1
nls and rbinom function: step factor 0.000488281 reduced below 'minFactor' of 0.000976562
I am trying to learn nls using a simple simulation. I assumed that the binomial prob varies linearly as 0.2 + 0.3*x in x {0,1},
and the objective is to recover the known parameters a=0.2, b=0.3
..data frame d has 1000 rows...
d$x<-runif(0,1)
d$y<-rbinom(1000,1,0.2+0.3*d$x)
table(d$y,cut(d$x,breaks=5));
(-0.000585,0.199] (0.199,0.399] (0.399,0.599] (0.599,0.799]
2011 Mar 28
1
error in nls, step factor reduced below minFactor
Hello,
I've seen various threads on people reporting:
step factor 0.000488281 reduced below `minFactor' of 0.000976563
While I know how to set the minFactor, what I'd like to have happen is for nls to return to me, the last or closest fitted parameters before it errors out. In other words, so I don't get convergence, I'd still like to acquire the values of the parameters
2004 Feb 16
0
error in nls, step factor reduced below minFactor
Hello,
I am trying to estimate 4 parameters of a non-linear
model using nls.
My model function is a Fourier integral and is very
expensive to
calculate. I get the following error:
> theta0 <- c(0.045, 1.02*10^(-4), 0.00169,
5.67*10^(-4))
> res <- nls(log(y) ~ log(model(theta,r,t)),
data=dataModel,
+ start=list(theta=theta0), trace=TRUE,
+ control=nls.control(tol=1e-2))
2006 Jan 09
1
trouble with extraction/interpretation of variance structure para meters from a model built using gnls and varConstPower
I have been using gnls with the weights argument (and varConstPower) to
specify a variance structure for curve fits. In attempting to extract the
parameters for the variance model I am seeing results I don't understand.
When I simply display the model (or use "summary" on the model), I get what
seem like reasonable values for both "power" and "const". When I
2006 Jun 15
2
Standard Deviation Distribution
I'm having trouble with the standard deviation distribution
as shown on http://mathworld.wolfram.com/StandardDeviationDistribution.html .
(Eric Weisstein references Kenney and Keeping 1951, which I can't check.)
I believe the graphs they show, but when I code the function in R, according to the listed formula,
I get very different graphs.
Would someone please point out my error or tell