PCA plotting?

316 views
Skip to first unread message

Candice Neves

unread,
Jan 21, 2018, 1:06:43 PM1/21/18
to geomorph R package
Hi All

Part of my MSc has me looking at shape change in molars using EFA and the analysis portion was conducted in Momocs, but I much prefer the visualisations for PCAs in geomorph. Is there possibly a way to create a PCA of my molar PC axes using geomorph? I'd like to maintain consistency in the way my figures look, so any help on this would be very much appreciated

Thanks
Candice Neves


Antigoni Kaliontzopoulou

unread,
Jan 22, 2018, 5:39:12 AM1/22/18
to geomorph-...@googlegroups.com
Dear Candice,

PCA and related graphical support is easily implemented in geomorph using the function plotTangentSpace. The function allows you to calculate the PCA, plot a scatterplot of PCA scores, and visualize associated shape changes. For a more detailed account, you can see the function´s help file, or this short tutorial:


All the best,
Antigoni

--
You received this message because you are subscribed to the Google Groups "geomorph R package" group.
To unsubscribe from this group and stop receiving emails from it, send an email to geomorph-r-package+unsub...@googlegroups.com.
To post to this group, send email to geomorph-r-package@googlegroups.com.
Visit this group at https://groups.google.com/group/geomorph-r-package.
To view this discussion on the web, visit https://groups.google.com/d/msgid/geomorph-r-package/4b37004d-8e42-44b0-826f-6716fc1115d9%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.



--


Antigoni Kaliontzopoulou

Assistant Researcher
CIBIO/InBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos
University of Porto
Campus Agrário de Vairão, 4485-661 Vairão
PORTUGAL

http://antigonik.wix.com/akaliontzopoulou
http://cibio.up.pt/people/details/akaliont
http://aphylogenetics.wix.com/phylogenetics

Candice Neves

unread,
Jan 22, 2018, 5:54:34 AM1/22/18
to geomorph-...@googlegroups.com, Antigoni Kaliontzopoulou
Hi Antigoni

I’ve been playing around with the plotTangentSpace function, but my problem is that I’m working with outline points. I can’t simply treat the outlines as landmarks, and if I did the PCAs come out looking completely different (see attached examples). I’m wondering if there is a way to use the PCs from my momocs analysis in a PCA visualisation in Geomorph?

Thanks
Candice


On 22 Jan 2018, at 12:39, Antigoni Kaliontzopoulou <ant...@gmail.com> wrote:


Antigoni Kaliontzopoulou

unread,
Jan 22, 2018, 6:24:37 AM1/22/18
to geomorph-...@googlegroups.com
Sorry, Candice, I misunderstood your question. 

geomorph works with landmark shape data. One possibility would be to do a subsampling of your outline data to produce semilandmarks, using the digit.curves() function of geomorph, and possibly complement these with some fixed landmarks, which could then be used in a landmark-based shape analysis. In that case, you would need to superimpose your data using gpagen() with the curves argument, and then do the PCA using plotTangentSpace. But this is an analytically completely different solution compared to outline analysis using EFA.

Alternatively, the PCA and scatter-plotting part of plotTangentSpace are both based on the prcomp() base function of R, so producing plots with your EFA-treated data should be straightforward by setting uniform plot() arguments in R.

Hope this helps,
Antigoni
Reply all
Reply to author
Forward
0 new messages