R-scripts/plot_simulation.R at master · barcaroli/R-scripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
est <- read.csv2("simul_results 1.csv",dec=".")
write.table(est,"simul_results 1.txt",row.names=F,sep="\t",quote=F)
est <- est[1:15,]

maxi = (max(c(est$n1),(est$n2),(est$n3))) + 100
plot(est$n1, col = "red", lty = 4, type = "b", pch = "*", cex = 1.3,cex.lab = 1.5,
     xlab = "Heteroscedasticity",
     ylab = "Sample size",
     ylim = c(0,maxi), axes = F)
points(est$n2, col = "black", lty = 4, type = "b", pch = "+", cex = 1.3)
points(est$n3, col = "blue", lty = 4, type = "b", pch = "o", cex = 1)

abline(0,0)
axis(1, at=1:nrow(est), labels = est$gamma,las=2)
axis(2, at=seq(0,maxi+10,50), labels=seq(0,maxi+10,50), las=1)
legend("topleft",
       legend = c("Linear model","Kriging","Spatial Linear Model"),
       pch = c("*","+","o"),
       col = c("red","black","blue"),
       ncol = 1, cex = 1.2, text.font = 1)

# title("Performance in terms of sampling size",cex.main=1.5)

maxi = (max(c(est$cv1),(est$cv2b),(est$cv3)))+0.02
plot(est$cv1, col = "red", lty = 4, type = "b", pch = "*", cex = 1.3,cex.lab = 1.3,
     xlab = "Heteroscedasticity",
     ylab = "CV",
     ylim = c(0,maxi+0.01), axes = F)
points(est$cv2, col = "black", lty = 4, type = "b", pch = "+", cex = 1.3)
points(est$cv3, col = "blue", lty = 4, type = "b", pch = "o", cex = 1)

abline(0,0)
abline(a=0.03, b=0.0,col="blue")
axis(1, at=1:nrow(est), labels = est$gamma,las=2)
axis(2, at=seq(0,maxi+10,0.01), labels=seq(0,maxi+10,0.01), las=1)
legend("topleft",
       legend = c("Linear model","Kriging","Spatial Linear Model"),
       pch = c("*","+","o"),
       col = c("red","black","blue"),
       ncol = 1, cex = 1.2, text.font = 1)

# title("Performance in terms of CV on target variable",cex.main=1.5)

maxi = (max(c(est$cv1),(est$cv2a),(est$cv3a))) + 0.01
plot(est$cv1, col = "red", lty = 4, type = "b", pch = "*", cex = 1.3,cex.lab = 1.3,
     xlab = "Heteroscedasticity",
     ylab = "Equalized CV",
     ylim = c(0,maxi+0.01), axes = F)
points(est$cv2a, col = "black", lty = 4, type = "b", pch = "+", cex = 1.3)
points(est$cv3a, col = "blue", lty = 4, type = "b", pch = "o", cex = 1)

abline(0,0)
abline(a=0.03, b=0.0,col="blue")
axis(1, at=1:nrow(est), labels = est$gamma,las=2)
axis(2, at=seq(0,maxi+10,0.01), labels=seq(0,maxi+10,0.01), las=1)
legend("topleft",
       legend = c("Linear model","Kriging","Spatial Linear Model"),
       pch = c("*","+","o"),
       col = c("red","black","blue"),
       ncol = 1, cex = 1.2, text.font = 1)

# title("Performance in terms of equalized CV on target variable",cex.main=1.2)