2015-07-03 40 views
2

我目前正在進行配電裝置。我使用了fitdistr函數,但在確定MLE的初始點時遇到了問題。例如,我想用伽瑪分佈擬合我的數據(降雨量13149乘以1矩陣)。如何確定最大似然法的初始點

fit.gamma = fitdistr(rainfall,dgamma,start=list(shape = ?, scale = ?),method="Nelder-Mead") 
+1

也許這可能幫助:http://stackoverflow.com/questions/14266354/how-can-i-estimate-the-shape-and-scale-of-a-gamma-dist-with-a -particular-mean-a –

+0

謝謝@Pascal的建議 –

回答

1

圖書館fitdistrplus是非常好的。如果您沒有起始值,它會爲您猜測gamma參數。此外,如果猜測失敗,您可以使用片刻方法。

x <- rgamma(100, 0.5, 0.5) 

library(fitdistrplus) 
(pars <- fitdist(x, "gamma")) 
# Fitting of the distribution ' gamma ' by maximum likelihood 
# Parameters: 
#  estimate Std. Error 
# shape 0.4443304 0.05131369 
# rate 0.5622472 0.10644511 
+0

謝謝Legalizelt ..我已經找到了幾個包,包括fitdistrplus,但困惑哪一個是最好的。你的意思是使用瞬間的方法來找到起始值? –

+0

@Legalizelt我必須做一些研究的初始值,我不知道是否有一些方法來找到最佳的出發點,可以給出最好的估計結果。 –

+0

好的,非常感謝有用的回覆@Legalizelt –