[RsR] robust PCA in robComposition

My data set includes lots of zeros and 100's. 
I thought that might be the problem, but I could find anything in the documentation about requirements for the data frame. 

Alice 


From: "Filzmoser Peter" <peter.filzmoser at tuwien.ac.at> 
To: "Alice Domalik" <adomalik at sfu.ca>, r-sig-robust at r-project.org 
Sent: Wednesday, September 6, 2017 8:26:40 PM 
Subject: Re: robust PCA in robComposition 

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 at tuwien.ac.at 
http://www.statistik.tuwien.ac.at/public/filz/ 
------------------------------------------------------------------------------

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 at r-project.org mailing list 
https://stat.ethz.ch/mailman/listinfo/r-sig-robust