我是R的新手。制作循环时遇到问题。我有一个错误消息,我想忽略。
我无法从“ try”或“ trycatch”中了解如何操作。我希望如果有错误,则在我的数据框“ Resum_RL_MM
”中将“ p值”标记为1。
for (IGE in C:D) {
row <- 1
print(IGE)
for(HLA in A:B){
print(HLA)
fit1 <- gee (formula = Egea9[[IGE]] ~ Egea9[[HLA]], id=Num_fam, data=Egea9, family = binomial, corstr = "exchangeable")
p_value <- (2 * pnorm(abs(coef(summary(fit1))[,5]), lower.tail = FALSE))[2]
p_value <- unname (p_value)
Resum_RL_MM [row, col] <- p_value
row <- row + 1
}
}