- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 79 个
- 通用积分
- 121.5301
- 学术水平
- 10 点
- 热心指数
- 14 点
- 信用等级
- 7 点
- 经验
- 107838 点
- 帖子
- 203
- 精华
- 0
- 在线时间
- 871 小时
- 注册时间
- 2014-2-7
- 最后登录
- 2021-6-15
|
- > f<-function(a=2,low=0,up=2,logic=TRUE){
- + A<-array();i<-n<-0
- + repeat{
- + rnorm(1)->x
- + if(x>low&x<up){
- + A[i+1]<-x;i=i+1}
- + n=n+1
- + if(i==a) break
- + }
- + r=(i/n);t<-(pnorm(up)-pnorm(low));c<-abs(r-t)
- + p<-data.frame(n=n,a=i,r=r,t=t,c=c)
- + list(A[logic],p)
- + }
- > f()
- [[1]]
- [1] 1.7379671 0.7411203
- [[2]]
- n a r t c
- 1 5 2 0.4 0.4772499 0.07724987
复制代码
|
-
总评分: 学术水平 + 1
热心指数 + 1
信用等级 + 1
查看全部评分
|