求解线性规划的Theano实例

Question

我有以下问题。假设我们有x_{jk}，它是一个样本k中基因j的表达值，它是s_{ij}细胞类型的平均表达水平，按各自的比例a_{ki} (i =1)加权。N，N是疾病类型)：

​

​

通常，这可以用矩阵形式表示。

​

​

我想要做的是解决这个方程

​

​

怎么才能用提亚诺来完成呢？

Answer 1

你可以在西亚诺做，也可以不在西亚诺做。西亚诺能帮你做的唯一一件事就是欧几里德范数的梯度，它可以为你计算，但也很容易用手写。解决这个问题的算法需要自己来实现。要么你写拉格朗日，然后用梯度上升和投影到约束来解决对偶问题，要么你直接用梯度下降和约束投影来解决原始问题。您需要自己编写这些优化步骤，这也是您在theano中所做的任何其他优化的情况。