Displaying 1 result from an estimated 1 matches for "countveccontrol".
2011 Aug 25
1
Possible Error in generic function rzigp in ZIGP Package
...l <-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],dispvectpoiss[i],0);
countveccontrol[i] <- rzigp(1,miuveccontrol[i],dispvectpoiss[i],0);
}
[[alt...