ブートストラップ法によるoptimism計算について

Question

前提・実現したいこと

モデルとブートストラップによるROC曲線の比較をしたくて(internal validation)、下記の関数入力をしています。
Optimismを調べるためにc統計値の差を調べるべく compare_aucを使っています。

データ名：X0618_2

ロジスティック回帰分析：BD~ Age65 + Bystander + witness +etiology

compare_auc <- function(formula1, formula2, X0618_2, response, alternative="two.sided"){

• olm1 <- glm(formula1, data=X0618_2, family="binomial")
• olm2 <- glm(formula2, data=X0618_2, family="binomial")
• ROC_1 <- pROC::roc(response=response, predictor=fitted(olm1))
• ROC_2 <- pROC::roc(response=response, predictor=fitted(olm2))
• print(pROC::roc.test(ROC_1, ROC_2, method="delong", alternative=alternative))
• print(pROC::roc.test(ROC_1, ROC_2, method = "bootstrap", boot.n = 2000, progress="none", alternative=alternative))}

compare_auc(BD~ Age65 + Bystander + witness +etiology,as.data.X0618_2,X0618_2$BD,"less")

発生している問題・エラーメッセージ　

エラーメッセージ
```eval(predvars, data, env) でエラー: 
   数値の 'envir' 引数の長さが 1 ではありません 

### 該当のソースコード

```R 
ソースコード

http://freeverb3vst.osdn.jp/sodan/tips/r3.html

試したこと

AUCが出ません。
"less"のところを何かに書き換えるべきかわかりません。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。ぶn

Comments

[kei344]

質問タグの「Bootstrap」はCSSフレームワークの Bootstrapなので、適切なタグに置き換えてください。