# HW 3 sim_basicstat.R JLCY, 19 Dec 2013

# ***calculate basic statistics****************

{
    for (j in 1:12) {
        armean[k, j] = mean(simdismon[, j])
        armax[k, j] = max(simdismon[, j])
        armin[k, j] = min(simdismon[, j])
        arvar[k, j] = var(simdismon[, j])
        arskw[k, j] = skew(simdismon[, j])
    }

    # correlation between one month to another..

    for (j in 2:12) {
        j1 = j - 1
        arcor[k, j] = cor(simdismon[, j], simdismon[, j1])
    }

    arcor[k, 1] = cor(simdismon[1:nyrs1, 12], simdismon[2:nyrs, 1])
}

## Error: 找不到物件 'simdismon'