# Statistics

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

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

**Page:** 20

---

## [How do you guys would plot and histogram this Distribution F code?](https://discourse.julialang.org/t/how-do-you-guys-would-plot-and-histogram-this-distribution-f-code/28675)

<div class="topic-metadata">

**Author:** [@Halor](https://discourse.julialang.org/u/Halor)\
**Replies:** 3\
**Last updated:** [September 12, 2019, 7:29pm UTC](https://discourse.julialang.org/t/how-do-you-guys-would-plot-and-histogram-this-distribution-f-code/28675 "2019-09-12T19:29:28Z")

</div>

Hi. I am a new user in programming and in statistics so I want to learn all the things that I can. I try to code a PDF of F Distribution which formula have taken from Wikipedia (Distribución F - Wikipedia, la encicloped…

---

## [Conditional logit regression -solved using cox regression](https://discourse.julialang.org/t/conditional-logit-regression-solved-using-cox-regression/28615)

<div class="topic-metadata">

**Author:** [@compleat](https://discourse.julialang.org/u/compleat)\
**Replies:** 11\
**Last updated:** [September 12, 2019, 11:52am UTC](https://discourse.julialang.org/t/conditional-logit-regression-solved-using-cox-regression/28615 "2019-09-12T11:52:47Z")

</div>

I noticed on here a thread which can’t now find about conditional/multinomial logit. Actually, there is no need to implement anything new - this is trivial if cox proportional hazards regression exists and allows for l…

---

## [F Distribution, Beta function](https://discourse.julialang.org/t/f-distribution-beta-function/28630)

<div class="topic-metadata">

**Author:** [@Halor](https://discourse.julialang.org/u/Halor)\
**Replies:** 7\
**Last updated:** [September 12, 2019, 2:26am UTC](https://discourse.julialang.org/t/f-distribution-beta-function/28630 "2019-09-12T02:26:43Z")

</div>

Hi. I am new in this programming world so I need help for coding a F Distribution, due to I do not have idea to do this. I am having trouble even programming beta function. I just need some ideas so I would be able to d…

---

## [Anything like Distributions.jl for GPU?](https://discourse.julialang.org/t/anything-like-distributions-jl-for-gpu/17147)

<div class="topic-metadata">

**Author:** [@dfdx](https://discourse.julialang.org/u/dfdx)\
**Replies:** 15\
**Last updated:** [September 8, 2019, 9:29pm UTC](https://discourse.julialang.org/t/anything-like-distributions-jl-for-gpu/17147 "2019-09-08T21:29:42Z")

</div>

Is there any package supporting multivariate distributions and basic operations on them like sampling and log pdf, but working on GPU?

---

## [How do I fit a distribution to these data sets?!](https://discourse.julialang.org/t/how-do-i-fit-a-distribution-to-these-data-sets/28243)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 5\
**Last updated:** [September 1, 2019, 5:59am UTC](https://discourse.julialang.org/t/how-do-i-fit-a-distribution-to-these-data-sets/28243 "2019-09-01T05:59:22Z")

</div>

Let me start by saying that I’m very inexperienced in this area, so go easy on me @Tamas\_Papp :wink:. I have two datasets that are similar in shape and I’m trying to fit a distribution to each of them with the goal of be…

---

## [Mixedmodels.jl branches: nlmm and dropmulαβ!](https://discourse.julialang.org/t/mixedmodels-jl-branches-nlmm-and-dropmul/28223)

<div class="topic-metadata">

**Author:** [@Juan](https://discourse.julialang.org/u/Juan)\
**Replies:** 2\
**Last updated:** [August 31, 2019, 11:48am UTC](https://discourse.julialang.org/t/mixedmodels-jl-branches-nlmm-and-dropmul/28223 "2019-08-31T11:48:43Z")

</div>

I’ve seen the Mixedmodels.jl package has thre differente branches: master nlmm dropmulαβ! How are they different? What’s new there?

---

## [GLM linear regression. How to extract the coefficients?](https://discourse.julialang.org/t/glm-linear-regression-how-to-extract-the-coefficients/28219)

<div class="topic-metadata">

**Author:** [@doomphoenix-qxz](https://discourse.julialang.org/u/doomphoenix-qxz)\
**Replies:** 3\
**Last updated:** [August 30, 2019, 9:56pm UTC](https://discourse.julialang.org/t/glm-linear-regression-how-to-extract-the-coefficients/28219 "2019-08-30T21:56:45Z")

</div>

I fit a simple linear regression using GLM. How do I extract the calculated slope and intercept? I looked through the documentation but didn’t see an example of it. If such an example exists, feel free to point there.

---

## [Sample size calculation + experimental design](https://discourse.julialang.org/t/sample-size-calculation-experimental-design/27433)

<div class="topic-metadata">

**Author:** [@Crghilardi](https://discourse.julialang.org/u/Crghilardi)\
**Replies:** 1\
**Last updated:** [August 28, 2019, 11:35pm UTC](https://discourse.julialang.org/t/sample-size-calculation-experimental-design/27433 "2019-08-28T23:35:32Z")

</div>

I started to write a small function that I use somewhat regularly to calculate sample size from coefficient of variation. I checked JuliaStats/StatsBase to see if any sample size functionality was included and I did not …

---

## [Creating graphs from data](https://discourse.julialang.org/t/creating-graphs-from-data/27199)

<div class="topic-metadata">

**Author:** [@iamsuddhasattwa](https://discourse.julialang.org/u/iamsuddhasattwa)\
**Replies:** 7\
**Last updated:** [August 7, 2019, 7:50pm UTC](https://discourse.julialang.org/t/creating-graphs-from-data/27199 "2019-08-07T19:50:01Z")

</div>

A common procedure in data science and mathematics is building graphs/networks from a large collection of data points. I am interested in the constructing the adjacency matrix for the k-nearest neighbor graph of a collec…

---

## [Balanced Repeated Replication in Julia](https://discourse.julialang.org/t/balanced-repeated-replication-in-julia/27003)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 5\
**Last updated:** [August 2, 2019, 4:55pm UTC](https://discourse.julialang.org/t/balanced-repeated-replication-in-julia/27003 "2019-08-02T16:55:33Z")

</div>

Hi All, I am working with a public data source and I’m having trouble understanding how to compute replicate weights that are needed for calculating the standard errors for my estimates. The data source is here, if you’…

---

## [Network structural balance analysis in Julia](https://discourse.julialang.org/t/network-structural-balance-analysis-in-julia/27045)

<div class="topic-metadata">

**Author:** [@kcyb](https://discourse.julialang.org/u/kcyb)\
**Replies:** 0\
**Last updated:** [July 31, 2019, 5:19pm UTC](https://discourse.julialang.org/t/network-structural-balance-analysis-in-julia/27045 "2019-07-31T17:19:31Z")

</div>

Dear Julia community, I have a question regarding network analysis in Julia. I apologize if this is the wrong section, I couldn’t decide for a better one. Are you aware of any packages for sound network structural balan…

---

## [Mean function unknown](https://discourse.julialang.org/t/mean-function-unknown/26844)

<div class="topic-metadata">

**Author:** [@Diego](https://discourse.julialang.org/u/Diego)\
**Replies:** 4\
**Last updated:** [July 30, 2019, 4:57pm UTC](https://discourse.julialang.org/t/mean-function-unknown/26844 "2019-07-30T16:57:15Z")

</div>

After adding and importing the Statistics package Julia still does not recognize the mean function, any ideas?

---

## [Distributions package invlogcdf() function returning the wrong answer?](https://discourse.julialang.org/t/distributions-package-invlogcdf-function-returning-the-wrong-answer/26867)

<div class="topic-metadata">

**Author:** [@dataSurfer](https://discourse.julialang.org/u/dataSurfer)\
**Replies:** 2\
**Last updated:** [July 27, 2019, 12:39am UTC](https://discourse.julialang.org/t/distributions-package-invlogcdf-function-returning-the-wrong-answer/26867 "2019-07-27T00:39:21Z")

</div>

Consider that the following answer: invlogcdf(Chi(11), log(0.975)) # 4.681885225101231 qchisq(0.975, 11) # 21.920049261021205 is wrong when comparing invlogcdf() to the equivalent R function, but the following is corre…

---

## [DynamicHMC API rewrite](https://discourse.julialang.org/t/dynamichmc-api-rewrite/18980)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 1\
**Last updated:** [July 25, 2019, 9:02am UTC](https://discourse.julialang.org/t/dynamichmc-api-rewrite/18980 "2019-07-25T09:02:01Z")

</div>

I will be rewriting the API of DynamicHMC.jl shortly. I started a discussion at https://github.com/tpapp/DynamicHMC.jl/issues/30 and invite users to participate there, or in this topic. In particular, if there is a fea…

---

## [Using RCall, R 3.5.3 -\> R 3.6.1](https://discourse.julialang.org/t/using-rcall-r-3-5-3-r-3-6-1/26675)

<div class="topic-metadata">

**Author:** [@Gunter\_Faes](https://discourse.julialang.org/u/Gunter_Faes)\
**Replies:** 1\
**Last updated:** [July 23, 2019, 10:08am UTC](https://discourse.julialang.org/t/using-rcall-r-3-5-3-r-3-6-1/26675 "2019-07-23T10:08:12Z")

</div>

First of all, this is just a field report! With R version 3.6.0, some major changes and enhancements have been implemented compared to the previous versions (e.g. R 3.5.3). This means that it became necessary to instal…

---

## [Is the following code using some regression technique?](https://discourse.julialang.org/t/is-the-following-code-using-some-regression-technique/26636)

<div class="topic-metadata">

**Author:** [@bsnyh](https://discourse.julialang.org/u/bsnyh)\
**Replies:** 1\
**Last updated:** [July 22, 2019, 3:57pm UTC](https://discourse.julialang.org/t/is-the-following-code-using-some-regression-technique/26636 "2019-07-22T15:57:32Z")

</div>

The following piece of code is from the OpenStreetMapX package. The logic is pretty straightforward to follow. I was just wondering, where is the constant 3.6 from? Does it come from some statistical fittings? Any commen…

---

## [Changing the binding of a variable in a closure](https://discourse.julialang.org/t/changing-the-binding-of-a-variable-in-a-closure/26297)

<div class="topic-metadata">

**Author:** [@dmbates](https://discourse.julialang.org/u/dmbates)\
**Replies:** 5\
**Last updated:** [July 13, 2019, 6:22pm UTC](https://discourse.julialang.org/t/changing-the-binding-of-a-variable-in-a-closure/26297 "2019-07-13T18:22:32Z")

</div>

I am labeling this question as “statistics” because the motivation comes from a statistical application. The tl: dr version is that I want to change the binding of a variable in a closure within a loop. The function it…

---

## [Google Summer of Code Proposal: Econometrics.jl](https://discourse.julialang.org/t/google-summer-of-code-proposal-econometrics-jl/2811)

<div class="topic-metadata">

**Author:** [@Nosferican](https://discourse.julialang.org/u/Nosferican)\
**Replies:** 20\
**Last updated:** [July 6, 2019, 8:29pm UTC](https://discourse.julialang.org/t/google-summer-of-code-proposal-econometrics-jl/2811 "2019-07-06T20:29:38Z")

</div>

Greetings y’all. I am applying for the Google Summer of Code with Julia. My proposal draft is here: Proposal. Any comments or suggestions are welcomed. Thanks! Bayoán PS: Latest Draft Version can be accessed here.

---

## [Odd result from Distributions quantile](https://discourse.julialang.org/t/odd-result-from-distributions-quantile/25820)

<div class="topic-metadata">

**Author:** [@abx](https://discourse.julialang.org/u/abx)\
**Replies:** 3\
**Last updated:** [July 4, 2019, 6:56pm UTC](https://discourse.julialang.org/t/odd-result-from-distributions-quantile/25820 "2019-07-04T18:56:34Z")

</div>

I am experiencing inconsistent results using the Distributions package. When I run the following code, I get two different results on different machines Distributions.quantile(Distributions.Normal(), 0.3631592166130919…

---

## [Gaussian Mixture Model estimation in Julia 1.0](https://discourse.julialang.org/t/gaussian-mixture-model-estimation-in-julia-1-0/17339)

<div class="topic-metadata">

**Author:** [@tanhevg](https://discourse.julialang.org/u/tanhevg)\
**Replies:** 8\
**Last updated:** [July 3, 2019, 12:16pm UTC](https://discourse.julialang.org/t/gaussian-mixture-model-estimation-in-julia-1-0/17339 "2019-07-03T12:16:53Z")

</div>

Hello, I wonder, if there are any up-to-date Julia packages for GMM estimation? I googled around, and came across https://github.com/davidavdav/GaussianMixtures.jl and https://github.com/lindahua/MixtureModels.jl. Both …

---

## [PSA: breaking changes in StatsModels v0.6.0 (Terms 2.0: Son of Terms)](https://discourse.julialang.org/t/psa-breaking-changes-in-statsmodels-v0-6-0-terms-2-0-son-of-terms/25626)

<div class="topic-metadata">

**Author:** [@dave.f.kleinschmidt](https://discourse.julialang.org/u/dave.f.kleinschmidt)\
**Replies:** 6\
**Last updated:** [July 2, 2019, 2:53pm UTC](https://discourse.julialang.org/t/psa-breaking-changes-in-statsmodels-v0-6-0-terms-2-0-son-of-terms/25626 "2019-07-02T14:53:28Z")

</div>

It’s finally time to release version 0.6.0 of StatsModels.jl, which includes the radically re-designed internals introduced in Terms 2.0: Son of Terms (announcement here). Most of these changes should not affect end use…

---

## [Turing - Syntax for Compositional Gibbs](https://discourse.julialang.org/t/turing-syntax-for-compositional-gibbs/25767)

<div class="topic-metadata">

**Author:** [@mrVeng](https://discourse.julialang.org/u/mrVeng)\
**Replies:** 4\
**Last updated:** [June 27, 2019, 4:47pm UTC](https://discourse.julialang.org/t/turing-syntax-for-compositional-gibbs/25767 "2019-06-27T16:47:58Z")

</div>

Hi everyone, I just finished reading the paper for the Julia Turing module and went through your tutorials – looks really great, congratulations! I did not find a specific forum for it, so I post my question here: In y…

---

## [MixedModels.jl confidence intervals and repeated measure, how to get?](https://discourse.julialang.org/t/mixedmodels-jl-confidence-intervals-and-repeated-measure-how-to-get/25534)

<div class="topic-metadata">

**Author:** [@PharmCat](https://discourse.julialang.org/u/PharmCat)\
**Replies:** 9\
**Last updated:** [June 26, 2019, 3:38pm UTC](https://discourse.julialang.org/t/mixedmodels-jl-confidence-intervals-and-repeated-measure-how-to-get/25534 "2019-06-26T15:38:17Z")

</div>

I have a data: julia\> df 298×6 DataFrame │ Row │ Subject │ Period │ Seq │ Formulation │ logData │ Data │ │ │ Categorical… │ Categorical… │ Categorical… │ Categorical… │ Float64 │ Float64 │ ├─…

---

## [Estimate most probable value from data](https://discourse.julialang.org/t/estimate-most-probable-value-from-data/25451)

<div class="topic-metadata">

**Author:** [@favba](https://discourse.julialang.org/u/favba)\
**Replies:** 7\
**Last updated:** [June 19, 2019, 9:33pm UTC](https://discourse.julialang.org/t/estimate-most-probable-value-from-data/25451 "2019-06-19T21:33:05Z")

</div>

I have discrete data from a continuous field. The data probability distribution is not known, but it has only one peak value, what I guess is called “mode” of the pdf. The normalized histograms from the data look like s…

---

## [Normalization of histograms](https://discourse.julialang.org/t/normalization-of-histograms/25007)

<div class="topic-metadata">

**Author:** [@woutersj](https://discourse.julialang.org/u/woutersj)\
**Replies:** 1\
**Last updated:** [June 6, 2019, 1:22pm UTC](https://discourse.julialang.org/t/normalization-of-histograms/25007 "2019-06-06T13:22:40Z")

</div>

This used to work in v0.6.3: using StatsBase normalize(fit(Histogram,rand(100),closed=:left)) In v1.1.0 I get UndefVarError: normalize not defined. Is this a bug or has the syntax changed?

---

## [Gaussian process in Julia](https://discourse.julialang.org/t/gaussian-process-in-julia/22900)

<div class="topic-metadata">

**Author:** [@yusri-dh](https://discourse.julialang.org/u/yusri-dh)\
**Replies:** 13\
**Last updated:** [May 29, 2019, 1:33pm UTC](https://discourse.julialang.org/t/gaussian-process-in-julia/22900 "2019-05-29T13:33:29Z")

</div>

Hello all, Julia newbie here. I am sorry I am not sure whether this is the right place to ask this question. If no I will delete it later. Is there any well-maintained package dedicated for gaussian process in Julia lik…

---

## [Find gradient of a posterior](https://discourse.julialang.org/t/find-gradient-of-a-posterior/24254)

<div class="topic-metadata">

**Author:** [@Uche](https://discourse.julialang.org/u/Uche)\
**Replies:** 12\
**Last updated:** [May 16, 2019, 5:25pm UTC](https://discourse.julialang.org/t/find-gradient-of-a-posterior/24254 "2019-05-16T17:25:45Z")

</div>

Please can anyone help me with this. I have been trying to use ForwardDiff or Flux.Tracker or Zygote to get the gradient of this posterior but has not had any success yet: using AltDistributions, Distributions using L…

---

## [Bayesian Neural Network (Multiclass Classification) in Turing.jl and Flux.jl](https://discourse.julialang.org/t/bayesian-neural-network-multiclass-classification-in-turing-jl-and-flux-jl/21134)

<div class="topic-metadata">

**Author:** [@alasaadstat](https://discourse.julialang.org/u/alasaadstat)\
**Replies:** 5\
**Last updated:** [May 9, 2019, 3:08pm UTC](https://discourse.julialang.org/t/bayesian-neural-network-multiclass-classification-in-turing-jl-and-flux-jl/21134 "2019-05-09T15:08:01Z")

</div>

I want to do a multiclass classification using Bayesian Neural Network (BNN) in Turing.jl and Flux.jl. There’s a good implementation already of the binary classification using BNN in Turing.jl, check it here. Hence, my g…

---

## [Rand() doesn't work for NormalInverseGaussian](https://discourse.julialang.org/t/rand-doesnt-work-for-normalinversegaussian/23948)

<div class="topic-metadata">

**Author:** [@BMval](https://discourse.julialang.org/u/BMval)\
**Replies:** 2\
**Last updated:** [May 7, 2019, 5:12pm UTC](https://discourse.julialang.org/t/rand-doesnt-work-for-normalinversegaussian/23948 "2019-05-07T17:12:33Z")

</div>

rand() doesn’t work for NormalInverseGaussian. x = rand(dnig, 100) Error: MethodError: no method matching iterate(::NormalInverseGaussian{Float64}) Would somebody help to generate random values?

---

## [SpecialFunctions ( Julia v1.0 )](https://discourse.julialang.org/t/specialfunctions-julia-v1-0/23920)

<div class="topic-metadata">

**Author:** [@BMval](https://discourse.julialang.org/u/BMval)\
**Replies:** 6\
**Last updated:** [May 6, 2019, 6:55pm UTC](https://discourse.julialang.org/t/specialfunctions-julia-v1-0/23920 "2019-05-06T18:55:11Z")

</div>

Good day. I need to use besselk function from SpecialFunctions package. Unfortunately, SpecialFunctions don’t work under v1.0 Julia version. At the same time Distributions.jl uses besselk function and work under v1.0 …

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

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