daewr
Design and Analysis of Experiments with R

Package index
Search the daewr package
Vignettes
Functions
141
Source code
36
Man pages
103
Home
/
R-Forge
/
daewr
/
R/DefScreen.R

R/DefScreen.R
In daewr: Design and Analysis of Experiments with R

Defines functions std.ord interleave DefScreen

Documented in DefScreen 

DefScreen <- function(m=0, c=0,center=0,randomize=FALSE) {
# Definitive Screening Designs for 3-level factors ##
if (m<4) {stop("Definitive Screening Designs only exist for 4-12 3-level factors") }
if (m>12) {stop("Definitive Screening Designs only exist for 4-12 3-level factors") }
if (c<0) {stop("Definitive Screening Designs only exist for 0-4 2-level categorical factors") }
if (c>4) {stop("Definitive Screening Designs only exist for 0-4 2-level categorical factors") }
if (c>0 && center>0) {stop("Cannnot add center points to design with 2-level categorical factors") }

  if (c==0)              {
# Design for 4 3-level factors 
if (m==4) {
  f1 <- matrix(c(0,1,-1,-1,-1,0,-1,1,-1,-1,0,-1,-1,1,1,0), ncol=m,byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D")
  des <- data.frame(des)
              }
# Design for 5 3-level factors 
if (m==5)  {
  f1 <- matrix(c(0,1,1,-1,-1, 1,0,-1,-1,1 ,1,-1,0,1,-1, 1,-1,1,0,1, 1,1,1,1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E")
  des <- data.frame(des)
}

if (m==6)     {
  f1 <- matrix(c(0,1,-1,-1,-1,-1,1,0,-1,1,1,-1,-1,-1,0,1,-1,-1,-1,1,1,0,1,-1,1,-1,1,-1,0,-1,1,1,1,1,-1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F")
  des <- data.frame(des)
              }
# Design for 7 3-level factors 
if (m==7)     {
  f1 <- matrix(c(0,1,-1,1,-1,1,-1, -1,0,1,-1,1,1,-1, 1,-1,0,1,1,1,1, 1,-1, -1,0,1,-1,-1, 
  -1,-1,1,1,0,-1,-1, -1,1,-1,1,1,0,1, 1,1,1,1,1,-1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F","G")
  des <- data.frame(des)
}
# Design for 8 3-level factors 
if (m==8)     {
  f1 <- matrix(c(0,-1,1,1,-1,1,1,1, -1,0,-1,1,1,1,1,-1, -1,-1,0,1,1,-1,-1,1,
  1,-1,1,0,1,1,-1,-1, -1,-1,1,-1,0,-1,1,-1, 1,-1,-1,-1,1,0,1,1, 
  -1,1,1,-1,1,1,0,1, 1,1,1,1,1,-1,1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F","G","H")
  des <- data.frame(des)
}
# Design for 9 3-level factors
if (m==9)     {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,1,  1,0,1,-1,1,-1,-1,1,-1,
                 -1,1,0,-1,1,-1,1,-1,-1,  -1,-1,1,0,1,-1,-1,-1,1,
                 1,-1,1,-1,0,1,1,-1,-1, -1,-1,-1,-1,1,0,1,1,1,
                 1,1,-1,-1,1,1,0,-1,1, -1,-1,-1,1,1,1,-1,0,-1,
                 -1,1,1,-1,-1,1,-1,1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F","G","H","J")
  des <- data.frame(des)
}
# Design for 10 3-level factors
if (m==10)     {
  f1 <- matrix(c(0,1,1,-1,1,1,1,1,-1,1, 1,0,-1,1,1,-1,1,1,-1,-1,
                 -1,1,0,-1,-1,-1,1,-1,-1,-1,  -1,1,1,0,1,-1,-1,1,1,-1,
                 -1,-1,-1,-1,0,1,1,1,1,-1, -1,1,-1,1,1,0,1,-1,1,1,
                 1,1,-1,-1,-1,-1,0,1,1,1, 1,1,1,1,-1,1,1,0,1,-1,
                 1,1,-1,-1,1,1,-1,-1,0,-1, 1,-1,1,-1,1,-1,1,-1,1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K")
  des <- data.frame(des)
}
# Design for 11 3-level factors
if (m==11)     {
  f1 <- matrix(c(0,-1,1,-1,-1,-1,-1,-1,1,-1,1, -1,0,-1,-1,1,-1,-1,-1,-1,1,1,
                 -1,-1,0,1,1,1,1,-1,-1,-1,1, -1,-1,-1,0,-1,1,1,-1,1,1,-1,
                 1,-1,-1,1,0,1,-1,1,1,1,1, -1,-1,1,1,-1,0,-1,1,-1,1,-1,
                 -1,-1,-1,1,1,-1,0,1,1,-1,-1, -1,1,1,1,-1,-1,1,0,1,1,1,
                 -1,1,-1,-1,-1,1,-1,1,0,-1,1, 1,-1,-1,-1,-1,-1,1,1,-1,0,1,
                 1,1,-1,1,-1,-1,-1,-1,-1,-1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L")
  des <- data.frame(des)
}
# Design for 12 3-level factors
if (m==12)     {
  f1 <- matrix(c(0,-1,-1,1,-1,1,-1,1,1,1,-1,1,
                 -1,0,1,1,1,1,1,1,1,-1,-1,-1,
                 1,1,0,-1,1,1,-1,1,1,-1,1,1,
                 1,-1,-1,0,1,-1,1,-1,1,-1,-1,1,
                 1,1,1,1,0,-1,1,1,1,1,1,1,
                 1,-1,1,-1,1,0,1,1,-1,1,-1,1,
                 1,1,1,1,-1,1,0,-1,-1,-1,-1,1,
                 -1,-1,1,1,1,-1,-1,0,-1,-1,1,1,
                 1,-1,1,1,1,1,-1,-1,0,1,1,-1,
                 1,1,-1,1,1,-1,-1,1,-1,0,-1,-1,
                 -1,1,-1,1,1,1,1,-1,-1,1,0,1,
                 1,-1,-1,1,-1,1,1,1,-1,-1,1,0), ncol=m, byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(m)
  des <- des[ord , ]
  des <- rbind(des,c(rep(0,m)))
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M")
  des <- data.frame(des)
}
                             }
# Definitive Screening Designs for 3-level factors  and one 2-level factor ##
if (c==1)                    {  
# Design for 4 3-level factors 
if (m==4)  {
  f1 <- matrix(c(0,1,1,1,1, 1,0,-1,1,1, 1,-1,0,-1,1, 1,1,-1,0,-1, 
                 1,1,1,-1,1, 1,-1,1,1,-1, 0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  #z <- rbind(c(rep(0,m),-1),c(rep(0,m),-1))
  #des <- rbind(des,z)
  colnames(des) <- c("A","B","C","D","E")
  des <- data.frame(des)
          }
if (m==5)  {
  f1 <- matrix(c(0,1,1,1,1,1, 1,0,-1,1,1,-1,
                 1,-1,0,-1,1,1, 1,1,-1,0,-1,1,
                 1,1,1,-1,0,-1, 1,-1,1,1,-1,1,
                 0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  #z <- rbind(c(rep(0,m),-1),c(rep(0,m),-1))
  #des <- rbind(des,z)
  colnames(des) <- c("A","B","C","D","E","F")
  des <- data.frame(des)
}
if (m==6)  {
  f1 <- matrix(c(0,1,1,1,1,1,1, 1,0,1,1,-1,1,-1,
                1,-1,0,1,1,-1,1, 1,-1,-1,0,1,1,-1,
                1,1,-1,-1,0,1,1, 1,-1,1,-1,-1,0,1,
                1,1,-1,1,-1,-1,1, 1,1,1,-1,1,-1,-1,
                0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G")
  des <- data.frame(des)
}
if (m==7)  {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,  1,0,1,1,-1,1,-1,-1,
                 1,-1,0,1,1,-1,1,-1, 1,1,-1,0,1,1,-1,1,
                 1,1,-1,-1,0,1,1,-1, 1,-1,1,-1,-1,0,1,1,
                 1,1,-1,1,-1,-1,0,1, 1,1,1,-1,1,-1,-1,1,
                 0,0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G","H")
  des <- data.frame(des)
}
if (m==8)  {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,
                 1,-1,0,-1,1,1,-1,-1,1, 1,-1,-1,0,1,1,1,1,-1,
                 1,-1,1,1,0,-1,-1,1,-1, 1,-1,1,1,-1,0,1,-1,1,
                 1,1,-1,1,-1,1,0,-1,-1, 1,1,-1,1,1,-1,-1,0,1,
                 1,1,1,-1,-1,1,-1,1,1, 1,1,1,-1,1,-1,1,-1,-1,
                 0,0,0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G","H","J")
  des <- data.frame(des)
}
if (m==9)  {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,1,
               1,-1,0,-1,1,1,-1,-1,1,1, 1,-1,-1,0,1,1,1,1,-1,-1,
               1,-1,1,1,0,-1,-1,1,-1,1, 1,-1,1,1,-1,0,1,-1,1,-1,
               1,1,-1,1,-1,1,0,-1,-1,1, 1,1,-1,1,1,-1,-1,0,1,-1,
               1,1,1,-1,-1,1,-1,1,0,-1, 1,1,1,-1,1,-1,1,-1,-1,1,
               0,0,0,0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K")
  des <- data.frame(des)
}
if (m==10)  {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,
               1,-1,0,1,-1,1,1,1,-1,-1,-1, 1,1,-1,0,1,-1,1,1,1,-1,-1,
               1,-1,1,-1,0,1,-1,1,1,1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,
               1,-1,-1,-1,1,-1,0,1,-1,1,1, 1,1,-1,-1,-1,1,-1,0,1,-1,1,
               1,1,1,-1,-1,-1,1,-1,0,1,-1, 1,1,1,1,-1,-1,-1,1,-1,0,1,
               1,-1,1,1,1,-1,-1,-1,1,-1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,
               0,0,0,0,0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L")
  des <- data.frame(des)
}
if (m==11)  {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,-1,
               1,-1,0,1,-1,1,1,1,-1,-1,-1,1, 1,1,-1,0,1,-1,1,1,1,-1,-1,-1,
               1,-1,1,-1,0,1,-1,1,1,1,-1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,-1,
               1,-1,-1,-1,1,-1,0,1,-1,1,1,1, 1,1,-1,-1,-1,1,-1,0,1,-1,1,1, 
               1,1,1,-1,-1,-1,1,-1,0,1,-1,1, 1,1,1,1,-1,-1,-1,1,-1,0,1,-1,
               1,-1,1,1,1,-1,-1,-1,1,-1,0,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,1,
               0,0,0,0,0,0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M")
  des <- data.frame(des)
}
if (m==12)  {
  f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1, 1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
               1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1, 1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
               1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1, 1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
               1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1, 1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
               1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1, 1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,
               1,-1,1,1,1,1,-1,-1,1,-1,0,-1,1, 1,-1,-1,1,1,1,1,-1,-1,1,-1,0,-1,
               1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1, 1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
               0,0,0,0,0,0,0,0,0,0,0,0,-1), ncol=(m+1), byrow=TRUE)
  f2 <- (-1)*f1
  des <- rbind(f1,f2)
  ord<-std.ord(nrow(f1))
  des <- des[ord , ]
  colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N")
  des <- data.frame(des)
}


                              }
if (c==2)                    {  
  # Design for 4 3-level factors 
  if (m==4)  {
    f1 <- matrix(c(0,1,1,1,1,1, 1,0,-1,1,1,-1, 1,-1,0,-1,1,1,
                   1,1,-1,0,-1,1, 1,1,1,-1,1,-1, 1,-1,1,1,-1,1,
                   0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F")
    des <- data.frame(des)
        }
  if (m==5)  {
    f1 <- matrix(c(0,1,1,1,1,1,1, 1,0,1,1,-1,1,-1,
                   1,-1,0,1,1,-1,1, 1,-1,-1,0,1,1,-1,
                   1,1,-1,-1,0,1,1, 1,-1,1,-1,-1,1,1,
                   1,1,-1,1,-1,-1,1, 1,1,1,-1,1,-1,-1,
                   0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G")
    des <- data.frame(des)
  }
  if (m==6)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1, 1,0,1,1,-1,1,-1,-1,
                 1,-1,0,1,1,-1,1,-1, 1,-1,-1,0,1,1,-1,1,
                 1,1,-1,-1,0,1,1,-1, 1,-1,1,-1,-1,0,1,1,
                 1,1,-1,1,-1,-1,1,1, 1,1,1,-1,1,-1,-1,1,
                 0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H")
    des <- data.frame(des)
  }
  if (m==7)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,
                   1,-1,0,-1,1,1,-1,-1,1, 1,-1,-1,0,1,1,1,1,-1,
                   1,-1,1,1,0,-1,-1,1,-1, 1,-1,1,1,-1,0,1,-1,1,
                   1,1,-1,1,-1,1,0,-1,-1, 1,1,-1,1,1,-1,-1,1,1,
                   1,1,1,-1,-1,1,-1,1,1, 1,1,1,-1,1,-1,1,-1,-1,
                   0,0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J")
    des <- data.frame(des)
  }
  if (m==8)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,1,
                 1,-1,0,-1,1,1,-1,-1,1,1, 1,-1,-1,0,1,1,1,1,-1,-1,
                 1,-1,1,1,0,-1,-1,1,-1,1, 1,-1,1,1,-1,0,1,-1,1,-1,
                 1,1,-1,1,-1,1,0,-1,-1,1, 1,1,-1,1,1,-1,-1,0,1,-1,
                 1,1,1,-1,-1,1,-1,1,1,-1, 1,1,1,-1,1,-1,1,-1,-1,1,
                 0,0,0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K")
    des <- data.frame(des)
  }
  if (m==9)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,
                 1,-1,0,1,-1,1,1,1,-1,-1,-1, 1,1,-1,0,1,-1,1,1,1,-1,-1,
                 1,-1,1,-1,0,1,-1,1,1,1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1, 
                 1,-1,-1,-1,1,-1,0,1,-1,1,1, 1,1,-1,-1,-1,1,-1,0,1,-1,1,
                 1,1,1,-1,-1,-1,1,-1,0,1,-1, 1,1,1,1,-1,-1,-1,1,-1,1,1,
                 1,-1,1,1,1,-1,-1,-1,1,-1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,
                 0,0,0,0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L")
    des <- data.frame(des)
  }
  if (m==10)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,-1,
                   1,-1,0,1,-1,1,1,1,-1,-1,-1,1, 1,1,-1,0,1,-1,1,1,1,-1,-1,-1,
                   1,-1,1,-1,0,1,-1,1,1,1,-1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,-1,
                   1,-1,-1,-1,1,-1,0,1,-1,1,1,1, 1,1,-1,-1,-1,1,-1,0,1,-1,1,1,
                   1,1,1,-1,-1,-1,1,-1,0,1,-1,1, 1,1,1,1,-1,-1,-1,1,-1,0,1,-1,
                   1,-1,1,1,1,-1,-1,-1,1,-1,1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,1,
                   0,0,0,0,0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M")
    des <- data.frame(des)
  }
  if (m==11)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1, 1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
                 1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1, 1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
                 1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1, 1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
                 1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1, 1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
                 1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1, 1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,
                 1,-1,1,1,1,1,-1,-1,1,-1,0,-1,1, 1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,
                 1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1, 1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
                 0,0,0,0,0,0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N")
    des <- data.frame(des)
  }
  if (m==12)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1,1, 
                   1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
                   1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
                   1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1,
                   1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
                   1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1,
                   1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
                   1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,
                   1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
                   1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,
                   1,-1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,
                   1,-1,-1,1,1,1,1,-1,-1,1,-1,0,-1,1,
                   1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,
                   1,-1,1,-1,-1,1,1,1,1,-1,-1,1,1,1,
                   0,0,0,0,0,0,0,0,0,0,0,0,-1,-1), ncol=(m+2), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N","O")
    des <- data.frame(des)
  }


                               }


if (c==3)                     {  
  # Design for 4 3-level factors 
  if (m==4)  {
    f1 <- matrix(c(0,1,1,1,1,1,1, 1,0,1,1,-1,1,-1,
                   1,-1,0,1,1,-1,1, 1,-1,-1,0,1,1,-1,
                   1,1,-1,-1,1,1,  1,-1,1,-1,-1,1,1,1,
                   1,1,-1,1,-1,-1,1,  1,1,1,-1,1,-1,-1,
                   0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G")
    des <- data.frame(des)
  }
  if (m==5)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1, 1,0,1,1,-1,1,-1,-1,
                 1,-1,0,1,1,-1,1,-1, 1,-1,-1,0,1,1,-1,1,
                 1,1,-1,-1,0,1,1,-1, 1,-1,1,-1,-1,1,1,1,
                 1,1,-1,1,-1,-1,1,1, 1,1,1,-1,1,-1,-1,1,
                 0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H")
    des <- data.frame(des)
  }
  if (m==6)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,
                   1,-1,0,-1,1,1,-1,-1,1, 1,-1,-1,0,1,1,1,1,-1,
                   1,-1,1,1,0,-1,-1,1,-1, 1,-1,1,1,-1,0,1,-1,1,
                   1,1,-1,1,-1,1,1,-1,-1, 1,1,-1,1,1,-1,-1,1,1,
                   1,1,1,-1,-1,1,-1,1,1, 1,1,1,-1,1,-1,1,-1,-1,
                   0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J")
    des <- data.frame(des)
  }
  if (m==7)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,1,
                   1,-1,0,-1,1,1,-1,-1,1,1, 1,-1,-1,0,1,1,1,1,-1,-1,
                   1,-1,1,1,0,-1,-1,1,-1,1, 1,-1,1,1,-1,0,1,-1,1,-1,
                   1,1,-1,1,-1,1,0,-1,-1,1, 1,1,-1,1,1,-1,-1,1,1,-1,
                   1,1,1,-1,-1,1,-1,1,1,-1, 1,1,1,-1,1,-1,1,-1,-1,1,
                   0,0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K")
    des <- data.frame(des)
  }
  if (m==8)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,  1,0,1,-1,1,1,1,-1,-1,-1,1,
                   1,-1,0,1,-1,1,1,1,-1,-1,-1, 1,1,-1,0,1,-1,1,1,1,-1,-1,
                   1,-1,1,-1,0,1,-1,1,1,1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,
                   1,-1,-1,-1,1,-1,0,1,-1,1,1,  1,1,-1,-1,-1,1,-1,0,1,-1,1,
                   1,1,1,-1,-1,-1,1,-1,1,1,-1, 1,1,1,1,-1,-1,-1,1,-1,1,1,
                   1,-1,1,1,1,-1,-1,-1,1,-1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,
                   0,0,0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L")
    des <- data.frame(des)
  }
  if (m==9)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,-1,
                 1,-1,0,1,-1,1,1,1,-1,-1,-1,1, 1,1,1,0,1,-1,1,1,1,-1,-1,-1,
                 1,-1,1,-1,0,1,-1,1,1,1,-1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,-1,
                 1,-1,-1,-1,1,-1,0,1,-1,1,1,1, 1,1,-1,-1,-1,1,-1,0,1,-1,1,1,
                 1,1,1,-1,-1,-1,1,-1,0,1,-1,1, 1,1,1,1,-1,-1,-1,1,-1,1,1,-1,
                 1,-1,1,1,1,-1,-1,-1,1,-1,1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,1,
                 0,0,0,0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M")
    des <- data.frame(des)
  }
  if (m==10)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1, 1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
                 1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1, 1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
                 1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1, 1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
                 1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1, 1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
                 1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1, 1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,
                 1,-1,1,1,1,1,-1,-1,1,-1,1,-1,1, 1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,
                 1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1, 1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
                 0,0,0,0,0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N")
    des <- data.frame(des)
  }
  if (m==11)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1,1,
                 1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
                 1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
                 1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1,
                 1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
                 1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1,
                 1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
                 1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,
                 1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
                 1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,
                 1,-1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,
                 1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,1,
                 1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,
                 1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1,
                 0,0,0,0,0,0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N","O")
    des <- data.frame(des)
  }
  if (m==12)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
                   1,0,1,1,-1,1,-1,-1,1,-1,1,1,-1,1,-1,
                   1,-1,0,1,1,-1,1,-1,1,1,1,-1,1,-1,-1,
                   1,-1,-1,0,1,1,-1,1,1,1,-1,1,-1,-1,-1,
                   1,1,-1,-1,0,1,1,-1,1,-1,1,-1,-1,-1,1,
                   1,-1,1,-1,-1,0,1,1,1,1,-1,-1,-1,1,1,
                   1,1,-1,1,-1,-1,0,1,1,-1,-1,-1,1,1,-1,
                   1,1,1,-1,1,-1,-1,0,1,-1,-1,1,1,-1,1,
                   -1,1,1,1,1,1,1,1,0,-1,-1,-1,-1,-1,-1,
                   -1,-1,1,1,-1,1,-1,-1,1,0,-1,-1,1,-1,1,
                   -1,1,1,-1,1,-1,-1,-1,1,1,0,-1,-1,1,-1,
                   -1,1,-1,1,-1,-1,-1,1,1,1,1,0,-1,-1,1,
                   -1,-1,1,-1,-1,-1,1,1,1,-1,1,1,1,-1,-1,
                   -1,1,-1,-1,-1,1,1,-1,1,1,-1,1,1,1,-1,
                   -1,-1,-1,-1,1,1,-1,1,1,-1,1,-1,1,1,1,
                    -1,-1,-1,1,1,-1,1,-1,1,-1,-1,1,-1,1,1,
                    0,0,0,0,0,0,0,0,0,0,0,0,-1,-1,-1), ncol=(m+3), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N","O","P")
    des <- data.frame(des)
  }


                               }


if (c==4)                     {  
  # Design for 4 3-level factors 
  if (m==4)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1, 1,0,1,1,-1,1,-1,-1,
                   1,-1,0,1,1,-1,1,-1, 1,-1,-1,0,1,1,-1,1,
                   1,1,-1,-1,1,1,1,-1, 1,-1,1,-1,-1,1,1,1,
                   1,1,-1,1,-1,-1,1,1, 1,1,1,-1,1,-1,-1,1,
                   0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H")
    des <- data.frame(des)
  }
  if (m==5)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,
                   1,-1,0,-1,1,1,-1,-1,1, 1,-1,-1,0,1,1,1,1,-1,
                   1,-1,1,1,0,-1,-1,1,-1, 1,-1,1,1,-1,1,1,-1,1,
                   1,1,-1,1,-1,1,1,-1,-1, 1,1,-1,1,1,-1,-1,1,1,
                   1,1,1,-1,-1,1,-1,1,1, 1,1,1,-1,1,-1,1,-1,-1,
                   0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J")
    des <- data.frame(des)
  }
  if (m==6)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1, 1,0,-1,-1,-1,-1,1,1,1,1,
                   1,-1,0,-1,1,1,-1,-1,1,1, 1,-1,-1,0,1,1,1,1,-1,-1,
                   1,-1,1,1,0,-1,-1,1,-1,1, 1,-1,1,1,-1,0,1,-1,1,-1,
                   1,1,-1,1,-1,1,1,-1,-1,1, 1,1,-1,1,1,-1,-1,1,1,-1,
                   1,1,1,-1,-1,1,-1,1,1,-1, 1,1,1,-1,1,-1,1,-1,-1,1,
                   0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K")
    des <- data.frame(des)
  }
  if (m==7)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,
                 1,-1,0,1,-1,1,1,1,-1,-1,-1, 1,1,-1,0,1,-1,1,1,1,-1,-1,
                 1,-1,1,-1,0,1,-1,1,1,1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,
                 1,-1,-1,-1,1,-1,0,1,-1,1,1, 1,1,-1,-1,-1,1,-1,1,1,-1,1,
                 1,1,1,-1,-1,-1,1,-1,1,1,-1, 1,1,1,1,-1,-1,-1,1,-1,1,1,
                 1,-1,1,1,1,-1,-1,-1,1,-1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,
                 0,0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L")
    des <- data.frame(des)
  }
  if (m==8)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1, 1,0,1,-1,1,1,1,-1,-1,-1,1,-1,
                 1,-1,0,1,-1,1,1,1,-1,-1,-1,1, 1,1,-1,0,1,-1,1,1,1,-1,-1,-1,
                 1,-1,1,-1,0,1,-1,1,1,1,-1,-1, 1,-1,-1,1,-1,0,1,-1,1,1,1,-1,
                 1,-1,-1,-1,1,-1,0,1,-1,1,1,1, 1,1,-1,-1,-1,1,-1,0,1,-1,1,1,
                 1,1,1,-1,-1,-1,1,-1,1,1,-1,1, 1,1,1,1,-1,-1,-1,1,-1,1,1,-1,
                 1,-1,1,1,1,-1,-1,-1,1,-1,1,1, 1,1,-1,1,1,1,-1,-1,-1,1,-1,1,
                 0,0,0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M")
    des <- data.frame(des)
  }
  if (m==9)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1, 1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
                   1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1, 1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
                   1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1, 1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
                   1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1, 1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
                   1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1, 1,1,1,1,1,-1,-1,1,-1,1,-1,1,-1,
                   1,-1,1,1,1,1,-1,-1,1,-1,1,-1,1, 1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,
                   1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1, 1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
                   0,0,0,0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N")
    des <- data.frame(des)
  }
  if (m==10)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1,1,
                   1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,
                   1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1,1,
                   1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,-1,
                   1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1,-1,
                   1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,1,
                   1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,1,
                   1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,1,
                   1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,1,
                   1,1,1,1,1,-1,-1,1,-1,0,-1,1,-1,-1,
                   1,-1,1,1,1,1,-1,-1,1,-1,1,-1,1,-1,
                   1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,1,
                   1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1,-1,
                   1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1,1, 
                   0,0,0,0,0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N","O")
    des <- data.frame(des)
  }
  if (m==11)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
                   1,0,1,1,-1,1,-1,-1,1,-1,1,1,-1,1,-1,
                   1,-1,0,1,1,-1,1,-1,1,1,1,-1,1,-1,-1,
                   1,-1,-1,0,1,1,-1,1,1,1,-1,1,-1,-1,-1,
                   1,1,-1,-1,0,1,1,-1,1,-1,1,-1,-1,-1,1,
                   1,-1,1,-1,-1,0,1,1,1,1,-1,-1,-1,1,1,
                   1,1,-1,1,-1,-1,0,1,1,-1,-1,-1,1,1,-1,
                   1,1,1,-1,1,-1,-1,0,1,-1,-1,1,1,-1,1,
                   -1,1,1,1,1,1,1,1,0,-1,-1,-1,-1,-1,-1,
                   -1,-1,1,1,-1,1,-1,-1,1,0,-1,-1,1,-1,1,
                   -1,1,1,-1,1,-1,-1,-1,1,1,0,-1,-1,1,-1,
                   -1,1,-1,1,-1,-1,-1,1,1,1,1,1,-1,-1,1,
                   -1,-1,1,-1,-1,-1,1,1,1,-1,1,1,1,-1,-1,
                   -1,1,-1,-1,-1,1,1,-1,1,1,-1,1,1,1,-1,
                    -1,-1,-1,-1,1,1,-1,1,1,-1,1,-1,1,1,1,
                    -1,-1,-1,1,1,-1,1,-1,1,-1,-1,1,-1,1,1,
                    0,0,0,0,0,0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N","O","P")
    des <- data.frame(des)
  }
  if (m==12)  {
    f1 <- matrix(c(0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
                   1,0,1,1,-1,1,-1,-1,1,-1,1,1,-1,1,-1,-1,
                   1,-1,0,1,1,-1,1,-1,1,1,1,-1,1,-1,-1,-1,
                   1,-1,-1,0,1,1,-1,1,1,1,-1,1,-1,-1,-1,1,
                   1,1,-1,-1,0,1,1,-1,1,-1,1,-1,-1,-1,1,1,
                   1,-1,1,-1,-1,0,1,1,1,1,-1,-1,-1,1,1,-1,
                   1,1,-1,1,-1,-1,0,1,1,-1,-1,-1,1,1,-1,1,
                   1,1,1,-1,1,-1,-1,0,1,-1,-1,1,1,-1,1,-1,
                   -1,1,1,1,1,1,1,1,0,-1,-1,-1,-1,-1,-1,-1,
                   -1,-1,1,1,-1,1,-1,-1,1,0,-1,-1,1,-1,1,1,
                   -1,1,1,-1,1,-1,-1,-1,1,1,0,-1,-1,1,-1,1,
                   -1,1,-1,1,-1,-1,-1,1,1,1,1,0,-1,-1,1,-1,
                   -1,-1,1,-1,-1,-1,1,1,1,-1,1,1,1,-1,-1,1,
                   -1,1,-1,-1,-1,1,1,-1,1,1,-1,1,1,1,-1,-1,
                   -1,-1,-1,-1,1,1,-1,1,1,-1,1,-1,1,1,1,-1,
                   -1,-1,-1,1,1,-1,1,-1,1,-1,-1,1,-1,1,1,1,
                    0,0,0,0,0,0,0,0,0,0,0,0,-1,-1,-1,-1), ncol=(m+4), byrow=TRUE)
    f2 <- (-1)*f1
    des <- rbind(f1,f2)
    ord<-std.ord(nrow(f1))
    des <- des[ord , ]
    colnames(des) <- c("A","B","C","D","E","F","G","H","J","K","L","M","N","O","P","Q")
    des <- data.frame(des)
  }

}
  
  nr <- nrow(des)
  ncd<-ncol(des)
  if(center>0){
    cpr<-c(rep(0,ncd))
    for (i in 1:center){
      des<-rbind(des,cpr)    
    }
   if(randomize==TRUE) {des <- des[sample(1:nr+center), ]}
  }
  if(center==0){
    if(randomize==TRUE) {des <- des[sample(1:nr), ]}
  }
    
return(des)
                                          }


####Interleve function
interleave <- function(v1,v2)
{
ord1 <- 2*(1:length(v1))-1
ord2 <- 2*(1:length(v2))
c(v1,v2)[order(c(ord1,ord2))]
}
### std.ord function
std.ord <-function(m) {
  v1 <- c(1:m)
  v2 <- c((m+1):(2*m))
  interleave(v1,v2)
                      }

Try the daewr package in your browser

Any scripts or data that you put into this service are public.

daewr documentation built on March 13, 2021, 3:01 a.m.

R Package Documentation

rdrr.io home R language documentation Run R code online

Browse R Packages

CRAN packages Bioconductor packages R-Forge packages GitHub packages

We want your feedback!

Note that we can't provide technical support on individual packages. You should contact the package authors for that.
GitHub issue tracker
ian@mutexlabs.com
Personal blog

 
Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.

Close