Displaying 1 result from an estimated 1 matches for "miuveccontrol".
2011 Aug 25
1
Possible Error in generic function rzigp in ZIGP Package
...P) # 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);
miuvectreat <- foldchange*miuveccontrol;
countvectreat <- matrix(0,nrow=nmiu,ncol=1);
c...