similar to: (no subject)

Dear R-users, I am trying to create a model using the NLS function, such that: Y = f(X) + q + e Where f is a nonlinear (Weibull: a*(1-exp(-b*X^c)) function of X and q is a covariate (continous variable) and e is an error term. I know that you can create multiple nonlinear regressions where x is polynomial for example, but is it possible to do this kind of thing when x is a function with unknown

Hi there, I am trying to compare nonlinear least squares regression with AIC and anova. The simplest model is one nonlinear curve, and in the more complex model I have a categorical variable (producing parameter estimates for four curves). Both models run fine, but when I try to produce an AIC value for the second model I get the error: > AIC(pow.nls1) [1] 114408.3 > AIC(pow.nls2) Error in

Dear All, A couple of questions about the nls package. 1. I'm trying to run a nonlinear least squares regression but the routine gives me the following error message: step factor 0.000488281 reduced below `minFactor' of 0.000976563 even though I previously wrote the following command: nls.control(minFactor = 1/4096), which should set the minFactor to a lower level than the default

  Hi, I fitted tree growth data with Chapman-Richards growth function using nls. summary(fit.nls) Formula: Parameters: Estimate Std. Error t value Pr Signif. codes: 0 ''***'' 0.001 ''**'' 0.01 ''*'' 0.05 ''.'' 0.1 '' '' 1 Residual standard error: 1.879 on 713 degrees of freedom Algorithm

Hi there, I am looking to compare nonlinear mixed effects models that have different nonlinear functions (different types of growth curve)embedded. Most of the literature I can find focuses on comparing nested models with likelihood ratios and AIC. Is there a way to compare model fits when models are not nested, i.e. when the nonlinear functions are not the same? Many thanks in advance! Lindsay

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:

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

Hello colleagues, I am attempting to determine the nonlinear least-squares estimates of the nonlinear model parameters using nls. I have come across a common problem that R users have reported when I attempt to fit a particular 3-parameter nonlinear function to my dataset: Error in nls(r ~ tlm(a, N.fix, k, theta), data = tlm.data, start = list(a = a.st, : step factor 0.000488281

When I run gnls I get the error: Error in nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xy, : step factor 0.000488281 reduced below 'minFactor' of 0.000976563 My first thought was to decrease minFactor but gnlsControl does not contain minFactor nor nlsMinFactor (see below). It does however contain nlsMaxIter and nlsTol which I assume are the analogs of

Dear Friends. I tried to use nls.control() to change the 'minFactor' in nls( ), but it does not seem to work. I used nls( ) function and encountered error message "step factor 0.000488281 reduced below 'minFactor' of 0.000976563". I then tried the following: 1) Put "nls.control(minFactor = 1/(4096*128))" inside the brackets of nls, but the same error message

Hello. I am trying to do a non-linear fit using the 'nls' command. The data that I'm using is as follows pH k 1 3.79 34.21 2 4.14 25.85 3 4.38 20.45 4 4.57 15.61 5 4.74 12.42 6 4.92 9.64 7 5.11 7.30 8 5.35 5.15 9 5.67 3.24 with a transformation of pH to H <- 10^-pH When using the nls command for a set of parameters - a, b and c, I receive two sets of errors: >

I'm trying to fit a Gompertz sigmoid as follows: x <- c(15, 16, 17, 18, 19) # arbitrary example data here; y <- c(0.1, 1.8, 2.2, 2.6, 2.9) # actual data is similar gm <- nls(y ~ a+b*exp(-exp(-c*(x-d))), start=c(a=?, b=?, c=?, d=?)) I have been unable to properly set the starting value '?'s. All of my guesses yield either a "singular gradient" error if they

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 :

Hi to All, I am fitting some models to a data using non linear least square, and whenever i run the command, parameters value have good convergence but I get the error in red as shown below. Kindly how can I fix this problem. Convergence of parameter values 0.2390121 : 0.1952981 0.9999975 1.0000000 0.03716107 : 0.1553976 0.9999910 1.0000000 0.009478433 : 0.2011017 0.9999798 1.0000000

Often it is useful to keep a "codebook" to document the contents of a dataset. (By "dataset" I mean a rectangular structure such as a dataframe.) The codebook has as many rows as the dataset has columns (variables, fields). The columns (fields) of the codebook may include: ? variable name ? type (character, factor, integer, etc) ? variable label

Dear R experts, I remember a similar function existed and have been mentioned in R-help before. I tried my best to search but I really can't find it out. suppose I have an data frame like this: > somedata <- data.frame(age.min = 1, age.max = 1.5, male = TRUE, l = -1.013, m=16.133, s=0.07656) In order to back up the data and I don't want to use write.table(), I would like to back

R-help, I'm using a function whose end result is a trellis plot. When I call the function I get sometimes the following message: "Note: The default device has been opened to honour attempt to modify trellis settings " leading up to any plot whatsoever. I call the package 'lattice' within the function and 'detach' after plotting. I have to close the graphics

Fellow R users, I am stumped on what would seem to be something fairly simple. I have a dataframe that has a variable named 'WEEK' that takes the numbers 1:26 (26 week time-period) with each number repeated five times consecutively (once for each weekday, Monday through Friday). Ex. 111112222233333.....2626262626. I would like to randomly extract two weekdays per five day week for

Hello All, I am trying to work on writing the following piece of (pseudo) code in an optimal fashion: ---------------------------------------------------- # Two data frames with some data a = data.frame(somedata) b = data.frame(somedata) for(i in 1:nrow(dt) { # Merge dates for a given date into a new data frame c = merge(a[a$dt==dt[i],),b[b$dt == dt[i],], by=c(some column)); }

Hi Eli, Slightly changed the example :- .section somedata .Ltmpa: .word 100 .Ltmpb: .word 60 .Ltmpc: .word 40 .Ltmpd: .word (.Ltmpa-.Ltmpb)-.Ltmpc The GNU assembler properly handles the above. $as -c x.s -o x.o $readelf -x 4 x.o Hex dump of section 'somedata': NOTE: This section has relocations against it, but these have NOT been applied to this dump. 0x00000000 64003c00