I am trying to compute the variance of a Dirichlet-multinomial within a for cycle.
dirmult <- dir_mult_GD(x=t(as.matrix(data_zeros)), k=1) alpha <- t(dirmul
