Re: [問題] R loop的問題

看板Statistics作者 (態度)時間17年前 (2008/12/06 19:39), 編輯推噓0(001)
留言1則, 1人參與, 最新討論串3/3 (看更多)
※ 引述《clover520 (clover)》之銘言: : 不好意思想請問各位高手 我的R code 是出了什麼問題> < 謝謝!! : 我要自己創一組資料 L1 T1 各有兩筆 : L1<T1 皆服從exp(1)分配 L1 T1獨立 : for (i in 1 : 2) { : tmp1 <- rexp(1, rate = exp(1)) : tmp2 <- rexp(1, rate = exp(1)) : if (tmp1 < tmp2) { : L1 <- c(L1,tmp1) : T1 <- c(T1,tmp2) : } : i<-i+1 : } : 我這個code寫時的想法是 R會一直重複做tmp1 tmp2 : 直到拿到前兩筆tmp1<tmp2 ==>然後就可以得到我要的兩筆L1 和兩筆T1 : ie (L11,T11) 和 (L12,T12) 其中L12<T12, L11<T11 : 但...我想的跟我寫出來的無法符合>_< : 所以想請問版上的高手們 我應該怎麼改?? 謝謝大家 不知道這是不是你要的,其實你只要設初始值給T,L存的空間。 R code L=T=rep(0,2) #2的地方是你要抓幾個要n改n個 for (i in 1 : 2) { #但這邊i也要改n,之後都一樣 tmp1 <- rexp(1, rate = exp(1)) tmp2 <- rexp(1, rate = exp(1)) if (tmp1 < tmp2) { L[i] <-tmp1 #存的方式也要改。 T[i]<- tmp2 } i<-i+1 } out put > L [1] 0.01565775 0.14340350 > T [1] 0.7275713 0.4400562 L1<T1, L2<T2 如果有錯,請邊小力點。 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 210.209.232.70

12/07 09:44, , 1F
謝謝!
12/07 09:44, 1F
文章代碼(AID): #19EcHhPQ (Statistics)
討論串 (同標題文章)
文章代碼(AID): #19EcHhPQ (Statistics)