similar to: nonlinear least squares trust region fitting ?

Dear Mr Graves, Thank you very much for your response. Nobody else from this mailing list ventured to reply to me for the two weeks since I posted my question. "nlminb" and "optim" are just optimization procedures. What I need is not just optimization, but a nonlinear CURVE FITTING procedure. If there is some way to perform nonlinear curve fitting with the

Hello! I am running R-2.3.1-i386-1 on Slackware Linux 10.2. I am a former matlab user, moving to R. In matlab, via the cftool, I performed nonlinear curve fitting using the method "nonlinear least squares" with the "Trust-Region" algorithm and not using robust fitting. Is it possible to perform the same analysis in R? I read quite a lot of R documentation, but I could not

Dear all, I would like to use the Levenberg-Marquardt algorithm for non-linear least-squares regression using function nls.lm. Can anybody help me to find a a way to compute confidence intervals on the fitted parameters as it is possible for nls (using confint.nls, which does not work for nls.lm)? Thank you for your help Michael

Hi, It seems to me that the most suitable method in R for curve-fitting is the use of nls, which uses a Gauss-Newton (GN) algorithm, while the use of the Levenberg-Marquardt (LM) algorithm does not seem to be very stressed in R. According to this [1] by Ripley, 'Levenberg-Marquardt is hardly competitive these days' which could imply the low emphasize on LM in R. The position of LM is, to

Hi All, Is the Levenberg-Marquardt algorithm available in R. This method combines the steepest descent algorithm and Newton's method. Thanks in Advance, Dermot MacSweeney. ************************************************************** Dermot MacSweeney NMRC, Email: dsweeney at nmrc.ucc.ie Lee Maltings, Tel: +353 21 904178 Prospect Row, Fax: +353 21 270271 Cork, WWW:

Hello, i have a problem with the function nls(). This are my data in "k": V1 V2 [1,] 0 0.367 [2,] 85 0.296 [3,] 122 0.260 [4,] 192 0.244 [5,] 275 0.175 [6,] 421 0.140 [7,] 603 0.093 [8,] 831 0.068 [9,] 1140 0.043 With the nls()-function i want to fit following formula whereas a,b, and c are variables: y~1/(a*x^2+b*x+c) With the standardalgorithm

Hello, I want to fit some non linear models with the Levenberg-Marquardt algorithm. It doesn''t seem to exist any function to do this in R ( well, maybe one does, but I''m a new user, and the only documentation I have is "An introduction to R"). I''d like to know if this function exists, maybe throught an additionnal package. I''d also like to know if if

I'm wondering what packages exist to implement nonlinear least squares regression in R other than 'nls'. Are there packages which implement methods to estimate the optimum values of the parameters which do not use the Gauss-Newton algorithm e.g. use Nelder Mead. In particular, I'd be interested where this is done where the methods of the plinear algorithm are also used (the initial

Dear All, In the past I have often used minpack (http://bit.ly/zXVls3) relying on the Levenberg-Marquardt algorithm to perform non-linear fittings. However, I have always dealt with a function of a single variable. Is there any difference if the function depends on two variables? To fix the ideas, please consider the function f(R,N)=(a/(log(2*N))+b)*R+c*N^d, where a,b,c,d are fit parameters. For

Dear All, I'd like to fit a "kind" of logistic model to small data-set using nonlinear least-squares regression. A transcript of R-script are reproduced below. Estimated B and T (the model's coeff, herein B=-8,50 and T=5,46) seem appropriate (at least visually) but are quite diff from those obtained w/ SPSS (Levenberg-Marquardt): B=-19,56 and T=2,37. Am I doing something wrong in

Hi, I have a problem fitting the following Weibull Model to a set of data. The model is this one: a-b*exp(-c*x^d) If I fitted the model with CurveExpert I can find a very nice set of coefficients which create a curve very close to my data, but when I use the nls.lm function in R I can't obtain the same result. My data are these: X Y 15 13 50 13 75 9 90 4 With the commercial

Hi! I'm translating some code from Matlab to R and I found a problem. I need to translate Matlab's function 'lsqnonlin' (http://www-ccs.ucsd.edu/matlab/toolbox/optim/lsqnonlin.html) into R, and at the beginning I thought it would be the same as R's 'optim'. But then I looked at the definition of 'lsqnonlin' and I don't quite see how to make

I'm trying to fit a model to solve a biological problem. There are multiple independent variables, and also there are multiple responses. Each response is a function of all the independent variables, plus a set of parameters. All the responses depend on the same variables and parameters - just the form of the function changes to define each seperate response. Any ideas how I can fit

I have used Mike Powell's optimization routine (VA05AD) from the Harwell Subroutine Library (HSL) for more than 20 years. It is no exaggeration to say that it has helped make my career (thanks Mike). I recently learned that I am not alone in this respect - apparently it still has a loyal following in all sorts of fields! It is an exceedingly fine piece of software - fast, reliable and easy to

In the course of revising a paper I have had occasion to attempt to maximize a rather complicated log likelihood using the function nlm(). This is at the demand of a referee who claims that this will work better than my proposed use of a home- grown implementation of the Levenberg-Marquardt algorithm. I have run into serious hiccups in attempting to apply nlm(). If I provide gradient and

Hello guys. I've built and uploaded to CRAN an R interface to MINPACK Fortran library, which solves non-linear least squares problem by modification of the Levenberg-Marquardt algorithm. The package includes one R function, which passes all the necessary control parameters to the appropriate Fortran functions. The package location is

Hello guys. I've built and uploaded to CRAN an R interface to MINPACK Fortran library, which solves non-linear least squares problem by modification of the Levenberg-Marquardt algorithm. The package includes one R function, which passes all the necessary control parameters to the appropriate Fortran functions. The package location is

Dear all, When doing nonlinear regression, we normally use nls if e are iid normal. i learned that if the form of the variance of e is not completely known, we can use the IRWLS (Iteratively Reweighted Least Squares ) algorithm: for example, var e*i =*g0+g1*x*1 1. Start with *w**i = *1 2. Use least squares to estimate b. 3. Use the residuals to estimate g, perhaps by regressing e^2 on

Hi guys I have data that contains the variances vt of the yields of 1, 2, 3, 4, 5,10, 20 year bonds. Assuming the Hull-White model for the yield of a t-year zero-coupon bond, I have to estimate the ? of the Hull-White model using nonlinear least squares and give a 95% con?dence interval for each parameter. Please can you guys tell how to find out ? using R. Any suggestion regarding what functions

Hi, I am new to R, and am trying to solve the following optimization problem: This is a nonlinear least squares problem. I have a set of 3D voxels. All I need is to find a least squares fit to this data. The data model actually represent a cube-like structure, consisting of seven straight lines. The lines have some intersections (and at this intersection both of the participating lines end).