All Questions
Tagged with distributions random
101 questions
0
votes
0
answers
21
views
Simulating data from a piecewise pdf occasionally returns an error [duplicate]
In the following code, I generate data from a piecewise PDF. Usually, I get a sample size of ns, but sometimes I get an error message:
Error: RandomVariate::unsdst:...
2
votes
1
answer
100
views
3
votes
1
answer
239
views
How to use Mathematica to find the maximum expected value of multiple normal distributions?
Is it possible to calculate the maximum expected value of $ n (n\geq 2) $ independent (identical distributed) random variables that follow a normal distribution? If possible, how to find it? If it ...
1
vote
1
answer
133
views
How to quickly sample from the following distributions
Consider the following distributions:
...
0
votes
0
answers
42
views
Combinations of distributions to be called with `RandomVariate`
I'd like to manipulate distributions that I'm going to retrieve using RandomVariate; e.g., I would like to take
...
1
vote
1
answer
98
views
How can I calculate the moments of a truncated random variable?
I am using the Moment function to calculate the moments of an exponential random variable.
as follows,
Moment[ExponentialDistribution[3], 9]
However, I need to ...
4
votes
1
answer
74
views
Speeding up RandomVariate with custom probability distribution
I have the defined a custom distribution, which is a Gaussian where the exponent can be any real n greater than 0 (set to actually be >0.5 in use), i.e.
...
0
votes
1
answer
94
views
How to generate random value under the given distribution fast?
Consider a distribution
distr[EN_, Enu_, mN_] =
Exp[-(Sqrt[EN^2 + mN^2]/Enu)]*Sqrt[(Sqrt[Enu^2-mN^2]/2)^2 - EN^2] Sqrt[
EN^4 - (Enu/4)^4]
where mN is a free ...
1
vote
0
answers
67
views
How to generate random points with weight distribution and specific criteria? [duplicate]
Consider some function $f(E_{1},E_{2})$, say
$$
\tag 1 f(E_{1},E_{2}) = \exp[-E_{1}^{3}+E_{2}^{2}]
$$
I want to generate random $n$ points $E_{1}$, $E_{2}$ using this function as weights, with an ...
5
votes
2
answers
693
views
How do you generate a random number with a pre-specified probability distribution?
I want to generate a random real number between 0 and 1 but not distributed uniformly but with a probability distribution say, 2 Sin[ pi x ]^2 which integrates to 1 ...
0
votes
2
answers
163
views
How to generate random dots with given distances?
I am interested in the following problem, and am curious if Mathematica could be used to solve it.
I start with a radius function $\rho:D\to \mathbf{R}_{\ge 0}$ on a region $D\subseteq \mathbf{R}^n$. ...
0
votes
1
answer
84
views
Is there a difference between generating and loading data?
EDIT-1
This is a minimal example where an error occurs after 44538 samples generated in d1 upon importing as d2. Note that if <...
1
vote
1
answer
162
views
A very basic question about weighted data
I have a dataset in the form
{{A1,B1,C1,...},{A2,B2,C2,...},...}
and the weights in the form
{w1,w2,w3,...}
Having this ...
2
votes
1
answer
61
views
Are there any issues/best practices/suggestions when working with a DataDistribution instead of a distribution?
Suppose there is a random variable $X$ that can take on values $-1,0,1$, with probabilities $P(X=1)=P(X=-1)=.25$ and $P(X=0)=.5$.
How can I work with such a random variable in mathematica? (Find the ...
3
votes
1
answer
194
views
How to generate random variate in custom domain for a distribution?
I have a distribution defined in a particular domain of the variable but now I want to generate a random variable not in the entire domain but only in a subset of the domain. Here is what I'm trying ...
2
votes
1
answer
103
views
Random Number Generation using LogNormalDistribution
I have a problem that I can't seem to figure out on my own due to limited knowledge of mathematics, and I was hoping I could seek help from the community. I would like to use ...
2
votes
1
answer
222
views
RandomVariate form a Specific Distribution with Restrictions
There are a lot of questions about generating random numbers form both uni-variate and bi-variate distributions. Mathematica has a lot of built in distributions for that purpose. However, I could not ...
1
vote
1
answer
70
views
distribution implement for RandomVariate
I have a scalar function of two random matrices of dimension $n$ which are in the Gaussian Unitary Ensemble
$$f(x,y) \in \mathbb R, \quad x^{\dagger}=x \: \operatorname{and} \: y^{\dagger} = y, \:\...
4
votes
2
answers
190
views
Maximum of all distances between any pairs of vertices of a random triangle
Given a set $S_n$ of $n$ points selected uniformly at random in a $d$-ball, I want to estimate the average length of the longest side of a triangle having as vertices any three distinct points of $S_n$...
-1
votes
1
answer
105
views
Can the WishartMatrixDistribution command be used for generating random density matrices?
I am interested in generating random members of two classes of $n \times n$ positive-definite matrices $A$ and $B$--the former symmetric in nature, the latter, Hermitian.
The standard (Ginibre-matrix-...
7
votes
1
answer
393
views
How to distribute points on a Sphere[] cap from a normal distribution?
I was wondering if one could distribute points on a "cap" of a sphere, following a normal distribution of the points instead of a uniform distribution. This normal could be centered at the ...
4
votes
2
answers
173
views
Optimize RandomVariate for a $1/x^4$ distribution
I'm working with the following code, where I define a probability distribution proportional to $1/x^4$.
...
6
votes
1
answer
584
views
RandomVariate cannot sample from a custom multivariate ProbabilityDistribution
Let us try to define 2D distribution through PDF:
...
21
votes
5
answers
1k
views
Sisyphus Random Walk
A Sisyphus Random Walk evolves as follows: With probability (r) you advance the position of the walker by +1. With probability (1-r) the walker resets at x0 = 0.
To simulate the probability of reset ...
0
votes
2
answers
3k
views
How to generate random numbers for Beta distribution in a range?
In Excel I can just write:
BETA.INV(RAND(),alpha,beta,lowerbound,upperbound)
To generate a beta distributed random between lower bound and upper bound
How can ...
2
votes
1
answer
151
views
Random points from an image using a BernoulliDistribution
I am a Mathematica newbie. I need to sample random points from an image which I have stored in a variable img as:
...
1
vote
1
answer
38
views
Can I sample from a distribution whose parameters depend on the current iteration?
Let's say I want to sample from 100 Beta distributions, and in each successive one, the first parameter increases by 1. So given a current iteration number n, I ...
0
votes
1
answer
183
views
Sampling two integer lists of given sizes from a distribution that have equal sums
(This question emerged from discussions in this post.)
Context and code sample:
I am trying to figure out if there is a way to generate two $a$ and $b,$ comprised of $n_a$ and $n_b$ integers which ...
51
votes
3
answers
6k
views
RandomVariate from 2-dimensional probability distribution
A probability distribution can be created in Mathematica (I am using 8.0.1) with e.g.
...
3
votes
1
answer
253
views
Generate a uniform distribution over an $n$-simplex employing quasi-random low-discrepancy points
This post is intended in a similar direction to an earlier one Can I use Compile to speed up InverseCDF? . I now wish to generate uniform points in an $n$-simplex (specifically $n=8$), making use of ...
2
votes
1
answer
372
views
RandomReal in combination with NormalDistribution [closed]
I am normally not working with Mathematica, but I have to understand a code written in i!
Can some body say what the arguments and parameters of code below say?
...
2
votes
1
answer
86
views
Drawing random numbers from TransformedDistribution
Suppose I have created the following TransformedDistribution:
...
7
votes
5
answers
1k
views
Survival Probability for Random Walks
The Survival Probability for a walker starting at the origin is defined as the probability that the walker stays positive through n steps. Thanks to the Sparre-Andersen Theorem I know this PDF is ...
6
votes
2
answers
903
views
Speed up RandomVariate on a custom probability distribution
I have created a custom probability distribution:
...
1
vote
2
answers
195
views
Sampling data from numerical values of density distribution of nucleons
How could i sample my data which follows distribution if i have only numerical values of density distribution as an example data attached as in picture
I tried but do not know how to do it.
The data:...
0
votes
1
answer
173
views
Negative values with RandomVariate
I am using the following code to generate some random number from a custom distribution:
...
8
votes
1
answer
898
views
How to implement Monte Carlo Importance Sample in 2D or in N dimensions?
So far i have managed to implement the Crude Monte Carlo (CMC) and the Importance Sample Monte Carlo (ISMC) for unidimensional problems.
Consider the function ...
4
votes
3
answers
876
views
Construct Distribution Histogram From Random Variable
Given a Beta Random Variable $X$ with parameters $\alpha, \beta$ and a positive constant $n$, suppose I am interested in the distribution of:
$$Y:=\lfloor nX\rfloor$$
Suppose I want a histogram ...
3
votes
1
answer
267
views
Generating the cumulative distribution function (CDF) of a tuple made by RandomFunction
I use the RandomFunction to generate a continuous random function (trying to model a rough surface):
...
7
votes
2
answers
2k
views
Can this code be changed to run faster?
I've worked for several days trying to get this code to run faster. The imported data is 15MB.
...
5
votes
1
answer
79
views
Random normalised variable generation
I want to generate 8 random variables (in reality to form 4 complex numbers) such that the sum of the 8 variables squared is equal to unity. The aim of generating such numbers is to perform a quantum ...
3
votes
2
answers
664
views
Distribution of Function of Random Sum of Random Variables
Let $ S = \sum_{i=1}^n X_i$ where:
Each $X_i$ is independently 3 or 9 (with equal probability), and
The sample size $n$ is itself an independent random variable where $N \sim \text{NegativeBinomial}(...
3
votes
3
answers
147
views
SkewNormalDistribution with a mean value
Since it got incredibly good answers the last time i posted a question, i try it again! ;-)
I programmed a chemical/physical simulation in mathematica, where i can use residence time distribution as ...
3
votes
1
answer
827
views
Simulate random sample from defined probability density function
Suppose I define a function in mathematica that has all of the properties of a valid probability density function. Would it be possible to have Mathematica simulate a random sample of a given size of ...
3
votes
1
answer
475
views
First passage time distribution in mathematica
I'm trying to calculate first passage time distribution for some processes in Mathematica. First I tried to calculate this for Wiener process with zero drift and $\sigma = 1$ and the absorbing ...
7
votes
3
answers
3k
views
How do I define a two point random variable?
How do I define a random variable $X$ such that $Pr[X=x_1]=p_1$ and $Pr[X=x_2]=p_2$ where the values $x_i$, $i=1,\,2$ are just symbols and the probabilities $p_i$, $i=1,\,2$ are symbolic expressions ...
1
vote
1
answer
246
views
Modeling Dependent Joint Distribution
I'm trying to model the joint pdf of (X,Y) from example 4.5.8 of the text Statistical Inference 2nd Edition, Casella, Berger (pg 199 of the pdf). Essentially,
X ~ Uniform[0,1]
Z ~ Unfiorm[0,1/10]
Y ...
10
votes
4
answers
845
views
Uniform distribution on z = 1 − √( x ^2 + y ^2) , z ≥ 0 [duplicate]
Let $S$ belong to ${\rm I\!R^3}$, $z = 1 − \sqrt{x^2 + y^2}$ , $z ≥ 0$.
I need to write a program that generates a random point on $S$ with uniform distribution. I've tried to use
...
4
votes
1
answer
274
views
Random Variate generating strange results when using ProbabilityDistribution
Bug introduced in 10.1 and fixed in 10.3
I have created a probability distribution that follows a general normal distribution, given by:
$$
\frac{b *e^{-\left(\frac{\sqrt{(x-u )^2}}{a }\right)^{b }}}...
1
vote
1
answer
93
views
How to generate number of points in a m/m/c/c process with Poisson arrival and exponential service time?
I am considering the Earling-B or m/m/c/c process.
Let the request arrival rate follows Poisson Process with average rate, $\lambda$
Let the service time follows Exponential Distribution, i.e., $\...