Displaying 1 result from an estimated 1 matches for "euv6qi".
Did you mean:
euv6qy
2012 Oct 04
2
Help with R Fitting an inverse Gamma
Dear all,
I am new in R and would like to ask for someone's help in understanding
where I go wrong with the following code:
rm(list=ls())
# Required packages
library(MCMCpack)
# Simulated data
set.seed(1)
data = rinvgamma(n=250, shape = 5, scale = 2) + 2
hist(data)
# log-likelihood
ll = function(par){
if(par[1]>0 & par[2]>0 & par[3]<min(data)) return(