R語言中壹組數據服從威布爾分布,怎麽判斷擬合的效果
首先可以利用fitdistr函數求得weibull分布的形狀參數和尺寸參數,假設數據為x:
library(MASS)?#fitdistr需要利用MASS包fitdistr(x,?densfun?=?"weibull",lower=0)
得到形狀參數shape與尺度參數scale
然後利用ks.test進行檢驗:
ks.test(jitter(x),"pweibull",shape,scale)上邊的jitter用來做小擾動,因為如果x中有重復數據的話ks.test會報錯,如果x中沒有重復數據則不需要jitter。shape是得到的形狀參數,scale是得到的尺度參數。
ks.test得到兩個結果,壹個是D,越小越好,壹個是p-value,這個值要大於0.05