similar to: nls.control()

I have data that looks like O.lengthO.age 176 1 179 1 182 1 ... 493 5 494 5 514 5 606 5 462 6 491 6 537 6 553 6 432 7 522 7 625 8 661 8 687 10 704 10 615 12 (truncated) with a simple VonB growth model from within nls(): plot(O.length~O.age, data=OS) Oto = nls(O.length~Linf*(1-exp(-k*(O.age-t0))), data=OS, start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE) mod <- seq(0, 12)

I apologize if this comes across as confusing. I will try to explain my situation as best I can. I have R bootstrapping my growth data for fish. It's resampling my database of age and length data and then produces several new datasets for me. In this case, it's resampling my data to create three new datasets of age and length data. Here is my code with my original data called

I'm tyring to fit a relatively simple nls model to some data, but keep coming up against the same error (code follows): Oto=nls(Otolith ~ Linf*(1-exp(-k(AGE-to))), data = ages, start = list(Linf=1000, k=0.1, to=0.1), trace = TRUE) The error message I keep getting is "Error in eval(expr, envir, enclos) : could not find function "k"". I've used this

Dear colleagues, Von Bertalanffy model is commonly adjust to data on fish length (TL) and age (AGE) TL= Linf*(1-exp(-K*(AGE-t0)). Linf, K and t0 are parameters of the model. One main goal of the growth study is the comparison of growth parameter estimates between sexes of the same species, or estimates from different populations. The realibility statistical tests normally applied are highly

Thanks to the help of people from this forum I was able to bootstrap my data and then apply a model to it. Thanks for all your help. Everything worked out well, but I am having a difficult time getting the new parameter values. I bootstrapped the data 300 times and I want to get the 300 sets of parameter estimates and plot them in Excel. Here is my code:

Dear R-list members, Some months ago I wrote a message on the usage of gnls (nlme library) and here I come again. Let me give an example: I have a 10 year length-at-age data set of 10 fishes (see growth.dat at the end of this message) and I want to fit a von Bertalanffy growth model, Li= Linf*(1-exp(-k*(ti-t0))) where Li = length at age i, Linf= asymptotic length, k= curvature parameter, ti=

Hi, I have great problems with my work in R. I look for to model the growth of fish. I have "Longitudinal data", a serie of repeated measures for each individual. Using the corresponding packages "nlme" in R. I treat to fit to the data different growth functions, wich were entered by me. Unfortunately for no it was arrived at the convergence, several error messages appeared. I

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

Full_Name: Ben Bolker Version: 2.2.1 OS: Windows XP and 2000 Submission from: (NULL) (128.227.60.124) The following code, using confint() to try to get confidence intervals on an nls object that has been fitted with algorithm="port" reliably crashes R 2.2.0 and 2.2.1 with the latest version of MASS on a Windows 2000 and a Windows XP machine here. I *think* earlier versions of MASS

All - I'm having an issue with trying to plot a model derived from nls() onto a simple plot.? I have included a sample data set and the code that I've been using. ?? year month day?????? date location mileage? cost gallon????? cpg ? mpg????????? x 2009???? 1?? 4?? 1/4/2009????? BZN? 124585 19.39? 14.37 1.349339 10.71677 2009-01-04 2009???? 1? 15? 1/15/2009????? BZN? 124888? 23.2? 16.12

Hello all, This is the error message that I get. > hyp.res <- nls(log(y)~log(pdf.hyperb(theta,X)), data=dataModel, + start=list(theta=thetaE0), + trace=TRUE) 45.54325 : 0.1000000 1.3862944 -4.5577142 0.0005503 3.728302 : 0.0583857346 0.4757772859 -4.9156128701 0.0005563154 1.584317 : 0.0194149477 0.3444648833 -4.9365149150 0.0004105426 1.569333 :

Howdy, Last week I got some great help on why I was getting an error code when trying to run this model, thanks everyone!  I was able to get the code up and running beautifully for several data sets.  Now I am getting different errors with this new data set.  I can't figure out why, I have more data points with this species, and it is ordered exactly the same as the other species I have been

Hi I?m trying to fit a nonlinear model to a derivative of the logistic function y = a/(1+exp((b-x)/c)) (this is the parametrization for the SSlogis function with nls) The derivative calculated with D function is: > logis<- expression(a/(1+exp((b-x)/c))) > D(logis, "x") a * (exp((b - x)/c) * (1/c))/(1 + exp((b - x)/c))^2 So I enter this expression in the nls function:

Hi, I' m having a hard time understanding the computation of degrees of freedom when runing nlme() on the following model: > formula(my data.gd) dLt ~ Lt | ID TasavB<- function(Lt, Linf, K) (K*(Linf-Lt)) my model.nlme <- nlme (dLt ~ TasavB(Lt, Linf, K), data = my data.gd, fixed = list(Linf ~ 1, K ~ 1), start = list(fixed = c(70, 0.4)), na.action= na.include,

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

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

Hello, I am attempting to fit a non-linear model (Von Bertalanffy growth model) to fish length-at-age data with the purpose of determining if any of the three parameters differ between male and female fish. I believe that I can successfully accomplish this goal assuming an additive error structure (illustrated in section 1 below). My trouble begins when I attempt this analysis using a model

Dear All., I've created R-code for which a user will be asked to choose between 2 analyses. I've written one function for each type of analysis. Within each function, the users is prompted to enter information. An example is: cat("Enter value for lower Linf :\n") L1<-scan(n=1) cat("Enter value for upper Linf :\n") L2<-scan(n=1)

Hi, I have some data I want to fit with a non-linear function using nls, but it won't solve. > regresjon<-nls(lcfu~lN0+log10(1-(1-10^(k*t))^m), data=cfu_data, > start=(list(lN0 = 7.6, k = -0.08, m = 2))) Error in nls(lcfu ~ lN0 + log10(1 - (1 - 10^(k * t))^m), data = cfu_data, : step factor 0.000488281 reduced below 'minFactor' of 0.000976562 Tried to increase minFactor

Dear all i want to estimated the parameter of the gamma density(a,b,d) f(x) = (1/gamma(b)*(a^b)) * ((x-d)^(b-1)) * exp{-(x-d)/a)} for x>d f(x) = Age specific fertility rate x = age when i run this in R by usling nls() gamma.asfr <- formula(asfr ~ (((age-d)^(b-1))/((gamma(b))*(a^b)))* exp(-((age-d)/a))) gamma.asfr1 <- nls(gamma.asfr, data= asfr.aus, start = list(b = 28, a = 1, d= 0.5),