similar to: Weights in nls function

> Date: Thu, 10 Aug 2000 15:33:39 -0700 > From: Jesus Maria Frias Celayeta <iosu at bureau.ucc.ie> > > hi R-help! > > I would like to know a simple and easy way (if posible) to do the following > data manipulation. I have a matrix of experimental data (with replicae in time) actually, they seem to be data frames. > >

I'm written a function that uses nls in a loop to fit a model to resampled data. Occasionally a sample is drawn such that the algorithm (apparently) can't converge and step factor reduced below minimum is printed to the screen. I tried a few different settings for the minFactor parameter in nls.control with no different result. When this happens, I would like my function to put NAs

Dear R-helpers, Would anybody have a compiled binary version (for Windows) of the Dopt library that is in the development site?. best regards, IOsu -------------------------------------------------------------- Jes?s Mar?a Fr?as Celayeta Dublin Institute of Technology School of Food Sci. and Env. Health. Cathal Bruha St., Dublin 1. Ireland Phone: +353 1 4024459 Fax: +353 1 4024495

Dear all, Has anyone written functions to compute (and ideally project / plot onto a plane) confidence ellipsoids for q linear model parameters, q > 2 ? Greetings and best wishes, Renaud -- Dr Renaud Lancelot, v?t?rinaire CIRAD, D?partement Elevage et M?decine V?t?rinaire (CIRAD-Emvt) Programme Productions Animales http://www.cirad.fr/presentation/programmes/prod-ani.shtml ISRA-LNERV

hi R-help! I would like to know a simple and easy way (if posible) to do the following data manipulation. I have a matrix of experimental data (with replicae in time) >experiment<-data.frame(times=c(0,0,10,10,20,20,30,30),expval=c(1,1,2,2,3,3,4,4)) > experiment times expval 1 0 1 2 0 1 3 10 2 4 10 2 5 20 3 6 20 3 7 30 4 8

I've been looking at the ellipse package I co-wrote with Ed Chow for S-Plus a few years ago, cleaning up the documentation now it's been ported to R (by Jesus Frias, thanks IOsu!) and I'm back as the maintainer. One of the things it did which we thought was kind of cute at the time was to use a few lines of Postscript code to define a font, where the glyphs were ellipses of various

Hi R-helpers! I have a problem with the xfig device and I don't know how to tackle. If anybody can give me some help I will be very thankful: *I have produced a series of figures using a function that makes a call to barplot (in fact I am making Pareto Charts using a function that was kindly provided through this list :-) *I have exported the figures to xfig in order to give them a

Last week (Wed 9/10/2003, "regression questions") I posted a question regarding the use of gnls( ) and its dissimilarity to the syntax that nls( ) will accept. No one replied, so I partly answered my own question by constructing indicator variables for use in gnls( ). The code I used to construct the indicators is at the end of this email. I do have a nagging, unanswered

Does the nls package actually allow for weighted nonlinear regression? (i.e., I have data with individual variances associated, I'd like to use 1/var to weight the points.) The "nls" function does have a weights argument, but it doesn't seem to do anything as far as I can tell ... thanks ... Ben Bolker -- 318 Carr Hall bolker at

Dear Colleagues, Our group is also working on implementing the use of R for pharmacokinetic compartmental analysis. Perhaps I have missed something, but > fit <- nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=0.5, k2=0.5)), + data=C1.lsoda, + start=list(K1=0.3, k2=0.7), + trace=T + ) Error in eval(as.name(varName), data) : Object

Forgive this naive question, but I'm not accustomed to using Win ports of Unix tools (actually I'm unfamiliar with Unix tools, period). Anyway, I've successfully compiled the R-1.1.1 sources into the full binary package. Everything seems to work. I'd like to remove all the "made" files and leave only the original source code and the compiled binary code. Is there

Dear R-helpers, I’ve got a summary of results from a by() call that I am making with a list of more than two of factors not very different from the example in the by() help page require(stats) by(warpbreaks[, 1], warpbreaks[, -1], summary) The result of the command gives a list of the form wool: A tension: L Min. 1st Qu. Median Mean 3rd Qu. Max. 25.00 26.00

Dear R-helpers, I would like to change the default value of a functions that is call by a package. The function deparse() is called by a set of functions from a library and because I don't think is appropriate to change the library, I am trying to change the default value in deparse(). At the moment the only solution I have come across is to place a function in my working environment like

Dear R-helpers, I recently upgraded to R 1.9.0 in my computer at work and at home: 1.-The computer at home has Windows XP and Office XP and it seems to work perfectly and I copy-paste graphics perfectly. 2.-The computer at work has Microsoft Windows 2000 (5.00.2195 Service Pack 2) and Word 2000 (9.0.4402 SR-1) I cannot copy-paste windows metafiles into Office applications in my computer at

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

Dear R-helpers, I need some help to produce a set of contour plots that I am trying to make in order to compare surfaces between the levels of a factor. For example: library(lattice) g <- expand.grid(x = 60:100, y = 1:25, ti = c("a","b","c")) g$z <-with(g, (-1e-4*x-1e-3*y-1e-5*x*y)*(ti=="a") + (1e-2*x-1e-3*y-1e-4*x*y)*(ti=="b") +

Hi there, I don't have much experience with fitting at all and I'd like to get some advice how to use the "weights"-argument with nls correctly. I have created some data with a sigmoidal curve shape. Each y-Value was generated by the mean of three values. A standard deviation was calculated too. Now, I'd like to weight the data points respective to its standard

good morning to the R-help list!, I would like to ask if there is an alternative to the S ploting parameter "cxy" in R. I am trying to port some S code to R that makes use of this "cxy" parameter to scale the graphics, and I don't find this parameter in help(par). thank you, IOsu ------------------------------------------------------------------------------- Jes?s

The following warning message occurs when running the nls on some data: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Warning message: In is.na(wts) : is.na() applied to non-(list or vector) of type 'NULL' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ and I would like to know what is causing it and how I can fix it. As an example, the following, from Venables and Ripley, which

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 3706 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010913/b4135232/attachment.bin