作者rockken (态度)
看板Statistics
标题Re: [问题] R loop的问题
时间Sat Dec 6 19:39:22 2008
※ 引述《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
1F：→ clover520:谢谢! 12/07 09:44