File: random_randbeta.html

package info (click to toggle)
freemat 4.2%2Bdfsg1-4
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 141,800 kB
  • ctags: 14,082
  • sloc: ansic: 126,788; cpp: 62,046; python: 2,080; perl: 1,255; sh: 1,146; yacc: 1,019; lex: 239; makefile: 100
file content (100 lines) | stat: -rw-r--r-- 4,516 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<title>FreeMat: RANDBETA Beta Deviate Random Number Generator</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<script type="text/javascript">
  $(document).ready(initResizable);
</script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  <td style="padding-left: 0.5em;">
   <div id="projectname">FreeMat
   </div>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.1.1 -->
  <div id="navrow1" class="tabs">
    <ul class="tablist">
      <li><a href="index.html"><span>Main&#160;Page</span></a></li>
      <li class="current"><a href="pages.html"><span>Related&#160;Pages</span></a></li>
    </ul>
  </div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
  <div id="nav-tree">
    <div id="nav-tree-contents">
    </div>
  </div>
  <div id="splitbar" style="-moz-user-select:none;" 
       class="ui-resizable-handle">
  </div>
</div>
<script type="text/javascript">
$(document).ready(function(){initNavTree('random_randbeta.html','');});
</script>
<div id="doc-content">
<div class="header">
  <div class="headertitle">
<div class="title">RANDBETA Beta Deviate Random Number Generator </div>  </div>
</div><!--header-->
<div class="contents">
<div class="textblock"><p>Section: <a class="el" href="sec_random.html">Random Number Generation</a> </p>
<h1><a class="anchor" id="Usage"></a>
Usage</h1>
<p>Creates an array of beta random deviates based on the supplied two parameters. The general syntax for <code>randbeta</code> is </p>
<pre class="fragment">   y = randbeta(alpha, beta)
</pre><p> where <code>alpha</code> and <code>beta</code> are the two parameters of the random deviate. There are three forms for calling <code>randbeta</code>. The first uses two vectors <code>alpha</code> and <code>beta</code> of the same size, in which case the output <code>y</code> is the same size as both inputs, and each deviate uses the corresponding values of <code>alpha</code> and <code>beta</code> from the arguments. In the other forms, either <code>alpha</code> or <code>beta</code> are scalars. </p>
<h1><a class="anchor" id="Function"></a>
Internals</h1>
<p>The probability density function (PDF) of a beta random variable is </p>
<p class="formulaDsp">
<img class="formulaDsp" alt="\[ f(x) = x^(a-1) * (1-x)^(b-1) / B(a,b) \]" src="form_138.png"/>
</p>
<p> for <code>x</code> between 0 and 1. The function <code>B(a,b)</code> is defined so that the integral of <code>f(x)</code> is 1. </p>
<h1><a class="anchor" id="Example"></a>
Example</h1>
<p>Here is a plot of the PDF of a beta random variable with <code>a=3</code>, <code>b=7</code>.</p>
<pre class="fragment">--&gt; a = 3; b = 7;
--&gt; x = (0:100)/100; t = x.^(a-1).*(1-x).^(b-1); 
--&gt; t = t/(sum(t)*.01);
--&gt; plot(x,t);
</pre><p>which is plotted as </p>
<div class="image">
<img src="betapdf.png" alt="betapdf.png"/>
</div>
  If we generate a few random deviates with these values, we see they are distributed around the peak of roughly <code>0.25</code>.</p>
<pre class="fragment">--&gt; randbeta(3*ones(1,5),7*ones(1,5))

ans = 
    0.4343    0.4220    0.3992    0.2727    0.2475 
</pre> </div></div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
  <ul>
    <li class="navelem"><a class="el" href="index.html">FreeMat Documentation</a></li><li class="navelem"><a class="el" href="sec_random.html">Random Number Generation</a></li>
    <li class="footer">Generated on Thu Jul 25 2013 17:17:45 for FreeMat by
    <a href="http://www.doxygen.org/index.html">
    <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.1.1 </li>
  </ul>
</div>
</body>
</html>