[RsR] robust PCA in robComposition

Filzmoser Peter peter@|||zmo@er @end|ng |rom tuw|en@@c@@t
Thu Sep 7 06:26:40 CEST 2017


Hard to say without no more details about the data frame. Maybe there are zeros included?
Best regards,
Peter

------------------------------------------------------------------------------
Prof. Dr. Peter Filzmoser
CSTAT - Computational Statistics
Institute of Statistics & Mathematical Methods in Economics
Vienna University of Technology
Wiedner Hauptstr. 8-10
1040 Vienna, Austria
(DVR: 0005886)
Tel. +43 1 58801 10560
FAX +43 1 58801 105699
e-mail: P.Filzmoser using tuwien.ac.at
http://www.statistik.tuwien.ac.at/public/filz/
------------------------------------------------------------------------------

________________________________________
From: R-SIG-Robust <r-sig-robust-bounces using r-project.org> on behalf of Alice Domalik <adomalik using sfu.ca>
Sent: Thursday, September 7, 2017 3:13 AM
To: r-sig-robust using r-project.org
Subject: [RsR] robust PCA in robComposition

I am trying to perform a robust PCA analysis on compositional data. My data set includes 400 diet samples from a species of predatory fish.
For each sample, I have the percent of the diet that is comprised of 7 species of prey item. Each row sums to 100.

sandlance rockfish salmon saury anchovy herring squid
100.00000 0.00000 0.00000 0.00000 0 0.000000 0.00000
100.00000 0.00000 0.00000 0.00000 0 0.000000 0.00000
75.00000 0.00000 0.00000 0.00000 0 25.000000 0.00000

I tried "pcaCoDa(diet.df)" (from the package robComposition), b ut I get the following error:



Error in x.ilr[, i] <- sqrt((i)/(i + 1)) * log(((apply(as.matrix(x[, 1:i]), : incorrect number of subscripts on matrix




What does this error mean and how can I resolve the problem?

        [[alternative HTML version deleted]]

_______________________________________________
R-SIG-Robust using r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-robust




More information about the R-SIG-Robust mailing list