全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2176 11
2012-11-19
未命名.jpg

请大家帮忙看看,谢谢了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2012-11-19 16:25:27
热心的高手们,帮帮忙吧
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-19 20:41:12
levyangel
翻了你之前的一个老帖https://bbs.pinggu.org/thread-1213353-2-1.html
for循环是可以用向量的方式来做的

m=10;c=5;a=2;b=4
k=rgamma(1,a,b)
z<-c()
for(i in 1:100000)

  {
                w=runif(m)
                S=sum(-log(w)/k)
                y=10-2*k
                y.star=1-(a+m)/(S+b)
                z=(y-y.star)^2
   
    }

测试了下, system.time(for...)很久。。。

m=10;c=5;a=2;b=4
k=rgamma(1,a,b)
z<-c()
zzz <- function(){
  w=runif(m)
  S=sum(-log(w)/k)
  y=10-2*k
  y.star=1-(a+m)/(S+b)
  z=(y-y.star)^2
  return(z)
}
c<-replicate(100000,zzz())#重复动作100000次,只要2.5s,快的不是一点点
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-19 20:47:32
我偏题了, 但又忍不住不说 呵呵
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-19 22:57:24
确实快多了,谢谢!今天这问题,特别是中括号里那部分该怎么表示呢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-20 00:32:51
那个连乘是说k从1到4,而且k不等于i,是吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群