[R] Solaris: --with-internal-tzcode

2014-11-09 Thread Vasantha Kumar Kesavan
Hi, what is this parameter on Solaris? when I am executing the below command, it is failed with the error message. configure ‘--with-internal-tzcode’ configure: WARNING: unrecognized options: --with-internal-tzcode Version output: platform x86_64-pc-solaris2.10 arch x86_64

Re: [R] Lme4 Package Help!

2014-11-09 Thread Søren Højsgaard
Actually, I believe that for a brief period (around 2006) lmer did return p-values (chisq-based), but then 2006 is a long time ago... If the question pertains p-values in the usual output of coef( summary( model ) ), I'll just mention that with doBy and pbkrtest you can do linest( model,

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread Daniel Nordlund
On 11/8/2014 8:40 PM, Muhuri, Pradip (SAMHSA/CBHSQ) wrote: Hello, The example data frame in the reproducible code below has 5 columns (1 column for id and 4 columns for dates), and there are 7 observations. I would like to insert the most recent date from those 4 date columns into a new

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Hi Dan, Thank you so much for sending me your code that provides me desired results. But, I don't understand why I am getting the follow warning message, In FUN(newX[, i], ...) : no non-missing arguments, returning NA. Any thoughts? Regards, Pradip data2x - within(data1, oidflag -

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread arun
You could try library(dplyr) data1 %% rowwise() %% mutate(oldflag=as.Date(max(mrjdate,cocdate, inhdate, haldate, na.rm=TRUE), origin='1970-01-01')) Source: local data frame [7 x 6] Groups: by row idmrjdatecocdateinhdatehaldate

[R] Compression variable in memory with gzcon

