-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathanalysis
75 lines (55 loc) · 3.64 KB
/
analysis
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
63
64
65
66
67
68
69
70
71
72
73
74
75
rm(list=ls())
library(ggplot2)
library(Rmisc)
setwd(dir = "Masters/DIS/Results/")
SeaInvaders <- read.table("jobs_to_run_ForR2.csv", sep=",", header=TRUE)
###########################################################################################################################
Pvol_Atl_Carib <- SeaInvaders[SeaInvaders$Model==1,]
ggplot(data=Pvol_Atl_Carib,aes(x=MostProbColonists,y=RPI,colour=as.factor(Theta))) +
geom_errorbar(aes(ymin=RPI-Lower50CI, ymax=RPI+Upper50CI,width=0)) +
ggtitle(label = "Pvol_Atl_Carib_RPIvMP") +
theme(legend.text=element_text(colour="#000000",size=18)) +
theme(legend.title=element_text(colour="#000000",size=18)) +
theme(axis.title.x=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.title.y=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.text.y=element_text(face="bold",colour="#000000",size=12)) +
theme(axis.text.x=element_text(face="bold",colour="#000000",size=12)) +
geom_point()
###########################################################################################################################
Pvol_Indo_Atl <- SeaInvaders[SeaInvaders$Model==2,]
ggplot(data=Pvol_Indo_Atl,aes(x=MostProbColonists,y=RPI,colour=as.factor(Theta))) +
geom_errorbar(aes(ymin=RPI-Lower50CI, ymax=RPI+Upper50CI,width=0)) +
ggtitle(label = "Pvol_Indo_Carib_RPIvMP") +
theme(legend.text=element_text(colour="#000000",size=18)) +
theme(legend.title=element_text(colour="#000000",size=18)) +
theme(axis.title.x=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.title.y=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.text.y=element_text(face="bold",colour="#000000",size=12)) +
theme(axis.text.x=element_text(face="bold",colour="#000000",size=12)) +
geom_point()
###########################################################################################################################
Lkas_Marq_Haw <- SeaInvaders[SeaInvaders$Model==3,]
ggplot(data=Lkas_Marq_Haw,aes(x=RPI,y=MostProbColonists,colour=as.factor(Theta))) +
geom_errorbar(aes(ymin=MostProbColonists-Lower50CI, ymax=MostProbColonists+Upper50CI),width=0) +
ggtitle(label = "Lkas_Marq_Haw_RPIvMP") +
theme(legend.text=element_text(colour="#000000",size=18)) +
theme(legend.title=element_text(colour="#000000",size=18)) +
theme(axis.title.x=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.title.y=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.text.y=element_text(face="bold",colour="#000000",size=12)) +
theme(axis.text.x=element_text(face="bold",colour="#000000",size=12)) +
geom_point() +
stat_smooth(method = 'nls', formula = 'y~a*x^b' , method.args = list(start= c(a = 200,b=-0.3)),se=FALSE)
###########################################################################################################################
Lkas_Soc_Haw <- SeaInvaders[SeaInvaders$Model==4,]
ggplot(data=Lkas_Soc_Haw,aes(x=RPI,y=MostProbColonists,colour=as.factor(Theta))) +
geom_errorbar(aes(ymin=MostProbColonists-Lower50CI, ymax=MostProbColonists+Upper50CI),width=0) +
ggtitle(label = "Lkas_Soc_Haw_RPIvMP") +
theme(legend.text=element_text(colour="#000000",size=18)) +
theme(legend.title=element_text(colour="#000000",size=18)) +
theme(axis.title.x=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.title.y=element_text(face="bold",colour="#000000",size=24)) +
theme(axis.text.y=element_text(face="bold",colour="#000000",size=12)) +
theme(axis.text.x=element_text(face="bold",colour="#000000",size=12)) +
geom_point() +
stat_smooth(method = 'nls', formula = 'y~a*x^b', method.args = list(start= c(a = 200,b=-0.3)),se=FALSE)