similar to: Results of nls

Hi, I would like to be able to access the calculated results from the nls package. Using the example in R, fm3DNase1 we can reurn certain parts of the calculations: > coef(fm3DNase1) Asym xmid scal 2.345179 1.483089 1.041454 > resid(fm3DNase1) [1] -0.0136806237 -0.0126806237 0.0089488569 0.0119488569 -0.0025803222 [6] 0.0064196778 0.0026723396 -0.0003276604

Hi! How to write in model's formula of type: nls(y~k/(x^n), data=data, start=list(k=1,n=1)) i.e the problem is on x^n, I(x^n) generate error thanks Jarek

Hi, I am trying to bootstrap the difference between each parameters among two non linear regression (distributed loss model) as following: # data.frame > Raies[1:10,] Tps SolA Solb 1 0 32.97 35.92 2 0 32.01 31.35 3 1 21.73 22.03 4 1 23.73 18.53 5 2 19.68 18.28 6 2 18.56 16.79 7 3 18.79 15.61 8 3 17.60 13.43 9 4 14.83 12.76 10 4 17.33 14.91 etc... # non

I have found some "issues" (bugs?) with nls confidence intervals ... some with the relatively new "port" algorithm, others more general (but possibly in the "well, don't do that" category). I have corresponded some with Prof. Ripley about them, but I thought I would just report how far I've gotten in case anyone else has thoughts. (I'm finding the code

Colleagues, I am using "nls" successfully (2.11.1, OS X) but I am having difficulties retrieving part of the output - residual sum of squares. I have assigned the output to FIT: > > FIT > Nonlinear regression model > model: NEWY ~ PMESOR + PAMPLITUDE * cos(2 * pi * (NEWX - POFFSET)/PERIOD) > data: parent.frame() > PMESOR PAMPLITUDE POFFSET >

Hello, In doing a series of non-linear estimations of a function which is a sum of a varying number of sinusoids, I would like to "autogenerate" the arguments needed by nls() depending on that number. For example, when there are two sinusoids: > nls( y ~ mu + A1 * cos(2*pi*f1*x - P1) + A2 * cos(2*pi*f2*x - P2), data = some.xy.data, start = list( mu=some.value0,

I've got a PhD student starting this year. She'll be working on data mining. She has asked for a reading list while she's still in her home country, which is a really good sign. Her cosupervisor and I don't (think we) have any problem with choosing things that are specifically about data mining, but there are some statistical ideas (sampling, exploratory-vs-confirmatory,

Hi all, I need to use the parameter estimates of nls() for further analysis. I know how to do in S+, e.g. nls(...)$parameters. In R, the attributes(nls(...)) does not have parameters, how would one get the parameter values out of nls()? Thanks in advance. Nancy

Hello, I've configured a squid proxy server to use the wb_auth module to authenticate NT Users . I'm aware that this is not the squid users mailing list but since the wb_auth modle uses the winbind application i thought this might be a good place to look for some help. Basically, the wb_auth module asks winbind to authenticate a users. Winbind then returns a '0' if the

You also need to reply-all so the mailing list stays in the loop. -- Sent from my phone. Please excuse my brevity. On December 19, 2017 4:00:29 PM PST, Timothy Axberg <axbergtimothy at gmail.com> wrote: >Sorry about that. Here is the code typed directly on the email. > >qe = (Qmax * Kl * ce) / (1 + Kl * ce) > >##The data >ce <- c(15.17, 42.15, 69.12, 237.7, 419.77)

We are trying to compile and run the ttest.c example that comes with R (in C:\Program Files\R\rw1081\src\library\windlgs\src\ttest.c). After compiling it with MS Visual C++ we load the DLL with dyn.load. So far it seems good, but when we try to call it from R (after running C:\Program Files\R\rw1081\src\library\windlgs\R\windlgs.R) R crashes. We have tried changing the exports from DLL but have

Dear Sir I am new user of R. I am interested in modeling hydrological extreme events. I found MSClaio2008 very interesting function. In this function four criterions for choosing distributions. Can we call these criterions as model selection techniques or goodness of fit techniques or both? Because goodness of fit techniques are usually performed after modle selection. Can I found

Saludos a todos, Quisiera saber como ocultar algunos componentes de una lista para que cuando crees una funcion solo se muestre una parte de esa lista, y que al resto tengas acceso cuando descompongas ese objeto, por ejemplo: Tengo la siguiente data: data=data.frame(Samples=1:10,Sobs=c(12.6,22.4,27.4,30.5,32.5,34.3,36,37.4,38.9,39.9)) Cuando uso la funcion nls: A <-

Hello I followed the example in page 59, chapter 11 of the 'Introduction to R' manual. I entered my own x,y data. I used the least squares. My function has 5 parameters: p[1], p[2], p[3], p[4], p[5]. I plotted the x-y data. Then I used lines(spline(xfit,yfit)) to overlay best curves on the data while changing the parameters. My question is how do I calculate the residual sum of squares.

Hola Javier, Creo que comprendio mal mi pregunta, quiero saber como ocultar algunos componentes de una lista producida por una funcion que uno crea, pero que solo esten ocultos y que pueda tener acceso a esos componentes cuando uso la funcion "str", la cual despliega todos los componentes de la lista. Además como hacer para que salga ordenado o de una forma preferida los componentes

Hola Eva, Gracias por la respuesta fue de mucha ayuda. Lo que me gustaria saber es por que solo algunas funciones que usan el metodo S3 son visibles en su estructura con la funcion "body". Por ejemplo si coloco body(print.ts) Me va a salir la estructura de esa funcion, sin embargo si coloco body(print.nls) Me sale el siguiente mensaje :Error in body(print.nls) :

Hello there, I would like to test the agreement amongst 6 raters for nominal data on a scale from 1-4, and conduct a latent class analysis in R. How should the data be formatted and what code should I use? Thank you very much Lisa Wang Princess Margaret Hospital Biostatistics tel:416 946 4501

Hi, I am a new user of R. I am trying to use R to do some nonlinear modeling. However, I can not find how to get the parameters's value back. Here is the code in S-plus: ycos1.nl <- nls(y ~ m +ac*cos(2*pi*f*t) +as*sin(2*pi*f*t), start=nl.st, trace=T) ## store model info after satisfactory model has been found cNL1 <- c(ycos1.nl$parameters[1:3],0) What is same code in R? Thanks a

> Henrik Bengtsson: > > I'm looking for a way to get the length of an object 'x' as given by > base data type without dispatching on class. The performance improvement you're looking for is implemented in the latest version of pqR (pqR-2016-10-24, see pqR-project.org), along with corresponding improvements in several other circumstances where unclass(x) does not

I'm looking for a way to get the length of an object 'x' as given by base data type without dispatching on class. Something analogous to how .subset()/.subset2(), e.g. a .length() function. I know that I can do length(unclass(x)), but that will trigger the creation of a new object unclass(x) which I want to avoid because 'x' might be very large. Here's a dummy example