Re: [問題] R loop的問題
※ 引述《clover520 (clover)》之銘言：
: 不好意思想請問各位高手 我的R code 是出了什麼問題> < 謝謝!!
: 我要自己創一組資料 L1 T1 各有兩筆
: L1<T1 皆服從exp(1)分配 L1 T1獨立
太久沒玩 R 了, 試寫一下比較笨的方法：
n <- 4 # 可換任意偶數, 依你的問題是 4
x <- rexp(n, rate=exp(1))
(T1 <- ifelse((x[1:(n/2)] - x[(n/2+1):n]) > 0, x[1:(n/2)], x[(n/2+1):n]))
(L1 <- ifelse((x[1:(n/2)] - x[(n/2+1):n]) <= 0, x[1:(n/2)], x[(n/2+1):n]))
all(T1 > L1) # 驗證
這樣就會是你要的結果了～！
不過我原本是想把 L1 那段程式換成 setdiff(x, T1), 這樣可以少重複計算一次,
可是這樣結果在驗證的時候偶爾會出現 FLASE 的情形,
就跟原 PO 原程式後面提出的疑問一樣,
所以我這算是比較笨的方法, 不知有啥高手可以幫我再更簡化一下～！
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.169.121.214
→
12/07 09:43, , 1F
12/07 09:43, 1F
討論串 (同標題文章)