Re: [問題] R的function問題
看板Statistics作者casella123 (煩煩煩煩...............)時間17年前 (2008/11/27 15:24)推噓1(1推 0噓 0→)留言1則, 1人參與討論串2/2 (看更多)
※ 引述《bonjones (原罪)》之銘言：
: 想請教一下
: R的function是否可以隨著迴圈改變而"變動"function的內容
: 因為程式的其中一個部分要解 "一元多次方程式" (我用uniroot解根)
: 要解的方程式會隨著每次執行迴圈而不一樣(但是有一個規則)
: 之前都沒遇過這樣的問題><
: 謝謝~
下面是我之前有寫過的程式(寫的很糟)
應該跟你的情況大同小異
真正參數是b0 b1 a但是函數另外跟非參數的w
我要在不同的w下，求出對大值
希望對你有幫助
y = ( rexp(n, rate = 1)-1 ) * mu^(k/2) + mu
gn_nlm = function(para,w)
{ b0 = para[1]
b1 = para[2]
a = para[3]
muu = x%*%c(b0,b1)
gn_nlm = -0.5*n*log(a) - 0.5*w*sum(log(muu)) -
(0.5/a)*sum((y-muu)^2 / muu^w)
gn_nlm = -gn_nlm
}
for (m in 1:7){
b = optim(c(3,2,1), gn_nlm, w=(m-1)/2., control=list(maxit=2000))
optimb[1,m] = b$par[1]
optimb[2,m] = b$par[2]
optimb[3,m] = b$par[3]
optimb[4,m] = b$convergence
}
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.45.119
推
11/28 00:19, , 1F
11/28 00:19, 1F
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 2 篇)：