Megh Dal
2008-Feb-27 10:24 UTC
[R] Calculating monthly var-cov matrix on non-overlapping rooling window basis
let create a 'zoo' object : library(zoo) date.data = seq(as.Date("01/01/01", format = "%m/%d/%y"), as.Date("06/25/02", format = "%m/%d/%y"), by = 1) len = length(date.data) data1 = zoo(matrix(rnorm(2*len), nrow = len), date.data ) head(data1) Now I want to create an 3 dimensional array (suppose name " var.cov") where, var.cov[,,i] gives the Variance-covariance matrix for i-th month of data1. That is I want to calculate monthly variance-covariance matrix on non-overlapping rolling window basis. Any suggestion? --------------------------------- [[alternative HTML version deleted]]
Henrique Dallazuanna
2008-Feb-27 11:15 UTC
[R] Calculating monthly var-cov matrix on non-overlapping rooling window basis
Perhaps something like this: lapply(split(data1, format(index(data1), "%m")), cov) On 27/02/2008, Megh Dal <megh700004 at yahoo.com> wrote:> let create a 'zoo' object : > > library(zoo) > date.data = seq(as.Date("01/01/01", format = "%m/%d/%y"), as.Date("06/25/02", format = "%m/%d/%y"), by = 1) > len = length(date.data) > data1 = zoo(matrix(rnorm(2*len), nrow = len), date.data ) > head(data1) > > Now I want to create an 3 dimensional array (suppose name " var.cov") where, var.cov[,,i] gives the Variance-covariance matrix for i-th month of data1. That is I want to calculate monthly variance-covariance matrix on non-overlapping rolling window basis. > > Any suggestion? > > > > --------------------------------- > > [[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. >-- Henrique Dallazuanna Curitiba-Paran?-Brasil 25? 25' 40" S 49? 16' 22" O