Dear R users, I want to subset a daily zoo series according to its month, find % of "NA" in each month. I am finding it difficult to subset the daily dataset into monthly for the given operation.I am planning to do this for a huge dataset. Thanks in advance. Regards Vikram [[alternative HTML version deleted]]
On Mon, 23 Apr 2012, Vikram Bahure wrote:> Dear R users, > > I want to subset a daily zoo series according to its month, find % of "NA" > in each month. > > I am finding it difficult to subset the daily dataset into monthly for the > given operation.I am planning to do this for a huge dataset.Use aggregate(z, as.yearmon, ...) for aggregation to a monthly series. A simple artificial example is: ## data with NAs set.seed(1) x <- rnorm(50) x[sample(1:50, 10)] <- NA ## as zoo z <- zoo(x, as.Date("2012-01-01") + 0:49) ## proportion of NAs aggregate(z, as.yearmon, function(x) mean(is.na(x))) hth, Z> Thanks in advance. > > Regards > Vikram > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >
Hi Vikram, maybe something like this? require("zoo") ## a daily zoo series z <- zoo(1:200, as.Date("2012-01-01") + 1:200) ## add some NAs z[ceiling(runif(10)*200)] <- NA ## aggregate to monthly aggregate(z, format(index(z), "%Y-%m"), length) aggregate(z, format(index(z), "%Y-%m"), function(x) sum(is.na(x))) aggregate(z, format(index(z), "%Y-%m"), function(x) sum(is.na(x))/length(x)) Regards, Enrico Am 23.04.2012 11:27, schrieb Vikram Bahure:> Dear R users, > > I want to subset a daily zoo series according to its month, find % of > "NA" in each month. > > I am finding it difficult to subset the daily dataset into monthly > for the given operation.I am planning to do this for a huge dataset. > > Thanks in advance. > > Regards Vikram > >-- Enrico Schumann Lucerne, Switzerland http://nmof.net
Apparently Analagous Threads
- How to do operations on zoo/xts objects with Monthly and Daily periodicities
- Using zoo() to aggregate daily data to monthly means
- confused with yearmon, xts and maybe zoo
- Time series: xts/zoo object at annual (yearly) frequency
- converting a zoo or an xts to a data frame