|
SSJ
3.3.1
Stochastic Simulation in Java
|
This class represents a gamma process [171] (page 82) \(\{ S(t) = G(t; \mu, \nu) : t \geq0 \}\) with mean parameter \(\mu\) and variance parameter \(\nu\). More...
Public Member Functions | |
| GammaProcess (double s0, double mu, double nu, RandomStream stream) | |
Constructs a new GammaProcess with parameters \(\mu= \mathtt{mu}\), \(\nu= \mathtt{nu}\) and initial value \(S(t_0) = \mathtt{s0}\). More... | |
| GammaProcess (double s0, double mu, double nu, GammaGen Ggen) | |
Constructs a new GammaProcess with parameters \(\mu= \mathtt{mu}\), \(\nu= \mathtt{nu}\) and initial value \(S(t_0) = \mathtt{s0}\). More... | |
| double | nextObservation () |
| double | nextObservation (double nextT) |
Generates and returns the next observation at time \(t_{j+1} = \mathtt{nextTime}\), using the previous observation time \(t_j\) defined earlier (either by this method or by setObservationTimes), as well as the value of the previous observation \(X(t_j)\). More... | |
| double [] | generatePath () |
| Generates, returns and saves the path \(\{X(t_0), X(t_1), …, X(t_d)\}\). More... | |
| double [] | generatePath (double[] uniform01) |
| Generates, returns and saves the path \( \{X(t_0), X(t_1), …, X(t_d)\}\). More... | |
| void | setParams (double s0, double mu, double nu) |
| Sets the parameters \(S(t_0) = \mathtt{s0}\), \(\mu= \mathtt{mu}\) and \(\nu= \mathtt{nu}\) of the process. More... | |
| double | getMu () |
| Returns the value of the parameter \(\mu\). | |
| double | getNu () |
| Returns the value of the parameter \(\nu\). | |
| void | setStream (RandomStream stream) |
Resets the umontreal.ssj.rng.RandomStream of the umontreal.ssj.randvar.GammaGen to stream. | |
| RandomStream | getStream () |
Returns the umontreal.ssj.rng.RandomStream stream. | |
Public Member Functions inherited from StochasticProcess | |
| void | setObservationTimes (double[] T, int d) |
Sets the observation times of the process to a copy of T, with. More... | |
| void | setObservationTimes (double delta, int d) |
| Sets equidistant observation times at \(t_j = j\delta\), for. More... | |
| double [] | getObservationTimes () |
| Returns a reference to the array that contains the observation times. More... | |
| int | getNumObservationTimes () |
| Returns the number \(d\) of observation times, excluding the time \(t_0\). | |
| abstract double [] | generatePath () |
| Generates, returns, and saves the sample path \(\{X(t_0), X(t_1), \dots, X(t_d)\}\). More... | |
| double [] | generatePath (RandomStream stream) |
Same as generatePath(), but first resets the stream to stream. | |
| double [] | getPath () |
| Returns a reference to the last generated sample path \(\{X(t_0), ... , X(t_d)\}\). More... | |
| void | getSubpath (double[] subpath, int[] pathIndices) |
Returns in subpath the values of the process at a subset of the observation times, specified as the times \(t_j\) whose indices. More... | |
| double | getObservation (int j) |
| Returns \(X(t_j)\) from the current sample path. More... | |
| void | resetStartProcess () |
| Resets the observation counter to its initial value \(j=0\), so that the current observation \(X(t_j)\) becomes \(X(t_0)\). More... | |
| boolean | hasNextObservation () |
Returns true if \(j<d\), where \(j\) is the number of observations of the current sample path generated since the last call to resetStartProcess. More... | |
| double | nextObservation () |
| Generates and returns the next observation \(X(t_j)\) of the stochastic process. More... | |
| int | getCurrentObservationIndex () |
| Returns the value of the index \(j\) corresponding to the time. More... | |
| double | getCurrentObservation () |
| Returns the value of the last generated observation \(X(t_j)\). | |
| double | getX0 () |
| Returns the initial value \(X(t_0)\) for this process. | |
| void | setX0 (double s0) |
Sets the initial value \(X(t_0)\) for this process to s0, and reinitializes. | |
| abstract void | setStream (RandomStream stream) |
Resets the random stream of the underlying generator to stream. | |
| abstract RandomStream | getStream () |
| Returns the random stream of the underlying generator. | |
| int [] | getArrayMappingCounterToIndex () |
| Returns a reference to an array that maps an integer \(k\) to \(i_k\), the index of the observation \(S(t_{i_k})\) corresponding to the \(k\)-th observation to be generated for a sample path of this process. More... | |
Protected Member Functions | |
| void | setLarger (double[] path, int left, int mid, int right) |
| double | setLarger (double[] path, int left, int right) |
| double | setLarger (double v) |
| void | init () |
Protected Member Functions inherited from StochasticProcess | |
| void | init () |
Protected Attributes | |
| boolean | usesAnti = false |
| RandomStream | stream |
| GammaGen | Ggen |
| double | mu |
| double | muOverNu |
| double [] | mu2dtOverNu |
Protected Attributes inherited from StochasticProcess | |
| boolean | observationTimesSet = false |
| double | x0 = 0.0 |
| int | d = -1 |
| int | observationIndex = 0 |
| int | observationCounter = 0 |
| double [] | t |
| double [] | path |
| int [] | observationIndexFromCounter |
Static Protected Attributes | |
| static final double | EPS = 1.0e-15 |
Package Attributes | |
| double | nu |
| double | mu2OverNu |
This class represents a gamma process [171] (page 82) \(\{ S(t) = G(t; \mu, \nu) : t \geq0 \}\) with mean parameter \(\mu\) and variance parameter \(\nu\).
It is a continuous-time process with stationary, independent gamma increments such that for any \(\Delta t > 0\),
\[ S(t + \Delta t) = S(t) + X,\tag{GammaEqn} \]
where \(X\) is a random variate from the gamma distribution Gamma \((\mu^2\Delta t / \nu, \mu/ \nu)\).
In this class, the gamma process is sampled sequentially using equation ( GammaEqn ).
| GammaProcess | ( | double | s0, |
| double | mu, | ||
| double | nu, | ||
| RandomStream | stream | ||
| ) |
Constructs a new GammaProcess with parameters \(\mu= \mathtt{mu}\), \(\nu= \mathtt{nu}\) and initial value \(S(t_0) = \mathtt{s0}\).
The gamma variates \(X\) in ( GammaEqn ) are generated by inversion using stream.
| GammaProcess | ( | double | s0, |
| double | mu, | ||
| double | nu, | ||
| GammaGen | Ggen | ||
| ) |
Constructs a new GammaProcess with parameters \(\mu= \mathtt{mu}\), \(\nu= \mathtt{nu}\) and initial value \(S(t_0) = \mathtt{s0}\).
The gamma variates \(X\) in ( GammaEqn ) are supplied by the gamma random variate generator Ggen. Note that the parameters of the umontreal.ssj.randvar.GammaGen object Ggen are not important since the implementation forces the generator to use the correct parameters (as defined above).
| double [] generatePath | ( | ) |
Generates, returns and saves the path \(\{X(t_0), X(t_1), …, X(t_d)\}\).
The gamma variates \(X\) in ( GammaEqn ) are generated using the umontreal.ssj.rng.RandomStream stream or the umontreal.ssj.rng.RandomStream included in the umontreal.ssj.randvar.GammaGen Ggen.
| double [] generatePath | ( | double [] | uniform01 | ) |
Generates, returns and saves the path \( \{X(t_0), X(t_1), …, X(t_d)\}\).
This method does not use the umontreal.ssj.rng.RandomStream stream nor the umontreal.ssj.randvar.GammaGen Ggen. It uses the vector of uniform random numbers \(U(0, 1)\) provided by the user and generates the path by inversion. The vector uniform01 must be of dimension \(d\).
| double nextObservation | ( | double | nextT | ) |
Generates and returns the next observation at time \(t_{j+1} = \mathtt{nextTime}\), using the previous observation time \(t_j\) defined earlier (either by this method or by setObservationTimes), as well as the value of the previous observation \(X(t_j)\).
Warning: This method will reset the observations time \(t_{j+1}\) for this process to nextT. The user must make sure that the \(t_{j+1}\) supplied is \(\geq t_j\).
| void setParams | ( | double | s0, |
| double | mu, | ||
| double | nu | ||
| ) |
Sets the parameters \(S(t_0) = \mathtt{s0}\), \(\mu= \mathtt{mu}\) and \(\nu= \mathtt{nu}\) of the process.
Warning: This method will recompute some quantities stored internally, which may be slow if called repeatedly.
1.8.14