Displaying 1 result from an estimated 1 matches for "gminvscale".
2011 Aug 25
1
Possible Error in generic function rzigp in ZIGP Package
...# 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, gmshape, gminvscale, gmscale)
foldchange <- runif(nmiu,min=...