report/pyp_clustering/acl09-short/code/cohnpred.m


1
2
3
4
5
6
7
8
9
10
11
12

function output = cohnpred(input,b)

uniqin = unique(input);
prediction = zeros(max(input),1);

p0=1/30114;
for i = 1:length(uniqin)
  prediction(uniqin(i)) = b*p0*(psi(b*p0+uniqin(i)) - psi(b*p0));
end

output = prediction(input);