w.r.t. problem 2. The curvature of shape space, itself, introduces some
purely geometric variation that is picked up by PCA. For 2D, you expect
to lose 4 dimensions of variation. In fact, you get three exactly zero
eigenvalues and one very, very small one. The latter (usually, but I can
construct counter examples) reflecting the curvature of GPA space.
-ds
morphmet wrote:
Original Message
Subject: Understanding Morphologika
Date: Wed, 3 Mar 2010 09:41:33 +
From:
To:
Dear Colleagues
I have two unrelated issues, for which I don’t understand what
Morphologika is calculating exactly.
Firstly, I am having difficulties with the Procrustes Superimposition
performed in Morphologika. I was under the impression that scaling,
rotation, reflection and translation could be switched on and off and
that the subsequent principal component analysis would be performed on
the Procrustes coordinates as you instructed Morphologika to calculate
them.
Now my problem is the following:
If I perform Procrustes superimposition without scaling and perform PCA
afterwards or Procrustes superimposition with scaling and PCA
afterwards, the Procrustes coordinates are different, as would be
expected. However, the PC scores are almost identical (a few small
differences 5 decimals behind the comma). This is not what I would
expect. Since some of my specimens are twice the size of some others, I
would expect the first PC to show size (both isometric and allometric),
however, it is showing the same signal as the PCA on the full Procrustes
superimposition coordinates.
I am confused about what Morphologika calculates exactly to come to
these results.
And I have not been able to reproduce either of the PCA plots of
Morphologika with SPSS, even though I am forcing SPSS to use the
covariance matrix instead of the correlation matrix.
Any help or suggestions of what might be going on would be greatly
appreciated.
My second issue relates to the number of principal components calculated
by Morphologika.
In the help file it is stated: Principal components analysis of
specimens with k landmarks in m dimensions results in km-m-m(m-1)-1
eigenvectors; the principal components of variation of shape.
In my dataset I have 15 landmarks and 3 dimensions, so I think that
should result in 15*3-3-3*(3-1)-1=35 principal components. However,
Morphologika is giving me 38 principal components in the output. I don’t
understand the discrepancy and would appreciate it if anybody could
explain where the extra three principal components come from.
Thanking you all in advance.
Best wishes,
Anneke van Heteren
Anneke H. van Heteren
School of Human and Life Sciences
Roehampton University
Whitelands College
Holybourne Avenue
London SW15 4JD
Tel: +44 (0) 20 8392 3728
E-Mail: a.vanhete...@roehampton.ac.uk
Consider the environment. Please don't print this e-mail unless you
really need to.
Consider the environment. Please don't print this e-mail unless you
really need to.
This email and any attachments are confidential and intended solely for
the addressee and may also be privileged or exempt from disclosure under
applicable law. If you are not the addressee, or have received this
e-mail in error, please notify the sender immediately, delete it from
your system and do not copy, disclose or otherwise act upon any part of
this email or its attachments.
Internet communications are not guaranteed to be secure or virus-free.
Roehampton University does not accept responsibility for any loss
arising from unauthorised access to, or interference with, any Internet
communications by any third party, or from the transmission of any viruses.
Any opinion or other information in this e-mail or its attachments that
does not relate to the business of Roehampton University is personal to
the sender and is not given or endorsed by Roehampton University.
Roehampton University is a company limited by guarantee incorporated in
England under number 5161359. Registered Office: Grove House, Roehampton
Lane, London SW15 5PJ. An exempt charity.
--
Dennis E. Slice
Associate Professor
Dept. of Scientific Computing
Florida State University
Dirac Science Library
Tallahassee, FL 32306-4120
-
Guest Professor
Department of Anthropology
University of Vienna
-
Software worth having/learning/using...
Linux (Operating System: Ubuntu, CentOS, openSUSE, etc.)
OpenOffice (Office Suite: http://www.openoffice.org/)
R package (Stats/Graphics environment: http://www.r-project.org/)
Eclipse (Java/C++/etc IDE: http://www.eclipse.org/)