powered by
dgam(y,mu,delta)
pgam
y<-0:10 mu<-5 delta<-1 out<-dgam(y,mu,delta) # out equals out1<-dgamma(y,shape=1/delta,rate=1/(delta*mu))
Run the code above in your browser using DataLab