[程式] 用R同時做多次迴圈的相關問題
------------------------------------------------------------------------
[軟體程式類別]:
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
→
08/27 07:49, , 1F
08/27 07:49, 1F
→
08/27 08:47, , 2F
08/27 08:47, 2F
→
08/27 08:48, , 3F
08/27 08:48, 3F
→
08/27 11:53, , 4F
08/27 11:53, 4F
推
08/28 11:16, , 5F
08/28 11:16, 5F
討論串 (同標題文章)
完整討論串 (本文為第 1 之 2 篇)：