CheckingUpdateFomular.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)
library(MASS)
devtools::load_all()


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


vecC    <- c(1,1,1)
nDesign <- Design + vecC %*% t(vecC)
test    <- ginv(Design)



InvOld  <- ginv(Design)
InvNew  <- ginv(nDesign)
xxx1    <- (InvOld %*% vecC %*% t(vecC) %*% InvOld)
xxx2    <- 1+t(vecC) %*% InvOld %*% vecC
upFor   <- InvOld - xxx1/xxx2[1]
#round(10^10*test %*%Design)/10^10
orgFor  <- ginv(nDesign)