Questions on the statistical functions of Mathematica.

learn more… | top users | synonyms

-1
votes
0answers
26 views

Experiment design in statistics [on hold]

Company X has 1000 employees and they are required to work 7 days a week. During one week,, a reporter collect attendance data that showed there were 200 employees who missed one day of work due ...
1
vote
1answer
68 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 ...
1
vote
1answer
52 views

Finding the true expectation

I have the density function and I want to determine its expected value ...
7
votes
1answer
115 views

How to derive SinglePredictionBands from a NonlinearModelFit model

I have found NonlinearModelFit to be extremely useful in my work. Here I have fit 4 points to a function and am being asked to explain how the 95% SinglePredictionBands associated with the plot below ...
4
votes
2answers
156 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
250 views

How to plot an implicit value funtion

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

PairwiseScatterPlot with histograms along the diagonal

The built-in function PairwiseScatterPlot has been discussed for example here: Scatterplots for subsets of multivariate data [Updated]. However, I find it not very ...
1
vote
1answer
149 views

Modify BinList to give elements indices in list instead of giving the list elements belonging to each bin

How could I use BinList to give elements indices in list instead of giving the list elements belonging to each bin ? Thanks.
27
votes
1answer
562 views

How to augment the realm of functions Mathematica thinks it knows how to integrate symbolically

My question involves extending the functionality of Integrate over specific integrals in the most generic manner. Specifically, is it possible to "hack into" ...
3
votes
1answer
81 views

How to make a diagnostic plot for the regression?

Context I am interested in LinearModelFit, and would like to make diagnostic plots using Mathematica. I cannot find them in Mathematica document center. For example, the following R program makes ...
2
votes
2answers
56 views

LogLikelihood[MultinomialDistribution[],{}] Not working?

LogLikelihood[MultinomialDistribution[100, {0.2, 0.3, 0.5}], {{10, 20, 70}}] -12.8338 No problem at all. ...
-2
votes
0answers
48 views

Statical Analyze a PDF (document file) for

I would like to import a set of PDF files to Mathematica and analyze it. And find out which words used most and find meaningful information about them. Finally, I would like to find out are they use a ...
2
votes
2answers
103 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 ...
0
votes
0answers
42 views

Robust (White) estimation with LinearModelFit?

You can make a robust estimate (White) by LinearModelFit or GeneralizedLinearModelFit? Thanks.
4
votes
2answers
104 views

Correlation test issues

No error in 10.0.2 When using Correlation[] function I sometimes get the strange warning: ...
1
vote
1answer
68 views

Method of moments using just the moments?

FindDistributionParameters (and EstimatedDistribution) can use "MethodOfMoments" as its ...
0
votes
0answers
32 views

Mixture / Hierarchical Distribution from Conditional

How does Mathematica handle Hierarchical distributions. For example, if I have X|Y ~ Poisson(Z) Y ~ Exp(lambda) Z ~ Normal(0,1) First, of all, the distribution of X is only given conditionally and ...
0
votes
1answer
70 views

probability calculation with two independent normal RVs

I don't understand why Probability[X < Y, {X \[Distributed] NormalDistribution[0, 1], Y \[Distributed] NormalDistribution[0, 1]}] gives 1/2 while we know ...
2
votes
0answers
58 views

The CorrelationMatrix is equal to corrcoef of Matlab?

I know that my question imply to know matlab, and this is not the right place, but i'm a physics student and i don't quite yet understand this concepts. When i make a LinearModelFit i can get a ...
1
vote
1answer
50 views

Find a probability with random variates from the same distribution

I'm trying to find Probability involving 3 random variates from the same distribution, but Mathematica does not seem to solve it. ...
0
votes
1answer
41 views

Central Limit Theorem - solve for a const, given probability

I hope this is not a very complicated question, but I have problems finding the answer. Quite often, I need to calculate something like this: $P(-a<\bar{x}<a) = p$, where $\bar{x}$ refers to the ...
2
votes
1answer
80 views

plot the numerical value of the average on a BoxWhiskerChart

Im using BoxWhiskerChart to get a boxplot of a matrix of two columns. I need to read the value of the sample average directly from the plot and not from the axes. ...
22
votes
3answers
673 views

Which Distributions can be Compiled using RandomVariate

Recently, Oleksandr kindly showed a list of Mathematica commands that can be compiled. RandomVariate was part of that list. However, whether this can be compiled depends upon the distribution that is ...
3
votes
2answers
151 views

How can I best plot timestamps by frequency-to-date?

I have a site called Sudomemo. I record when I first see a user in the database; and I've exported it in JSON format, located here. http://www.sudomemo.net/statistics/firstSeenDump.php I've done ...
4
votes
2answers
171 views

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

I calculate the Kolmogrov-samirnov test statistic for the following data: ...
5
votes
2answers
299 views

