所以我为我的论文创建了一个数据库,我想比较不同情况下的物质来获得P值。为了做到这一点,我在R中创建了一个不自动的脚本。我对如何使用for
循环创建一个脚本来执行此操作非常感兴趣。我想为每一行选择组,但也要为每三列选择组。
如果可能,请告诉我。提前致谢!
现在我这样做。首先选择第1行前三个值,然后在第一行选择接下来的三个值,以便在下一步中比较它们。
row1_1 <- my_data[1, 2:4]
row1_2 <- my_data[1, 6:8]
根据问题,以下将使用for循环遍历每一行
for (i in 1:nrow(my_data))
{
case1<- my_data[i, 2:4]
case2<- my_data[i, 6:8]
compare_cases <- t.test(case1, case2, alternative = "two.sided", var.equal = FALSE)
print(case1)
print(case2)
print(compare_cases)
}
根据注释修改了错误处理的变体
pvalues <- c()
for (i in 1:nrow(my_data))
{
case1<- my_data[i, 2:4]
case2<- my_data[i, 6:8]
pval <- NA
test <- tryCatch({ t.test(case1, case2, alternative = "two.sided", var.equal = FALSE) },
warning = function(war) { NA },
error = function(err) { NA },
finally = { })
print(paste0("Test for row i=",i))
print(test)
if (identical(test,NA)) pval <- NA
else pval <- test$p.value
pvalues <- c(pvalues,pval)
}
print("All my pvalues for all rows are:")
print(pvalues)