Questions on the statistical functions of Mathematica.

learn more… | top users | synonyms

1
vote
1answer
90 views

How do you graph the confidence intervals that results from a GARCH over a time series?

Please, does someone know how to draw a confidence intervals (from the results of a GARCH (1,1)) in a way similar to these: *Asume an ARMA(1,1) fit the time series. Thank you very much
5
votes
1answer
78 views

How to formally tell if one time series affects another?

In time series analysis, a correlogram, also known as an autocorrelation plot, is a plot of the sample autocorrelations. I'm rolling my own now and I'm not a statistician, but I think this sort of ...
0
votes
4answers
65 views

Mathematica & Statistics? [closed]

So, let's say you flip a coin 10 times, and call that 1 "event". If you run, 1,000,000 of these "events", what is the proportion of events that have heads between 0.4 and 0.6? Binomial probability ...
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 ...
1
vote
1answer
70 views

Finding the true expectation

I have the density function and I want to determine its expected value ...
0
votes
0answers
42 views

How do I add custom deviation when using LinearModelFit?

I have a set of data in the form a list S1={{x,y},{x,y}} etc. I also have a function that calculates a corrected standard deviation in the y direction for ...
7
votes
1answer
139 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
178 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
291 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
193 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 ...
3
votes
1answer
104 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
59 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
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 ...
0
votes
0answers
43 views

Robust (White) estimation with LinearModelFit?

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

Correlation test issues

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

Method of moments using just the moments?

FindDistributionParameters (and EstimatedDistribution) can use "MethodOfMoments" as its ...
0
votes
0answers
36 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
77 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
97 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
69 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
50 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
93 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. ...
3
votes
2answers
179 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 ...
5
votes
2answers
324 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
304 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
211 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
69 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
179 views

on Kolmogorov-Smirnov test

I found a problem to calculate Kolmogorov-Smirnov test statistic or p-value the data is ...
0
votes
1answer
68 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
126 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
137 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
67 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
472 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
171 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
167 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
140 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
154 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
158 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
141 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
77 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
106 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
88 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 ...
0
votes
1answer
638 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
107 views

How to include confidence intervals into ProbabilityScalePlot?

I use ProbabilityScalePlot[data, "Weibull"] to plot measured data and ...
4
votes
1answer
172 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 ...
3
votes
1answer
59 views

TransformedDistribution with respect to new variable

Given a custom pdf $f_x(x)$, I'm trying to find it's transformation $f_y(y)$ where $$y=x^2$$ and $$f_x(x)=30*x^2 (1 - x)^2, 0<x<1$$ I tried to using the following commands: ...
2
votes
0answers
59 views

k-means clustering using squared correlation in Mathematica

I am trying to cluster data using the squared correlation coefficient since in the case of my data, -1 correlation is equivalent to 1 correlation (polarities do not matter). How can this be done in ...