#include <stdio.h>
#include <stdlib.h>
#include <time.h>

#define MATHLIB_STANDALONE
#include <Rmath.h>

#include "beta.h"

int main(int argc, char **argv)
{
    int i, j, m, n, nrep, mhat, nhat, dsize, k = 100;
    double *x;
    if (argc != 6) {
	printf("\nUsage: %s m n dsize nrep k\n\n", argv[0]);
	return 1;
    }
    set_seed((unsigned int) time(NULL), 20120210);
    m = atoi(argv[1]);
    n = atoi(argv[2]);
    dsize = atoi(argv[3]);
    nrep = atoi(argv[4]);
    k = atoi(argv[5]);
    if (dsize > 0) {
	x = (double *) malloc(dsize * sizeof(double));
	for (i = 0; i < nrep; i++) {
	    for (j = 0; j < dsize; j++) x[j] = rbeta((double) m, (double) n);
	    beta_mle(x, dsize, 2 * k, &mhat, &nhat);
	    printf("%d %d\n", mhat, nhat);
	}
    }
    free(x);
    return 0;
}