Hi,
We don't know anything about your data. Providing a small reproducible
example would go a long way toward making it possible to answer your
question. The help file for s.chull gives some fake data; perhaps you
could use that to more clearly describe what you are doing and what is
going wrong.
Or you could use dput() to give us a small portion of your actual data.
Sarah
On Wed, May 9, 2012 at 4:43 AM, Julian WITTISCHE <jwittis...@hotmail.com> wrote:
>
> Dear all,
> I have a problem to plot a convex hulls on PCA results. The PCA has reduced
> the 19 BIOCLIM variables to two axes and I want to plot a convex hull to be
> able to illustrate the environmental space overlap of two species.
> This is my script:
> library(Hmisc);library(ade4);
> #Import data
> batrabio<-read.table(file="Bioclimfull.txt", sep="\t",
> header=TRUE);batrabio<-batrabio[,c(1,4:22)]; #To select only species name
> and the 19 bioclim variables columns
> batrapair<-batrabio[batrabio$SCIENTIFIC_NAME=="Batrachoseps minor"|
> batrabio$SCIENTIFIC_NAME=="Batrachoseps incognitus",]; #To select two species
> to compare and illustrate their environmental space
> batrapair<-unique(batrapair[,]);#To select only unique localities
> #Run PCAs based on columns with climate data
> pca.ade4<-dudi.pca(batrapair[,2:20], center=T, scale=T, scannf=T); 2;#number
> of axes
> #Create convex hull using PCA axes, and colour-code points according to
> species name
> convex.hull<-s.chull(pca.ade4$li,as.factor(batrapair$SCIENTIFIC_NAME),add.plot=F,
> grid=F, clabel=0, origin=c(0,0),addaxes=T,col=c("red","blue"), cpoint=.65)
>
> I get the axes, the grid, the species name label but nothing else. I would
> greatly appreciate any advice as I have not been able to figure myself how to
> solve this problem.
> Thank you very much,
> Julian
>
--
Sarah Goslee
http://www.functionaldiversity.org
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
