For questions about distributions built-in in Mathematica, and functions that operate on them. Also includes questions about defining your own distributions.

learn more… | top users | synonyms (1)

5
votes
2answers
513 views

How to fit one distribution to another?

I have a custom distribution created to model some experimental observations. While too complicated to include in this question, I can provide an example and some illustrations to convey a sense of ...
2
votes
2answers
51 views

FindDistributionParameters gives an error for a mixture of user defined sinh-arcsinh distributions

I have some data (1593 values). I first binned it using the Knuth rule and fitted a function that is a linear combination of two sinh-arcsinh functions. I then ...
0
votes
0answers
100 views

Heuristic Method to create an uniform distribution of points

Introduction I am now trying to program an algorithm to achieve an uniform distribution of points within a space of dimension n. I started working in 2D. The idea is to achieve an uniform ...
1
vote
0answers
33 views

FisherHypergeometricDistribution and WalleniusHypergeometricDistribution, how to specify boundary conditions

I'm trying to write a function that consists of a noncentral hypergeometric distribution. When I'm passing a set of parameters to the function, the following boundary problems occur: ...
2
votes
0answers
25 views

Marginal distribution with weight

Given a multi-variate distribution $P(\mathbf{x},\mathbf{y})$, I would like to obtain the marginal distribution $$P_m(\mathbf{x})=\sum\int_{\mathbf{y}}P(\mathbf{x},\mathbf{y})W(\mathbf{y})$$ In my ...
1
vote
0answers
68 views

Truncated Distribution Parameter Search Produces Singular Jacobian Using FindRoot [closed]

I'm trying to find the 2 parameters of a truncated Weibull distribution using version 10 and FindRoot. I would like to truncate the distribution between t =0 and t = 3600 with the probability of ...
2
votes
2answers
44 views

Cannot use the MultinormalDistribution code

I am trying to use the MultinormalDistribution code with the Quantile function to get the 0.01 quartile of the multinormal ...
10
votes
2answers
475 views

Generating uniform random points over a binary image

So I recently read this excellent post about generating random points over arbitrary regions. I hope Wolfram eventually builds that functionality into Mathematica's Random functions. I want to do this ...
22
votes
7answers
846 views

How to generate random points in a region?

The Mathematica 10 documentation was updated for FindInstance adding support for regions. In my use case, I'm trying to sample points in a set of disks: ...
1
vote
1answer
64 views

Mean of inverse probability distribution does not agree with explicit integral

I tried to work with the inverse exponential distribution and encountered the following, which I think is an inconsistency in Mathematica (version 10.0.2): ...
0
votes
1answer
214 views

Fitting a theoretical density function with log-return sample

I am working on a exercise to fit a theoretical density function,presented below, with log-return series of a stock. $\frac{1}{2}a e^{a(x-\mu)} \quad \text{if x } < \mu $ or $ ...
2
votes
0answers
33 views

Custom graph distribution

I am working with random graphs a lot. It's quite easy to describe random graph using, adjacency matrix or just a list of edges. However I am thinking that I might benefit from describing random ...
-2
votes
1answer
55 views

How to generate normally distributed random number [closed]

I want to sample or return random number from a distribution which is normally distributed with mean of zero and standard deviation of 1. I am using RandomReal[NormalDistribution[0, 1]] I s this ...
2
votes
1answer
38 views
4
votes
2answers
113 views

Iterating distributions

I want to calculate the distribution that comes from summing up random numbers between -0.5 and 0.5. A computationally intensive version of doing this is to simply do the following: ...
6
votes
2answers
172 views

Nearest neighbor distribution

I need to construct nearest neighbor distribution for a given finite set of real numbers. I have no clue how this can be realized. Any suggestion would be very much appreciated.
1
vote
1answer
64 views

calculate variance of conditional sum of Normals

Let $\quad W \sim N(\mu_0, \sigma_0^2), \quad X_1 \sim N(\mu_1, \sigma_1^2), \quad X_2 \sim N(\mu_2, \sigma_2^2)$ denote independent random variables, and let $c$ denote a constant. I want to find ...
4
votes
2answers
143 views

Deriving a Rayleigh distribution

When I tried to find this Rayleigh distribution Mathematica couldn't come with an answer: ...
1
vote
1answer
93 views

Accumulating elements of lists taking the criteria out of another list

In order to find the critical values for the significance of a discrete distribution, for example BinomialDistribution[20,0.6], i need a function which works on the ...
9
votes
2answers
117 views

Why are Binormal and Multinormal giving me different results when I use NExpectation?

fixed in 10.1 (windows) Binormal and Multinormal can both be used to represent normal random variables: ...
0
votes
2answers
71 views

Plot Bivariate Data with Marginals

I would like to plot bivariate data along with the marginal distributions. I've arranged the plots in a GraphicsGrid as below but the marginals are centered away from the bivariate plot and I wonder ...
1
vote
1answer
45 views

Generate a sample from a multivariate Cauchy distribution

I want to generate a random sample from a multivariate Cauchy distribution, however I couldn't find a function for the multivariate Cauchy in Mathematica. I know how to define a distribution in 1D and ...
4
votes
1answer
58 views

How to get probabilities for multinomial & hypergeometric distribution ranges more quickly?

I often need to calculate the probability of a result from a multinomial distribution falling within some range, e.g., the maximum of the categories having an upper bound, or a lower bound for the ...
1
vote
0answers
84 views

