similar to: Computing an ordering on subsets of a data frame

Full_Name: Lixin Han Version: 2.4.0 OS: Windows 2000 Submission from: (NULL) (155.94.110.222) A character vector c('a','b') is supplied to rm(). As a result, 'c' is deleted unintentionally. > a <- 1:5 > b <- 'abc' > c <- letters > ls() [1] "a" "b" "c" > rm(c('a','b')) > ls() character(0)

Dears Sirs During my computational work I encountered unexpected behavior when calling "ar" function, namely # time series x<-ts(c(-0.2052083,-0.3764986,-0.3762448,0.3740089,0.2737568,2.8235722,- 1.7783313,0.2728676,-0.3273164),start=c(1978,3),frequency=4,end=c(1980,3)) # ar function res.ar<-ar(x,aic=TRUE,demean=F) # call "ar" again and ............

Dear R gurus, On the CRAN website, it says that a 64bit version for Mac OS Tiger would be release shortly. Do we know what are the expected dates? Will the packages be also compiled for 64bit? We are running large microarray analysis and we keep hitting the 3Gb memory limit. I saw that there is a version available on the development mirrors, but I am not too excited to replace our very stable

Running R.app on Mac OS X 10.4 > version _ platform powerpc-apple-darwin8.6.0 arch powerpc os darwin8.6.0 system powerpc, darwin8.6.0 status major 2 minor 3.1 year 2006

Running R.app on Mac OS X 10.4 > version _ platform powerpc-apple-darwin8.6.0 arch powerpc os darwin8.6.0 system powerpc, darwin8.6.0 status major 2 minor 3.1 year 2006

Hi all, Can anyone explain why the following use of the subset() function produces a different outcome than the use of the "[" extractor? The subset() function as used in density(subset(mydf, ht >= 150.0 & wt <= 150.0, select = c(age))) appears to me from documentation to be equivalent to density(mydf[mydf$ht >= 150.0 & mydf$wt <= 150.0, "age"])

Hi all, What are current methods people use in R to identify mis-spelled column names when selecting columns from a data frame? Alice Johnson recently tackled this issue (see [BioC] posting below). Due to a mis-spelled column name ("FileName" instead of "Filename") which produced no warning, Alice spent a fair amount of time tracking down this bug. With my fumbling fingers

Full_Name: Mike Ulrich Version: 2.9 OS: Mac OSX Submission from: (NULL) (69.169.178.34) The help documentation for loadings() lists more then one parameter. The function call only expects one parameter. The digits, cutoff, and sort parameters are not used in the function. ## S3 method for class 'loadings': print(x, digits = 3, cutoff = 0.1, sort = FALSE, ...) ## S3 method for class

I all I am trying to combine columns from two dataframes to make a completely new dataframe consisting of one column of dates (ie a combination of dates from the two dataframes). >From the following dataframes a b 1 2008-07-27 1 2 2008-10-01 2 3 2008-08-15 3 4 2008-08-14 4 5 2008-08-14 5 6 2008-09-20 6 c d 1 2008-07-27 1 2 2008-10-01 2 I would like to get: z

Hi useRs, A recent coding infelicity along these lines yielded a corrupt data frame. foo <- matrix(1:12, nrow = 3) bar <- data.frame(foo) bar$NewCol <- foo[foo[, 1] == 4, 4] bar lapply(bar, length) > foo <- matrix(1:12, nrow = 3) > bar <- data.frame(foo) > bar$NewCol <- foo[foo[, 1] == 4, 4] > bar X1 X2 X3 X4 NewCol 1 1 4 7 10 <NA> 2 2 5 8 11

Dear list, I have a list of number sequences. Each number sequence has different numbers of elements. Is there a quick way (other than to iterate through the entire list) way to coerce list to matrix with NAs filling in the short sequences? An example of what I mean is this: A <- list(c(3,2,3),c(6,5)) I'd like to get A so that it is 3 2 3 6 5 NA Best, Ken

I recently received a Sigmaplot file (*.jnb) from a customer and would like to know if I can input it to a data frame and then manipulate the data in R.    I did a search on Google and on RSeek (www.rseek.org), but did not get any good hits.  Thank for any feedback and insight you can provide.   P.S. Love the flexibility of R and would love to keep using it.   Just wanting to know if this is

I'm trying to redirect where temporary files go under R to D:\temp\somerandomname from its default (C:\Documents and Settings\username\Temp\somerandomname) -- how do I go about doing this? --j -- Jonathan A. Greenberg, PhD Postdoctoral Scholar Center for Spatial Technologies and Remote Sensing (CSTARS) University of California, Davis One Shields Avenue The Barn, Room 250N Davis, CA 95616

On 04/08/2009 7:33 PM, Steven McKinney wrote: >> -----Original Message----- >> From: Duncan Murdoch [mailto:murdoch at stats.uwo.ca] >> Sent: Tuesday, August 04, 2009 8:03 AM >> To: Steven McKinney >> Cc: r-devel at stat.math.ethz.ch; R-bugs at r-project.org >> Subject: Re: [Rd] Wishlist: Navigate to "Index" page of help when no >> topic

Hi all, Hi Gavin, thx for your help i finally found out what i want to do and how to fix it. just needed to get some more level my cut level was too small... two question remain... a) can i somehow scale the twigs after cutting ? b) how can i label the nodes and how to label which one... thx !! -m.

I would like to figure out the equations for calculating "AIC" in both "step() function" and "AIC () function". They are different. Then I just type "step" in the R console, and found the "AIC" used in "step() function" is "extractAIC". I went to the R help, and found: "The criterion used is AIC = - 2*log L + k *

Does anyone know if the lmer function of lme4 works fine for unbalanced designs? I have the examination results of 1000 pupils on three subjects, one score every term. So, I have three scores for English (one for every term), three scores for maths etc. However, not everybody was examined in maths, not everybody was examined in English etc, but everybody was in effect examined on four subjects. I

I'd like to know why I cannot get a plot and the QQnorm in the same sheet. The commands are simple but: library(nlme) glmod1 <- gls(upfmla,correlation=corAR1(),method="ML") summary(glmod1) par(mfrow = c(2,1)) plot(glmod1, main="GLS Residuals vs. GLS Fitted") qqnorm(glmod1) No matter what (I tried different permutations of the plotting commands) the second drawing

I noticed taht R cannot understand certain Fortran real constant formats. For instance: c14 <- as.double( 7.785205408500864D-02) Error: unexpected symbol in " c14 <- as.double( 7.785205408500864D" The above "D" is used in Fortran language to indicate the memory starage mode. That is for instructing Fortran compiler to store such a REAL constant in DOUBLE

I am trying to replicate the repeated measures example from Dr.Faraway's book (Extending the linear model with R) as follows: data(vision) vision$npower <- rep(1:4,14) mmod <-lmer(acuity~power+(1|subject)+(1|subject:eye),vision) When I look at the fixed effects p-value, it is missing. Am I missing something here? Thanks. Fixed effects: Estimate Std. Error t value