I tried to use below R code to model inverse gamma distribution (alpha=1,beta=1). However, the resulting histogram is not alike the one plotted in the wiki. Could anyone provide any hint about this? Thank for you help.
v <- rgamma(1,1,100000)
v <- 1 / v
dv <- as.data.frame(v)
ggplot(data=dv,aes(x=v)) + geom_histogram()