Questions on the statistical functions of Mathematica.

learn more… | top users | synonyms

0
votes
2answers
30 views

Simulate from an AR(1) process where the parameters vary over time

I would like to simulate from an AR(1) process, where the $\rho$ parameter in the process of the form: $X_t = \rho X_{t-1} + \epsilon_t$ varies over time. The path of $\rho$ and $var(\epsilon_t)$ ...
2
votes
2answers
111 views

How can I make an animated demonstration of the Wigner semicircle distribution?

With Mathematica 10, it's possible to draw the semicircle distribution of R. Wigner: ...
1
vote
2answers
61 views

Creating histograms with specified mean and intervals

I have a set of data that I want to create a histogram for. The thing is, I would like to set the center of the histogram to the mean of the dataset and the intervals to be the standard deviation. How ...
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 ...
3
votes
2answers
68 views

PredictionBands and the presence of {0,0} point

We know from this post, that the default fitting method cannot handle data with a {0,0} point, and that there are alternative methods. But, with these alternative methods, the prediction bands ...
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 ...
1
vote
2answers
73 views

Performance in calculating maximum-likelihood- based estimates

My probability density function is a complicated one for which numerical estimation is necessary. Here my pdf: ...
2
votes
0answers
54 views

Computational Bayesian analysis in Mathematica: Any plans to develop MCMC?

Does anyone know if there are any plans to develop an MCMC capability in Mathematica? My reasoning for asking is that as it stands, I can't seem to find any 'out-of-the-box' functions/capabilities ...
1
vote
1answer
27 views

LogitModelFit bias term when using design matrix

I'm using the LogitModelFit function with a 32-columns design matrix and 0 or 1 output vector: LogitModelFit[{datax, datay}]. ...
0
votes
0answers
63 views

Autocorrelation of 3D Points

I'm trying to compute the autocorrelation of a list of points (x,y,z) in 3D. Can the CorrelationFunction or maybe ListCorrelate ...
0
votes
1answer
34 views

ANOVA on data with different lengths

I have data in the following form: ...
0
votes
1answer
68 views

how to scale the plot of PDF (t-distribution)

The directions in my book say: "To sketch the t distribution in Figure 3.11, simply multiply the abscissa t value by the scale factor and plot this against the ordinate of t at that point." So ...
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}\\ ...
1
vote
1answer
59 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
1
vote
0answers
72 views

How can I model the volatility, using a GARCH(1,1), of a time series of returns (and plot it) using Mathematica?

I know that if I use this code: tsm = TimeSeriesModelFit[ret, {"GARCH", {1, 1}}]; I can get the parameters using this: ...
5
votes
1answer
66 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
60 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
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 ...
1
vote
1answer
56 views

Finding the true expectation

I have the density function and I want to determine its expected value ...
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. ...
3
votes
1answer
92 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 ...
3
votes
0answers
186 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 ...
2
votes
2answers
57 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
105 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.
0
votes
0answers
35 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
73 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 ...
1
vote
1answer
71 views

Method of moments using just the moments?

FindDistributionParameters (and EstimatedDistribution) can use "MethodOfMoments" as its ...
2
votes
0answers
79 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
55 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
45 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 ...
1
vote
1answer
88 views

How to group box-and-whisker in parallel for comparison in BoxWhiskerChart?

I have two collections of data: xdata and ydata ...
2
votes
1answer
83 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
159 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
193 views

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

I calculate the Kolmogrov-samirnov test statistic for the following data: ...
3
votes
2answers
166 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
54 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 ...
0
votes
1answer
57 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 ...
2
votes
1answer
136 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
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 ...
3
votes
2answers
125 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. ...
7
votes
1answer
118 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 ...
5
votes
2answers
142 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 ...
0
votes
1answer
113 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. ...
1
vote
1answer
117 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
115 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 ...