plot {RobustGaSP} | R Documentation |
Function to make plots on Robust GaSP models after the Robust GaSP model has been constructed.
## S4 method for signature 'rgasp' plot(x,y, ...)
x |
an object of class |
y |
not used. |
... |
additional arguments not implemented yet. |
Three plots: the leave-one-out fitted values versus exact values, standardized residuals and QQ plot.
Mengyang Gu [aut, cre], Jesus Palomo [aut], James Berger [aut]
Maintainer: Mengyang Gu <mgu6@jhu.edu>
M. Gu. (2016). Robust uncertainty quantification and scalable computation for computer models with massive output. Ph.D. thesis. Duke University.
library(RobustGaSP) #------------------------ # a 3 dimensional example #------------------------ # dimensional of the inputs dim_inputs <- 3 # number of the inputs num_obs <- 30 # uniform samples of design input <- matrix(runif(num_obs*dim_inputs), num_obs,dim_inputs) # Following codes use maximin Latin Hypercube Design, which is typically better than uniform # library(lhs) # input <- maximinLHS(n=num_obs, k=dim_inputs) ##maximin lhd sample # outputs from the 3 dim dettepepel.3.data function output = matrix(0,num_obs,1) for(i in 1:num_obs){ output[i]<-dettepepel.3.data (input[i,]) } # use constant mean basis, with no constraint on optimization m1<- rgasp(design = input, response = output, lower_bound=FALSE) # plot plot(m1)