plotnonp: Plotting Nonparametric Function Estimates In BayesX: R Utilities Accompanying the Software Package BayesX

Description

Plots nonparametric function estimates obtained from BayesX

Usage

 ```1 2 3``` ```plotnonp(data, x = 2, y = c(3, 4, 5, 7, 8), ylim = NULL, lty = c(1, 2, 3, 2, 3), cols = rep(1, length(y)), month, year, step=12, xlab, ylab, ...) ```

Arguments

 `data` Either the name of a file or a data frame containing the estimation results. `x` Defines the x-axis in the plot. Either the name of a variable in `data` or the index of the corresponding column. `y` Defines the variables to be plotted against `x`. May be either a vector of names of variables in `data` or the corresponding indices. The default choice corresponds to the point estimate plus two confidence bands. `ylim` Since `plotnonp` plots multiple y-variables, it automatically determines the appropriate `ylim` to make all curves visible. Argument `ylim` allows to override this default behaviour with fixed values. `lty` Vector of line types used for plotting (must have the same length as `y`). The default corresponds to solid lines for the point estimate and dashed and dotted lines for the confidence bands. `cols` Vector of colors used for plotting (must have the same length as `y`). Default are black lines. `month, year, step` Provide specific annotation for plotting estimation results for temporal variables. `month` and `year` define the minimum time point whereas step specifies the type of temporal data with `step=4`, `step=2` and `step=1` corresponding to quartely, half yearly and yearly data. `xlab, ylab` `plotnonp` constructs default labels that can be overwritten by these arguments `...` Further arguments to be passed to the interval call of `plot` such as `type`, etc.

Author(s)

Felix Heinzl, Andreas Brezger and Thomas Kneib

`drawmap`,`plotautocor`,`plotsample`,`plotsurf`
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19``` ```res <- read.table(system.file("examples/nonparametric_f_x_pspline.res", package="BayesX"), header=TRUE) plotnonp(res) plotnonp(res, x="x") plotnonp(res, x="x", y="pmean") plotnonp(res, x="x", y="pmed") plotnonp(res, x="x", y="pmed", ylim=c(-2,2)) plotnonp(res, x="x", y=c("pmean", "pqu10", "pqu90"), lty=c(1,1,1), col=c("red","blue","blue")) plotnonp(res, xlab="some variable", ylab="f(some variable)", main="Nonlinear effect of some variable", sub="penalised spline") res <- read.table(system.file("examples/nonparametric2_f_time_pspline.res", package="BayesX"), header=TRUE) plotnonp(res) plotnonp(res, month=1, year=1980, step=12) res <- res[1:18,] plotnonp(res, month=1, year=1980, step=12) ```