randvector.m

The BNL toolbox is a set of Matlab functions for defining and estimating the parameters of a Bayesi

function R=randvector(r,c,alpha,beta)
% R=randvector(r,c,alpha,beta) genereert r random prob vectoren van lengte c
%alpha en beta zijn de parameters van de gamma verdeling (default alpha=1, beta=1/2)
if nargin==1
    error('aantal en lengte van vectoren opgeven');
elseif nargin==2
    alpha=1;
    beta=1/2;
elseif nargin==3
    beta=1/2;
end
 G=gamrnd(alpha, beta, r,c);
 R=G./repmat(sum(G')',1,c);