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...
Hi! It has been two years I haven't updated this blog ever since Azizi was admitted to the hospital for schizophrenia. Now, he's back to school and is determined more than ever to complete his PhD in Fluid Dynamics (hooray!). If you're reading this Zee, I am above the clouds for this great news and I wish you good luck! No.2 on the agenda, is to update all the R codings and lab protocols that I used during my PhD into this humble blog. I was really busy with thesis writing and attending conferences. Now, I can dedicate my time to contribute my knowledge and experience here and hopefully be of use. Watch this space guys! Love lots, Nas Denis