#ifndef __BETA_H__
#define __BETA_H__ 1

double logbeta(int m, int n, int tsize);
double beta_loglik(double sum1, double sum2, double dsize, int m, int n, int k);
void beta_mle(double *x, int dsize, int k, int *mest, int *nest);

#endif