search for: geopca

...the following script to obtain a PCA plot but the end result is rather disappointing as the points are very very small and there are no titles etc geochemdata<-read.csv(file.choose(),header=TRUE) names(geochemdata) library(vegan) bstick<-function(n, tot.var=1) rev(cumsum(tot.var/n:1)/n) geopca<-rda(geochemdata, scale=TRUE) geopca geopca$CA$eig/geopca$tot.chi summary(geopca, scaling=2) apply(scores(geopca, choices=1.6, display= "species"), + 2,function(x) x/sd(x)) plot(geopca$CA$eig, type= "o", col= "red", xlab= "PCA Axis", ylab= "Var...