source("met.r") cat("\nDATA1:\n\n") data1 <- as.matrix(read.table("data1")) s1 <- c(0.4,0.8,0.1,0.1,0.1) set.seed(1) res1.0 <- met (20000, s1, data1, rep(0,5)) set.seed(2) res1.1 <- met (20000, s1, data1, rnorm(5)) res1.2 <- met (20000, s1, data1, rnorm(5)) res1.3 <- met (20000, s1, data1, rnorm(5)) thin1 <- seq(2001,20000,length=1000) pdf("plots1.pdf") plot(as.data.frame(res1.0[thin1,]),pch=".") dev.off() cat("\nDATA2:\n\n") data2 <- as.matrix(read.table("data2")) s2 <- c(0.5,0.6,0.2,0.1,0.1) set.seed(1) res2.0 <- met (500000, s2, data2, rep(0,5)) set.seed(2) res2.1 <- met (500000, s2, data2, rnorm(5)) res2.2 <- met (500000, s2, data2, rnorm(5)) res2.3 <- met (500000, s2, data2, rnorm(5)) thin2 <- seq(50001,500000,length=1000) pdf("plots2.pdf") plot(as.data.frame(res2.0[thin2,]),pch=".") dev.off() cat("\nDATA3:\n\n") data3 <- as.matrix(read.table("data3")) s3 <- c(0.5,0.5,0.4,0.3,0.3) set.seed(1) res3.0 <- met (20000, s3, data3, rep(0,5)) set.seed(2) res3.1 <- met (20000, s3, data3, rnorm(5)) res3.2 <- met (20000, s3, data3, rnorm(5)) res3.3 <- met (20000, s3, data3, rnorm(5)) thin3 <- seq(2001,20000,length=1000) pdf("plots3.pdf") plot(as.data.frame(res3.0[thin1,]),pch=".") dev.off()