TestingCodes.R

In nhavt/Block3AugDesigns:

#' Stuttgart 06 April 2017
#' @author Nha Vo-Thanh
#' reading data from text files

rm(list=ls())
setwd("D:/user/vthanh/Google Drive/Project Stuttgart/R_programing/UpdateFormular/Block3AugDesigns/")
library(devtools)
library(roxygen2)
devtools::load_all()


data    <- split(read.table("./Designs/TestDesign.txt"), gl(1, 3))
repli   <- 3
Design <- data$`1`

#adding observation, and deleting observation

nTreatments <- 3
nRows       <- 3
nCols       <- 3
irow        <- 3 
icol        <- 1
itrt        <- 1

# information matrix of an orginal design

infoMat <- infoMatrix(Design,nTreatments, nRows, nCols)
Cn      <- infoMat$Cn
AnBn    <- infoMat$AnBn
InvCO   <- infoMat$InvCn

#deleting a treatment#
dlist  <- delInfoMatrix(Design,nTreatments, nRows, nCols, irow, icol)
#adding a treatment 
alist <- addInfoMatrix(Design,nTreatments, nRows, nCols, irow, icol, itrt)

Cn  <- alist$Cn
Cn1 <- alist$Cn1
c2  <- alist$c2
Cnx <- Cn1 + c2 %*% t(c2)

InvCn1 <- ginv(Cn1)
InvCn  <- ginv(Cnx)
zzz     <- c2 %*% t(c2)
zzzstar <- InvCn1%*%Cn1%*%zzz



xxx1   <- (InvCn1 %*% c2 %*% t(c2) %*% InvCn1)
xxx2   <- 1+t(c2) %*% InvCn1 %*% c2
InvCn1 - xxx1/xxx2[1]

#---------------------------------------------#
# replace a new treatment
#Design[irow,icol] <- itrt

#infoMat <- infoMatrix(Design,nTreatments, nRows, nCols)
# comparing 
#Cnew    <- infoMat$invCn
#Cupdate <- alist$invCn
cat("\014")