R中不同處理的相關矩陣

Question

用下面的代碼我創建了一個相關矩陣。下面的代碼只是爲所有數據創建一個矩陣，而不考慮治療。但是，我的數據中的一列是治療。我想製作兩個不同的矩陣（每個治療類型一個）。我的治療在第6欄「治療」中是一個明確的價值。第10至第44列我想創建一個矩陣。

correlations <- cor(Plants[,c(10:44)], use="pairwise.complete.obs", method="pearson") 
correlations<-as.matrix(correlations) 

Answer 1

如果您只有兩個治療類別，那麼您可以對每個級別進行兩次單獨分析（僅針對特定級別選擇數據）。只需將名稱Treatment1和Treatment2替換爲治療的實際名稱即可。

kor1<-cor(Plants[Plants$Treat=="Treatment1",c(10:44)], use="pairwise.complete.obs", method="pearson") 
kor2<-cor(Plants[Plants$Treat=="Treatment2",c(10:44)], use="pairwise.complete.obs", method="pearson") 

Answer 2

您可以使用by將函數應用於不同的數據子集。

by(Plants[10:44], Plants["Treat"], 
    cor, use = "pairwise.complete.obs", method = "pearson")