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
... 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 ..
CRAN (and crantastic) updates this week
New packages
* anim.plots (0.1)
Maintainer: David Hugh-Jones
Author(s): David Hugh-Jones
License: GPL-2
http://crantastic.org/packages/anim-plots
Simple animated versions of basic R plots, using the 'animation'
package. Includes a
> 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(Da
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) 4
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
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
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 sam
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?
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 wrote:
>
> Dear member list,
>
> In many experimental sciences, there is a lower detection limit (LDL) when a
> dosage of
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
> 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 (
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
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 to
Erica Cseko Nolasco 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 would be
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 prop
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 em
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:
idmrjdatecocdateinhdatehaldateo
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 <- apply(
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 col
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, diag(1
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
o
22 matches
Mail list logo