2014-11-09 Thread Wush Wu
Hello, I cannot figure out how to compress a raw vector in memory with gzip format that is compatible with `gzfile`. (`memCompress` does not compatible with `gzfile`.) I modified the example in `rawConnection`, but it does not work: ``` zz - gzcon(rawConnection(raw(0), r+b)) # start with empty

Re: [R] Solaris: --with-internal-tzcode

2014-11-09 Thread Prof Brian Ripley
On 09/11/2014 08:50, Vasantha Kumar Kesavan wrote: Hi, what is this parameter on Solaris? You are using a current option (for any platform, not just for Solaris) on a long-obsolete version of R. Please do study the posting guide and update (to R 3.1.2) *before posting*. (And also send

Re: [R] Distance matrix for Spatial Filtering {spdep}

2014-11-09 Thread Roger Bivand
Erica Cseko Nolasco ecnolasco at gmail.com writes: Hi all, I´m working on a species distribution modeling, and I want to build eigenvectors to represent my spatial varuable. I´m trying to use SpatialFiltering(spdep) for it, but I´m having problems to create the nb object. My weights

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Dear Arun, Thank you so much for sending me the dplyr/mutate() solution to my code. But, I am getting the following warning message. Any suggestions on how to avoid this message? Pradip Warning message: In max(13081, NA_real_, NA_real_, 15282, na.rm = TRUE) : no non-missing arguments

[R] Is it possible to define another kind of NA

2014-11-09 Thread Marc Girondot
Dear member list, In many experimental sciences, there is a lower detection limit (LDL) when a dosage of a product is done. Then some samples are evaluated to be below this limit. I search for the best way to indicate in a data.frame that some values are such LDL. Ideally, an equivalent of NA

Re: [R] Is it possible to define another kind of NA

2014-11-09 Thread Lee, Chel Hee
LDL - NA_real_ is.LDL - is.na v - c(0.2, 0.28, LDL, 0.9) v [1] 0.20 0.28 NA 0.90 is.LDL(v) [1] FALSE FALSE TRUE FALSE Hope this helps. Chel Hee Lee On 11/9/2014 12:07 PM, Marc Girondot wrote: Dear member list, In many experimental sciences, there is a lower detection limit (LDL)

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread arun
Dear Pradip, From the documentation of ?max: The minimum and maximum of a numeric empty set are ‘+Inf’ and ‘-Inf’ One of the rows in your dataset is all `NAs.` I am not sure you want to keep that row with all NAs. You could remove it and run the code or keep it and run with

Re: [R] Is it possible to define another kind of NA

2014-11-09 Thread Dustin Tran
Hi Marc, This may be a helpful link: http://stackoverflow.com/questions/16074384/specify-different-types-of-missing-values-nas On Nov 9, 2014, at 1:07 PM, Marc Girondot marc_...@yahoo.fr wrote: Dear member list, In many experimental sciences, there is a lower detection limit (LDL) when a

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread Daniel Nordlund
On 11/9/2014 3:05 AM, Muhuri, Pradip (SAMHSA/CBHSQ) wrote: Hi Dan, Thank you so much for sending me your code that provides me desired results. But, I don't understand why I am getting the follow warning message, In FUN(newX[, i], ...) : no non-missing arguments, returning NA. Any thoughts?

Re: [R] Getting the most recent dates in a new column from dates in four columns using the dplyr package (mutate verb)

2014-11-09 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Hi Arun and Dennis, This is just an FYI. You're right - In one row, there are all NA's in the four date columns. I have tested below the TRUEness of the condition Arun has set. is.logical(data1[rowSums(is.na(data1[,-1]))!=4,]) [1] FALSE All these 3 approaches below provide the exact same

[R] Working with data-frame

2014-11-09 Thread Christofer Bogaso
Hi again, Let say, I have following data frame: Dat - structure(list(A1 = structure(c(3L, 3L, 1L, 3L, 3L, 3L, 3L, 2L, 3L, 3L, 1L, 2L, 3L, 2L, 1L, 1L, 3L, 3L, 2L, 3L, 2L, 2L, 3L, 3L, 3L, 2L, 3L, 1L, 1L, 3L), .Label = c(a, b, c), class = factor), A2 = c(2, 3, 2, 1, 3, 3, 2, 2, 3, 1, 3, 1, 3,

Re: [R] Working with data-frame

2014-11-09 Thread Bert Gunter
Christopher: If I understand correctly, see ?ave, or ?tapply, depending on what form you want to be returned. The trick is to first paste the columns together on whose unique you want to split to form a singtl factor. e.g. of the form lapply(split(yourcolumn,paste0(...)),FUN= sum) However, the

Re: [R] Is it possible to define another kind of NA

2014-11-09 Thread Bert Gunter
Ouch! The values are **NOT** missing -- they are (left) censored, and need to be handled by appropriate censored data methods. I suggest you (all!) either read up on this or consult someone locally who has knowledge of such methods. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650)

Re: [R] Working with data-frame

2014-11-09 Thread William Dunlap
I tried with spilt() function. However it looks to me that, it can split a data-frame w.r.t. only one column. (I assume you you meant 'split', not 'spilt'.) You did not show what you tried, but the following splits Dat by its A1 and A2 columns (creating a list of data.frames): split(Dat,

[R] CRAN (and crantastic) updates this week

2014-11-09 Thread Crantastic
CRAN (and crantastic) updates this week New packages * anim.plots (0.1) Maintainer: David Hugh-Jones Author(s): David Hugh-Jones davidhughjo...@gmail.com License: GPL-2 http://crantastic.org/packages/anim-plots Simple animated versions of basic R plots, using the

Re: [R] Working with data-frame

2014-11-09 Thread Jeff Newmiller
... so... #1 ... flexible syntax for split-apply-combine, not very efficient for large data library(plyr) ddply(Dat,c(A1, A2), function(DF){data.frame(C1=sum(DF$C1))}) #2 ... compatible with large data on disk library(sqldf) sqldf(select A1,A2,sum(C1) as C1 from Dat group by A1, A2) #3 ...

[R] the bug of function base::order

2014-11-09 Thread 岳�S
Hi, all: I find a bug of the function base::order. For example, y = rep(9, 9:1); rbind(y,order(y, decreasing=FALSE)); the result is: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] y998765432 1 1098765431