remove.projections {demu} | R Documentation |
remove.projections()
is a helper function to identify all lower-dimensional marginal projection points of the existing design points indexed by curpts
. This function can be used to remove a subset of points from the candidate set in order to enforce non-collapsingness of when sequentially adding design points.
remove.projections(curpts,X)
curpts |
Indices of points currently in the design |
X |
An |
A list containing the vector curpts
, the vector projpts
which contains the identified projection points of the current design, and allpts
.
demu-package
sim.dpp.modal.seq
library(demu)
n1=3
n2=3
n3=3
rho=rep(1e-10,2)
ngrid=10
x=seq(0,1,length=ngrid)
X=as.matrix(expand.grid(x,x))
l.d=makedistlist(X)
# Initial design
R=rhomat(l.d,rho)$R
pts.1=sim.dpp.modal(R,n1)
pts.1.proj=remove.projections(pts.1,X)
# Plot - design points in black, design+projection points in grey.
#plot(X,xlim=c(0,1),ylim=c(0,1))
#points(X[pts.1.proj$projpts,],pch=20,cex=2,col="grey")
#points(X[pts.1,],pch=20,cex=2)