similar to: Robust Non-linear Regression

Dear R-wizards, I'm trying to plot "binned scatterplots", or 2d histograms, if you wish, for a number of groups by using the lattice functionality it works fine for one group at a time, and probably I could find a work-around, but I prefer to do it the elegant way here's an example of what I want, what I tried and where it goes wrong: require(gregmisc) require(lattice) #toy

Hello useRs (please don't kill me), I've fairly new to R having only a few months of playing around with R. What little I've learned has been extremely useful. If someone could point me as to how to replicate the Boltzman Sigmodial curve fit as provided by Graphpad software I'd be eternally grateful. Where we currently use Graphpad for only this one function,its seems highly

Hi friends I just started using wine on my Linux-machine. Mostly because I want to switch finally from Windows to Linux on my computer at work. There is just one single application I want to "take with me" and that is "GraphPad Prism" (http://www.graphpad.com/prism/Prism.htm). So I installed wine0.9.20 on my Mandriva2007. Then I used wine to install Prism in its most recent

I'm having problems getting nls to agree that convergence has occurred in a toy problem. nls.out never gets defined when there is an error in nls. Reaching the maximum number of iterations is alway an error, so nls.out never gets defined when the maximum number of iterations is reched. >From ?nls.control: tol: A positive numeric value specifying the tolerance level for the

Dear R-help subscribers Would you tell me how to apply SSfpl with binary data as below? Unfortunately, there is not the EXAMPLE in help(SSfpl) for binary data but for quantitative data(Chick). V1: dose V2: log-transformed dose V3: response (rate) V1 V2 V3 1 0.775 -0.2548922 0.1666667 2 5.000 1.6094379 0.8148148 3 10.000 2.3025851 0.5000000 4 20.000 2.9957323

Hi, I was hoping to get some advice on how to derive estimates of slopes from four parameter logistic models fit with SSfpl. I fit the model using: model<-nls(temp~SSfpl(time,a,b,c,d)) summary(model) I am interested in the values of the lower and upper asymptotes (parameters a and b), but also in the gradient of the line at the inflection point (c) which I assume tells me my rate of

Hi there, I have installed MS Office 2007 in wine 1.2.2 running in LinuxMint KDE 9, together with other softwares (including GraphPad Prism 4). After install, I set riched20.dll as native (windows), as stipulated here: http://appdb.winehq.org/objectManager.php?sClass=version&iId=4992 > Post Installation Instructions? > > Once installed, one override is necessary. Without it,

Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to calculate the kriging parameters using REML, including the parameters of a sigmoidal trend function

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

Dear R Community, I am using the package DRC ( to fit a boltzman model to my data. I can fit the model and extract the lower limit, upper limit, and ED50 (aka V50), but I cannot figure out how to get the slope of the curve at ED50. Is there a simple way to do this? I've searched the mailing list and looked through the package documentation, but could not find anything. I am new to r, and

HI I'm a student in chemical engineering, and i have to implement an algoritm about FIVE PARAMETERS INTERPOLATION for a calibration curve (dose, optical density) y = a + (c - a) /(1+ e[-b(x-m]) where x = ln(analyte dose + 1) y = the optical absorbance data a = the curves top asymptote b = the slope of the curve c = the curves bottom asymptote m = the curve X intercept Have you never seen

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

First, thanks to all who helped me with my question about rescaling axes on the fly. Using unlist() and range() to set the axis ranges in advance worked well. I've since plotted about 300 datasets with relative ease. Now I'm trying to fit a lossy oscillator resonance to (the square root of) a lorentzian (testframe$y is oscillator amplitude, testframe$x is drive frequency): lorentz

Hello, We use drc to fit dose-response curves, recently we discovered that there are quite different standard error values returned for the same dataset depending on the drc-version / R-version that was used (not clear which factor is important) On R 2.9.0 using drc_1.6-3 we get an IC50 of 1.27447 and a standard error on the IC50 of 0.43540 Whereas on R 2.7.0 using drc_1.4-2 the IC50 is

Your understanding isn't similar to mine. Mine says robust/resistant methods are for data with heavy tails, not heteroscedasticity. The common ways to approach heteroscedasticity are transformation and weighting. The first is easy and usually quite effective for dose-response data. The second is not much harder. Both can be done in R with nls(). Andy From: Quin Wills > > I am

Dear all I am struggling to fit data which form something like CDF by lognorm. Here are my data: proc <- c(0.9, 0.84, 0.5, 0.16, 0.1) size <- c(0.144, 0.172, 0.272, 0.481, 0.583) plot(size, proc, xlim=c(0,1), ylim=c(0,1)) fit<-nls(proc~SSfpl(size, 1, 0, xmid, scal), start=list(xmid=0.2, scal=.1)) lines(seq(0,1,.01), predict(fit, newdata=data.frame(sito=seq(0,1,.01))), col=2) I tried

Howdee, My question appears at #6 below: 1. I want to model the growth of each of a large number of individuals using a 4-parameter logistic growth curve. 2. nlme does not converge with the random structure that I want to use. 3. nlsList does not converge for some individuals. 4. I decided to go around nlsList using: t(sapply(split(data, list(data$id)), function(subd){coef(nls(mass ~

I am using "nls" to fit dose-response curves but am not sure how to approach more robust regression in R to get around the problem of the my error showing increased variance with increasing dose. My understanding is that "rlm" or "lqs" would not be a good idea here. 'Fairly new to regression work, so apologies if I'm missing something obvious.

Dear all I try to refine my nlme models and with partial success. The model is refined and fitted (using Pinheiro/Bates book as a tutorial) but when I try to plot plot(augPred(fit4)) I obtain Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels (0,3.5],(3.5,5],(5,7],(7,Inf] not allowed for vykon.fac > Is it due to the fact that I have unbalanced

Dear all, I am using the nlsLM function to fit a Lorentzian function to my experimental data. The LM algorithm should allow to specify limits, but the upper limit appears not to work as expected in my code. The parameter 'w', which is peak width at half maximuim always hits the upper limit if the limit is specified. I would expect the value to be in-between the upper and lower limit with