[RsR] robust PCA in robComposition

Alice Domalik @dom@||k @end|ng |rom @|u@c@
Thu Sep 7 03:13:36 CEST 2017


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]]




More information about the R-SIG-Robust mailing list