Examples.R

In nhavt/Block3AugDesigns:

#' Stuttgart 05 April 2017
#' @author Nha Vo-Thanh
#' reading data from text files
#'
rm(list=ls())
setwd("D:/user/vthanh/Google Drive/Project Stuttgart/R_programing/")
library("Block3AugDesigns")
require(Matrix)
library(MASS)
#data <- split(read.table("D:/user/vthanh/Google Drive/Project Stuttgart/R_programing/3blocks/Example3_1.txt"), gl(1, 22))#
#data <- split(read.table("D:/user/vthanh/Google Drive/Project Stuttgart/R_programing/3blocks_Search1_4/Example3_S4.txt"), gl(1, 22))
data  <- split(read.table("UpdateFormular/Block3AugDesigns/Designs/Example_with_3_blocks.txt"), gl(1, 8))


pDesign <- data$`1`

nTreatments <- 144
rowgroup    <- 1
colgroup    <- 3
nRows       <- 8
nCols       <- 24
repli       <- -1
ncheck      <- 24
flag        <- 0
subrgroup   <- c(1) # beta
subcgroup   <- c(8,8,8) # alpha
modelD      <- model3Matrix(nTreatments, nRows, nCols, rowgroup, colgroup, subrgroup, subcgroup, pDesign)
eff         <- A3BlkDesign(pDesign, nTreatments, nRows, nCols, rowgroup, colgroup, subrgroup, subcgroup, repli, ncheck, flag, cRandom=0, reps=1)
cat("\f")
eff