我正在写一个剧本,我还没有最初创作。这是核心记录。我有一个数据框,显示给定深度(第一列)中列中的不同元素组成。我想用它来执行PCA,我对我必须选择的标准化方法感到困惑。
您是否有人使用clr()
为prcomp()
准备数据?或者它是否掺假我的解决方案。除了使用clr()
中的prcomp()
属性之外,我还尝试在使用scale
函数之前对数据使用prcomp()
。
data_f_clr<- clr(data_f)
data_pca <- prcomp(data_f, center = TRUE, scale. = TRUE)
https://stat.ethz.ch/R-manual/R-devel/library/stats/html/prcomp.html
scale
用于扩展数据,因此它们具有单位差异。我认为,由于我的数据具有非常不同的规模。问题是,当我使用上面的代码或者当我跳过clr()
(这会产生更多想要的结果)时,我会收到一个不同的解决方案。但我想知道为什么clr()
在这种情况下令人不安?
如果你们中的任何人能够帮助我,那就太棒了。