similar to: plot via xyplot not being saved

Hi everyone, I have 2 environments (2 different R sessions) as described below: Session 1: Name of the environment: "CrlmmInfo" Objects in the environment: index1: logical index - length 238304 index2: logical index - length 238304 priors: list of 4 - (matrix 6x6, 2 vectors of length 6, vector of length 2) - all num params: list of 4: centers [238304 x 3 x

Hi everyone, I have 2 environments (2 different R sessions) as described below: Session 1: Name of the environment: "CrlmmInfo" Objects in the environment: index1: logical index - length 238304 index2: logical index - length 238304 priors: list of 4 - (matrix 6x6, 2 vectors of length 6, vector of length 2) - all num params: list of 4: centers [238304 x 3 x

Hi, Until now, I thought that the results of glm() and bigglm() would coincide. Probably a naive assumption? Anyways, I've been using bigglm() on some datasets I have available. One of the sets has >15M observations. I have 3 continuous predictors (A, B, C) and a binary outcome (Y). And tried the following: m1 <- bigglm(Y~A+B+C, family=binomial(), data=dataset1, chunksize=10e6)

Hi everyone, say I have a function called 'foo', which takes the argument arg1. Is there any mechanism that I can use to "learn" about the variable where foo(arg1) is going to be stored? For example: x <- foo(arg1) so, inside foo() I'd like to be able to get the string "x". if, foo(arg1) was used insted, I'd like to get NA. thank you very much, b

Hi everyone, Firstly, I'm sorry for the cross-post. I re-read the posting guide and it appears my question is more related to r-devel.... (i think :-) ) I downloaded the source code available at: http://cran.fhcrc.org/src/base/R-2/R-2.4.0.tar.gz to a linux machine (Linux 2.6.11-1.1369_FC4smp #1 SMP). I successfully configured and compiled it, which means I'm able to run it from the

I'm bringing this issue to r-devel as my idea of sending it to r-help appears to be wrong. The following is reproducible in R-patched and R-devel (also in older versions). An outlier is inserted in group E and the last plot by plot.lm suggests that the point belongs to D. It appears to me that: * (facval, rsp) and * x$xlevels[[1L]][order(sapply(split(yh,mf[,1]), mean))] are sorted

Hi, I'm observing the following warning with R-2.10.0 on a machine I have remote access to. Like what happened to Roger ( http://www.mail-archive.com/r-devel at r-project.org/msg08669.html ), two extra instances of X11 are open after dev.off(). Although I'm reporting this on R-2.10.0, it's reproducible on R-devel r50541. Any thoughts about this? Thanks a lot for any info, b >

Hi, first of, my apologies for the possibly off-topic message regarding a problem when installing the XML package (I meant to post at omegahat- help, but I get 404 at http://www.omegahat.org/mailman/listinfo) When installing the XML package using: install.packages("XML") everything seems to go fine... but then I try to load the package and I get: > library(XML) Error in

Hi everyone, sorry if this was discussed before (and in this situation, could you please point me to the discussion in the archive? My search didn't seem to be effective). Is there a way of getting the names of objects in a .rda file without having to load it? Thank you very much, benilton -- PhD Candidate Department of Biostatistics Bloomberg School of Public Health Johns Hopkins

Hi Everyone, I want to generate some random numbers according to some empirical distribution. Therefore I am looking for the inverse of an empirical cumulative distribution function. I haven't found any in R. Can anyone give a pointer? Thanks, Geoffrey _______________________________________________________=0A= =0A= =0A= The information in this email or in any file attached

I think that you misunderstood me. As far as I know, RMA does three things: background correction, quantile normalization, and summary from probes to probesets. I want the probe values after background correction and quantile normalization but before the summary. On Sat, Dec 26, 2009 at 12:07 PM, Benilton Carvalho <bcarvalh at jhsph.edu> wrote: > pm(data) > > b > > On Dec

I want to remove the column names from a data frame. I do it the long way, can any body show me a better way ? df= data.frame(chrN= c(“chr1”, “chr2”, “chr3”), start= c(1, 2, 3), end= c(4, 5, 6), score= c(7, 8, 9)) df #I write a txt file without row or column names write.table(df,"df1.txt",sep='\t',quote=FALSE,row.names=F,col.names=F) #then I read it with the header = F

For some time I have been using rsync to backup data between different machines and my NAS drive. But within the last 2 days I have found out that rsync is not removing old files from destination - directories renamed or removed from source. In this example I want to sync files on my Debian Linux server to a NAS drive. NAS drive is mounted using smbmount. Using "-avhzn" dry-run it

Hi, I'm trying the brlr function in a penalized logistic regression function. However, I am not sure why I am encountering errors. I hope to seek your advice here. (output below) Thank you! Your help is truly appreciated. Min-Han #No error here, the glm seems to work fine >

I find one row in my large dataset. But when I use the "rownames" for the data on the 100,000 row, the result show nothing. I try it by the following example, it still likes that. > tst[rownames(tst)==100000,] [1] x y <0 rows> (or 0-length row.names) ############## > tst=data.frame(x=1:200000,y=200000:1) > tst[rownames(tst)==1,] x y 1 1 200000 >

Simplest of template files, tst.rb, contains: PROJECT_NAME = File.basename(RAILS_ROOT) run ''echo '' + PROJECT_NAME Running command "rails tst -m ~/Rails_Tools/tst.rb" produces: create create app/controllers . . create log/development.log create log/test.log applying template: /Users/rick/RailsTools/tst.rb Anonymous modules

I am moving this from r-help to r-devel. Based on offline communications with Jim, suppose dat is defined as follows: set.seed(123) dat <- data.frame(ID= c(rep(1,2),rep(2,3), rep(3,3), rep(4,4), rep(5,5)), var1 =rnorm(17, 35,2), var2=runif(17,0,1)) # Then this ave call works as expected: ave(dat$ID, dat$ID, FUN = function(x) seq_along(x)) # but this apparently identical calculation

Hi, All: Is there a simple way to detect complex conjugates in the roots returned by 'polyroot'? The obvious comparison of each root with the complex conjugate of the next sometimes produces roundoff error, and I don't know how to bound its magnitude: (tst <- polyroot(c(1, -.6, .4))) tst[-1]-Conj(tst[-2]) [1] 3.108624e-15+2.22045e-16i

centOs5.5/samba4/named here is a short guide setting it up to work. First of all do not install the bind package coming with centos 5.5!! Install needs for samba yum install libacl* gnutls* readline* python* gdb* autoconf* Named installation: Here is a description on what to do: http://jason.roysdon.net/2009/10/16/building-bind-9-6-on-rhel5-centos5-for-d nssec-nsec3-support/ The steps, yum

Dear useRs, How come the first attempt to sort a POSIXt vector fails (Error: non-atomic type in greater), while the second succeeds? (Code inserted below.) The documentation says that POSIXt is used to allow operations such as subtraction, so I'd expect sorting to work. Is this perhaps an OS issue? (I run R 2.0.1 on Win xp.) Thank you, b. #------------code test <- c("2005-02-08