nsprcomp（）：為R中的主成分分析強制執行非負性

Question

我的目標是對數據使用nsprcopm() （ 版本0.5.1.2 ）進行非負 PCA並提取PC的值。 我在R版本3.5.0中按如下方式進行操作：

nnpca <- nsprcomp(mydata[,10:16], center = T, scale. = T, nneg = T)

關鍵是我需要強制執行非負性運算 （因為在下一步中，我將使用glm(..., family = Gamma()) ，而glm不會接受Gamma的負值）

問題在於nsprcopm()為PC產生負值。 這是nnpca$x的快照：

任何想法？？！ 這是nsprcomp()的錯誤還是什么？ 我錯過了什么嗎？

Answer 1

我找到了這個問題，如果有人遇到類似的問題，請把它放在這里以備將來參考：

問題來自於論點center = T 如果為True ，則nsprcomp嘗試移動所有PC，使其位於坐標中心附近。 這將不可避免地使某些PC變成負數（即使原始值是正數）。

可以通過center = F克服