--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/rDiff/src/tests/calcl_nbin_pdf.m	Thu Feb 14 23:38:36 2013 -0500
@@ -0,0 +1,11 @@
+function [RET]=calcl_nbin_pdf(VAL,R,P)
+%uses formula used by matlab for compatability
+
+RET=zeros(size(VAL,1),size(VAL,2));
+IDX=VAL==round(VAL);
+%RET(IDX) = exp(gammaln(R+VAL(IDX))-gammaln(VAL(IDX)+1)-gammaln(R)+R*log(P)+VAL(IDX)*log1p(-P));
+
+
+nk = (gammaln(R + VAL(IDX)) - gammaln(VAL(IDX) + 1) - gammaln(R));
+RET(IDX) = nk + (R.*log(P)) + (VAL(IDX).*log(1 - P));
+RET(IDX) = exp(RET(IDX));