作者KnightX (........)
看板Statistics
标题Re: [问题] R loop的问题
时间Sat Dec 6 15:02:08 2008
※ 引述《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
1F：→ clover520:谢谢! 12/07 09:43