Stuttgart_April_24_2017.R

#' check cnonected designs
rm(list = ls())
library("rJava")
library(xlsx)

library(devtools)
library(roxygen2)
library(MASS)
library(Matrix)
setwd(
  "D:/user/vthanh/Google Drive/Project Stuttgart/R_programing/UpdateFormular/Block3AugDesigns/"
)
devtools::load_all()
data    <- split(read.table("D:/user/vthanh/Google Drive/Project Stuttgart/R_programing/2BlocksSearch/Example_two_blocking_factor_search1.txt"), gl(1, 10))
nRows   <- 10
nCols   <- 11
pDesign <- data$`1`

arcopt <- aRCDesign(Design= pDesign, nTreatments=91, nRows=nRows, nCols=nCols, repli=-1, ncheck=3, flag=0, cRandom=0, reps=1)

cat("\f")
arcopt
#ll <- setdiff(c(1:580),c(1,2,3,4))
#for(ii in ll){
#  pDesign[pDesign==ii] <- NA
#}
#write.xlsx(pDesign, "pDesign.xlsx")
nhavt/Block3AugDesigns documentation built on May 7, 2019, 11:15 a.m.