оценке максимального правдоподобия для плотности гамма в R

Question

я просто имитируемых 100 Randoms наблюдения с помощью гамма-плотности с альфа (параметр формы) = 5 и лямбда (параметр скорости) = 5:

x=rgamma(100,shape=5,rate=5)

В настоящее время, Я хочу запланировать оценки максимального правдоподобия альфа и лямбда с функцией, которая возвращает оба параметра и которые используют эти наблюдения.

Любые советы были бы оценены. Спасибо.

Answer 1

Вы можете использовать fitdistr(...) для этого в пакете MASS.

set.seed(1) # for reproducible example 
x <- rgamma(100,shape=5,rate=5) 

library(MASS) 
fitdistr(x, "gamma", start=list(shape=1, rate=1))$estimate 
# shape  rate 
# 6.603328 6.697338 

Обратите внимание, что с небольшим образцом, подобным этому, вы не получите больших оценок.

x <- rgamma(10000,shape=5,rate=5) 
library(MASS) # may be loaded by default 
fitdistr(x, "gamma", start=list(shape=1, rate=1))$estimate 
# shape  rate 
# 4.984220 4.971021 

fitdistr(...) также возвращает стандартную ошибку оценки и лог-правдоподобия.