Re: [R] Computing sd across an array with missing values

2009-02-25 Thread Matt Oliver
c(5, 5, 3)) >> >> # First >> res<-apply(p,3,function(X) >>c(scols=apply(X,2,sd,na.rm=TRUE),srows=apply(X,2,sd,na.rm=TRUE)) >>) >> res >> >> # Second >> res2<-apply(p,3,function(X) >>list(scols=apply(X,2,sd,na.rm

[R] Computing sd across an array with missing values

2009-02-25 Thread Matt Oliver
Dear help, suppose I have this array and want to compute sd aross rows and columns. p <- array(c(1:5, rep(NA, times = 3)), dim = c(5, 5, 3)) apply(p, 1:2, sd) fails because sd requires at least 2 numbers to compute sd apply(p, 1:2, sd, na.rm = TRUE) fails for the same reason I crafted my own fu

Re: [R] [R-sig-Geo] LIDAR Problem in R (THANKS for HELP)

2008-08-05 Thread Matt Oliver
you can try memory.limit(size=4000) only if you have 4GB of memory on the system This is not guaranteed to solve your problem though With big datasets like lidar, you are much better off getting access to a 64bit system with a ton of RAM (>64GB). Cheers Matt On Tue, Aug 5, 2008 at 1:47 PM,