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)

2
votes
2answers
38 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
97 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 ...
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 ...
10
votes
2answers
474 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 ...
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): ...
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

Is this behavior of DistributionFitTest together with StableDistribution a bug?

I have been using DistributionFitTest to fit the following data: ...
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
142 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 ...
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 ...
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: ...
2
votes
2answers
43 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 ...
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, ...
1
vote
1answer
59 views

Find maximum Log likelihood for repeated random variates

I have the following code ...
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 ...
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 ...
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)}$, ...
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 ...
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
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 ...
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
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 ...
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.
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
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 ...
1
vote
1answer
56 views

eigenvalues distribution

How does one plot eigenvalues distribution of a large matrix? I want to plot the distribution to show that the set of eigenvalues of the matrix of current interest is symmetric with respect to the ...
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 ...
1
vote
1answer
140 views

multivariate lognormal distribution in mathematica

Is the LogMultinormalDistribution function really the multivariate lognormal distribution? Because I get different results with the same parameters in ...
0
votes
0answers
21 views

Simulating from CopulaDistribution

I am currently using CopulaDistribution and have been generating RandomVariates from this distribution. Say my code looks sort of like this (this is much simpler) ...
1
vote
0answers
46 views

Why is Plot3D applied to TransformedDistribution in 2 random variables so slow?

I am trying to plot 3D, the pdf of a TransformedDistribution made up of random variables. So I can see the pdf in 3D, similar to this answer here Visualizing a 2-dimensional PDF But it seems to take ...
2
votes
0answers
89 views

Copula - Correlation or Covariance

Can someone explain to me why the CopulaDistribution function with a "Multinormal" kernal asks for the covariance matrix instead of the correlation matrix. Just by looking at the formula for copula, ...
0
votes
0answers
72 views

Creating a Distribution

Say I want to create a multivariate distribution that perhaps, as an example, looks something like this. f(x1,x2,x3)=f1(x1)*f2(x2)*f3(x3)InverseStandardNormal[F1(x1)]... Obviously the function I am ...
1
vote
1answer
96 views

Mixture coefficients and Parameter estimation for moving Normal distributions

I am searching for a way to do the following steps below. I am stumped and I have stumped a few profs too. Any help would be appreciated. I would like to import data from, say, excel or a .txt file ...