search for: miupareto

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