twd <- getwd()
source('BuildEmulator/BuildEmulator.R')
source('HistoryMatching/HistoryMatching.R')
load("Demonstrations/globalCanAM4W1.RData")
head(tData)
cands <- names(tData)[2:14]
Noise <- rnorm(dim(tData)[1],0,0.5)
tData <- cbind(tData,Noise)
cands <- c(cands, "Noise")
myem.lm.canTOA <- EMULATE.lm(Response="BALX", tData=tData, tcands=cands,tcanfacs=NULL,TryFouriers=TRUE,maxOrder=2,maxdf = 20)
myem.gp.canTOA = EMULATE.gpstan(meanResponse=myem.lm.canTOA, tData=tData, additionalVariables=NULL, FastVersion = TRUE)
tLOOs <- LOO.plot(StanEmulator = myem.gp.canTOA, ParamNames = myem.gp.canTOA$Names)
myem.lm.canPCP <- EMULATE.lm(Response="PCP", tData=tData, tcands=cands,tcanfacs=NULL,TryFouriers=TRUE,maxOrder=2,maxdf = 20)
myem.gp.canPCP = EMULATE.gpstan(meanResponse=myem.lm.canPCP, tData=tData, additionalVariables=NULL, FastVersion = TRUE)
tLOOs <- LOO.plot(StanEmulator = myem.gp.canPCP, ParamNames = myem.gp.canPCP$Names)