a=c(5,2,3)K=length(a)rSBM(K=K,p_small=0.5,p_large=0.0,eta=1.0e3,gam=1.5,del=1.5)gamdel=shape_Dir2genDir(a)# tapered parameters to mimic Dirichlet special caserSBM(K=K,p_small=0.5,p_large=0.0,eta=1.0e3,gam=gamdel[,1],del=gamdel[,2])# tapered to mimic Dirichlet special case