search for: ssasympoff

Hello: I am new on R, this is the first time I work with it, I am trying to run some R example to get learn of them, but, I have a problem when running the example below: http://stat.ethz.ch/R-manual/R-patched/library/datasets/html/zCO2.html CO2.Qn1 <- CO2[CO2$Plant == "Qn1", ] SSasympOff(CO2.Qn1$conc, 32, -4, 43) # response only Asym <- 32; lrc <- -4; c0 <- 43 SSasympOff(CO2.Qn1$conc, Asym, lrc, c0) # response and gradient getInitial(uptake ~ SSasympOff(conc, Asym, lrc, c0), data = CO2.Qn1) ## Initial values are in fact the converged values fm1 <- nls(uptake ~ SSasympO...

...e one major difference is that in my case the 'conc' steps are not the same for each 'Plant'. I have replicated the problem using the CO2 data in nlme (based off of the Ch08.R script). This works (when 'conc' is the same for each 'Plant': (fm1CO2.lis <- nlsList(SSasympOff, CO2)) (fm1CO2.nlme <- nlme(fm1CO2.lis, control = list(tolerance = 1e-2))) (fm2CO2.nlme <- update(fm1CO2.nlme, random = Asym + lrc ~ 1)) CO2.nlme.var <- update(fm2CO2.nlme, fixed = list(Asym ~ Type * Treatment, lrc + c0 ~ 1), start = c(32.412, 0, 0, 0, -4.5603, 49.344), weights=varCons...

....1-65. on a WINDOWS 2000 machine I am trying to run the scripts from the Mixed Models book and am running into some difficulty - could any one tell me why this doesn't work (from pages 369 - 370) CO2 #FINE plot(CO2,outer= ~Treatment*Type,layout=c(4,1)) #FINE fm1CO2.lis <- nlsList( SSasympOff, CO2 ) #FINE fm1CO2.lis #FINE fm1CO2.nlme <- nlme( fm1CO2.lis ) #COPIED DIRECTLY FROM BOOK!!?? Gives the error message:- Error in nlme.formula(model = uptake ~ SSasympOff(conc, Asym, lrc, c0), : Maximum number of iterations reached without convergence Any assistance would be...

Hi there are some predefined self-start functions, like SSmicmen, SSbiexp, SSasymp, SSasympOff, SSasympOrig, SSgompertz, SSflp, SSlogis, SSweibull, Quadratic, Qubic, SSexp (nlrwr) Btw, do you know graphic examples for this functions? The SSexpDecay (exponential decay) for y = (y0 - plateau)*exp(-k*x) + plateau from http://stats.stackexchange.com/questions/17126/how-can-i-make-my-r-nls-mode...

...wo treatments and two types (84 observations in total), but den df for each parameter of the model is 64. Isn’t it too high? Your help is greatly appreciated, Julia Summary of the CO2 example: > summary(fm4CO2.nlme) Nonlinear mixed-effects model fit by maximum likelihood Model: uptake ~ SSasympOff(conc, Asym, lrc, c0) Data: CO2 AIC BIC logLik 388.4185 420.0191 -181.2092 Random effects: Formula: list(Asym ~ 1, lrc ~ 1) Level: Plant Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr Asym.(Intercept) 2.349640 As.(I) lrc.(...

...extranames, : variable lengths differ I cannot find out why the grouping factor cannot be used (it has the same length as the other variables...) Another strange thing occurs: in the example given in the help of nlsList.selfstart, the following command works well: fm1 <- nlsList(SSasympOff, CO2) However its seemingly equivalent applied to the case above fails: mymod4<-nlsList(SSfol,data=mydata) Error in eval(expr, envir, enclos) : object "input" not found Error in eval(expr, envir, enclos) : object "input" not found Error in eval(expr, envir, enclos) : object...

...1.00 1 AksToKa 7 3. Summer foothills Half 1.25 1 AksToKa 8 4. Fall foothills Half 1.50 1 AksToKa 9 1. Fall-Winter foothills Meat 0.60 1 AksToKa 10 2. Spring foothills Meat 1.00 1 AksToKa > kew11.nlme$call nlme.formula(model = lw ~ SSasympOff(ageyr2, mw, lgr, age0), data = kew, fixed = list(mw + lgr + age0 ~ Season * MBreed + ecoreg + ecoreg:Season), random = list(farm = list(mw ~ 1, lgr ~ 1), sheep = list(mw ~ 1)), start = c(fixef (kew8.nlme)[1:15], 0, 0, 0, 0, 0, 0, 0, 0, 0, fixef(kew8.nlme)[16:30], 0,...

...uot;coefficients" and "parameters"? The have the same "Estimate", different se (therefore t value), but same p values. R.2.8.0 on winxp with nlme_3.1-89 Thanks, ...Tao +++++++++++++++++++++++++++++++++++++++++++++++++ > library(nlme) > fm1 <- nlsList(uptake ~ SSasympOff(conc, Asym, lrc, c0), data = CO2, start = c(Asym = 30, lrc = -4.5, c0 = 52)) > summary(fm1)$para[,,1] Estimate Std. Error t value Pr(>|t|) Qn1 38.13977 0.9911148 38.48169 1.991990e-06 Qn2 42.87169 1.0932089 39.21638 2.583953e-06 Qn3 44.22800 1.0241029 43.18706 1.809264e-07 Qc1...

I am trying to compare Vbert growth curves from several years of fish data. I am following the code provided by: http://www.ncfaculty.net/dogle/fishR/gnrlex/VonBertalanffy/VonBertalanffy.pdf. Specifically the section on VBGM Comparisons between groups. ? This code is pretty cut and dry. I am able to run it perfectly with the "fake" data that is provided. But when I run it with my own

Hi there: I'm trying to run a large number of nonlinear regressions on a time series dataset, where the data will be formatted something to the effect of: ObservationID,time,dependentvar I'll have a number of time and dependentvars for each observation, and I want to apply a nonlinear regression to one ObservationID at a time, and I want to have a dataset that is the parameter

...sestX.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/NLSstLfAsymptote.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/NLSstRtAsymptote.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/SSasymp.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/SSasympOff.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/SSasympOrig.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/SSbiexp.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/SSfol.R action: create D: file: /usr/local/lib/R/library/nls/R-ex/SSfpl.R action: creat...

...html latex example NLSstRtAsymptote text html latex example Orange text html latex example Puromycin text html latex example SSasymp text html latex example SSasympOff text html latex example SSasympOrig text html latex example SSbiexp text html latex example SSfol text html latex example SSfpl text...