search for: gmscale

...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=0.5,max=3); miuvec...