[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
入門(寫過其他程式，只是對語法不熟悉)
[問題敘述]: 矩陣的範圍不能是負值
[程式範例]:
想找出一個scalar function的optimal value所以用了optim這個function~
目標函式是個向量化函式~
也就是說給一個n x k的matrix會吐出n x 1的vector
有丟過一些值進去算過~
都有給出想要的結果~
但不知道為什麼~
一丟進optim裡馬上就掛了~
跳出"錯誤在matrix(rep(mu_x, length(x)), nrow = s) : 矩陣的範圍不能是負值"~
附上程式碼如下：
http://nopaste.info/e864a2b39d.html
本來這是跑得動的~
但是因為後來加了幾個變數(sig_v, sig_x跟sig_epo)必須是正的~
在負數上的函數值要為0(因為函數統一回傳log值，所以在code中是指定為-Inf)~
ind_vs那幾行程式碼之後就出問題了....
謝謝大家
[關鍵字]:
矩陣 範圍不能是負值