# changes for Poisson model
set.seed(1)
x1<-c(rpois(125,1),rpois(100,5),rpois(50,1),rpois(75,5),rpois(50,1))
x2<-c(rpois(100,1),rpois(100,3),rpois(75,2),rpois(50,0.5),rpois(75,3))
M<-rbind(x1,x2)
E<-EBSProfiles(M)
C1 <- GetCondition(E,1)
Run the code above in your browser using DataLab