→ obarisk:相同值不是好事嗎？ 08/27 07:49
→ simonwu:但是不同的參數組合 應該是要得到不同的結果 08/27 08:47
→ simonwu:相同的值 表示 每次都用了相同的參數去計算 08/27 08:48
→ clickhere:for(l ...) 一開始就錯了. 見下篇. 08/27 11:53
推 obarisk:恩！抱歉沒有認真看 08/28 11:16
------------------------------------------------------------------------
[軟體程式類別]:
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