Plot an rpart object with all splits labeled with p-values and all nodes labeled with confidence intervals.

Essentially a wrapper function for rpart.plot() from the rpart.plot package, but additional arguments allow user to add p-values and confidence intervals to plots.

treeval.plot(
  tree,
  sigma_y = NULL,
  nn = TRUE,
  printn = TRUE,
  inferenceType = 2,
  digits = 3,
  alpha = 0.05,
  permute = FALSE,
  ...
)

Arguments

tree	An rpart tree. The tree must have been build with parameter model=TRUE.
sigma_y	Provide the standard deviation of y, if known. If not provided, the sample standard deviation of y will be used as a conservative estimate.
nn	boolean- would you like node numbers to be printed? Nodes are numbered using the same methodology as the rpart package. If node n has children, its children are numbered 2n and 2n+1.
printn	boolean - would you like the number of observations to be printed in each node?
inferenceType	An integer specifying which pieces of inference information should be added to the plot. The options currently available are (0) No confidence intervals, p-values, or "fitted mean" label. Just calls rpart.plot(). (1) No confidence intervals. Each split labeled with a p-value. (2) Label each internal node with a confidence interval and label each split with a p-value. This is the default, but can also be a little messy/hard to read. Options 3 and 4 print the same information but with small formatting tweaks.
digits	Integer- how many digits would you like the text in the plot rounded to.
alpha	If inferenceType is such that confidence intervals will be printed, (1-alpha) confidence intervals will be printed.
permute	If inferenceType is such that confidence intervals will be printed, should the conditioning set for the confidence intervals include all permutations of the relevant branch? Setting this to TRUE will lead to slightly narrower confidence intervals, but will make computations more expensive. See paper for more details.
...	Additional arguments are passed on to rpart.plot(). Examples include "cex".

Examples

bls.tree <-rpart::rpart(kcal24h0~hunger+disinhibition+resteating,
    model = TRUE, data = blsdata, maxdepth=1)
treeval.plot(bls.tree, inferenceType=0)
treeval.plot(bls.tree, inferenceType=1)
treeval.plot(bls.tree, inferenceType=2)