Displaying 20 results from an estimated 800 matches similar to: "using nls to fit a four parameter logistic model"
2011 Nov 17
3
Obtaining a derivative of nls() SSlogis function
Hello, I am wondering if someone can help me. I have the following function
that I derived using nls() SSlogis. I would like to find its derivative. I
thought I had done this using deriv(), but for some reason this isn't
working out for me.
Here is the function:
asym <- 84.951
xmid <- 66.90742
scal <- -6.3
x.seq <- seq(1, 153,, 153)
nls.fn <- asym/((1+exp((xmid-x.seq)/scal)))
2008 Jan 04
3
nls (with SSlogis model and upper limit) never returns (PR#10544)
Full_Name: Hendrik Weisser
Version: 2.6.1
OS: Linux
Submission from: (NULL) (139.19.102.218)
The following computation never finishes and locks R up:
> values <- list(x=10:30, y=c(23.85, 28.805, 28.195, 26.23, 25.005, 20.475,
17.33, 14.97, 11.765, 8.857, 5.3725, 5.16, 4.2105, 2.929, 2.174, 1.25, 1.0255,
0.612, 0.556, 0.4025, 0.173))
> y.max <- max(values$y)
> model <- nls(y ~
2001 Oct 07
1
Bug in Deriv? (PR#1119)
deriv seems to have problems with a minus-sign before a bracket.
Below are four examples of the same function, the top one
is wrong, all others are correct (hopefully).
Rest of expression not shown, it is the same for all versions.
_
platform i386-pc-mingw32
arch x86
os Win32
system x86, Win32
status
major 1
minor 3.0
year 2001
month 06
day 22
language R
2011 Aug 09
1
nls, how to determine function?
Hi R help,
I am trying to determine how nls() generates a function based on the
self-starting SSlogis and what the formula for the function would be.
I've scoured the help site, and other literature to try and figure
this out but I still am unsure if I am correct in what I am coming up
with.
**************************************************************************
dat <-
2000 Oct 14
2
Access to calculations in nls
Hi,
I would like to be able to access the calculated results from the nls package.
Using the example in R, fm3DNase1 we can reurn certain parts of the
calculations:
> coef(fm3DNase1)
Asym xmid scal
2.345179 1.483089 1.041454
> resid(fm3DNase1)
[1] -0.0136806237 -0.0126806237 0.0089488569 0.0119488569 -0.0025803222
[6] 0.0064196778 0.0026723396 -0.0003276604
2006 May 17
1
nlme model specification
Hi folks,
I am tearing my hair out on this one.
I am using an example from Pinheiro and Bates.
### this works
data(Orange)
mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange )
### This works
mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange,
fixed = Asymp + xmid + scal ~ 1,
start =
2006 Jan 23
1
nlme in R v.2.2.1 and S-Plus v. 7.0
Dear R-Users,
I am comparing the nlme package in S-Plus (v. 7.0) and R (v. 2.2.1, nlme
package version 3.1-68.1; the lattice, Matrix, and lme4 have also just
been updated today, Jan. 23, 2006) on a PC (2.40 GHz Pentium 4 processor
and 1 GHz RAM) operating on Windows XP. I am using a real data set with
1,191 units with at most 4 repeated measures per unit (data are
incomplete, unbalanced). I
2005 May 10
2
predict nlme syntax
Dear all
Please help me with correct syntax of predict.nlme.
I would like to predict from nlme object for new data.
I used predict(fit.nlme6, data=newdata) but I have always got
fitted values, no matter how I changed newdata.
I have
> summary(fit.nlme6)
Nonlinear mixed-effects model fit by maximum likelihood
Model: konverze ~ SSfpl(tepl, A, B, xmid, scal)
Data: limity.gr
AIC
2004 Feb 20
1
nlme and multiple comparisons
This is only partly a question about R, as I am not quite sure about the
underlying statistical theory either.
I have fitted a non-linear mixed-effects model with nlme. In the fixed
part of the model I have a factor with three levels as explanatory
variable. I would like to use Tukey HSD or a similar test to test for
differences between these three levels.
I have two grouping factors:
2006 Sep 11
4
syntax of nlme
Hello,
How do I specify the formula and random effects without a startup object
? I thought it would be a mixture of nls and lme.
after trying very hard, I ask for help on using nlme.
Can someone hint me to some examples?
I constructed a try using the example from nls:
#variables are density, conc and Run
#all works fine with nls
DNase1 <- subset(DNase, Run == 1 )
fm2DNase1 <- nls(
2004 Oct 01
4
gnls or nlme : how to obtain confidence intervals of fitted values
Hi
I use gnls to fit non linear models of the form y = alpha * x**beta
(alpha and beta being linear functions of a 2nd regressor z i.e.
alpha=a1+a2*z and beta=b1+b2*z) with variance function
varPower(fitted(.)) which sounds correct for the data set I use.
My purpose is to use the fitted models for predictions with other sets
of regressors x, z than those used in fitting. I therefore need to
2009 Nov 09
1
Parameter info from nls object
Hi!
When checking validity of a model for a large number
of experimental data I thought it to be interesting
to check the information provided by
the summary method programmatically.
Still I could not find out which method to
use to get to those data.
Example (not my real world data, but to show the point):
[BEGIN]
> DNase1 <- subset(DNase, Run == 1)
> fm1DNase1 <- nls(density ~
2001 Aug 08
1
NLME augPred error
Could someone explain the meaming of this error message from augPred:
> augPred(area3.pen.nlme, primary=~day)
Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop =
FALSE], :
Levels 1,2,3 not allowed for block
>
predict.nlme(area3.pen.nlme) does not produce an error.
area3.pen.nlme was created with:
> area3.pen.nlme <- nlme(area ~ SSlogis(day, Asym, xmid, scal),
2008 Apr 14
3
Logistic regression
Dear all,
I am trying to fit a non linear regression model to time series data.
If I do this:
reg.logis = nls(myVar~SSlogis(myTime,Asym,xmid,scal))
I get this error message (translated to English from French):
Erreur in nls(y ~ 1/(1 + exp((xmid - x)/scal)), data = xy, start =
list(xmid = aux[1], :
le pas 0.000488281 became inferior to 'minFactor' of 0.000976562
I then tried to set
2009 Oct 02
1
nls not accepting control parameter?
Hi
I want to change a control parameter for an nls () as I am getting an error
message "step factor 0.000488281 reduced below 'minFactor' of 0.000976562".
Despite all tries, it seems that the control parameter of the nls, does not
seem to get handed down to the function itself, or the error message is
using a different one.
Below system info and an example highlighting the
2009 May 04
1
how to change nlme() contrast parametrization?
How to set the nlme() function to return the answer without the intercept parametrization?
#=========================================================================================
library(nlme)
Soybean[1:3, ]
(fm1Soy.lis <- nlsList(weight ~ SSlogis(Time, Asym, xmid, scal),
data = Soybean))
(fm1Soy.nlme <- nlme(fm1Soy.lis))
fm2Soy.nlme <- update(fm1Soy.nlme,
2008 Sep 27
1
seg.fault from nlme::gnls() {was "[R-sig-ME] GNLS Crash"}
>>>>> "VW" == Viechtbauer Wolfgang (STAT) <Wolfgang.Viechtbauer at STAT.unimaas.nl>
>>>>> on Fri, 26 Sep 2008 18:00:19 +0200 writes:
VW> Hi all, I'm trying to fit a marginal (longitudinal)
VW> model with an exponential serial correlation function to
VW> the Orange tree data set. However, R crashes frequently
VW>
2003 Feb 22
2
4-parameter logistic model
Dear R users
I'm a new user of R and I have a basic question about the 4-parameter
logistic model. According to the information from Pinheiro & Bates the model
is:
y(x)=theta1+(theta2-theta1)/(1+exp((theta3-x)/theta4)) ==
y(x)=A+(B-A)/(1+exp((xmid-input)/scal))
from the graph in page 518 of the book of the same authors (mixed models in
S) theta 1 corresponds to the horizontal asymptote
2007 Mar 03
2
Sigmoidal fitting
I am trying to write a function that fits a sigmoid given a X and Y vector guessing the start parameters.
I use nls. What I did (enclosed) seems to work well with many data points but if I want to fit small
vectors like :
pressure <- c(5,15,9,35,45)
gas <- c(1000,2000,3000,4000,5000)
it do not work. The help page says that it do no not work on zero residual data.
Massimo Cressoni
2009 Mar 27
3
nls, convergence and starting values
"in non linear modelling finding appropriate starting values is
something like an art"... (maybe from somewhere in Crawley , 2007) Here
a colleague and I just want to compare different response models to a
null model. This has worked OK for almost all the other data sets except
that one (dumped below). Whatever our trials and algorithms, even
subsetting data (to check if some singular