产生随机数出错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

n=30
f<-function(x,sigma=0.84,theta=1.12){
  sqrt(2/pi)*(theta/x)*(x*sigma)^theta*exp(-(x*sigma)^(2*theta)*(1/2))
}
rf<-function(f,n,sigma,theta,N=500){
  df<-seq(1,50,length.out=N)
  p<-f(df,sigma,theta)
  p<-p/sum(p)#正则化
  sample(df,size=n,replace=F,prob=p)
}
x=matrix(rep(1,s*n),nrow=s,ncol=n)
xcensor=matrix(rep(1,s*r),nrow=s,ncol=r)
for (k in 1:s){
    x[k,]=rf(f,n,sigma,theta)
    xcensor[k,]= sort(x[k,])[1:r]
}
为什么产生的随机数每一行都一样

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：随机数