# evoTS v.1.0.2 on CRAN was used in all analyses. library(evoTS) ### Simulating time series using an multivariate URW model simulation_symmetric.R_n10<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(10, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n10[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n20<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(20, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n20[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n30<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(30, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n30[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n40<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(40, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n40[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n50<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(50, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n50[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n60<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(60, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n60[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n70<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(70, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n70[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } simulation_symmetric.R_n80<-list() for(i in 1:1000) { tryCatch({ sim.data<-sim.multi.URW(80, R = matrix(c(0.5,0.25,0.25,0.5), nrow=2, byrow=2), anc=c(0,0), vp=0.0001) simulation_symmetric.R_n80[[i]]<-fit.multivariate.URW(sim.data, R = "symmetric", trace = TRUE) print(i) }, error=function(e){cat("ERROR :",conditionMessage(e), "\n")}) } ### prepare and print output R_1_10<-rep(NA, 1000) R_2_10<-rep(NA, 1000) R_1_2_10<-rep(NA, 1000) R_1_20<-rep(NA, 1000) R_2_20<-rep(NA, 1000) R_1_2_20<-rep(NA, 1000) R_1_30<-rep(NA, 1000) R_2_30<-rep(NA, 1000) R_1_2_30<-rep(NA, 1000) R_1_40<-rep(NA, 1000) R_2_40<-rep(NA, 1000) R_1_2_40<-rep(NA, 1000) R_1_50<-rep(NA, 1000) R_2_50<-rep(NA, 1000) R_1_2_50<-rep(NA, 1000) R_1_60<-rep(NA, 1000) R_2_60<-rep(NA, 1000) R_1_2_60<-rep(NA, 1000) R_1_70<-rep(NA, 1000) R_2_70<-rep(NA, 1000) R_1_2_70<-rep(NA, 1000) R_1_80<-rep(NA, 1000) R_2_80<-rep(NA, 1000) R_1_2_80<-rep(NA, 1000) for (i in 1:1000){ R_1_10[i]<-simulation_symmetric.R_n10[[i]]$R[1,1] R_2_10[i]<-simulation_symmetric.R_n10[[i]]$R[2,2] R_1_2_10[i]<-simulation_symmetric.R_n10[[i]]$R[1,2] R_1_20[i]<-simulation_symmetric.R_n20[[i]]$R[1,1] R_2_20[i]<-simulation_symmetric.R_n20[[i]]$R[2,2] R_1_2_20[i]<-simulation_symmetric.R_n20[[i]]$R[1,2] R_1_30[i]<-simulation_symmetric.R_n30[[i]]$R[1,1] R_2_30[i]<-simulation_symmetric.R_n30[[i]]$R[2,2] R_1_2_30[i]<-simulation_symmetric.R_n30[[i]]$R[1,2] R_1_40[i]<-simulation_symmetric.R_n40[[i]]$R[1,1] R_2_40[i]<-simulation_symmetric.R_n40[[i]]$R[2,2] R_1_2_40[i]<-simulation_symmetric.R_n40[[i]]$R[1,2] R_1_50[i]<-simulation_symmetric.R_n50[[i]]$R[1,1] R_2_50[i]<-simulation_symmetric.R_n50[[i]]$R[2,2] R_1_2_50[i]<-simulation_symmetric.R_n50[[i]]$R[1,2] R_1_60[i]<-simulation_symmetric.R_n60[[i]]$R[1,1] R_2_60[i]<-simulation_symmetric.R_n60[[i]]$R[2,2] R_1_2_60[i]<-simulation_symmetric.R_n60[[i]]$R[1,2] R_1_70[i]<-simulation_symmetric.R_n70[[i]]$R[1,1] R_2_70[i]<-simulation_symmetric.R_n70[[i]]$R[2,2] R_1_2_70[i]<-simulation_symmetric.R_n70[[i]]$R[1,2] R_1_80[i]<-simulation_symmetric.R_n80[[i]]$R[1,1] R_2_80[i]<-simulation_symmetric.R_n80[[i]]$R[2,2] R_1_2_80[i]<-simulation_symmetric.R_n80[[i]]$R[1,2] } setEPS() postscript("multivariate_URW_1.eps") par(mfrow=c(4,3)) hist(R_1_10, breaks = 40, col="gray93", main="10", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_10), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_10), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.46"),bty="n") mean(R_1_10) median(R_1_10) hist(R_2_10, breaks = 40, col="gray93", main="10", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_10), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_10), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.51", "0.45"),bty="n") mean(R_2_10) median(R_2_10) hist(R_1_2_10, breaks = 40, col="gray93", main="10", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_10), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_10), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.23"),bty="n") mean(R_1_2_10) median(R_1_2_10) hist(R_1_20, breaks = 40, col="gray93", main="20", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_20), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_20), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.48"),bty="n") mean(R_1_20) median(R_1_20) hist(R_2_20, breaks = 40, col="gray93", main="20", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_20), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_20), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50, 0.50", "0.48"),bty="n") mean(R_2_20) median(R_2_20) hist(R_1_2_20, breaks = 40, col="gray93", main="20", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_20), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_20), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.24"),bty="n") mean(R_1_2_20) median(R_1_2_20) hist(R_1_30, breaks = 40, col="gray93", main="30", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_30), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_30), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_1_30) median(R_1_30) hist(R_2_30, breaks = 40, col="gray93", main="30", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_30), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_30), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_2_30) median(R_2_30) hist(R_1_2_30, breaks = 40, col="gray93", main="30", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_30), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_30), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.25"),bty="n") mean(R_1_2_30) median(R_1_2_30) hist(R_1_40, breaks = 40, col="gray93", main="40", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_40), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_40), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_1_40) median(R_1_40) hist(R_2_40, breaks = 40, col="gray93", main="40", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_40), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_40), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_2_40) median(R_2_40) hist(R_1_2_40, breaks = 40, col="gray93", main="40", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_40), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_40), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.25"),bty="n") mean(R_1_2_40) median(R_1_2_40) dev.off() setEPS() postscript("multivariate_URW_2.eps") par(mfrow=c(4,3)) hist(R_1_50, breaks = 40, col="gray93", main="50", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_50), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_50), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50","0.50","0.49"),bty="n") mean(R_1_50) median(R_1_50) hist(R_2_50, breaks = 40, col="gray93", main="50", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_50), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_50), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.50"),bty="n") mean(R_2_50) median(R_2_50) hist(R_1_2_50, breaks = 40, col="gray93", main="50", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_50), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_50), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.24"),bty="n") mean(R_1_2_50) median(R_1_2_50) hist(R_1_60, breaks = 40, col="gray93", main="60", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_60), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_60), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_1_60) median(R_1_60) hist(R_2_60, breaks = 40, col="gray93", main="60", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_60), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_60), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_2_60) median(R_2_60) hist(R_1_2_60, breaks = 40, col="gray93", main="60", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_60), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_60), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.24"),bty="n") mean(R_1_2_60) median(R_1_2_60) hist(R_1_70, breaks = 40, col="gray93", main="70", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_70), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_70), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.51", "0.50"),bty="n") mean(R_1_70) median(R_1_70) hist(R_2_70, breaks = 40, col="gray93", main="70", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_70), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_70), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.50"),bty="n") mean(R_2_70) median(R_2_70) hist(R_1_2_70, breaks = 40, col="gray93", main="70", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_70), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_70), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.26", "0.26"),bty="n") mean(R_1_2_70) median(R_1_2_70) hist(R_1_80, breaks = 40, col="gray93", main="80", xlab="R1") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_80), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_80), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.49"),bty="n") mean(R_1_80) median(R_1_80) hist(R_2_80, breaks = 40, col="gray93", main="80", xlab="R2") abline(v=0.5, col="Black", lty=1, lwd=2.5) abline(v=mean(R_2_80), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_2_80), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.50", "0.50", "0.50"),bty="n") mean(R_2_80) median(R_2_80) hist(R_1_2_80, breaks = 40, col="gray93", main="80", xlab="R1.2") abline(v=0.25, col="Black", lty=1, lwd=2.5) abline(v=mean(R_1_2_80), col="darkgreen", lty=2, lwd=2.5) abline(v=median(R_1_2_80), col="royalblue1", lty=3, lwd=2.5) legend("topright", legend = c("0.25", "0.25", "0.25"),bty="n") mean(R_1_2_80) median(R_1_2_80) dev.off()