similar to: Cleaning up after compiling R source

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 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

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

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

Thanks to those who suggested I use make clean. I *had* done that before asking the question, but it would never run. I kept getting: nothing to make. After getting a couple of replies suggesting "make clean" again, I blew away the entire directory, reinstalled the sources, reset the paths, and first did a normal "make". Then I followed that with "make

I need to completely remove row.names from a dataframe. Are there other ways to remove them (and not anything else) besides: mydataframe<-data.frame(mydataframe, row.names=NULL) I realize that this doesn't really remove the row.names; it merely replaces the current row.names vector with the numbers 1..nrow (in quotes). ===================== Dr. Marc R. Feldesman Professor and

I need to completely remove row.names from a dataframe. Are there other ways to remove them (and not anything else) besides: mydataframe<-data.frame(mydataframe, row.names=NULL) I realize that this doesn't really remove the row.names; it merely replaces the current row.names vector with the numbers 1..nrow (in quotes). ===================== Dr. Marc R. Feldesman Professor and

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 using the package rpart with R 1.3.0 for Windows to produce simple classification trees for some measurement data from paleontological specimens. Both the rpart documentation and the output confirm that the program produces splits on continuous data that leave "holes" in the data. It is probably of little practical importance, but is there a reason why the binary

I'm running a series of simple bivariate linear regressions on grouped data. I want to test the slopes to see if they are parallel. I normally use analysis of covariance to do so, looking at interaction between the covariate and the factor to make this determination. VR3 pp.149 - 154 has a very nice example of an ANOCOVA, ending with a discussion of this very operation. My question has

Sorry for the uninformative subject. The following piece of code gives different output in S-Plus (2K & 6) vs R 1.2.1 (Win) apply(apehum[,6:15], 2, function(x) which(is.na(x))) S-Plus results: > apply(apehum[,6:15], 2, function(x) which(is.na(x))) $latsupri: numeric(0) $medepico: numeric(0) $pdhtcapi: numeric(0) $mlhtcapi: numeric(0) $aphttroc: [1] 151 152 $mlhttroc: numeric(0)

I have a matrix (1000 x 224) that contains the output from the function boot.array. For each row, there are 224 columns that represent the case numbers of the selected individuals from a bootstrap. The cases are listed in selection order, not numerical order. I'd like to rearrange the cells so that every row has the column entries listed ascending numeric order. I've been able to

I'm doing kmeans partitioning on a small (n=26) dataset that has 5 variables. I noticed that if I repeatedly run the same command, the cluster centers change and the cluster membership changes. Using RW1022 under Windows NT & Windows 2000 >kmeans(pottery[,1:5], 4, 20) [...snip] $size [1] 7 3 9 7 [...snip] $size [1] 7 10 4 5 [...snip] $size [1] 6 10 5 5 yields a different

I just downloaded the binaries for the released 1.2.0 for Windows (from Guido Masarotti's site). In using the rwinst.exe installer, I noted two significant problems. First, the installer seems extraordinarily slow compared to the the 1.1.1 installer. It took almost 20 minutes to perform a full base extraction from binaries to installation (all base pages, html help, windows help, pdf

My earlier message was not as clear as it could have been. I already have Emacs for my NT system, but prefer not to use it or vi. It is really a case of not needing a sledgehammer to swat a gnat. The advantage of EditPadPro (for me anyway) is that it has a footprint of 1MB including helpfiles, whereas the versions of Emacs I've seen for NT require about 30MB of disk space to store.

I'm trying to adapt some S-Plus scripts to run in R (1.2.2, Windows). In one of these scripts, I've bootstrapped the prediction success rate under the discriminant function (lda). The bootstrap() functions are proprietary to S-Plus and there aren't exact equivalents in R. The closest is Canty's library boot based on the Davidson and Hinkley book. Unfortunately, I

I want to produce some boxplots (and qqnorm plots) for each of a series of groups. It is easy to produce the boxplots or qqplots by groups using either by() or tapply(), but the moment that I want unique labels for each boxplot, I seem to be stuck using a for loop to achieve the effect. For example: tapply(fstat$femur, fstat$race, boxplot) by(fstat$femur, fstat$race, boxplot) gives nicely

I've run into another postscript/rpart problem unrelated to the issues I've mentioned in a previous query. I'm using 1.3.1 on a Win2K box. >plotcp(some.rpart.object) draws a very nice plot on the windows graphic device. If I save this as a postscript file, either by opening a postscript device before calling plotcp, or by saving the graphics window as a postscript file, the

I'm trying to develop a systematic way of placing legends on (among other things) scatterplots of grouped data. The journals in my field can be very fussy about the placement of figure legends. Typically they want them to be in either of two places - the upper left corner of the plotting area, not obscuring any data points, or outside the plotting area in the lower right hand corner of

I'm not really sure this is the place to raise this, but here goes anyway. I've downloaded the 1.4.0 recommended package sources R-1.4.0-recommended.tgz multiple times from the main CRAN site. The sources have arrived corrupted each time (from multiple computers BTW). If I download the same file from the UCLA or University of Bristol mirror, I also get corrupt archive files. On the