Simulation/simulation5_2/1/sigma_DDE/05/main.R

In shijiaw/smcDE: Adaptive semiparametric Bayesian differential equations via sequential Monte Carlo methods

rm(list=ls())
library("fda")
library("truncnorm")
library(numDeriv)
library(fda)
library(deSolve)
library(NLRoot)
library(Hmisc)

NITER <- 12
times <- seq(0, 100, 0.5)
knots = seq(min(times), max(times),by=2)
sigma1 <- 0.5
sigma1_simulator <- sigma1
seeds <- 1:NITER


alambda <- 1
blambda <- 1
sigmac <- 0.3
CESSthresholds <- 0.9
NP <- 500

filename_nu = './nu'
filename_P = './P'
filename_tau = './tau'
filename_W = './W'
filename_sigma = './sigma'
filename_MSE = './MSE'
filename_lambda = './lambda'

filename_nuCI = './nuCI'
filename_PCI = './PCI'
filename_tauCI = './tauCI'
filename_WCI = './WCI'
filename_sigmaCI = './sigmaCI'

filename_iter = './iterations'

ifAppend=FALSE;
for(iter in 1:NITER){
  if(iter == 1){
    ifAppend=FALSE
  }else{
    ifAppend=TRUE
  }
  set.seed(seeds[iter])
  source('mainDDE.R')
  
  write.table(matrix(round(mean_r, digits = 5),nr=1), file = paste(filename_nu,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(mean_P, digits = 5),nr=1), file = paste(filename_P,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(mean_tau, digits = 5),nr=1), file = paste(filename_tau,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(meansigma, digits = 5),nr=1), file = paste(filename_sigma,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(lambda, digits = 5),nr=1), file = paste(filename_lambda,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(sqrt(sum((log(output[,2]) - fitted1)^2)/length(fitted1)), digits = 5),nr=1), file = paste(filename_MSE,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(cMean[1], digits = 5),nr=1), file = paste(filename_W,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  
  write.table(matrix(round(r_Credible, digits = 5),nr=1), file = paste(filename_nuCI,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(P_Credible, digits = 5),nr=1), file = paste(filename_PCI,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(tau_Credible, digits = 5),nr=1), file = paste(filename_tauCI,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(sigma_Credible, digits = 5),nr=1), file = paste(filename_sigmaCI,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(W_Credible, digits = 5),nr=1), file = paste(filename_WCI,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);
  write.table(matrix(round(iter, digits = 5),nr=1), file = paste(filename_iter,".txt",sep=""), append = ifAppend, row.names = FALSE,
              col.names = FALSE);

}