しょっちゅう分からなくなってはTwitterで強い人に教えてもらうので，ここにメモ． CartesianIndexは.Iで添え字を取り出せる．これを'collect'すると所望のvecになる. ただ，これだとちょっと気持ち悪いと思う人もいるようなので， Tuple(CartesianIndex(2,3…

前回，補助行列が一つだけのNMMFの実装を公開した． 今回は同様に，これらの論文にのっとり，補助行列が２つある場合のNMMFを実装した． 行列X,Y,Zを，WH, AH,WBに同時に分解するというタスクである．コスト関数はKL情報量．（LSの場合の更新式ってどこかに…