質問編集履歴
4
数値の修正
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -104,7 +104,7 @@
|
|
|
104
104
|
|
|
105
105
|
FF = stats.norm.cdf(x=XX, loc=mu, scale=seg)
|
|
106
106
|
|
|
107
|
-
S
|
|
107
|
+
S += (2 * i - 1)*(math.exp(F0) + math.exp(1-FF))
|
|
108
108
|
|
|
109
109
|
MAD = (- n - S/n) *(1 + 0.75 / n + 2.25 / n**2)
|
|
110
110
|
|
3
コードの修正
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -106,17 +106,17 @@
|
|
|
106
106
|
|
|
107
107
|
S -= (2 * i - 1)*(math.exp(F0) + math.exp(1-FF))
|
|
108
108
|
|
|
109
|
-
|
|
109
|
+
MAD = (- n - S/n) *(1 + 0.75 / n + 2.25 / n**2)
|
|
110
110
|
|
|
111
|
-
|
|
111
|
+
print(MAD)#統計量
|
|
112
112
|
|
|
113
|
-
|
|
113
|
+
if MAD >= 0.752:
|
|
114
114
|
|
|
115
|
-
|
|
115
|
+
count += 1
|
|
116
116
|
|
|
117
|
-
|
|
117
|
+
else:
|
|
118
118
|
|
|
119
|
-
|
|
119
|
+
pass
|
|
120
120
|
|
|
121
121
|
power = count / sim
|
|
122
122
|
|
2
説明の追加
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,4 +1,10 @@
|
|
|
1
1
|
Rのコードをpythonに書き換えたいと思っています。
|
|
2
|
+
|
|
3
|
+
実装したい内容は、以下で、ziは標準正規分布の累積分布関数をあらわします。
|
|
4
|
+
|
|
5
|
+
アンダーソンダーリング検定を1から実行して、最終的にA*>0.752 がシミュレーション回数のうち何個あるかで、検出力を出力したいと考えています。
|
|
6
|
+
|
|
7
|
+
|
|
2
8
|
|
|
3
9
|
|
|
4
10
|
|
|
@@ -22,11 +28,11 @@
|
|
|
22
28
|
|
|
23
29
|
#rdata<-rlogis(n,0,1)
|
|
24
30
|
|
|
25
|
-
rdata<-rt(n,1)
|
|
31
|
+
rdata<-rt(n,1)#累積分布関数に使用する分布から乱数をとる。
|
|
26
32
|
|
|
27
|
-
data<-sort(rdata)
|
|
33
|
+
data<-sort(rdata)#順序統計量(小さい順に並べかえる)
|
|
28
34
|
|
|
29
|
-
cdf<-pnorm(data, mean(data), sd(data))
|
|
35
|
+
cdf<-pnorm(data, mean(data), sd(data))(b11f778f7ed7ac3a11db5dfa9fb94c49.png)
|
|
30
36
|
|
|
31
37
|
S <- 0
|
|
32
38
|
|
1
数値の修正
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -10,7 +10,7 @@
|
|
|
10
10
|
|
|
11
11
|
n = 50
|
|
12
12
|
|
|
13
|
-
sim<-100
|
|
13
|
+
sim<-100
|
|
14
14
|
|
|
15
15
|
count<-0
|
|
16
16
|
|