# Statistics

**URL:** https://discourse.julialang.org/c/domain/stats/12.md?page=8

[Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md)

**Page:** 9

---

## [LsqFit.jl Residuals Interpretation](https://discourse.julialang.org/t/lsqfit-jl-residuals-interpretation/88011)

<div class="topic-metadata">

**Author:** [@wiseac](https://discourse.julialang.org/u/wiseac)\
**Replies:** 4\
**Last updated:** [September 30, 2022, 8:03am UTC](https://discourse.julialang.org/t/lsqfit-jl-residuals-interpretation/88011 "2022-09-30T08:03:16Z")

</div>

I have a question regarding interpreting the results of a residual plot. Basically light is decaying inside a cell and we are doing an exponential fit to get the decay time (In layman terms) (Cavity RingDown). There are …

---

## [Quadratic (paraboloid) model fitting](https://discourse.julialang.org/t/quadratic-paraboloid-model-fitting/88007)

<div class="topic-metadata">

**Author:** [@bertulli](https://discourse.julialang.org/u/bertulli)\
**Replies:** 0\
**Last updated:** [September 29, 2022, 7:37pm UTC](https://discourse.julialang.org/t/quadratic-paraboloid-model-fitting/88007 "2022-09-29T19:37:47Z")

</div>

Hi, this is a off-branch topic from this other thread of mine about multiple linear regression. This is not the solution to the original question, but I thought it’s useful, so I’m making a recap here. Turned out, linea…

---

## [MultivariateStats.jl: predict for MulticlassLDA is throwing error](https://discourse.julialang.org/t/multivariatestats-jl-predict-for-multiclasslda-is-throwing-error/87799)

<div class="topic-metadata">

**Author:** [@Sourish](https://discourse.julialang.org/u/Sourish)\
**Replies:** 2\
**Last updated:** [September 27, 2022, 2:56am UTC](https://discourse.julialang.org/t/multivariatestats-jl-predict-for-multiclasslda-is-throwing-error/87799 "2022-09-27T02:56:16Z")

</div>

Hi I am trying to fit MulticlassLDA on the iris dataset. The predict is throwing an error. I was wondering if you can help me to figure out what mistakes I am making. using MultivariateStats, RDatasets; iris = dataset…

---

## [How can I tell Julia what are the parameters in a linear model?](https://discourse.julialang.org/t/how-can-i-tell-julia-what-are-the-parameters-in-a-linear-model/87692)

<div class="topic-metadata">

**Author:** [@bertulli](https://discourse.julialang.org/u/bertulli)\
**Replies:** 14\
**Last updated:** [September 26, 2022, 5:24pm UTC](https://discourse.julialang.org/t/how-can-i-tell-julia-what-are-the-parameters-in-a-linear-model/87692 "2022-09-26T17:24:11Z")

</div>

Hi all! So, I’m experimenting with the linear regression ecosystem (GLM, Plots, DataFrames, StatsPlots and CSV). I understood how to produce a simple linear regression model and a multivariate one. I know that OLS work…

---

## [How do I fit TGARCH-X model using ARCHModels.jl](https://discourse.julialang.org/t/how-do-i-fit-tgarch-x-model-using-archmodels-jl/87558)

<div class="topic-metadata">

**Author:** [@ultratrader](https://discourse.julialang.org/u/ultratrader)\
**Replies:** 2\
**Last updated:** [September 21, 2022, 12:48pm UTC](https://discourse.julialang.org/t/how-do-i-fit-tgarch-x-model-using-archmodels-jl/87558 "2022-09-21T12:48:00Z")

</div>

How do I fit a GARCH model with leverage effect as well as additional explanatory variables using ARCHModels.jl library? I can’t find any documentation about this in: https://s-broda.github.io/ARCHModels.jl/stable/

---

## [Smoothing empirical CDFs](https://discourse.julialang.org/t/smoothing-empirical-cdfs/87544)

<div class="topic-metadata">

**Author:** [@amrods](https://discourse.julialang.org/u/amrods)\
**Replies:** 6\
**Last updated:** [September 21, 2022, 12:20pm UTC](https://discourse.julialang.org/t/smoothing-empirical-cdfs/87544 "2022-09-21T12:20:51Z")

</div>

I have 5 empirical CDFs coming from a model. I’d like to average them and plot that average against the empirical CDF of the data. The issue I’m facing is that the points of support of those 5 CDFs are not the same, so I…

---

## [Φ(−600) = 6.546588 205807 692852 105927 713888 108782 119412 831853 177211 16943e − 78177](https://discourse.julialang.org/t/600-6-546588-205807-692852-105927-713888-108782-119412-831853-177211-16943e-78177/87309)

<div class="topic-metadata">

**Author:** [@blackeneth](https://discourse.julialang.org/u/blackeneth)\
**Replies:** 3\
**Last updated:** [September 16, 2022, 5:50am UTC](https://discourse.julialang.org/t/600-6-546588-205807-692852-105927-713888-108782-119412-831853-177211-16943e-78177/87309 "2022-09-16T05:50:19Z")

</div>

The paper Approximations to the Normal Distribution Function and An Extended Table for the Mean Range of the Normal Variables by Kiani, M and Panaretos, J and Psarakis, S and Saleem, M (2008) reviews algorithms for compu…

---

## [More flexible Distributions package](https://discourse.julialang.org/t/more-flexible-distributions-package/87016)

<div class="topic-metadata">

**Author:** [@pedro.rafae](https://discourse.julialang.org/u/pedro.rafae)\
**Replies:** 10\
**Last updated:** [September 14, 2022, 8:09pm UTC](https://discourse.julialang.org/t/more-flexible-distributions-package/87016 "2022-09-14T20:09:13Z")

</div>

Dear friends, I am a statistician and I have all my training in this area. In my PhD I worked with probability distribution generators and recently I became interested in the Julia language. I come from the R language an…

---

## [Suggestions for the design of Survey.jl?](https://discourse.julialang.org/t/suggestions-for-the-design-of-survey-jl/86381)

<div class="topic-metadata">

**Author:** [@iuliadmtru](https://discourse.julialang.org/u/iuliadmtru)\
**Replies:** 3\
**Last updated:** [August 28, 2022, 8:18pm UTC](https://discourse.julialang.org/t/suggestions-for-the-design-of-survey-jl/86381 "2022-08-28T20:18:05Z")

</div>

Hello, We are building Survey.jl, a package for statistical estimates from surveys. The inspiration for the package comes from the corresponding R survey package created by Thomas Lumley. We started out imitating the R…

---

## [Automated distribution fitting](https://discourse.julialang.org/t/automated-distribution-fitting/30161)

<div class="topic-metadata">

**Author:** [@mike\_k](https://discourse.julialang.org/u/mike_k)\
**Replies:** 11\
**Last updated:** [August 26, 2022, 9:43pm UTC](https://discourse.julialang.org/t/automated-distribution-fitting/30161 "2022-08-26T21:43:33Z")

</div>

Dear community, I am not an expert in statistics, so maybe this question is a bit naive. I have a data vector with ~10^6 entries (between \[0,1\]) and want to identify a distribution which fits best to that data. I wonder…

---

## [Counting weeks in a timeseries](https://discourse.julialang.org/t/counting-weeks-in-a-timeseries/86181)

<div class="topic-metadata">

**Author:** [@chiraganand](https://discourse.julialang.org/u/chiraganand)\
**Replies:** 7\
**Last updated:** [August 25, 2022, 6:50pm UTC](https://discourse.julialang.org/t/counting-weeks-in-a-timeseries/86181 "2022-08-25T18:50:34Z")

</div>

Hello, For one of the codes I was working on I had to count the number of unique weeks in a timeseries. I ended up using the Dates.week() function and found that the result can be incorrect in cases where the data spans…

---

## [How to do a OneSampleZTest for the proportion instead of the mean?](https://discourse.julialang.org/t/how-to-do-a-onesampleztest-for-the-proportion-instead-of-the-mean/85779)

<div class="topic-metadata">

**Author:** [@floxo115](https://discourse.julialang.org/u/floxo115)\
**Replies:** 1\
**Last updated:** [August 22, 2022, 10:40am UTC](https://discourse.julialang.org/t/how-to-do-a-onesampleztest-for-the-proportion-instead-of-the-mean/85779 "2022-08-22T10:40:10Z")

</div>

I am playing around with the HypothesisTests package and wanted to do a OneSampleZTest but with the population proportion instead of the populatin mean as the parameter of interest. But I really don’t find anyting in the…

---

## [Query CDF and quantiles for MCMC Chains](https://discourse.julialang.org/t/query-cdf-and-quantiles-for-mcmc-chains/86088)

<div class="topic-metadata">

**Author:** [@TheLateKronos](https://discourse.julialang.org/u/TheLateKronos)\
**Replies:** 4\
**Last updated:** [August 21, 2022, 3:12pm UTC](https://discourse.julialang.org/t/query-cdf-and-quantiles-for-mcmc-chains/86088 "2022-08-21T15:12:09Z")

</div>

I have a Chains object chain\_normal: julia\> @show chain\_normal chain\_normal = MCMC chain (1000×14×4 Array{Float64, 3}) Chains MCMC chain (1000×14×4 Array{Float64, 3}): Iterations = 501:1:1500 Number of chains =…

---

## [CDF of a copula from Copulas.jl?](https://discourse.julialang.org/t/cdf-of-a-copula-from-copulas-jl/85786)

<div class="topic-metadata">

**Author:** [@HJW019](https://discourse.julialang.org/u/HJW019)\
**Replies:** 24\
**Last updated:** [August 19, 2022, 2:41pm UTC](https://discourse.julialang.org/t/cdf-of-a-copula-from-copulas-jl/85786 "2022-08-19T14:41:15Z")

</div>

I am experimenting with Copulas.jl to see if I could use it for my work. I am in general impressed. However, I found I could not get the CDF value from a joint distribution constructed by the copula. Getting the PDF valu…

---

## [\[FixedEffectModels.jl\] Switching from Dummy to Categorical Variables](https://discourse.julialang.org/t/fixedeffectmodels-jl-switching-from-dummy-to-categorical-variables/85800)

<div class="topic-metadata">

**Author:** [@teamtelefunken](https://discourse.julialang.org/u/teamtelefunken)\
**Replies:** 2\
**Last updated:** [August 16, 2022, 6:01pm UTC](https://discourse.julialang.org/t/fixedeffectmodels-jl-switching-from-dummy-to-categorical-variables/85800 "2022-08-16T18:01:13Z")

</div>

I am trying to use FixedEffectModels.jl with a dataset with a continuous explanatory variable, and want to run regressions on a binned version of this variable. I have two approaches: a) I can either, via cut, produce ca…

---

## [Error with Bayesian Gaussian Process "PosDefException: matrix is not positive definite; Cholesky factorization failed."](https://discourse.julialang.org/t/error-with-bayesian-gaussian-process-posdefexception-matrix-is-not-positive-definite-cholesky-factorization-failed/85683)

<div class="topic-metadata">

**Author:** [@vickiluv](https://discourse.julialang.org/u/vickiluv)\
**Replies:** 2\
**Last updated:** [August 13, 2022, 3:57am UTC](https://discourse.julialang.org/t/error-with-bayesian-gaussian-process-posdefexception-matrix-is-not-positive-definite-cholesky-factorization-failed/85683 "2022-08-13T03:57:04Z")

</div>

I got the error of “PosDefException” when trying Bayesian Gaussian Process in Turing.jl. The Bayesian framework looks as follows @model function GP\_μ(μ, D; cov\_fn::Function=exp\_cov\_fn) n\_station = length(μ) # w…

---

## [ANOVA Tests in Julia?](https://discourse.julialang.org/t/anova-tests-in-julia/1421)

<div class="topic-metadata">

**Author:** [@ibadr](https://discourse.julialang.org/u/ibadr)\
**Replies:** 76\
**Last updated:** [August 11, 2022, 7:23pm UTC](https://discourse.julialang.org/t/anova-tests-in-julia/1421 "2022-08-11T19:23:02Z")

</div>

What is the status of ANOVA tests in Julia? I looked at various packages under the JuliaStats organization, including HypothesisTests.jl, but I can’t seem to find this functionality. The closest I was able to find was in…

---

## [Simple multi-parameter bayesian linear regression returning un-expected results!](https://discourse.julialang.org/t/simple-multi-parameter-bayesian-linear-regression-returning-un-expected-results/81285)

<div class="topic-metadata">

**Author:** [@charperflow](https://discourse.julialang.org/u/charperflow)\
**Replies:** 14\
**Last updated:** [August 9, 2022, 11:20am UTC](https://discourse.julialang.org/t/simple-multi-parameter-bayesian-linear-regression-returning-un-expected-results/81285 "2022-08-09T11:20:05Z")

</div>

The problem is from chapter 10 question 7 of “Regression and Other Stories”, and the general idea is to use data which includes a continuous parameter “beauty” and then several binary predictors. I have solved the probl…

---

## [Pushing Julia/statistics development](https://discourse.julialang.org/t/pushing-julia-statistics-development/80111)

<div class="topic-metadata">

**Author:** [@Ajay\_Shah](https://discourse.julialang.org/u/Ajay_Shah)\
**Replies:** 14\
**Last updated:** [August 8, 2022, 6:52pm UTC](https://discourse.julialang.org/t/pushing-julia-statistics-development/80111 "2022-08-08T18:52:09Z")

</div>

Julia is great, but there are many limitations in putting it to work for doing statistics. With the help of Viral, Alan, and @bkamins, a bunch of us got going on a systematic attack on this problem. We have gone after th…

---

## [GaussianMixtures.jl - Equivalent of sklearn GaussianMixture.predict](https://discourse.julialang.org/t/gaussianmixtures-jl-equivalent-of-sklearn-gaussianmixture-predict/85480)

<div class="topic-metadata">

**Author:** [@vikram-s-narayan](https://discourse.julialang.org/u/vikram-s-narayan)\
**Replies:** 0\
**Last updated:** [August 8, 2022, 2:30pm UTC](https://discourse.julialang.org/t/gaussianmixtures-jl-equivalent-of-sklearn-gaussianmixture-predict/85480 "2022-08-08T14:30:49Z")

</div>

I’m trying to recreate the equivalent of SKLearn’s documentation example with GaussianMixtures.jl. The SKLearn example is below: import numpy as np from sklearn.mixture import GaussianMixture X = np.array(\[\[1, 2\], \[1, …

---

## [Copulas in Julia?](https://discourse.julialang.org/t/copulas-in-julia/10285)

<div class="topic-metadata">

**Author:** [@floswald](https://discourse.julialang.org/u/floswald)\
**Replies:** 17\
**Last updated:** [August 6, 2022, 11:24am UTC](https://discourse.julialang.org/t/copulas-in-julia/10285 "2022-08-06T11:24:47Z")

</div>

hi all, I am looking for a native julia Copula package. I have come across https://github.com/JuliaFinMetriX/Copulas.jl but it seems quite something to install - plus I can’t figure out the most basic things from the d…

---

## [Polynomial fitting with all fitting feedback (chi-squared, etc.)](https://discourse.julialang.org/t/polynomial-fitting-with-all-fitting-feedback-chi-squared-etc/85122)

<div class="topic-metadata">

**Author:** [@CosmoProf](https://discourse.julialang.org/u/CosmoProf)\
**Replies:** 1\
**Last updated:** [August 5, 2022, 9:00pm UTC](https://discourse.julialang.org/t/polynomial-fitting-with-all-fitting-feedback-chi-squared-etc/85122 "2022-08-05T21:00:33Z")

</div>

Hi, I’m familiar with ‘’‘PolyFit’‘’ for fitting polynomials to data, but that algorithm just returns the best-fit polynomial coefficients. Is there any easy way to automatically get back more of the useful fitting result…

---

## [Should model matrix for nested factors be full-rank?](https://discourse.julialang.org/t/should-model-matrix-for-nested-factors-be-full-rank/84733)

<div class="topic-metadata">

**Author:** [@PharmCat](https://discourse.julialang.org/u/PharmCat)\
**Replies:** 15\
**Last updated:** [August 1, 2022, 8:45pm UTC](https://discourse.julialang.org/t/should-model-matrix-for-nested-factors-be-full-rank/84733 "2022-08-01T20:45:00Z")

</div>

In some cases when levels of one factor all “inside” levels of another we have not full-rank model matrix from StatsModels. For example - calssical bioequivlence crossover study (2X2) with factors: Var ~ Formulation + …

---

## [Iterative proportional fitting package?](https://discourse.julialang.org/t/iterative-proportional-fitting-package/60580)

<div class="topic-metadata">

**Author:** [@amrods](https://discourse.julialang.org/u/amrods)\
**Replies:** 4\
**Last updated:** [July 27, 2022, 10:49am UTC](https://discourse.julialang.org/t/iterative-proportional-fitting-package/60580 "2022-07-27T10:49:25Z")

</div>

Is there any package that implements the iterative proportional fitting procedure?

---

## [How to estimate uncertainty in a LinearMixedModel?](https://discourse.julialang.org/t/how-to-estimate-uncertainty-in-a-linearmixedmodel/84124)

<div class="topic-metadata">

**Author:** [@jar1](https://discourse.julialang.org/u/jar1)\
**Replies:** 8\
**Last updated:** [July 14, 2022, 2:14pm UTC](https://discourse.julialang.org/t/how-to-estimate-uncertainty-in-a-linearmixedmodel/84124 "2022-07-14T14:14:37Z")

</div>

In this example, I want to estimate the uncertainty around the beta value for each Cut & Clarity pair. I’m not sure how to interpret the allpars data. Is it producing a single sigma shared across all values of Cut & Clar…

---

## [GeneralizedLinearMixedModel no method matching \_ranef\_refs](https://discourse.julialang.org/t/generalizedlinearmixedmodel-no-method-matching-ranef-refs/82859)

<div class="topic-metadata">

**Author:** [@jar1](https://discourse.julialang.org/u/jar1)\
**Replies:** 5\
**Last updated:** [June 26, 2022, 5:44pm UTC](https://discourse.julialang.org/t/generalizedlinearmixedmodel-no-method-matching-ranef-refs/82859 "2022-06-26T17:44:53Z")

</div>

How can I fix this error? form = @formula(y ~ (1 | a \* b \* c)) n = 1000 df = DataFrame( y=rand(Bool, n), a=CategoricalArray(rand('a':'z', n)), b=CategoricalArray(rand('a':'z', n)), c=CategoricalArray(…

---

## [Morris Vs Sobol Method Global Sensitivity analysis](https://discourse.julialang.org/t/morris-vs-sobol-method-global-sensitivity-analysis/83571)

<div class="topic-metadata">

**Author:** [@xztt](https://discourse.julialang.org/u/xztt)\
**Replies:** 25\
**Last updated:** [July 7, 2022, 7:46pm UTC](https://discourse.julialang.org/t/morris-vs-sobol-method-global-sensitivity-analysis/83571 "2022-07-07T19:46:23Z")

</div>

I am comparing a sensitivity analysis using both the Morris and Sobol method for the most part the methods agree. However 2 parameters appear to have low influence according to the Sobol method. But, when using the Morr…

---

## [LinearRegressionKit results correctness](https://discourse.julialang.org/t/linearregressionkit-results-correctness/83579)

<div class="topic-metadata">

**Author:** [@Eric](https://discourse.julialang.org/u/Eric)\
**Replies:** 0\
**Last updated:** [June 30, 2022, 7:52pm UTC](https://discourse.julialang.org/t/linearregressionkit-results-correctness/83579 "2022-06-30T19:52:11Z")

</div>

I haven’t posted for a while, but I noticed some discussion on the forum about results correctness, and that got me thinking that maybe I should do more on my little package. Hence I used the tests that everyone seems t…

---

## [Add errorbars with dodge in barplot](https://discourse.julialang.org/t/add-errorbars-with-dodge-in-barplot/83504)

<div class="topic-metadata">

**Author:** [@Aesth\_Estheticism](https://discourse.julialang.org/u/Aesth_Estheticism)\
**Replies:** 1\
**Last updated:** [June 29, 2022, 9:02am UTC](https://discourse.julialang.org/t/add-errorbars-with-dodge-in-barplot/83504 "2022-06-29T09:02:40Z")

</div>

How to add errorbars with dodge in barplot.

---

## [Implementing an empircal r-th q-quantile](https://discourse.julialang.org/t/implementing-an-empircal-r-th-q-quantile/83453)

<div class="topic-metadata">

**Author:** [@rikh](https://discourse.julialang.org/u/rikh)\
**Replies:** 1\
**Last updated:** [June 28, 2022, 12:28pm UTC](https://discourse.julialang.org/t/implementing-an-empircal-r-th-q-quantile/83453 "2022-06-28T12:28:55Z")

</div>

I’m trying to implement an algorithm from a paper which makes use of empirical r-th q-quantiles of the marginals X^{(1)}, ..., X^{(p)} for a dataset X with p features (Bénard et al., 2021). If I understand this correctl…

[Previous page](https://discourse.julialang.org/c/domain/stats/12.md?page=7)

[Next page](https://discourse.julialang.org/c/domain/stats/12.md?page=9)