簡體   English   中英

使用R的Kolmogorov-Smirnov

[英]Kolmogorov-Smirnov using R

長話短說,我想手動編寫Kolmogorov-Smirnov單樣本統計信息的代碼,而不是在R中使用ks.test()。據我了解,KS測試可以分解為分子與分子之間的比率。分母。 我對寫分子很感興趣,據我所知,這是觀察樣本與理論假設之間的最大絕對差。 讓我們以以下情況為例:

         Data    Expected
1  0.01052632 0.008864266
2  0.02105263 0.010969529
13 0.05263158 0.018282548
20 0.06315789 0.031689751
22 0.09473684 0.046315789
24 0.26315789 0.210526316
26 0.27368421 0.220387812
27 0.29473684 0.236232687
28 0.30526316 0.252520776
3  0.42105263 0.365650970
4  0.42105263 0.372299169
5  0.45263158 0.398781163
6  0.49473684 0.452853186
7  0.50526316 0.460277008
8  0.73684211 0.656842105
9  0.74736842 0.665484765
10 0.75789474 0.691523546
11 0.77894737 0.718005540
12 0.80000000 0.735955679
14 0.84210526 0.791135734
15 0.86315789 0.809972299
16 0.88421053 0.838559557
17 0.89473684 0.857950139
18 0.96842105 0.958337950
19 0.97894737 0.968642659
21 0.97894737 0.979058172
23 0.98947368 0.989473684
25 1.00000000 1.000000000

在這里,我想獲得最大絕對差(數據-期望值)。

有人有主意嗎? 如有必要,我可以改寫這個問題。 謝謝!

我一直在尋找與此代碼類似的答案:

> A <- with(df, max(abs(Data-Expected)))

,其中df是數據幀。

在這里,我獲得每個數據與期望值之間的差異,將值轉換為絕對值,然后從絕對差的向量中選擇最大值。 因此,答案是:

> A
0.082

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM