mosses<-c(53, 85, 69, 98, 130, 90)
liverworts<-c(47, 67, 85, 113, 190, 239)
plot(mosses)

plot(mosses, ylab="Number of records", xlab="Site number", pch=15,col="blue", main ="Altitudinal distribution of bryophytes")
points(liverworts,pch=16,col="brown")

plot(mosses, ylab="Number of records", xlab="Site number", pch=15,col="blue", main ="Altitudinal distribution of bryophytes",ylim=c(0,250))
points(liverworts,pch=16,col="brown")

plot(mosses, ylab="Number of records", xlab="Site number", pch=15,col="blue", main ="Altitudinal distribution of bryophytes", ylim=c(min(mosses,liverworts), max(mosses,liverworts)))
points (liverworts,pch=16,col="brown")

plot(c(1:6,1:6), c(mosses,liverworts), pch=c(rep(15,6), rep(16,6)), col=c(rep("blue",6) ,rep("brown",6)), ylab="Number of records", xlab="Site number",main ="Altitudinal distribution of bryophytes")
lines(c(1:6), mosses, col="blue", lty=1) 
lines(c(1:6), liverworts, col="brown", lty=2)


altitude<-c(400,600,800,1000,1200, 1300)

plot(c(altitude,altitude), c(mosses,liverworts), pch=c(rep(15,6), rep(16,6)), col=c(rep("blue",6) ,rep("brown",6)), ylab="Number of records", xlab="Altitude a.s.l. (m)",main ="Altitudinal distribution of bryophytes")
lines(altitude, mosses, col="blue", lty=1) 
lines(altitude, liverworts, col="brown", lty=2)
A<-"Liverworts"
B<-"Mosses"
legend("topright",pch=c(15,16),col=c("blue","brown"),c(A,B))

#repeat plot
plot(c(altitude,altitude), c(mosses,liverworts), pch=c(rep(15,6), rep(16,6)), col=c(rep("blue",6) ,rep("brown",6)), ylab="Number of records", xlab="Altitude a.s.l. (m)",main ="Altitudinal distribution of bryophytes")
lines(altitude, mosses, col="blue", lty=1) 
lines(altitude, liverworts, col="brown", lty=2)
legend("topleft",pch=c(15,16),col=c("blue","brown"),c(A,B))

#repeat plot
plot(c(altitude,altitude), c(mosses,liverworts), pch=c(rep(15,6), rep(16,6)), col=c(rep("blue",6) ,rep("brown",6)), ylab="Number of records", xlab="Altitude a.s.l. (m)",main ="Altitudinal distribution of bryophytes")
lines(altitude, mosses, col="blue", lty=1) 
lines(altitude, liverworts, col="brown", lty=2)
legend("topleft",pch=c(15,16), col=c("blue","brown"), c(A,B), cex=0.8,inset=0.05)


######
mosses<-c(53, 85, 69, NA, 130,  90)
liverworts<-c(47, 67, 85, 113, 190, 239)
plot(altitude,mosses,pch=15,col="blue",ylab="Number of bryophytes", xlab="Altitude a.s.l. (m)", main ="Altitudinal distribution of bryophytes",ylim=c(min(na.omit(mosses,liverworts)), max(mosses,liverworts, na.rm=TRUE)))
points(altitude,liverworts,pch=16,col="brown")
lines(altitude,mosses,col="blue",lty=1)
lines(altitude,liverworts,col="brown",lty=2)

######
#going back to our original data
mosses[4]<-98
moss_species<- c(12, 19, 20, 24, 27, 19)
liverwort_species<-c(17, 19, 23, 31, 42, 46)
bryos<-c(mosses,liverworts)

alt<-c(altitude,altitude) #  WHY THIS?
use_colours<-c(rep("blue",6),rep("brown",6))
use_symbols<-c(rep(15,6),rep(16,6))
plot(alt,bryos,pch=use_symbols,col=use_colours,ylab="Number of records", xlab="Altitude a.s.l. (m)",main ="Altitudinal distribution of bryophytes",ylim=c(min(bryos), max(bryos)))



plot(NULL,axes=F,ylab="",xlab="",main="", xlim=c(400,1300), ylim= c(min(c(moss_species, liverwort_species)), max(c(moss_species,liverwort_species))))
points(altitude,moss_species,pch=3,col="blue")
lines(altitude,moss_species,lty=3,col="blue") # using a different dotted line style
points(altitude,liverwort_species,pch=4,col="red")
lines(altitude,liverwort_species,lty=4,col="red") # using a different dotted line style
axis(4,c(10,20,30,40,50),at=c(10,20,30,40,50))
axis(4,las=0) 
mtext(side=4,line=2.5,"Numbers of species")