質問編集履歴
2
式の訂正
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -8,31 +8,31 @@
|
|
|
8
8
|
|
|
9
9
|
|
|
10
10
|
|
|
11
|
+
データ名：X0618_2
|
|
12
|
+
|
|
13
|
+
|
|
14
|
+
|
|
15
|
+
ロジスティック回帰分析：BD~ Age65 + Bystander + witness +etiology
|
|
16
|
+
|
|
17
|
+
|
|
18
|
+
|
|
11
19
|
compare_auc <- function(formula1, formula2, X0618_2, response, alternative="two.sided"){
|
|
12
20
|
|
|
13
|
-
+
|
|
21
|
+
+ olm1 <- glm(formula1, data=X0618_2, family="binomial")
|
|
14
22
|
|
|
15
|
-
+
|
|
23
|
+
+ olm2 <- glm(formula2, data=X0618_2, family="binomial")
|
|
16
24
|
|
|
17
|
-
+
|
|
25
|
+
+ ROC_1 <- pROC::roc(response=response, predictor=fitted(olm1))
|
|
18
26
|
|
|
19
|
-
+
|
|
27
|
+
+ ROC_2 <- pROC::roc(response=response, predictor=fitted(olm2))
|
|
20
28
|
|
|
21
|
-
+
|
|
29
|
+
+ print(pROC::roc.test(ROC_1, ROC_2, method="delong", alternative=alternative))
|
|
22
30
|
|
|
23
|
-
+
|
|
31
|
+
+ print(pROC::roc.test(ROC_1, ROC_2, method = "bootstrap", boot.n = 2000, progress="none", alternative=alternative))}
|
|
24
32
|
|
|
25
33
|
|
|
26
34
|
|
|
27
35
|
> compare_auc(BD~ Age65 + Bystander + witness +etiology,as.data.X0618_2,X0618_2$BD,"less")
|
|
28
|
-
|
|
29
|
-
|
|
30
|
-
|
|
31
|
-
|
|
32
|
-
|
|
33
|
-
データ名：X0618_2
|
|
34
|
-
|
|
35
|
-
ロジスティック回帰分析：BD~ Age65 + Bystander + witness +etiology
|
|
36
36
|
|
|
37
37
|
|
|
38
38
|
|
|
@@ -64,7 +64,7 @@
|
|
|
64
64
|
|
|
65
65
|
|
|
66
66
|
|
|
67
|
-
|
|
67
|
+
AUCが出ません。
|
|
68
68
|
|
|
69
69
|
"less"のところを何かに書き換えるべきかわかりません。
|
|
70
70
|
|
1
質問タグの変更をしました。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,8 +1,6 @@
|
|
|
1
1
|
### 前提・実現したいこと
|
|
2
2
|
|
|
3
3
|
|
|
4
|
-
|
|
5
|
-
ここに質問の内容を詳しく書いてください。
|
|
6
4
|
|
|
7
5
|
モデルとブートストラップによるROC曲線の比較をしたくて(internal validation)、下記の関数入力をしています。
|
|
8
6
|
|