Re: [R] count NAs with aggregate

2013-03-25 Thread Rui Barradas
Hello, The problem seems to be with the formula interface. # Make up some data var1 <- rnorm(100) var2 <- rnorm(100) var3 <- rnorm(100) subject <- sample(4, 100, replace = TRUE) time <- sample(10, 100, replace = TRUE) var1[sample(100, 10)] <- NA var2[sample(100, 10)] <- NA var3[sample(100, 10)]

[R] count NAs with aggregate

2013-03-25 Thread Timo Stolz
Dear members of this list, I'd like to count missing values using the aggregate function. Something like this: count_nas <- function(arg1) { return(sum(is.na(arg1))) } aggregate(cbind(var1, var2, var3) ~ subject + time, data = mydataset, count_nas) It's not working: I end up with a matri