I have morphological data from two different determined groups (It and Nd), where the variables are heterogeneous (continuous, semi-qualitative, binomial). I want to know if the groups differ morphologically and if so, which variables account most for the difference between them.

For this I thought using gowdis() {FD} and then performing a Principal Coordinate analysis (PCoA, which is equivalent to a MDS) using cmdscale() {MASS} for my analysis. The use of a Gower distance should allow me to use the complete data set, with all the different kind of variables. Most of the examples found, so also the R-help, concerned ecological data (like dune and dune.env).

However, I have following questions:

  1. Is my basic thinking correct? Somehow, I would have expected to use some kind of constrained multivariate analysis.

  2. I have quite a lot of missing values, which are not only linked to a specific sample or variable (a part of my data set is from a previous study). How can I handle it? From my reading, I thought that Gower distance can handle missing values but see 3)

  3. How can I manage the warning missing species scores? Because of this, I cannot plot the centroids of my variables nor create an arrow showing the impact of variables. As far as I know, this issue is linked to the missing values.

  • 1
    $\begingroup$ Although nobody can prohibit you doing PCoA with Gower metric, it is not warranted from the view of geometrical correctness. Gower isn't euclidean or even metric distance. PCoA is basically a metric MDS based on PCA; it is appropriate mainly to euclidean distance. I may recommend you to use nonmetric MDS with Gower: you are likely to get better, more interpretable results. $\endgroup$ – ttnphns Dec 5 '13 at 14:05
  • 1
    $\begingroup$ Gower "handles" missing data in a manner similar to "pairwise deletion" (as far as I know). If you have many missing values that will distort the matrix of distances. Consider doing some imputation (for example hot-deck imputation) of the missing. $\endgroup$ – ttnphns Dec 5 '13 at 14:09
  • $\begingroup$ But I was not questionning your intention of doing MDS. It might be that MDS is not what you really need. $\endgroup$ – ttnphns Dec 5 '13 at 14:11
  • $\begingroup$ Hello all! Thanks for these answers. So maybe you can help me...I'm pretty new in R and far from beeing an expert with multivariate analysis. $\endgroup$ – Bettina Dec 9 '13 at 9:18
  • $\begingroup$ 1) I can understand your remark about the compatibility between Gower and the metric MDS. I will try to use a nMDS instead. 2) for the missing variables, I was thinking of replacing the NA with the group mean in cases where less than 25% of values are missing. I will take a closer look at the variables having a lot of missing values and maybe reject some (>50% of NA). $\endgroup$ – Bettina Dec 9 '13 at 9:24

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.