R-处理功能中的错误，请通知然后跳过

Question

我在处理一个可以输出回归摘要的函数和过程回归分析。因此，代码看起来像这样：

我有三个预测变量：Age1（连续），性别1（分类0/1），flushot（分类0/1）

在文件中，前100列是我要测试的响应变量（所有分类0/1）。

目标是对每个响应变量（1：100）进行回归分析，并且仅输出p值或和ci。

。

所以我拥有的代码看起来像这样：

fun1<-function(x){
  res<-c(paste(as.character(summary(x)$call),collapse = " "),
         summary(x)$coefficients[4,4],
         exp(coef(x))[4],
         exp(confint(x))[4,1:2],"n")
  names(res)<-c("call","p-value","OR","LCI","UCI","")
  return(res)}

res2=NULL
lms=list()
for(i in 1:100)
{
  lms[[i]]=glm(A[,i]~age1+gender1+as.factor(FLUSHOT),family="binomial",data=A)               
  res2<-rbind(res2,fun1(lms[[i]]))
}
write.csv(res2,"A_attempt1.csv",row.names=F)

例如，我们在每个类别中都有足够的样本量，或者边际频率看起来像这样：

table(variable1,FLUESHOT)
   0    1
0  15   3
1  11   19

此代码运行良好，但是如果我们有类似的内容：

 table(variable15,FLUESHOT)
       0    1
    0  15   0
    1  11   19

代码陷入错误，报告和停止。

我尝试了使用try()和tryCatch()的多种方法，但似乎对我不起作用。

Answer 1

您看到什么错误消息？您可以尝试使用RMS软件包中的LRM来估计逻辑回归模型。和TexReg将其输出到CSV。