# Commands to try out:

# To see list of files in working directory, type: 
dir() 

# To run the KS test to compare observed Jan maxT and model 3 of Jan MaxT:
ks.test(jan_maxT,jan_maxT_model3)

# To run the Student t-test to compare observed Jan maxT and model 3 of Jan MaxT:
t.test(jan_maxT,jan_maxT_model3)

# Some other plotting commands, try:
hist(jan_maxT)
boxplot(jan_maxT)
qqnorm(jan_maxT)

# For a boxplot including all 12 months, try (remember to change the names for different variables):
boxplot(jan_maxT,feb_maxT,mar_maxT,apr_maxT,may_maxT,jun_maxT,jul_maxT,aug_maxT,sep_maxT,oct_maxT,nov_maxT,dec_maxT,xlab="month",ylab="maxT",main="Plot name",names=c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"))

# To plot the empirical cumulative distribution function (CDF), try:

jan_maxT_c<-c(jan_maxT)
plot(ecdf(jan_maxT_c),do.points=FALSE,verticals=TRUE)

# To plot two CDFs together (e.g. jan_MaxT and feb_maxT), try:

jan_maxT_c<-c(jan_maxT)
feb_maxT_c<-c(feb_maxT)
plot(ecdf(jan_maxT_c),do.points=FALSE,verticals=TRUE)
plot(ecdf(feb_maxT_c),do.points=FALSE,verticals=TRUE,add=TRUE)