# \donttest{
da<-"https://github.com/pasraia/RRmorph_example_data/raw/refs/heads/main/RRmorphdata.rda"
download.file(url=da,destfile = paste0(tempdir(),"/RRmorphdata.rda"))
load(paste0(tempdir(),"/RRmorphdata.rda"))
require(Morpho)
pca<-procSym(endo.set)
ldm_homo<-endo.set[,,"Homo_sapiens"]
sur_homo<-endo.sur[["Homo_sapiens"]]
ldm_macaca<-endo.set[,,"Macaca_fuscata"]
sur_macaca<-endo.sur[["Macaca_fuscata"]]
# Convergence within group plotted on reconstructed surfaces
cm1<-conv.map(x1=c("Pan_troglodytes","Gorilla_gorilla","Pongo_abelii"),
scores=pca$PCscores,pcs=pca$PCs,mshape=pca$mshape,
focal=c("Pan_troglodytes","Gorilla_gorilla"))
# Convergence between group plotted on reconstructed surfaces
cm2<-conv.map(x1=c("Pongo_abelii"),x2=c("Alouatta_caraya"),
scores=pca$PCscores,pcs=pca$PCs,mshape=pca$mshape,
focal="Alouatta_caraya")
# Convergence within group plotted on real surfaces
cm3<-conv.map(x1=c("Homo_sapiens","Gorilla_gorilla","Pongo_abelii"),
scores=pca$PCscores,pcs=pca$PCs,mshape=pca$mshape,
refsur=list("Homo_sapiens"=sur_homo),
refmat=list("Homo_sapiens"=ldm_homo))
# Convergence between group plotted on real surfaces
cm3<-conv.map(x1=c("Homo_sapiens","Pongo_abelii"),x2=c("Macaca_fuscata"),
scores=pca$PCscores,pcs=pca$PCs,mshape=pca$mshape,
refsur=list("Homo_sapiens"=sur_homo,"Macaca_fuscata"=sur_macaca),
refmat=list("Homo_sapiens"=ldm_homo,"Macaca_fuscata"=ldm_macaca))
# }
Run the code above in your browser using DataLab