2015-10-27 62 views
1

我遇到以下問題。假設我們有x_{jk}它是樣本k中基因的表達值j。它是表達水平的平均值在整個細胞類型S_ {IJ},通過各自的比例加權A_ {き}I = 1 ... N,N是該疾病類型):解決線性編程的Theano示例

enter image description here

一般來說,這可以表示爲矩陣形式

enter image description here

我想要做的就是要解決這個方程

enter image description here

如何使用它來Theano做?

回答

3

你可以在theano中做到這一點,或不在theano中做到這一點。 theano在這裏可以幫助你的唯一東西是歐幾里得標準的梯度,它可以爲你計算,但也可以通過手工編寫。解決問題的算法需要由您自己來實現。您可以編寫拉格朗日函數,然後通過梯度上升和投影到約束上來解決對偶問題,或者直接通過梯度下降和投影到約束上來解決原始問題。您需要自己編寫這些優化步驟,對於您在theano中進行的任何其他優化,情況也是如此。

+0

你能舉個簡單的例子嗎? – neversaint

+3

對不起,這會花費太多時間。博伊德的凸優化書中解釋了很多。另見他們關於近端操作者的技術報告,這在這裏很有用。 – eickenberg

+1

這是[QP](https://en.wikipedia.org/wiki/Quadratic_programming)問題,不是[LP](https://en.wikipedia.org/wiki/Linear_programming)。我不認爲在theano中實施是明智的。一般來說,除非是你的研究領域,否則自己實現任何經典優化問題並不是一個好主意。選擇一個求解器[在這裏](https://en.wikipedia.org/wiki/Quadratic_programming#Solvers_and_scripting_.28programming.29_languages)! – BGabor