作者simonwu (人生就像乘船)
看板Statistics
标题[程式] 用R同时做多次回圈的相关问题
时间Sat Aug 27 02:32:00 2011
------------------------------------------------------------------------
[软体程式类别]:
R
[程式问题]:
optim 与 回圈
[软体熟悉度]:
中(3个月到1年)
[问题叙述]:
在分析data时碰到了 需要将三个已知参数
放进去optim 里面求另一个未知参数的值 这样的情形
假设三个参数分别有 5个可能值
所以一共需要求出 5 * 5 * 5 一共 125 个解
但是我在跑以下的回圈时
却总是只能跑出 125个 相同的值
希望各位可以帮我看看 我的回圈指令的顺序是否有误 感谢大家 <(_ _)>
[程式范例]:
Toptim <- function (T,a,Tm,T0){0.25-(a*T*(T-T0)*sqrt(Tm-T))}
Tm <- seq(31.70682,33.89119,length=5)
T0 <- seq(6.044388,9.184368,length=5)
a <- seq(1.355866e-04-1e-04,1.355866e-04+1e-04,length=5)
Topt <- rep(NA,125)
for (l in 1:125){
for (i in 1:5){
for (j in 1:5){
for (k in 1:5){
Topt[l] <- optim(25,Toptim,a=a[i],T0=T0[j],Tm=Tm[k])$par
}}}}
-----------------------------------------------------------------------------
--
如果我还有一天的时间,那天我要陪你看日出日落
我还有一天的时间吗?有.
所以,很自然.任何时间我都随传随到!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 124.8.76.127
1F：→ obarisk:相同值不是好事吗？ 08/27 07:49
2F：→ simonwu:但是不同的参数组合 应该是要得到不同的结果 08/27 08:47
3F：→ simonwu:相同的值 表示 每次都用了相同的参数去计算 08/27 08:48
4F：→ clickhere:for(l ...) 一开始就错了. 见下篇. 08/27 11:53
5F：推 obarisk:恩！抱歉没有认真看 08/28 11:16