Re: [問題] optim function
※ 引述《buzz20374 ()》之銘言：
:
: - 問題: 當你想要問問題時，請使用這個類別
:
: [問題類型]:
:
: 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
:
: [軟體熟悉度]:
: 請把以下不需要的部份刪除
: 新手(沒寫過程式，R 是我的第一次)
: [問題敘述]:
: 請簡略描述你所要做的事情，或是這個程式的目的
: 想請教optim function的問題。?optim裡面的一個例子:
: ## Rosenbrock Banana function
: fr<-function(x){
: x1<-x[1]
: x2<-x[2]
: 100*(x2-x1*x1)^2+(1-x1)^2
: }
: optim(c(-1.2,2),fr)
: 這裡我就想了很久，不太懂，
: 第一個argument是par:initial values for the parameters to be optimized over
你要先把fr看懂，他的x是向量
從裡面可以看得出來為長度二的向量，
x[1], x[2]分別assign到x1, x2
c(-1.2,2)是x1跟x2的初始值
: 但我不太懂 c(-1.2,2)代表要找值的區間嗎?
: 跑出來的結果也看不太懂:
: >optim(c(-1.2,2),fr)
: $par
: [1] 1.000260 1.000506
最小值的解
: $value
: [1] 8.825241e-08
最小值
: $counts
: function gradient
: 195 NA
number of function evaluation (fr的執行次數)
gradient 你沒提供所以沒有值
: $convergence
: [1] 0
收斂情況編碼，詳細的在optim下方
: $message
: NULL
如果非0，通常都有一些相關訊息
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.73.190
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1445821446.A.AAD.html
推
10/26 10:49, , 1F
10/26 10:49, 1F
→
10/26 10:49, , 2F
10/26 10:49, 2F
推
10/26 11:03, , 3F
10/26 11:03, 3F
195次就只是他在optimization過程中，fr被執行了195次，通常代表執行cost~~
好的演算法跟初始值通常有較低的cost
→
10/26 11:04, , 4F
10/26 11:04, 4F
→
10/26 11:05, , 5F
10/26 11:05, 5F
→
10/26 11:05, , 6F
10/26 11:05, 6F
→
10/26 11:06, , 7F
10/26 11:06, 7F
→
10/26 11:08, , 8F
10/26 11:08, 8F
→
10/26 11:08, , 9F
10/26 11:08, 9F
→
10/26 11:09, , 10F
10/26 11:09, 10F
推
10/26 11:11, , 11F
10/26 11:11, 11F
另外，題外話，你說的搜尋範圍在lower跟upper這兩個地方設定
這兩個用的演算法會稍微不同，你可以看method的解釋
推
10/26 11:29, , 12F
10/26 11:29, 12F
→
10/26 11:30, , 13F
10/26 11:30, 13F
→
10/26 11:30, , 14F
10/26 11:30, 14F
我以前用過其他最佳化套件， maxLik, bbmle 這兩個我覺得都不錯用
R內建的，我用沒幾次就放棄了~"~...
E大的這篇好豐富@@
※ 編輯: celestialgod (140.109.73.190), 10/26/2015 11:35:51
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 2 篇)：
R_Language 近期熱門文章
PTT數位生活區 即時熱門文章
