Description Details Author(s) References Examples
The 'minimaxdesign' package provides functions for generating minimax designs and minimax projection designs.
Package: | minimaxdesign |
Type: | Package |
Version: | 0.1.4 |
Date: | 2019-05-03 |
License: | GPL (>= 2) |
Provides two main functions, minimax()
and miniMaxPro()
, for computing minimax and minimax projection designs using the minimax clustering algorithm in Mak and
Joseph (2018) <DOI:10.1080/10618600.2017.1302881>. Current design region options
include the unit hypercube ("hypercube"), the unit simplex ("simplex"), the unit ball
("ball"), as well as user-defined constraints on the unit hypercube ("custom"). Minimax
designs can also be computed on user-provided images using the function minimax.map()
. Design quality can be assessed using the function mMdist()
, which computes the minimax (fill) distance of a design.
Simon Mak
Maintainer: Simon Mak <smak6@gatech.edu>
Mak, S. and Joseph, V. R. (2018). Minimax and minimax projection designs using clustering. Journal of Computational and Graphical Statistics, 27(1):166-178.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ## Not run:
#20-point minimax design on the hypercube [0,1]^2
D <- minimax(N=20,p=2)
plot(NULL,xlim=c(0,1),ylim=c(0,1),xlab="x1",ylab="x2") #set up plot
polygon(c(0,0,1,1),c(0,1,1,0),col="gray") #design space
points(D,xlim=c(0,1),ylim=c(0,1),xlab="x1",ylab="x2",pch=16) #design points
mM <- mMdist(D)
mM$dist #minimax (fill) distance
lines(rbind(mM$far.pt,mM$far.despt),col="red",lty=2,lwd=2) #plot farthest point
#20-point minimax design on custom design space (inequalities on [0,1]^2)
ineqs <- function(xx){ #user-defined inequalities
bool.vec <- rep(TRUE,length(xx))
bool.vec[1] <- (xx[2]<=2-2*xx[1]) #inequality 1: x2 <= 2 - 2*x1
bool.vec[2] <- (xx[1]>=xx[2]) #inequality 2: x1 >= x2
return(all(bool.vec))
}
D <- minimax(N=20,p=2,region="custom",const=ineqs)
plot(NULL,xlim=c(0,1),ylim=c(0,1),xlab="x1",ylab="x2") #set up plot
polygon(c(0,2/3,1),c(0,2/3,0),col="gray") #design space
points(D,pch=16) #design points
mM <- mMdist(D,region="custom",const=ineqs)
mM$dist #minimax (fill) distance
lines(rbind(mM$far.pt,mM$far.despt),col="red",lty=2,lwd=2) #plot farthest point
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.