交替泊松回归拟合稀疏计数数据CP
内容
设一个例题
我们遵循E. C. Chi和T. G. Kolda在《关于张量、稀疏性和非负因子分解》一书中概述的一般程序，arXiv:1112.2414[数学。2011年12月(http://arxiv.org/abs/1112.2414）.
挑选尺寸和等级Sz = [100 80 60];R = 5;生成因子矩阵，每列有几个大的条目;这将是我们销售的基础。=细胞(3,1);为n = 1:长度(sz) A{n} = rand(sz(n)， R);为r = 1: r p = randperm(sz(n));nbig = round((1/R)*sz(n));A{n}(p(1:nbig)，r) = 100 * A{n}(p(1:nbig)，r);结束结束λ=兰德(R, 1);S = ktensor(， A);=正常化(年代,“排序”1);根据提供的解决方案创建稀疏测试问题。Nz = prod(sz) * .05;信息= create_problem (“溶液”年代,“Sparse_Generation”、新西兰);%提取数据和解决方案X = info.Data;M_true = info.Soln;
叫CP-APR
%计算一个解M = cp_apr(X, R，“printitn”10);%为溶液评分= (M, M_true，“贪婪”,真正的)
CP_PQNR(拟牛顿交替泊松回归)预计算稀疏索引集…完成10。20. Ttl Inner Its: 648, KKT viol = 2.32e-02, obj = -8.28190514e+03, nz: 301Ttl Inner Its: 240, KKT viol = 9.81e-05, obj = -8.28188887e+03, nz: 302 ===========================================最终log-likelihood = -8.281889e+03最终最小二乘拟合= 5.784574e-01最终KKT违例= 9.8085093e-05总内部迭代= 20238总执行时间= 5.85秒factor_match_score = 0.9609