Numerical Integration with InverseErfc

I am trying to numerically integrate an equation that involves InverseErfc (embedded in the copula defined). The equation looks like the following: $$ \int_0^T \int_0^\infty \int_0^\infty ...
3
votes
2answers
154 views

How do I display the plane of best fit in 3D?

I'm trying to display the plane of best fit in the same 3D-box. Currently, I have the following code which spits out a plot of the points in 3D. I want to have a plane of best fit, can Mathematica do ...
0
votes
0answers
45 views

Comparing heavy tailed and non-heavy tailed distribution using Mathematica

I have data from a heavy tailed distribution (see picture 1) and a non-heavy tailed distribution (see picture 2) and want to compare it using a histogram on log-scale. However the picture 3 does not ...
2
votes
1answer
128 views

on Kolmogorov-Smirnov test

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

Distribution of block maxima of gaussian random variables

I try to verify the following result of extreme value statistics that the distribution of the block maxima of Gaussian random variables converges to a Gumbel distribution. To this end, I first ...
1
vote
1answer
107 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 ...
3
votes
2answers
121 views

RandomFunction and Histograms

RandomFunction[QueueingProcess[3, 5], {0, 15}] Histogram[data, Automatic, "PDF"] gives a very nice way to see the histogram based on the specified process. ...
2
votes
1answer
58 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 ...
7
votes
1answer
396 views

Estimation of vector autoregressive (VAR) process

To estimate 1-dimensional VAR process (which is AR process) you can easily use functions EstimatedProcess or ...
14
votes
2answers
1k views

Obtaining joint distributions and conditional distributions using Mathematica

I have two multivariate Gaussian distributions $p(x)$ and $p(z)$ with mean vectors $m_x$ and $m_z$, and covariance matrices $\Sigma_x$ and $\Sigma_z$. my model is a simple linear model $x = W z+n$ ...
5
votes
2answers
137 views

Getting the Cumulative Distribution Function of f(x) from a dataset of x

I've a data set $\mathcal{X}$ which consists of randomly generated numbers. Let's say for example data = RandomVariate[ExponentialDistribution[1], 10^4]; My aim ...
1
vote
3answers
158 views

Plotting a frequency function

I have a list of numbers and I want to plot the fraction of numbers >= k against k which will naturally be a decreasing curve. How do I do that?
0
votes
1answer
99 views

On plotting the Weibull distribution

First: How can I plot the Weibull distribution with parameters a = 2 and b = 3.5 on the range of ...
4
votes
2answers
152 views

QuantilePlot silliness

When I do a quantile plot of two large data sets (in my case, each of them of size around 3 million), the resulting notebook cell is huge, and if you save it as PDF, the resulting PDF is 13MB in size. ...
2
votes
1answer
105 views

Choosing The Best Fitted Probability Distribution Model

The following Example 20 Obs and i used Goodness of fit test with Exponential distribution, another time with Weibull distribution and third time with Pareto distribution .All of them are significant ...
1
vote
1answer
105 views

Estimating Confidence Interval Variance covariance matrix from weibull distribution

In The following Example 50 obs . I estimated distribution parameters (by point) but *i want to estimate parameters by confidence interval and calculate covariance matrix i need help * ...
2
votes
1answer
61 views

Meaning output ANOVA post hoc tests

I am carrying out an ANOVA in Mma, and I wanted to apply a post-hoc test, but I do not understand the output that is given. It says e.g. Model-> Duncan, {4, 5} ...
0
votes
1answer
77 views

Estimating Parameters of a Custom Distribution

The problem is I can not get estimation of parameters of any custom probability distribution (distributions are not built in Mathematica), here is an example: First: I entered random numbers of ...
2
votes
1answer
86 views

Probability evaluation fails with equals condition

The following evaluation fails (the result is just the copy of the input) $$\text{Probability}[0<x\leq a | x-y=t,\{x\ {ExponentialDistribution}[\lambda ],y\ {ExponentialDistribution}[\lambda ...
23
votes
5answers
4k views

Showing the correlation of two variables using a plot

I have a list of $\{x_i,y_i\}$ pairs and I want to show that $x_i$'s and $y_i$'s have positive correlation. The result of ...
0
votes
1answer
584 views

PDF of binomial distribution and mixed binomial distribution

I want to plot the PDF of binomial distributed function and a mixed binomial distributed function, therefore I entered for the binomial distribution function the following: ...
2
votes
0answers
80 views

How to include confidence intervals into ProbabilityScalePlot?

I use ProbabilityScalePlot[data, "Weibull"] to plot measured data and ...
4
votes
1answer
155 views

Calculate the covariance of a large matrix

I have to compute the covariance of 50 very large integer matrices (2500x2000 elements). However, according to my estimation this will take around 10 days. Do you have any ideas how to speed things ...