search for: gminvscal

...# 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...