機能のエラー処理、通知、スキップ

回帰サマリをcsvファイルに出力し、回帰分析を処理できる機能を処理するのに問題があります。だから、コードは次のようになります。 AGE1（連続）、gender1（カテゴリ0/1）、FLUSHOTファイル（カテゴリ0/1）機能のエラー処理、通知、スキップ

、最初の100列：

私は3つの予測の変数を持っています私がテストしたいレスポンス変数（すべてのカテゴリ0/1）です。

目標は、応答変数（1：100）ごとに回帰分析を行い、p値、OR、および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()の複数の方法を試しましたが、私にとってはうまくいかないようでした。

出典

2017-02-17 Yayi Zhao

どのようなエラーメッセージが表示されますか？ rmsパッケージからlrmを使用して、ロジスティック回帰モデルを推定することができます。そしてtexregをcsvに出力します。

出典

2017-02-17 01:00:29 katharina

キャット、返信ありがとうございます。私が得たエラーコードはこうです。約（sp $ y、sp $ x、xout = cutoff）のエラー：補間するために少なくとも2つの非NA値が必要このエラーを回避する方法はありますか？ –

こんにちはYayi、これはあなたのデータが時には完全に欠落していることを意味します。また、try-catchは、削除する必要のあるデータのエラーを見つけることができないようにします。データをループして、試行された各モデル（DVと機能のクロス集計など）が正常であることを確認することをおすすめします。 – katharina

はい、それはまさに私がやったことです、実行する前にそれをきれいにして、それは動作します。どうもありがとうございました。しかし、私はまだエラーを出力する方法があるとホッピングしているし、プロセス全体を止めることなく移動します。うまくいけば私はそれを把握します。 –

機能のエラー処理、通知、スキップ

答えて

関連する問題