R语言:批量计算COD与其他指标的相关性和显著性
要计算R语言中df数据框中COD与其他指标的相关性和显著性,可以使用cor.test()函数进行计算,并根据p值的阈值判断显著性。
下面是一个示例代码,可以帮助你批量计算相关性和显著性,并将结果整理成表格:
# 创建一个空的数据框,用于存储结果
result <- data.frame(Indicator = character(), Correlation = numeric(), Significance = character(), stringsAsFactors = FALSE)
# 循环遍历df中的每个指标
for (col in colnames(df)) {
if (col != 'COD') { # 排除COD列本身
# 计算相关性和显著性
cor_test <- cor.test(df$COD, df[[col]])
# 提取相关系数和p值
correlation <- cor_test$estimate
p_value <- cor_test$p_value
# 根据p值的阈值判断显著性
if (p_value < 0.001) {
significance <- '***'
} else if (p_value < 0.01) {
significance <- '**'
} else if (p_value < 0.05) {
significance <- '*'
} else {
significance <- ''
}
# 将结果添加到数据框中
result <- rbind(result, data.frame(Indicator = col, Correlation = correlation, Significance = significance, stringsAsFactors = FALSE))
}
}
# 打印结果
print(result)
该代码将计算df数据框中COD与其他指标的相关性,并根据p值的阈值判断显著性。最后,将结果以表格形式打印出来。
请注意,上述代码中的阈值为0.001、0.01和0.05,你可以根据需要进行修改。
原文地址: https://www.cveoy.top/t/topic/pgk0 著作权归作者所有。请勿转载和采集!