bingen.m

计量工具箱

function s = bingen(p0,p1,m)
% PURPOSE: generate binomial probability
% ----------------------------------------
% USAGE:  s = bingen(p0,p1,m)
% where: p0 = probability 1
%        p1 = probability 2
%         m = number of outputs
% ----------------------------------------
% RETURNS: an (m x 1) vector of 0,1 values
%          based on unif_rnd(m) >= p0/(p0+p1)        

% written by:
% James P. LeSage, Dept of Economics
% University of Toledo
% 2801 W. Bancroft St,
% Toledo, OH 43606
% [e-mail]


pr0 = p0./(p0+p1);
u = rand(m,1);
s = (u >= pr0);