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")

R中不同处理的相关矩阵

2 个答案: