title: R qbeta with vectors slow speed
body: <p>I have 3 vectors called <code>vp</code>, <code>va</code> and <code>vb</code>:  </p>
  
  <p><code>vp</code> contains 1.5 million random probabilities  </p>
  
  <p><code>va</code> contains 1.5 million alpha values  </p>
  
  <p><code>vb</code> contains 1.5 million beta values  </p>
  
  <p>I am trying to create an output that is a vector that looks like:  </p>
  
  <pre><code>vanswer &lt;- c(qbeta(vp,va,vb)) 
  </code></pre>
  
  <p>I know that this works, however the speed is tremndously slow with such large amounts of data and I am trying to find a way to speed it up. I have also tried doing the computation in a matrix with <code>cbind()</code> as well as using <code>sapply()</code>, but I cannot find the right way. </p>
  
  <p>Any help would be helpful!</p>