Questions on the statistical functions of Mathematica.

learn more… | top users | synonyms

5
votes
2answers
520 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 ...
5
votes
3answers
259 views

How to find max difference with the same x between two functions?

For example, I have function TF[y_] = 2 / Pi * ArcSin[y] And I have CDF: ...
4
votes
1answer
74 views

ANOVA test on dataset

I have a dataset with 3 groups and two variables. I want to do ANOVA test on the three groups for each variable. Does anyone know how to do this? ...
0
votes
0answers
6 views

How to find the mode when none of the terms repeat? [migrated]

What is the mode of a series / distribution of values when none of the values of the distribution are repeated , for example : What is the mode of 1 , 3 , 45, 67 , 89 , 100
2
votes
0answers
25 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 ...
0
votes
0answers
39 views

Better chances on a gambling game [closed]

I'm working on a program that I hope it will help me to get better chances on a gambling game. I have a very large database of numbers how occur in this game and my program pass through those numbers ...
3
votes
4answers
282 views

How to plot frequency polygon?

I have data sample. Built and plotted histogram based on the sample. How can I plot frequency polygon now? That's example of the desired result. Here is code: ...
1
vote
0answers
68 views

Truncated Distribution Parameter Search Produces Singular Jacobian Using FindRoot [closed]

I'm trying to find the 2 parameters of a truncated Weibull distribution using version 10 and FindRoot. I would like to truncate the distribution between t =0 and t = 3600 with the probability of ...
0
votes
0answers
34 views

Attributes of statistics functions

Do statistics functions have some kind of attribute which ignore non-number (Indeterminate, ComplexInfinity, calculation "errors") elements of a list?
2
votes
2answers
44 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 ...
1
vote
1answer
52 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}]. ...
1
vote
1answer
73 views

How do I create a histogram of the binomial distribution for, say, p=.5 after 10 trials?

I apologize if this is question is too simple, if so it's because I'm a beginner. Anyway, it was easy enough to create a list plot, but I couldn't find out how to create an output that displays the ...
1
vote
2answers
42 views

BoxWhiskerChart Vertical Line

I have made a box plot to show the spread of past paper results for my A2 exams so I can easily see which subject need more work! ...
10
votes
2answers
262 views

How to remove outliers from data

I have a data with noise which some times includes significant outliers. The position of the outliers are random. For example: ...
0
votes
1answer
46 views

Binomial Regression where n > 1

I am not well versed in Generalized Linear Models, but am trying to use this approach to model some discrete data. I have binomial data where the sample size for each observation is the same size, ...
0
votes
0answers
17 views

Derived distribution

I'm trying to solve a parametric DE in the unknown z(t,r), where the distribution of r is known, and obtain the derived ...
0
votes
1answer
52 views

Interpreting P-Value on ParameterTable

When using LinearModelFit[] we can generate a table of parameter values and statistics for the corresponding model using the property "ParameterTable". Looking in the documentation I can't find an ...
2
votes
1answer
144 views

How to draw confidence ellipse from a covariance matrix?

I am studying a two-dimensional dataset, whose mean vector and covariance matrix are the following: ...
26
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 ...
1
vote
1answer
74 views

how to calculate the covariance of sample central moments

$\bar{x}= \sum_{i=1} ^{n} \frac{1}{n}x_{i}$ $m_1=\sum _{i=1} ^{n} \frac{1}{n} (x_{i}-\bar{x})$ $m_2=\sum _{i=1} ^{n} \frac{1}{n} (x_{i}-\bar{x})^2 $ $m_3=\sum _{i=1} ^{n} ...
28
votes
1answer
609 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" ...
2
votes
1answer
38 views
4
votes
2answers
113 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: ...
8
votes
1answer
271 views

How can I improve my prediction for the end time of the Reddit button game?

There is an April Fools game on Reddit where users with accounts can press a button (once) to reset a 60 second countdown timer: https://www.reddit.com/r/thebutton ...
0
votes
0answers
42 views

Statistical significance test of FindCluster

I have a list of vectors which could be roughly divided into two clusters. The list is huge so I have to use RandomChoice[list, 20000] to speed up the calculation. ...
0
votes
0answers
46 views

How to Fit a Generalized Pareto Distribution to data in Mathematica

I'd like to know how to fit a Generalized Pareto distribution in Mathematica to implement the Peak Over Threshold method for tail estimation of operational risk losses. Does anybody have a toy example ...
0
votes
0answers
37 views

Pdf for sum of unknown named distribution and normal random variable

Is it possible to apply convolution theorem or software like Mathematica to find a closed form expression for the pdf of Z = R + X where $f_R(r;k,d) = kdr^{d-1}{(1-r^d)}^{k-1}$ and $X$ is zero mean ...
2
votes
1answer
70 views

Find confidence interval for negative binomial using Mathematica

I would like to find the exact 95% confidence intervals on the mean of a negative binomial random variable using Mathematica. It is possible to approximate the distribution using a gamma/normal, but ...
4
votes
0answers
69 views

RandomVariate, MSE and Cramer-Rao bound calculation

I'm using the following code to plot the variance of a moment estimate of a parameter of Rice distribution and its' Cramer-Rao bound. First I generate $10^3$ vectors $\texttt{R}$ ($10^3$ elements ...
2
votes
0answers
58 views

ANCOVA (ANalysis of COVAriance)

I have two time series of the form {time, value} I have linear fits of each set individually, and I am interested in the statistical differences between the ...
0
votes
1answer
35 views

Import command gives different solution to equation than cut & paste

I would like to import values from a text file to run a notebook. The notebook works fine when I cut and paste in all the input values manually. Here is a working version. But when I import the same ...
5
votes
1answer
285 views

BoxWhiskerChart - 10-90th percentile instead of min - max

I am trying to draw a custom BoxWhiskerChart with alternative values determining fences. I have data that has a lot of outliers and I have to compare it to another chart that has 10-90th percentile ...
0
votes
1answer
184 views

Random numbers in C++ and Mathematica gives wrong results? [closed]

I have the following strange phenomenon which puzzles me!: I have a piecewise constant probability density given as ...
2
votes
3answers
203 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.
0
votes
2answers
37 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)$ ...
1
vote
2answers
83 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 ...
2
votes
2answers
155 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: ...
3
votes
2answers
72 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
34 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
76 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
94 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 ...
0
votes
0answers
70 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
44 views
0
votes
1answer
105 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
70 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
0answers
105 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: ...
1
vote
1answer
89 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 ...