Imbrication of two experimental design sets

Description

Procedure to nest two experimental design sets.

Usage

SubstDesign(PX2 = NULL,PX1 = NULL)

Arguments

Details

The procedure does not change the experimental design set PX2 which must have a number of points smaller than the one of PX1. During the procedure, the points of PX1 the closest to those of PX2 are removed and they are replaced by the points of PX2. Thus, the length of the final PX1 could be larger than the one of the initial PX1.

Value

Author(s)

Loic Le Gratiet

See Also

MuFicokm, NestedDesignBuild, NestedDesign

Examples

	dimension <- 2
	nD1 <- 100
	nD2 <- 50		
	set.seed(1);D1 <- matrix(runif(n=nD1*dimension, 0,1),ncol=dimension)
	set.seed(2);D2 <- matrix(runif(n=nD2*dimension, 0,1),ncol=dimension)
	subDes <- SubstDesign(PX2 = D2, PX1 = D1)

	op <- par(mfrow=c(2,1))
	plot(rbind(D1,D2),col=c(rep(1,nD1),rep(2,nD2)),
		pch=c(rep(1,nD1),rep(2,nD2)),xlab="x1",ylab="x2")
	plot(rbind(subDes$PX,D2),col=c(rep(1,dim(subDes$PX)[1]),rep(2,nD2)),
		pch=c(rep(1,dim(subDes$PX)[1]),rep(2,nD2)),xlab="x1",ylab="x2")

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(MuFiCokriging)
Loading required package: DiceKriging
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/MuFiCokriging/SubstDesign.Rd_%03d_medium.png", width=480, height=480)
> ### Name: SubstDesign
> ### Title: Imbrication of two experimental design sets
> ### Aliases: SubstDesign
> 
> ### ** Examples
> 
> 	dimension <- 2
> 	nD1 <- 100
> 	nD2 <- 50		
> 	set.seed(1);D1 <- matrix(runif(n=nD1*dimension, 0,1),ncol=dimension)
> 	set.seed(2);D2 <- matrix(runif(n=nD2*dimension, 0,1),ncol=dimension)
> 	subDes <- SubstDesign(PX2 = D2, PX1 = D1)
> 
> 	op <- par(mfrow=c(2,1))
> 	plot(rbind(D1,D2),col=c(rep(1,nD1),rep(2,nD2)),
+ 		pch=c(rep(1,nD1),rep(2,nD2)),xlab="x1",ylab="x2")
> 	plot(rbind(subDes$PX,D2),col=c(rep(1,dim(subDes$PX)[1]),rep(2,nD2)),
+ 		pch=c(rep(1,dim(subDes$PX)[1]),rep(2,nD2)),xlab="x1",ylab="x2")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>