The lab notebook

Drawing completed by using SketchUp http://www.sketchup.com/ Maludam National Park in its glory, 3-D perception. Acknowledgement: Thank you Raini!

And to view multiple plots in one page, I referred to the below: > old.par <- par(mfrow=c(1, 2)) > plot(faithful, main="Faithful eruptions") > plot(large.islands, main="Islands", ylab="Area") > par(old.par) (Source: http://www.dummies.com/programming/r/how-to-put-multiple-plots-on-a-single-page-in-r/ ( inserted with my data: chns.par<-par(mfrow=c( 2 ,2)) barplot(as.matrix(Maludam.December.2014.Station.1CHNS), main="Elemental Analysis for Station 1", ylab= "Carbon Ratio", beside=TRUE, col=rainbow(5)) + legend("topright", c("Top","Middle","Bottom"), cex=0.6,bty="n", fill=rainbow(5)) barplot(as.matrix(Maludam.December.2014.Station.2CHNS), main="Elemental Analysis for Station 2", ylab= "Carbon Ratio", beside=TRUE, col=rainbow(5)) + legend("topright", c("Top","Middle","Bottom"), cex=0.6,bty="n", f...

Correlation.R # ============================================================ # Tutorial on drawing a correlation map using ggplot2 # by Umer Zeeshan Ijaz (http://userweb.eng.gla.ac.uk/umer.ijaz) # =============================================================   abund_table<- read.csv ( "SPE_pitlatrine.csv" , row.names = 1 , check.names= FALSE ) #Transpose the data to have sample names on rows abund_table<- t ( abund_table ) meta_table<- read.csv ( "ENV_pitlatrine.csv" , row.names = 1 , check.names= FALSE )   #Filter out samples with fewer counts abund_table<-abund_table [ rowSums ( abund_table ) > 200 , ]   #Extract the corresponding meta_table for the samples in abund_table meta_table<-meta_table [ rownames ( abund_table ) , ]   I changed the parameters #You can use sel_env to specify the variables you want to use and sel_env_label to specify the labes for the pannel sel_env<- c ( "pH" , "Temp" , "TS...

So this works for my data! I had started with the commands under the section NMDS from http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/ecological.html the sample and data samples were prepared according to the above. and proceed with the following: library("vegan", lib.loc="~/R/win-library/3.2") library("lattice", lib.loc="C:/Program Files/R/R-3.2.5/library") set.seed(2) abund_table<-read.csv("Sample.csv",row.names=1,check.names=FALSE) #Transpose the data to have sample names on rows abund_table<-t(abund_table) abund_table<-t(abund_table) example_NMDS=metaMDS(abund_table, k=2) plot(example_NMDS) ordiplot(example_NMDS,type="n") orditorp(example_NMDS,display="species",col="red",air=0.01) orditorp(example_NMDS,display="sites",cex=1.25,air=0.01) And I changed "Treatment 1" and "Treatment 2" to: treat=c(rep("pH",5),rep("Cond...