Extrapolation of area for a 2D integration

I have a set of points distributed almost uniformly in a certain area of a 2D plane as follow: An example of data points (orange) data = {{0.919443, 1.68921*10^-22}, {0.262277, ...
0
votes
2answers
121 views

SplicedDistribution - different definition across Mathematica 9?

I am having troubles with my notebook. I succesfully defined Asymmetric Laplace Distribution at school using SplicedDistribution function, and Mathematica was able ...
0
votes
0answers
32 views

How to fix the PlotRange on LogPlot? [duplicate]

I'm plotting families of Boltzmann'esque distribution functions $\propto\exp\left(-{\dfrac{\varepsilon+{\mathcal O}(\varepsilon^2)}{k_\mathrm{B} T}}\right)$. But it won't let me fix the ...
1
vote
1answer
67 views

Cannot evaluate TransformedDistribution for vectorial sum

I'm trying to obtain the distribution of the magnitude (length) of the vectorial sum of two vectors under a random (uniformly distributed) angle $\phi$. The sum is $\sqrt{x^2+y^2+2xy\cos(\phi)}$, ...
2
votes
1answer
56 views

Why is RandomVariate of special probability distribution doesn't working in all cases

I would like to generate a random variates in Mathematica 9 with this distribution: ...
0
votes
1answer
33 views

Fit process parameters to a transformed AR(1) process

I would like to fit the parameters of an exponented AR(1) process using Mathematica's EstimatedProcess, however, the function does not seem to evaluate to anything. First of all I create the ...
4
votes
2answers
197 views

Calculate PDF and CDF of a product of independent exponentially distributed random variables

Let x[1], x[2], ..., x[n] be independent random variables each of which is distributed exponentially, i.e. with the PDF ...
0
votes
0answers
76 views

Von Neumann's method to generate points uniformly distributed on a region of a n-sphere surface

I am trying to generate points uniformly distributed on a region of a single n-sphere surface (all angles in hyperspheric coordinates are between 0 and Pi/2). I decided to use von Neumann's method ...
2
votes
1answer
81 views

Time inhomogeneous Markov Chain in Mathematica

I would like to create a discrete 2-state Markov process, where the switching probabilities in the transition matrix vary with time. I can currently do the following, which creates a process with ...
1
vote
0answers
31 views

Fit a custom process to a data: inhomogeneous 2-state Markov chain

I would like to fit a custom process - a time inhomogeneous 2-state Markov chain, to data. The time inhomogeneity is a result of the transition probabilities varying sinusoidally through time with a ...
1
vote
1answer
55 views

Using Inner to calculate a log likelihood

I am trying to use the function Inner to calculate the sum of log likelihoods of a custom probability distribution 'custom' applied to a vector of data 'vData'. The custom distribution changes ...
0
votes
1answer
70 views

Sketch Probability density function [closed]

My question is : Suppose that the PDF of a random variable X is as follows: \begin{equation*} f(x) = \begin{cases} ce^{-2x} & \text {$x > 0$}\\ 0 & \text {otherwise}\\ ...
0
votes
1answer
75 views

Custom distribution [closed]

I am trying to define a custom distribution using the code below. Unfortunately, it did not produce a satisfactory result ...
1
vote
1answer
92 views

Mean Preserving PDF Spreading

I have a histogram representing the PDF (probability density function) of an unknown discrete RV. The histogram is asymmetrical. Is there a known way in Mathematica to increase/decrease the variance ...
4
votes
2answers
177 views

When analytical and numerical methods do not agree - Case study with Maximum Likelihoods methods

Here is the probability distribution I am interested in: $$P(q)=C e^{4 n s q} q^{4 n \nu - 1} (1 - q)^{4 n \mu - 1}$$ , where $e$ is the constant of Euler and $C$ is constant so that the whole thing ...
4
votes
2answers
290 views

How to plot an implicit value funtion

Previous question has been solved (Gap in a continuous plot). Here is a new related question. ...
2
votes
2answers
65 views

RandomVariate from Dirichlet into Multinomial Not Summing to 1

I have a draw from a Dirichlet distribution that I would like to use as input into a MultinomialDistribution. For instance: ...
-5
votes
1answer
71 views

Plot 2 student T distributions with different center values

Plot 2 student T distributions with different center values µ, but same DOF. The center values,µ has to be 45 and 48, DOF shall be 9.
2
votes
2answers
107 views

Interpretation of CopulaDistribution [closed]

I was just wondering, when I call the CopulaDistribution function in Mathematica, am I calling its cumulative function or its density function? I have looked up ...
1
vote
2answers
71 views

Random Variate from a transformed Distribution

I would like to simulate random variates from a transformed distribution of a joint distribution and a constant. i.e. ...
1
vote
1answer
65 views

Joint Distribution with a Constant

I would like to simulate from a joint distribution multiplied by a constant. e.g. 3*NormalDistribution[0,1]*NormalDistribution[1,1] I tried using something like ProductDistribution but that doesn't ...
0
votes
1answer
41 views

Product Distribution with Copula

I am trying to create a joint probability function that is of the form: f(a,b,c)=f1(a).f2(b).f3(c).copula where ...
9
votes
2answers
214 views

Sampling phase space uniformly for a given energy

For example, for the 1D harmonic oscillator, the energy is given by: $$ E = \frac{1}{2 m} p^2 + \frac{m \omega^2}{2} x^2 $$ where x is position and p is the momentum. I would like to sample phase ...