Displaying 1 result from an estimated 1 matches for "miupareto".
2011 Aug 25
1
Possible Error in generic function rzigp in ZIGP Package
...:
rm(list=ls(all=TRUE))
set.seed(123);
library(MCMCpack)
# You can use rzigp() in the ZIGP package. rzigp(n,mu,phi,omega=0) for generalized poisson.
library(ZIGP) # for generalized poisson
library(VGAM) # for pareto ditribution
# parameters for pareto to generate mius
nmiu <- 500; miupareto <- 3 ;
disppareto <- 7;
gmshape <- 1;
gminvscale <- 0.5;
gmscale <- 1/gminvscale;
#generate means for control from Pareto distribution
miuveccontrol <-rpareto(nmiu,miupareto,disppareto);
# generate disperssions with gamma distribution
dispvectpoiss <- rgamma(nmiu, gm...