####关于在R里定义三角分布函数#####

RT
现在用的三角分布随机抽取函数是这样的：
TriAng_Pre <- function(vecabc)
{
  min <- vecabc[1]
  mode <- vecabc[2]
  max <- vecabc[3]

  u <- runif(1,0,1)
  if (u <= (mode-min)/(max-min) ){
     r=min+sqrt(u*(max-min)*(mode-min))  }else{
   r=max-sqrt((1-u)*(max-min)*(max-mode))
  }
r
}
   a=50
   b=70
   c=100
M <- matrix(c(rep(a,10000),rep(b,10000),rep(c,10000)),10000)
Result <- apply(M,1,TriAng_Pre)


我想把它定义成一个简单的函数，只要给出最小，中间，最大值 就可以得出结果，不用每次都把这些Code运行一遍。我想到的做法是这样的。分两步：
TriAng_Pre <- function(vecabc)
{
  min <- vecabc[1]
  mode <- vecabc[2]
  max <- vecabc[3]

  u <- runif(1,0,1)
  if (u <= (mode-min)/(max-min) ){
     r=min+sqrt(u*(max-min)*(mode-min))  }else{
   r=max-sqrt((1-u)*(max-min)*(max-mode))
  }
r
}


TriAng_Dis <- function(A,B,C)
{
   a=A
   b=B
   c=C
M <- matrix(c(rep(a,80000),rep(b,80000),rep(c,80000)),80000)
result <- apply(M,1,TriAng_Pre)
return(result)
}


得到的TriAng_Dis（A, B, C）就是我的三角分布函数。但是这个是分了两步定义的，有没有办法把TriAng_Pre 部分直接嵌入到TriAng_Dis的定义中，一次性地定义好这个函数呢？


谢谢！

你在TriAng_Dis里面加个循环

1. TriAng_Pre <- function(vecabc, n)
   
2. {
   
3.     min <- vecabc[1]
   
4.     mode <- vecabc[2]
   
5.     max <- vecabc[3]
   
6.    
   
7.     out <- rep(NA, n)
   
8.     for (i in 1:n) {
   
9.         u <- runif(1,0,1)
   
10.         if (u <= (mode-min)/(max-min) ){
    
11.             r=min+sqrt(u*(max-min)*(mode-min))  }else{
    
12.                 r=max-sqrt((1-u)*(max-min)*(max-mode))
    
13.             }
    
14.         out[i] <- r
    
15.     }
    
16.     return(out)
    
17. }

复制代码

qoiqpwqr 发表于 2016-1-24 22:29
你在TriAng_Dis里面加个循环

赞！！！
谢谢！