Visualize {GaSP}

Visualize a GaSPModel object.

Description

Carry out a functional analysis of variance (ANOVA) of a GaSPModel object and generate plotting coordinates for its estimated main and 2-input joint effects.

Usage

Visualize(GaSP_model, x_description, main_percent = 0, interaction_percent = 0)

Arguments

GaSP_model

Object of class GaSPModel.

x_description

A data frame describing the input variables. See DescribeX.

main_percent

An optional minimum percentage of variation explained by an input's main effect to return the effect's plotting coordinates; the default of zero gives plotting coordinates for all inputs.

interaction_percent

An optional minimum percentage of variation explained by the interaction effect of a pair of inputs to return the plotting coordinates for their joint effect (main effects plus interaction effect); the default of zero gives plotting coordinates for all pairs of inputs.

Details

If there are many inputs, to avoid excessive plotting of many trivial joint effects set interaction_percent = 1 say.

Value

A list with the following elements:

anova_percent

A data frame containing the ANOVA percentages for all main effects and 2-input interaction effects.

main_effect

A data frame with plotting coordinates for the estimated main effects.

joint_effect

A data frame with plotting coordinates for the estimated 2-input joint effects.

total_percent

Total percentage of the prediction variation accounted for by all main effects and 2-input interaction effects.

average

Overall average of the prediction function, averaged with respect to all inputs.

SE_average

Standard error of the overall average.

References

Schonlau, M. and Welch, W.J. (2006), "Screening the Input Variables to a Computer Model Via Analysis of Variance and Visualization", in Screening: Methods for Experimentation in Industry, Drug Discovery, and Genetics, Dean. A. and Lewis, S., eds., pp. 308-327, Springer, New York, doi:10.1007/0-387-28014-6_14.

Examples

borehole_vis <- Visualize(borehole_fit, borehole_x_desc)