similar to: Suggestion on how to improve efficiency when using MASS:::hubers on high-dimensional arrays

Hello: I found the "hubers" function in MASS library is NOT working on the following data: > a <-

Hi, I have a long list of lists from which I want to efficiently extract and rbind elements. So I'm using the approach below: f <- function(i){ out <- replicate(5, list(matrix(rnorm(80), nc=20))) names(out) <- letters[1:5] out } set.seed(1) lst <- lapply(1:1.5e6, f) (t0 <- system.time(tmp <- do.call(rbind, lapply(lst, '[[', 'b')))) Is there

Hi Everyone, After searching the subject and not being successful, I was wondering if any you could explain me the idea behind the following fact: all(NULL == 2) ## TRUE any(NULL == 2) ## FALSE Thanks a lot, Benilton -- Benilton Carvalho PhD Candidate Department of Biostatistics Johns Hopkins University

Hi everybody, I've been looking for information on "how to create administrative shares", but I didn't find... Can you help me on this matter? Regards, -- Benilton Carvalho DE / IMECC / UNICAMP Red Hat Linux i18n Team

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 all, working on the problem that I cannot log in from my NT4 box I noticed that: - Test 8 from DIAGNOSIS.TXT fails: from the NT machine I try: "net view \\server" and I have: "System error 5 has occurred. Access denied" I don't know how to handle this... Could you please help me? Regards, -- Benilton Carvalho DE / IMECC / UNICAMP Red Hat Linux i18n Team

Hi everyone, For reasons beyond the scope of this message, I'd like to append a NULL element to the end of a list. tmp0 <- list(a=1, b=NULL, c=3) append(tmp0, c(d=4)) ## works as expected append(tmp0, c(d=NULL)) ## list with a/b/c only Given that I could use tmp0$a <- NULL to remove 'a', I seem to understand why appending NULL returns me the original list... But how should I

I need to read in csv files, created by 3rd party, with fields containing single quotes (as shown below). "header1","header2","header3","header4" "field1r1","field2r1","field3r1","field4r1" "field1r2","field2r2","field3r2PartA), field3r2PartB Very" Long","field4r2"

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

How to get robust M-estimators of multivariate scatter using Huber's psi? Which package/function should I look into? Ideally, I hope I can self-define thresholds of Huber's psi function. Thanks a lot!!! -- View this message in context: http://www.nabble.com/How-to-get-robust-M-estimator-of-multivariate-scatter-using-Huber%27s-psi--tp20585755p20585755.html Sent from the R help mailing

Hi, I have a BioC infra-structure package that works fine (I can build, check and load it successfully) on revision r53950. The very same package fails on r54591 with the error below: Error in loadNamespace(package, c(which.lib.loc, lib.loc), keep.source = keep.source) : cyclic name space dependency detected when loading ?oligoClasses?, already loading ?oligoClasses? I don't see anything

Hi, I'm working on a package and observed a behaviour that was suggested to be incompatible with the documentation. The code below: x <- matrix(1:4, 2) dimnames(x) <- list() works just fine. But checking the documentation for dimnames(), it states that if value is a list "its components are either 'NULL' or a character vector with positive length of the appropriate

Hi, First of all, I apologize for sending out this message, as I'm sure the answer is on the archives, but I just can't find it. Not long ago, there was a discussion about building Windows packages from the source code and someone posted a link to a website to which we could submit the source and get the Windows version a while later. Could someone point me to that website or to

Hi, what is the proper of of "passing a missing value" so I can extract the entire i-th row of a matrix (in a list of lists) without pre-computing the number of cols? For example, if I know that the matrices have 2 columns, I can do the following: set.seed(1) x0 <- lapply(1:10, function(i) replicate(4, list(matrix(rnorm(10), nc=2)))) lapply(lapply(x0, '[[', 3), '[',

Hi, I'm working on an S4 class that is expected to behave like an array. I have some difficulties when defining '[' and I wonder if someone could point me to the right direction: 1) Call the S4 object "obj" 2) Assume dim(obj) = c(10, 4, 2) 3) Suppose someone calls: obj[1:3,] , which is a mistake, given dim(obj); how do I detect such situations? Thank you very much for

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 Everyone, when building (say) R-2.4.0 from the source, is it still the recommendation to use GCC 3.4? Thank you very much, Benilton Carvalho PhD Candidate Department of Biostatistics Johns Hopkins University

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

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)