[程式] 關於Gauss code
------------------------------------------------------------------------
TITLE請著名軟體類別
按ctrl+v可預覽 發文前請把灰色的註解刪除
可使用ctrl+y刪除一整行
[軟體程式類別]:
請填入軟體程式類別 例如SAS、SPSS、R、EVIEWS...等
gauss
[程式問題]:
資料處理、迴歸、敘述統計、logistic、probit...等
線性規劃
[軟體熟悉度]:
請把以下不需要的部份刪除
新手(不到1個月)
根據不同熟悉度
板上高手可以給予你不同指導
新手或許需要的還是語法的熟悉
而其他部分則可能往往是盲點上的點破
[問題敘述]:
請詳盡敘述遭遇到的問題，可能的話，分點敘述你要處理的流程
能夠分點敘述處理流程，基本上在撰寫程式或者跑統計時是相當有幫助的
甚至有經驗的人可以藉此告知你步驟是否錯誤
[程式範例]:
雖然張貼程式很可怕，但基本上有些程式還是要張貼才能解決
x1 = ones(n1,1)~dat1[.,4:24]; @ x matrix = constraint coefs @
x2 = ones(n2,1)~dat2[.,4:24];
x3 = ones(n3,1)~dat3[.,4:24];
x4 = ones(n4,1)~dat4[.,4:24];
x5 = ones(n5,1)~dat5[.,4:24];
x = x1|x2|x3|x4|x5;
yhat1 = x1*parm[.,1]; @ right-hand-sides of constraints@
yhat2 = x2*parm[.,2];
yhat3 = x3*parm[.,3];
yhat4 = x4*parm[.,4];
yhat5 = x5*parm[.,5];
yhat = yhat1|yhat2|yhat3|yhat4|yhat5;
n = rows(x); @ no. of observations @
q = zeros(22,22); r = -x'yhat; @ LP objective function @
@q = x'x; r = x'yhat;@ @ QP objective function @
{b,u1,u2,u3,u4,ret} = QProg(parm[.,6],q,r,0,0,x,yhat,0);
@ minimise objective function @
請問如果要改他的限制式方向要怎麼改呢?就是此篇code的限制式是
X'(b*) >= X'(b) <不等式右邊已知>
現在我想改成X x'(b*) <= X '(b)
也就是限制式理面的不等式 符號方向相反....
感謝.......
-----------------------------------------------------------------------------
--
And as I draw near
The scene becomes clear
Like watching my life on a screen
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.59.238.241
討論串 (同標題文章)
完整討論串 (本文為第 2 之 2 篇)：