# HW 3 annflow.R JLCY, 19 Dec 2013

############################################################ 

# boxplots of the Annual PDF....  Obtain the PDF of the observed Annual
# flows

{

    zz = sm.density(Annual, eval.points = annxeval, display = "none")
    xdensityorig = zz$estimate

    par(mfrow = c(1, 1))
    xs = 1:length(annxeval)
    zz = boxplot(split(t(annsimpdf), xs), plot = F, cex = 1)
    zz$names = rep("", length(zz$names))
    z1 = bxp(zz, ylim = range(annsimpdf, xdensityorig), xlab = "Annual flow MAF", 
        ylab = "PDF", cex = 1.25)

    z2 = 1:6
    n1 = 1:6
    z2[1] = z1[1]
    z2[2] = z1[20]
    z2[3] = z1[40]
    z2[4] = z1[60]
    z2[5] = z1[80]
    z2[6] = z1[100]

    n1[1] = xeval[1]
    n1[2] = xeval[20]
    n1[3] = xeval[40]
    n1[4] = xeval[60]
    n1[5] = xeval[80]
    n1[6] = xeval[100]
    n1 = round(n1, dig = 0)
    n1 = as.character(n1)

    axis(1, at = z2, labels = n1, cex = 1)
    lines(z1, xdensityorig, lty = 2, lwd = 2, col = "red")

}

## Error: 沒有這個函數 "sm.density"