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

0
votes
0answers
37 views

Tricky Poisson distribution question [migrated]

Am new to probability and have been struggling to solve the following question. Assume X1,...,X4 are IID with Xi ∼ Po(λ). Let Y = 1/4 (X1 + ··· + X4). Find Pr (Y < 2) in terms of lambda. Can ...
0
votes
2answers
117 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
28 views

Need help for Introducing model missmatch [closed]

I intend do introduce the model mismatch issue from the basic level. The above figure is a very nice example and I found it from a mathematica code which is complicated to explain with simple ...
0
votes
0answers
3 views

Expected value in probability? [migrated]

Isn't the definition of expected value little non mathematical. I mean average is defined when experiment is done, so how can we before hand tell the expected value. Let say for toss, if I get ₹1 for ...
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
64 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
44 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
23 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
167 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
56 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
64 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
19 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
47 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
53 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
70 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
84 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
165 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
261 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
61 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
58 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.
1
vote
2answers
62 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
56 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 ...
9
votes
2answers
200 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
51 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 ...
1
vote
1answer
111 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 ...
1
vote
0answers
42 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 ...
0
votes
0answers
70 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
62 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 ...
2
votes
1answer
104 views

How to get a more compact form of this probability calculation?

Inspired by the probability calculation here, I am trying to solve a little general one: $$\mathbb{P} (\sum_{i=1}^{m-1} A_i + \sum_{i=1}^{m} S_i < L < \sum_{i=1}^{m} A_i + \sum_{i=1}^{m+1} ...
4
votes
2answers
193 views

what is The Formula behind the Kolmogorov-Smirnov test statistics?

I calculate the Kolmogrov-samirnov test statistic for the following data: ...
5
votes
1answer
292 views

How to solve this probability symbolically or numerically?

I am trying to calculate the following probability $$\mathbb{P} \big(\sum_{i=1}^{m} (A_i + S_i) \le L < \sum_{i=1}^{m+1} (A_i + S_i) \big)$$ where, $$A_i \sim \exp(\lambda), \quad S_i \sim ...
1
vote
1answer
121 views

Country density plot of lines (highways)

I love solution from question (histogram as well as distrubution) http://mathematica.stackexchange.com/questions/31948/show-density-plot-on-us-map#=. But both cases are a bit more difficult in my ...
2
votes
1answer
48 views

TransformedDistribution using n iid random variables

How do I create a TransformedDistribution that uses $n$ independent identically distributed (i.i.d.) random variables? For example, I can derive a chi-squared ...
0
votes
0answers
40 views

Using List/Histogram Data As A Weighted Function?

I have a list of data points that I would like to use as a weighted fucntion. The data represents the number of particles within a bin of a particle distribution. Ultimately I would like to multiply ...
2
votes
1answer
137 views

on Kolmogorov-Smirnov test

I found a problem to calculate Kolmogorov-Smirnov test statistic or p-value the data is ...
1
vote
1answer
110 views

In there a way to derive the inverse value of CDF format MultinormalDistribution?

Is there a way to produce the inverse value of multinormal distribution function in CDF format with 0.05 significance level? I found the function called ...
2
votes
0answers
93 views

Distance between vectors does not match distance between probability distribution [closed]

I have vectors that describe shapes. I am doing Euclidian distance between vectors to get the similar shapes and I get a good results of all classes. Theoretically, two similar shapes (shapes in a ...
0
votes
0answers
49 views
1
vote
0answers
32 views

Dependent probabilities distributions?

What's the best way to model situations with dependencies--I pick a world state, then pick a distribution to draw samples from based on that state. Here's a concrete version: If I choose a 6-sided ...
2
votes
1answer
62 views

Statistical Analyzis of Observed Data (Queuing Process)

Basically I have a process that belongs to a certain state at any point of time. However, the time intervals are not fixed (in fact they are exponentially distributed). Essentially, the data points ...
0
votes
1answer
40 views

Understanding `TransformedDistribution` behavior on two similar lines of code

The two following lines of code does not give the same output. Mathematica manages to calculate the first one but fails to calculate the second one. ...
0
votes
0answers
21 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: ...
1
vote
0answers
41 views

Why does integration of a radical times HeavisideTheta give a conditional expression?

If I want to compute an integral over x^p HeavisideTheta[y-x] where p is some fraction, I do not get the desired result, but ...
1
vote
1answer
121 views

Where can I find a good tutorial for stochastic modeling using Mathematica?

I am a newbie in Mathematica and I aim to use Mathematica mostly for stochastic modeling (to develop and analyze models in evolutionary biology). Note that I am not necessarily a great mathematician ...
1
vote
1answer
135 views

What is the PDF of a variable where a parameter is itself a random variable?

Question How do we analyze in Mathematica cases where a parameter in a Probability Density Function (PDF) varies according to another PDF? Consider the following example... Example The text below ...
-1
votes
1answer
142 views

Uniform distribution on unit circle

I need a uniform distribution on the unit circle (x; y) : x^2 + y^2 <= 1 with the density f(x) = (2/Pi) Sqrt[1-x^2] ...