#' 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")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.