감마, 로그 정규 분포 및와 이블의 로그 가능성을 찾는 방법은 무엇입니까?

Question

R에서 감마,와 이블 및 로그 정규 분포가 주어진 데이터의 로그 우도를 찾고 싶습니다. 각각의 분포의 매개 변수를 이미 추정 했으므로 어떻게 진행합니까?

Answer 1

다음은 감마의 예입니다. Weibull과 log-Normal은 똑같은 절차를 따릅니다.

set.seed(101) 
x <- rgamma(20,shape=3,rate=2.5) 

library(MASS) 
(ff <- fitdistr(x,"gamma")) 
##  shape  rate 
## 4.452775 4.175653 
## (1.358630) (1.348722) 

fitdistr

이 방법 로그-가능성을 접근 있습니다

logLik(ff) 
## 'log Lik.' -13.14535 (df=2) 

을 또는 당신이 손으로 작업을 수행 할 수 있습니다

sum(dgamma(x,shape=coef(ff)["shape"],rate=coef(ff)["rate"],log=TRUE)) 
## [1] -13.14535 

설탕/R-약간의 마법 :

with(as.list(coef(ff)), 
     sum(dgamma(x,shape=shape,rate=rate,log=TRUE))) 

분포

• densfun="weibull" ->dweibull()
• densfun="lognormal" - 두 경우>dlnorm()

파라미터 화/fitdistr 간의 파라미터가 일치하고 대응하는 밀도 함수의 이름.