preso

git-svn-id: https://ws10smt.googlecode.com/svn/trunk@435 ec762483-ff6d-05da-a07a-a48fb63a330f

1 files changed, 12 insertions, 0 deletions

diff --git a/report/pyp_clustering/acl09-short/code/crppred.m b/report/pyp_clustering/acl09-short/code/crppred.m
new file mode 100644
index 00000000..17f22652
--- /dev/null
+++ b/report/pyp_clustering/acl09-short/code/crppred.m
@@ -0,0 +1,12 @@
+function output = crppred(input,b)
+
+uniqin = unique(input);
+prediction = zeros(max(input),1);
+
+p0=1/30114;
+for i = 1:length(uniqin)
+  prediction(uniqin(i)) = b*p0*sum(1./((1:uniqin(i))+b*p0-1));
+end
+
+output = prediction(input);
+