param <- vector('list', 3)
param[[1]] <- list(prop = 1/3, # component 1 : alpha > beta
alpha = .04,
beta = .02,
deltaIn = 100,
deltaOut = 100,
R = 500
)
param[[2]] <- list(prop = 1/3, # component 2 : just permute alpha and beta ;
alpha = .01,
beta = .02,
deltaIn = 100,
deltaOut = .1,
R = 1000
)
param[[3]] <- list(prop = 1/3, # component 3 : alpha=beta
alpha = .015,
beta = .015,
deltaIn = .1,
deltaOut = .1,
R = 1000
)
obs <- sampleDPAMixture(M=20, param)
Run the code above in your browser using DataLab