data <- array(rnorm(45), dim = c(member = 7, sdate = 4, time = 3))
Dates <- c(seq(as.Date("2000-11-01", "%Y-%m-%d", tz = "UTC"),
as.Date("2001-01-01", "%Y-%m-%d", tz = "UTC"), by = "month"),
seq(as.Date("2001-11-01", "%Y-%m-%d", tz = "UTC"),
as.Date("2002-01-01", "%Y-%m-%d", tz = "UTC"), by = "month"),
seq(as.Date("2002-11-01", "%Y-%m-%d", tz = "UTC"),
as.Date("2003-01-01", "%Y-%m-%d", tz = "UTC"), by = "month"),
seq(as.Date("2003-11-01", "%Y-%m-%d", tz = "UTC"),
as.Date("2004-01-01", "%Y-%m-%d", tz = "UTC"), by = "month"))
dim(Dates) <- c(sdate = 4, time = 3)
res <- PeriodVariance(data, dates = Dates, start = list(01, 12), end = list(01, 01))
Run the code above in your browser using DataLab