此統計量是經由直接比較樣本分配函數(empirical distribution function)與其最接近的具單峰機率之分配函數間的差距,所得到的一個非參數檢定統計量,它的優點之一是其不須要倚賴非參數曲線估計,但是,它跟其他檢定峰數的非參數統計量有同樣的困難,就是不易以它建造一個具有準確水準的檢定。其難處在於如何對此統計量之值刻度 (亦即決定判斷點(critical point) 的問題)。因為虛無假設下之所有單峰分配的集合構成一無限維度的函數空間。
[R] 先安裝diptest package
載入library: library("diptest")
dip(x, full.result = FALSE, debug = FALSE)
Ref:
1. The Dip test of unimodality, by Hartigan (1985)
Ann. Statist. Volume 13, Number 1 (1985), 70-84
2. Package ‘diptest’ (2009)
http://cran.r-project.org/web/packages/diptest/diptest.pdf
3.鄭明燕-檢定機率密度函數峰數的非參數方法
http://nr.stpi.org.tw/ejournal/natscinews/v12n3/97-98.pdf
沒有留言:
張貼留言