PDF: Probability density function and histogram
In brasmus/preciparea: Precipitation area ananlysis.

Description Usage Arguments Author(s) See Also Examples

A graph showing the marginal distribution of the time series

1	PDF(x, is = 1, plot=TRUE, ...)

`x`	a time series
`is`	index for which series to show
`plot`	true: plot the results
`...`	arguments passed on to `plot`

Rasmus.Benestad@met.no

AC, timeseries

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, is = 1, ...) 
{
    z <- coredata(trend(x[, is], result = "residual"))
    hist(z, col = "grey", freq = FALSE, main = "", ylab = "", 
        xlab = paste(attr(x, "longname")[is], " (", attr(x, "unit")[is], 
            ")", sep = ""), ...)
    x <- seq(-max(abs(z)), max(abs(x)), length = 100)
    lines(x, dnorm(x, mean = mean(z), sd = sd(z)), lwd = 3, col = rgb(1, 
        0, 0, 0.4))
    grid()
    par(new = TRUE, fig = c(0.75, 0.98, 0.75, 0.98), mar = rep(0, 
        4), cex.axis = 0.7)
    qqnorm(z, main = "")
    qqline(z, col = "red")
  }