请教，此概率如何表示呢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请大家帮忙看看，谢谢了

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：如何

热心的高手们，帮帮忙吧

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,快的不是一点点

我偏题了， 但又忍不住不说 呵呵

确实快多了，谢谢！今天这问题，特别是中括号里那部分该怎么表示呢

那个连乘是说k从1到4，而且k不等于i，是吗？

不好意思才看到，是这样的

昨晚12点后 还是花了30min了解了 什么是连乘， 是不是 i还不能等于5？

这里的连乘只是对 k ，i 是前面求和的下标，i 的值是从1到m ，而 m 的值由那个条件确定

c <- function(i){
  (5-i)/4
}

p <-function(x){
  for(i in 6:m){
    k <-1:4
1-sum(c((c(1)-x)^4/(c(1)*prod(c(1)-c(k[-1]))),
         (c(2)-x)^4/(c(2)*prod(c(2)-c(k[-2]))),
         (c(3)-x)^4/(c(3)*prod(c(3)-c(k[-3]))),
         (c(4)-x)^4/(c(4)*prod(c(4)-c(k[-4]))),
         sum((c(i)-x)^4/(c(i)*prod(c(i)-c(k))))))
}
}
i是不能等于5的 不然 c(5)=0, 1/c(5)= inf

程序有问题 改改ing！！！！！！！！！！！！！！！！！！