作者casella123 (烦烦烦烦...............)
看板Statistics
标题Re: [问题] R的function问题
时间Thu Nov 27 15:24:28 2008
※ 引述《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
1F：推 bonjones:感谢~ 有一点idea了^^ 11/28 00:19