plot.spv  R Documentation 
Produce Variance Dispersion Graphs and/or Fraction of Design Space plots for
experimental designs. There are methods for the S3 classes spv
,
spvlist
, spvforlist
and spvlistforlist
– see
spv
.
## S3 method for class 'spv' plot( x, which = c("fds", "vdgsim", "vdgquantile", "vdgboth", "boxplots"), np = 50, alpha = 7/sqrt(length(x$spv)), points.colour = "#39BEB1", points.size = 2, tau = c(0.05, 0.95), radii = 21, hexbin = FALSE, bins = 80, df = 10, lines.size = 1, origin = rep(0, ncol(x$sample)), method, arrange = FALSE, ... ) ## S3 method for class 'spvforlist' plot( x, which = c("fds", "vdgsim", "vdgquantile", "vdgboth", "boxplots"), np = 50, alpha = 7/sqrt(length(x[[1]]$spv)), points.colour = "#39BEB1", points.size = 2, tau = c(0.05, 0.95), radii = 21, hexbin = FALSE, bins = 80, df = 10, lines.size = 1, origin = rep(0, ncol(x[[1]]$sample)), method, arrange = FALSE, ... ) ## S3 method for class 'spvlist' plot( x, which = c("fds", "vdgsim", "vdgquantile", "vdgboth", "boxplots"), np = 50, alpha = 7/sqrt(length(x[[1]]$spv)), points.colour = "#39BEB1", points.size = 2, tau = c(0.05, 0.95), radii = 21, hexbin = FALSE, bins = 80, VRFDS = FALSE, df = 10, lines.size = 1, origin = rep(0, ncol(x[[1]]$sample)), method, arrange = FALSE, ... ) ## S3 method for class 'spvlistforlist' plot( x, which = c("fds", "vdgsim", "vdgquantile", "vdgboth", "boxplots"), np = 50, alpha = 7/sqrt(length(x[[1]][[1]]$spv)), points.colour = "#39BEB1", points.size = 2, tau = c(0.05, 0.95), radii = 21, hexbin = FALSE, bins = 80, df = 10, lines.size = 1, origin = rep(0, ncol(x[[1]][[1]]$sample)), method, arrange = FALSE, ... )
x 
an object of type 
which 
either a numeric vector of integers or a character vector indicating which plots to produce. The possible plots are:

np 
scalar; the number of points to use for calculating the fraction of design space criterion. 
alpha 
the alpha transparency coefficient for the plots 
points.colour 
colour for points in scatterplot of SPV against the radius 
points.size 
size for points in scatterplot of SPV against the radius 
tau 
the tau parameter for 
radii 
either a numeric vector containing the radii to use for calculating the mean spherical SPV over the spherical design space, or an integer (length one vector) giving the number of radii to use for calculationg the mean spherical SPV. If missing, the mean spherical SPV is not used. 
hexbin 
logical indicating whether hexagonal binning should be used to display density instead of colour transparency 
bins 
argument passed to 
df 
degreesoffreedom parameter passed to 
lines.size 
line size passed to 
origin 
numeric vector specifying the origin of the design space 
method 
optional; passed to 
arrange 
Logical indicating whether to return a single graphical object arranging the
resulting plots in a single plot window via 
... 
additional arguments passed to 
VRFDS 
logical indicating whether to construct a variance ratio FDS plot or not (only for class 
Returns a list of ggplot
graphical objects (or grobs) with names corresponding
to the character version of which
. These plot objects can be manipulated by changing plot
aesthetics and theme elements.
Pieter C. Schoonees
Pieter C. Schoonees, Niel J. le Roux, Roelof L.J. Coetzer (2016). Flexible Graphical Assessment of Experimental Designs in R: The vdg Package. Journal of Statistical Software, 74(3), 122. doi: 10.18637/jss.v074.i03.
# Single design (class 'spv') # Larger n should be used in actual cases library(rsm) bbd3 < as.data.frame(bbd(3)[,3:5]) colnames(bbd3) < paste0("x", 1:3) quad.3f < formula(~ x1*x2*x3  x1:x2:x3 + I(x1^2) + I(x2^2) + I(x3^2)) set.seed(1234) out < spv(n = 1000, design = bbd3, type = "spherical", formula = quad.3f) out plot(out) # List of designs (class 'spvlist') ## Not run: data(SCDH5); data(SCDDL5) des.list < list(SCDH5 = SCDH5, SCDDL5 = SCDDL5) quad.5f < formula(~ x1 + x2 + x3 + x4 + x5 + x1:x2 + x1:x3 + x1:x4 + x1:x5 + x2:x3 + x2:x4 + x2:x5 + x3:x4 + x3:x5 + x4:x5 + I(x1^2) + I(x2^2) + I(x3^2) + I(x4^2) + I(x5^2)) out2 < spv(n = 500, design = des.list, type = "spherical", formula = quad.5f) out2 plot(out2) ## End(Not run) # List of formulae (class 'spvforlist') ## Not run: fact3 < expand.grid(x1 = c(1,1), x2 = c(1, 1), x3 = c(1,1)) lin.3f < formula(~ x1 + x2 + x3) int.3f < formula(~ (x1+x2+x3)^2) set.seed(4312) out3 < spv(n = 500, design = fact3, type = "cuboidal", formula = list(linear = lin.3f, interaction = int.3f)) out3 plot(out3) ## End(Not run) # List of formulae and designs (class 'spvlistforlist') ## Not run: fact3.n < rbind(fact3, 0, 0, 0) set.seed(4312) out4 < spv(n = 200, design = list(factorial = fact3, factorial.with.cntr = fact3.n), type = "cuboidal", formula = list(linear = lin.3f, interaction = int.3f)) out4 plot(out4) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.