数据引自:http://www.ruanyifeng.com/blog/2013/03/cosine_similarity.html
关键步骤:
对应代码:
note:这样只是为了和公式更好对应。事实上可以用a×b(转置)的方法写,更加简洁。
proc iml;
a={1 2 2 1 1 1 0};
b={1 2 2 1 1 2 1};
ab=a#b;
cos=ab[+]/(sqrt(a*a`)*sqrt(b*b`));
aa=a*a`;
bb=b*b`;
print cos aa bb;
quit;