r语言 - 数据帧中的输出回路柯尔莫果洛夫斯米尔诺夫测试 (ks.test)

Question

我想将循环中记录的 ks.test 的输出合并到数据帧或文件中，而不是在控制台中打印 1155 个测试的输出... ：-(。

column_equality_stats = function(data, lab_stats1, lab_stats2, min_count=100) {
for(i in 1:length(lab_stats1)) {
lab_testcodes_1 = lab_stats1[i]
lab_testcodes_2 = lab_stats2[i]
equal_columns <- filter(data, lab_testcode==lab_testcodes_1 | lab_testcode==lab_testcodes_2)
col1 <- equal_columns[equal_columns$lab_testcode==lab_testcodes_1, 'lab_result']
col2 <- equal_columns[equal_columns$lab_testcode==lab_testcodes_2, 'lab_result']
if(sum(!is.na(col1))>min_count && sum(!is.na(col2))>min_count){
stats <- ks.test(col1, col2)
print(stats)
}
}
}

我想要一个包含以下列的数据帧：col1 和 col2 的名称(方程值(、p 值和 D 值。

乌托邦数据框

提前非常感谢！！

Answer 1

这是我能做到的。请编辑您的问题并添加可重现的示例。

我需要知道您的数据是如何定义的。我无法运行您的代码！:-)

无论如何，只需为每次循环运行创建一个数据帧并将它们绑定在一起即可。

map_dfr咕噜咕噜包为您做到了。

library(purrr)
.column_equality_stats <- function(i, data, lab_stats1, lab_stats2, min_count = 100){
lab_testcodes_1 <- lab_stats1[i]
lab_testcodes_2 <- lab_stats2[i]
equal_columns <- filter(data, lab_testcode==lab_testcodes_1 | lab_testcode==lab_testcodes_2)
col1 <- equal_columns[equal_columns$lab_testcode==lab_testcodes_1, 'lab_result']
col2 <- equal_columns[equal_columns$lab_testcode==lab_testcodes_2, 'lab_result']
if(sum(!is.na(col1))>min_count && sum(!is.na(col2))>min_count){
stats <- ks.test(col1, col2)
res <- data.frame(col1 = lab_testcodes_1,
col2 = lab_testcodes_2,,
pvalue = stats$p.value,
dvalue = stats$statistics)
} else {res <- data.frame()}
res
}
column_equality_stats <- function(data, lab_stats1, lab_stats2, min_count=100) {
map_dfr(seq_along(lab_stats1), 
.column_equality_stats,
data       = data,
lab_stats1 = lab_stats1,
lab_stats2 = lab_stats2,
min_count  = min_count)
}