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)

0
votes
0answers
149 views

Binomial Distribution (binomial model inspired by Merton)

I want to do the following and I hope that anyone of you can help me out. In short, I have a portfolio consisting of 10 companies, whereas their asset values are driven by a common economic ...
16
votes
2answers
966 views

How to compute the inverse CDF of HyperbolicDistribution properly?

Fixed in version 9. I want to compute the CDF and inverse CDF of the hyperbolic distribution: ...
4
votes
1answer
70 views
6
votes
2answers
209 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 ...
5
votes
2answers
444 views
9
votes
3answers
586 views

The default bandwidth of the SmoothKernelDistribution function

I am researching the Silverman rule of thumb for bandwidth selection. Below is a sample of my code, from which it seems that BW=63 is close to the default parameter given by Mathematica (2 images are ...
2
votes
1answer
44 views

Documentation for the Bounded option of SmoothKernelDistribution

I have found the Bounded option useful in the SmoothKernelDistribution function, but I can't find any documentation for exactly ...
0
votes
1answer
33 views

Error message when exporting plot as XLS file

I imported data from excel which has 300 data (which are change in CDS spread). I want to compute corresponding CDF for the data using kernel density estimation. I get the CDF plot but cannt get the ...
0
votes
1answer
52 views

kernel density estimation

I am trying to compute CDF of historical data ( weekly change in CDS spread) using kernel density estimation. The original data has 300 observation. I uploaded data from excel ...
1
vote
0answers
39 views

How should I iteratively refine a MixtureDistribution?

I am trying to find a stopping time for a randomized simulation such that say 95% of trials that will be successful will have already finished. I would like to do this by creating and then refining ...
10
votes
1answer
306 views

How does LocationTest choose its “AutomaticTest”?

Is there a way of determining how LocationTest chooses its "AutomaticTest"? Sometimes it's clear from the ...
0
votes
0answers
78 views
1
vote
1answer
80 views

How to simplify and evaluate Integral with lots of DiracDelta dists

I have a question about an Integrals which I can do by hand but I want to implement this Integral in Mathematica. Actually it is the definition of a 2 Particle Phase Space. The Expression is the ...
2
votes
3answers
168 views

Fitting and plotting Weibull distribution

How can I find the best fitted 2-parameter Weibull's distribution (PDF) to describe my data? How can i plot my data and estimated distribution all in one graph? How can I show/compare results of ...
5
votes
2answers
90 views

Dirichlet marginals - why does Mathematica not give the *Kth* marginal?

Observe: MarginalDistribution[DirichletDistribution[{a, b, c, d}], #] & /@ Range@4 {BetaDistribution[a, b + c + d], BetaDistribution[b, a + c + d], ...
27
votes
8answers
1k 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
2answers
86 views

TargetFunctions option in FindDistribution yields error

I have some data set that I am trying to find the maximums . I figured I could do this by fitting the many peaks of the dataset to Gaussian distributions and then finding the peaks of this data from ...
0
votes
2answers
142 views

Random non-overlapping disks in a square

I've found a nice code snippet by @belisarius from this question, that I'll reproduce here for reference: ...
0
votes
1answer
46 views

How to get DistributionChart to show normalised densities?

I am having some issue getting DistributionChart to show bin widths as expected, or perhaps it is in my interpretation of the plots. I have four data sets, all of the same length, and plot them in ...
1
vote
1answer
111 views

How can I draw a box with moving particles inside and a hole in one side? [closed]

The idea is to create a simulation for ideal gas motion in a container that have a hole and then find the speed distribution for the gas particles that leave the box through the hole. I have a ...
2
votes
1answer
159 views

Confusion with SplineFit (Angle calculation)

This question is continuation to this other one I am creating polymers (where each monomer is of equal length) using this method: ...
3
votes
2answers
79 views

Efficient Sampling of A (Countable) Discrete Distribution

I have a sequence of positive terms $(a_n)$, for which $\sum_{n=1}^\infty a_n = A <\infty$, and wish to take a psuedo-random sample from the discrete probability distribution \begin{equation} ...
2
votes
1answer
86 views

P-Values are inconsistent when using DistributionFitTest repeatedly

I am attempting to become more familiar with the DistributionFitTest, specifically the distributions it yields and their fittedness. What I have found is that when I generate random data, fit it, and ...
1
vote
0answers
103 views

Summing the probability distribution to 1 to convince myself

I have worked out a probability distribution and want to check its sum which is necessarily 1. First we write $$ r \triangleq \frac{(2 \lambda + \mu)^2}{2(\mu + \lambda)^2}, \quad s \triangleq ...
0
votes
1answer
153 views

Prime number The Ulam spiral [closed]

I want a simple method of visualizing the prime numbers that reveals the apparent tendency of certain quadratic polynomials to generate unusually large numbers of primes. I was able to display the ...
2
votes
1answer
112 views

Mathematica equivalent of `qnorm` or `invNorm`

In R and TI-89 you can do qnorm(prob, mu, sigma) or invNorm(prob, mu, sigma) to get a value of ...
11
votes
4answers
244 views

Frequency of punctuation characters used in Mathematica languages

I would like to determine the computer language character frequency in Mathematica code. I want to remap my keyboard to the most used character such as [], {}, -, +. There are the frequency of ...
0
votes
0answers
44 views

Multinormal Distribution with an affine transformation

The affine property of Gaussians says: if $x \sim N_x(\mu_x, \Sigma_x)$, one can apply an affine transformation to $x$, $y = Ax + b$, which gives the distribution --> $y \sim N_y(A\mu_x+b, A\Sigma_x ...
1
vote
1answer
148 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 ...
0
votes
1answer
230 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 $ ...
4
votes
1answer
214 views

Heuristic Method to create low discrepancy sequences

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 ...
0
votes
1answer
92 views

Solve ${\rm BinomialDistribution}(x;\,n_1,\,p_1) = {\rm BinomialDistribution}(a-x;\,n_2,\,p_2)$ in terms of $x$ [closed]

I'm totally a novice with Wolfram software, so please forgive me if the question is nonsense. I want to solve the following equation in terms of $x$: $\quad \quad {\rm ...
24
votes
3answers
3k views

Standard errors for maximum likelihood estimates in FindDistributionParameters

I am trying to fit a non-standard PDF to data and FindDistributionParameters works great, and gives me the parameters of the distribution back using maximum ...
16
votes
1answer
460 views

Speed up adding Poisson noise to an image

I'm simulating some images corrupted with Poisson noise, but I'm encountering a few problems with performance. According to the documentation on ImageEffect, one ...
5
votes
2answers
551 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
72 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 ...
1
vote
0answers
36 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
41 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 ...
2
votes
2answers
55 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
529 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
70 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
38 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
73 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
44 views
4
votes
2answers
116 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
198 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.