# HW 3 mayflow.R JLCY, 19 Dec 2013

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

# plot mayflow

{

    par(mfrow = c(1, 1))

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

    par(mfrow = c(1, 1))

    xs = 1:length(xeval)
    zz = boxplot(split(t(simpdf), xs), plot = F, cex = 1)
    zz$names = rep("", length(zz$names))
    z1 = bxp(zz, ylim = range(simpdf, xdensityorig), xlab = "May 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)

    # original May PDF as red line
    lines(z1, xdensityorig, lty = 2, lwd = 2, col = "red")

}

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