similar to: Ref: Variable scope or function behaviour or array reassign

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 >

Hi thank you for talking the time to help me with this. I have a sequence of numbers in a file and an equal sequence of various character, say(a b c d) each occurs more than once. I need to plot the numbers so that numbers corresponding to a in the other sequence would have green dots, those corresponding to b a red dot, nothing on c and blue square for d. i.e 2 a show a green dot 4 b show a

Dear All, I'd like to know if it is possible to create animations with R. To be specific, I attach a code I am using for my research to plot some analytical results in 3D using the lattice package. It is not necessary to go through the code. Simply, it plots some 3D density profiles at two different times selected by the user. I wonder if it is possible to use the data generated for different

Hi, I have a data.frame that I need to construct iteratively. At the moment, I'm doing: d<-data.frame(x=c(),y=c(),z=()); # {and, within some loop} d<-rbind(d,data.frame(x=newx,y=newy,z=newz); While this works, it is horribly verbose and probably not efficient, either. My real data.frame has, of course, many more columns, which can be of different modes. I vaguely recall that

Good afternoon, Hope you all have a wonderful day. I am glad to be here. Hope you could help me with the following errors that i have been trying to figure it all out since last week. I am using Splus from Insightful, and as i read, R and Splus are very similar. So hope you could help me. I have been continously received these error messages after i ran my small program for couple of

Thanks for the reproducable example. I can confirm that it fails on my machine using survival 2-37.5, the next soon-to-be-released version, The issue is with NextMethod, and my assumption that the called routine inherited everything from the parent, including the environment chain. A simple test this AM showed me that the assumption is false. It might have been true for Splus. Working this

I'm trying to convert an S-Plus program to R.  Since I'm a SAS programmer I'm not facile is either S-Plus or R, so I need some help.  All I did was convert the underscores in S-Plus to the assignment operator <-.  Here are the first few lines of the S-Plus file:   sshc _ function(rc, nc, d, method, alpha=0.05, power=0.8,              tol=0.01, tol1=.0001, tol2=.005, cc=c(.1,2),

Hello R-users, I'm currently facing a pretty hard problem which I'm hopping you'll be able to help me with. I'm using R to create images. That alone is not the problem, the problem is that I'm using R to create 168 000 images... My code (which is given below) use different package (raster and rgdal) to import a image (size 20gig) and divide it into 168 000 pictures that are

Hi, I am new to R. and even though I've made my code to run and do what it needs to . It is taking forever and I can't use it like this. I was wondering if you could help me find ways to fix the code to run faster. Here are my codes.. the data set is a bunch of 0s and 1s in a data.frame. What I am doing is this. I pick a column and make up a new column Y with values associated with that

Hi, I am trying to simulate animal movement in a gridded landscape made up of cells. At each time step (iteration), the animal moves from one cell to another in a random fashion. This is how I am simulating movement, where a and b are the x,y co-ordinates of the animal at the previous time step: for (i in 1:no.of.steps){ direction <- sample(1:8, 1) if(direction == 1){ a <- a b <- b -

Is there a way to make a shading interpolation on an image plot? Something similar to matlab 'shading interp', I think it is called Gouraud shading. What I want is to make a image plot look nicer. with image() it looks very facetted, and I would like to make it look smoother. I've tried with interp.surface() in fields package but it (obviously) makes nan values at the borders and

Hallo, I have a set of S4 and S3 classes together in one script. While running this script I create a lot of new functions and objects An example for S3 and S4 classes: ## S3 classes pt <- list(x=1,y=2) class(pt) <- "xypoint" xpos <- function(x, ...) UseMethod("xpos") xpos.xypoint <- function(x) x$x ypos <- function(x, ...) UseMethod("ypos")

Dear all, As you can see from the attachment I'm using R to automatically annotate peptide fragmentation mass spectra, which are represented by impulse plots. I'd like to poll you on approaches of how to deal as generally as possible with the two biggest annotation issues I run into: 1) very close annotated masses (impulses) with similar y-axis dimensions - resulting in overlapping labels

Hi - When using the randomForest function for regression, I get different results for mean-squared error of the predictions depending on whether or not I specify to calculate variable importance. There is an example below. I looked briefly at the source code, but couldn't find anything that would indicate why calculating variable importance would (or should) change predictions. I'm

Are the results dramatically different? The result would be expected to be somewhat different, as setting importance=TRUE would make many calls to the random number generator (for permuting OOB data in each variable), making all but the first tree in the forest different than if importance=FALSE. Cheers, Andy > From: Scott Gilpin > > Hi - > > When using the randomForest

Dear R users I need to auto scale the left y axis in the code below, so that when I scroll left or right the left y-axis scale changes to accumulate the range of the displayed data with in the max hight of the y-axis. also how can I make the crosshair horizontal since it is only vertical in this code. this code with a kind help from "Gregory (Greg) L. Snow Ph.D." just

Please see the code below. When I try to run the variogram - vg.deft<-variog(rd,uvec=10) I keep getting this error- variog: computing omnidirectional variogram Error in vector("double", length) : vector size specified is too large Also, when I try to define distance-based neighborhood - nb.tr=dist.neighbors(tr.locs,2) I get this error - Error in vector("double", length) :

An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/r-help/attachments/20041028/dbf887a9/attachment.pl -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: random1.txt Url: https://stat.ethz.ch/pipermail/r-help/attachments/20041028/dbf887a9/random1.txt

Hi all, I am new to R, hence this question. I have a set of points with X and Y coordinates. I would like to build an R code to calculate the distances of points within a specified neighborhood (circular range) for each point. I would like the code to be a function, so that I can call the function, specifying parameters (data file, distance), from another piece of code. Also, I would like the

Dear all, I have the problem as the title shows. It is time series of measurements. I want to estimate the value at time point when we have no actual measured value (Only in the time range, not for prediction in future time). Hope you could give me some suggestion or hint. Thank you so much. Wang Zhipeng [[alternative HTML version deleted]]