Displaying 1 result from an estimated 1 matches for "miuvectreat".
2011 Aug 25
1
Possible Error in generic function rzigp in ZIGP Package
...mscale <- 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);
countveccontrol <- matrix(0,nrow=nmiu,ncol=1);
#### rzigp function gives error
lgh3 <- length(miuvectreat);
for (i in 1:lgh3)
{
countvectreat[i] <- rzigp(1,miuvectreat[i],dispvectpois...