## Examples for verification
w<-2 # Number of completion modes
Rs<-2 # Number of resources
ratio<-1.0 # Percentage 0-1
# Specification of Logic Domain
LD<-matrix(c(1,1,0,1,0,0,
0,0.8415,0,0,0,0,
0,0,1,1,0,0,
0,0,0,0.6478,0,0,
0,0,0,0,0,0,
0,0,0,0,0,0),nrow=6, byrow=TRUE)
# Specification of Time Domain
TD<-matrix(c(1.9755,2.3408,
2.6564,2.7002,
0.9018,1.1077,
1.0237,1.2117,
0,0,
0,0),nrow=6, byrow=TRUE)
# Specification of Cost Domain
CD<-matrix(c(2.0768,2.3008,
0.4819,0.5279,
3.4315,3.7682,
3.5902,3.8245,
0,0,
0,0),nrow=6, byrow=TRUE)
# Specification of Cost Domain
QD<-matrix(c(0.8723,0.9102,
0.9017,0.8215,
0.7250,0.9911,
1.0000,0.8245,
0,0,
0,0),nrow=6, byrow=TRUE)
# Specification of Resource Domain
RD<-matrix(c(2.2130,2.4636,1.7948,2.0603,
1.9709,2.1952,0.8216,0.8890,
1.6255,1.9432,1.8715,2.2341,
1.1377,1.3528,1.6768,1.8769,
0,0,0,0,
0,0,0,0),nrow=6, byrow=TRUE)
# PDM without QD
PDM<-cbind(LD,TD,CD,RD)
# Percentiles without QD
CONST<-percent(PDM,type=c("c","r","s","t"),w,Rs,ratio)
CONST
# PDM with QD
PDM<-cbind(LD,TD,CD,QD,RD)
# Percentiles with QD
CONST<-percent(PDM,type=c("c","q","r","s","t"),w,Rs,ratio)
CONST
Run the code above in your browser using DataLab