# Statistics

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

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

**Page:** 15

---

## [Inverse prediction](https://discourse.julialang.org/t/inverse-prediction/59535)

<div class="topic-metadata">

**Author:** [@riot](https://discourse.julialang.org/u/riot)\
**Replies:** 0\
**Last updated:** [April 18, 2021, 12:19pm UTC](https://discourse.julialang.org/t/inverse-prediction/59535 "2021-04-18T12:19:36Z")

</div>

I am looking for a way to do what I used to do with the DRC-package in R. Especially the prediction of x given y (inverse prediction e.g. in probit dose-response models). In R I used the functin ED() for this. Is there a…

---

## [Adding related work to a multilevel Turing model](https://discourse.julialang.org/t/adding-related-work-to-a-multilevel-turing-model/59182)

<div class="topic-metadata">

**Author:** [@rikh](https://discourse.julialang.org/u/rikh)\
**Replies:** 4\
**Last updated:** [April 14, 2021, 12:43pm UTC](https://discourse.julialang.org/t/adding-related-work-to-a-multilevel-turing-model/59182 "2021-04-14T12:43:50Z")

</div>

Suppose that we have fitted a multilevel model over data from multiple groups called A, B, C and D and that we’re trying to estimate the population means for these 4 groups. For example, this could look something like th…

---

## [ARFIMA estimation in Julia](https://discourse.julialang.org/t/arfima-estimation-in-julia/59041)

<div class="topic-metadata">

**Author:** [@askvorts](https://discourse.julialang.org/u/askvorts)\
**Replies:** 4\
**Last updated:** [April 11, 2021, 9:52pm UTC](https://discourse.julialang.org/t/arfima-estimation-in-julia/59041 "2021-04-11T21:52:16Z")

</div>

Is there a package able to estimate an ARFIMA model in Julia (or just ARIMA for a start)? I believe that the ARFIMA.jl package is just for simulating stochastic time-series that follow an ARFIMA process. Thank you in ad…

---

## [Joinpoint Trend Analysis](https://discourse.julialang.org/t/joinpoint-trend-analysis/59035)

<div class="topic-metadata">

**Author:** [@Rafael\_Brus](https://discourse.julialang.org/u/Rafael_Brus)\
**Replies:** 1\
**Last updated:** [April 11, 2021, 4:00pm UTC](https://discourse.julialang.org/t/joinpoint-trend-analysis/59035 "2021-04-11T16:00:15Z")

</div>

Does anyone know a Pkg capable of fits the “joinpoint model” like this software: Joinpoint Trend Analysis Software Thanks

---

## [Distribution cdf calculation](https://discourse.julialang.org/t/distribution-cdf-calculation/58252)

<div class="topic-metadata">

**Author:** [@statspy](https://discourse.julialang.org/u/statspy)\
**Replies:** 2\
**Last updated:** [March 30, 2021, 6:20pm UTC](https://discourse.julialang.org/t/distribution-cdf-calculation/58252 "2021-03-30T18:20:15Z")

</div>

i have the following code that i am using and its working: using Distributions d=Exponential(1000) for i in collect(100:100:1000) print(i,"-\>",cdf(d,i) end The code it working properly. What i need now…

---

## [Any GAM or Spline packages like mgcv and splines in R?](https://discourse.julialang.org/t/any-gam-or-spline-packages-like-mgcv-and-splines-in-r/56714)

<div class="topic-metadata">

**Author:** [@NightRider](https://discourse.julialang.org/u/NightRider)\
**Replies:** 0\
**Last updated:** [March 8, 2021, 5:11am UTC](https://discourse.julialang.org/t/any-gam-or-spline-packages-like-mgcv-and-splines-in-r/56714 "2021-03-08T05:11:15Z")

</div>

So I come from a stat background and primarily use R. Julia has been amazing but I’m wondering is there a GAM package at all in the ecosystem? I have used GLM.jl so I’m wondering if there is something that is kind of lik…

---

## [Does a Terror Attack Cause Another Terror Attack? A look at Hawkes Processes in Julia](https://discourse.julialang.org/t/does-a-terror-attack-cause-another-terror-attack-a-look-at-hawkes-processes-in-julia/56578)

<div class="topic-metadata">

**Author:** [@dm13450](https://discourse.julialang.org/u/dm13450)\
**Replies:** 1\
**Last updated:** [March 5, 2021, 6:25pm UTC](https://discourse.julialang.org/t/does-a-terror-attack-cause-another-terror-attack-a-look-at-hawkes-processes-in-julia/56578 "2021-03-05T18:25:14Z")

</div>

I’ve already posted this i n the Julia subreddit but thought I would share here. This is a blog post I’ve written based on a chapter from my PhD. http://dm13450.github.io/2021/03/02/Terror-Attacks-Hawkes-Processes.html …

---

## [\[ANN\] EarlyStopping.jl 0.1.0: For terminating iterative algorithms](https://discourse.julialang.org/t/ann-earlystopping-jl-0-1-0-for-terminating-iterative-algorithms/55740)

<div class="topic-metadata">

**Author:** [@ablaom](https://discourse.julialang.org/u/ablaom)\
**Replies:** 0\
**Last updated:** [February 21, 2021, 8:51pm UTC](https://discourse.julialang.org/t/ann-earlystopping-jl-0-1-0-for-terminating-iterative-algorithms/55740 "2021-02-21T20:51:58Z")

</div>

The EarlyStopper objects defined in EarlyStopping.jl consume a sequence of numbers called losses generated by some external algorithm - generally the training loss or out-of-sample loss of some iterative statistical mode…

---

## [Creating a new distribution: 3-parameter Weibull](https://discourse.julialang.org/t/creating-a-new-distribution-3-parameter-weibull/55538)

<div class="topic-metadata">

**Author:** [@Peter\_Staab](https://discourse.julialang.org/u/Peter_Staab)\
**Replies:** 2\
**Last updated:** [February 18, 2021, 7:36pm UTC](https://discourse.julialang.org/t/creating-a-new-distribution-3-parameter-weibull/55538 "2021-02-18T19:36:14Z")

</div>

I have been working with a 3-parameter Weibull, which can be found at the Weibull wikipedia site. I’m working to get a pull request to add to the Distributions.jl library. I have two somewhat unrelated questions: shou…

---

## [Hyper Priors in Turing](https://discourse.julialang.org/t/hyper-priors-in-turing/20752)

<div class="topic-metadata">

**Author:** [@krenova](https://discourse.julialang.org/u/krenova)\
**Replies:** 5\
**Last updated:** [February 17, 2021, 9:34pm UTC](https://discourse.julialang.org/t/hyper-priors-in-turing/20752 "2021-02-17T21:34:16Z")

</div>

Hi, I’m wondering if hyperpriors work on Turing. I’m playing around with Turing and trying to use it to estimate a Hidden Markov Model. In my HMM model, I’ve used a multivariate normal prior on parameters of my transiti…

---

## [Standard error of predictions?](https://discourse.julialang.org/t/standard-error-of-predictions/55398)

<div class="topic-metadata">

**Author:** [@tomtom](https://discourse.julialang.org/u/tomtom)\
**Replies:** 2\
**Last updated:** [February 16, 2021, 4:32pm UTC](https://discourse.julialang.org/t/standard-error-of-predictions/55398 "2021-02-16T16:32:23Z")

</div>

seems like StatsBase.predict() gives only the point estimate (the mean) of prediction. How could I get more information like confidence interval or standard error of the prediction (not the coefficients)? thanks.

---

## [Is there a logspline density estimation toolbox in Julia?](https://discourse.julialang.org/t/is-there-a-logspline-density-estimation-toolbox-in-julia/55214)

<div class="topic-metadata">

**Author:** [@drbenvincent](https://discourse.julialang.org/u/drbenvincent)\
**Replies:** 1\
**Last updated:** [February 15, 2021, 4:25pm UTC](https://discourse.julialang.org/t/is-there-a-logspline-density-estimation-toolbox-in-julia/55214 "2021-02-15T16:25:48Z")

</div>

For density estimation, I know that there are various kernel density estimation packages (KernelDensity.jl, KernelDensityEstimate.jl, KernelEstimator.jl), but I wonder if there are any spline-based packages? The goal he…

---

## [Correlated Random Effects in MixedModels.jl](https://discourse.julialang.org/t/correlated-random-effects-in-mixedmodels-jl/54612)

<div class="topic-metadata">

**Author:** [@danielw2904](https://discourse.julialang.org/u/danielw2904)\
**Replies:** 1\
**Last updated:** [February 10, 2021, 8:40am UTC](https://discourse.julialang.org/t/correlated-random-effects-in-mixedmodels-jl/54612 "2021-02-10T08:40:57Z")

</div>

I recently read up on Fixed Effects vs. Random Effects and learned about the Correlated Random Effects model in Wooldridge’s Introductory Econometrics. The model looks as follows: y\_{it} = \\beta x\_{it} + \\gamma \\bar x…

---

## [Proper use of scaling on axes from ordinations (like MDS and PCA)?](https://discourse.julialang.org/t/proper-use-of-scaling-on-axes-from-ordinations-like-mds-and-pca/24606)

<div class="topic-metadata">

**Author:** [@kevbonham](https://discourse.julialang.org/u/kevbonham)\
**Replies:** 16\
**Last updated:** [February 8, 2021, 3:51pm UTC](https://discourse.julialang.org/t/proper-use-of-scaling-on-axes-from-ordinations-like-mds-and-pca/24606 "2021-02-08T15:51:27Z")

</div>

Moving a conversation from slack over to discourse because there’s a bunch of useful stuff that I’d like preserved, but also because we haven’t actually reached a consensus answer. The original question (by me): Quest…

---

## [Seeking advice on contributing to Julia: A faster Lasso solver](https://discourse.julialang.org/t/seeking-advice-on-contributing-to-julia-a-faster-lasso-solver/42023)

<div class="topic-metadata">

**Author:** [@barankarakus](https://discourse.julialang.org/u/barankarakus)\
**Replies:** 9\
**Last updated:** [February 4, 2021, 7:02am UTC](https://discourse.julialang.org/t/seeking-advice-on-contributing-to-julia-a-faster-lasso-solver/42023 "2021-02-04T07:02:31Z")

</div>

Hi there, I’m a Master’s student in Statistics. I’ve recently been researching ‘safe screening rules’ for accelerating Lasso solvers. These can improve computation times by orders of magnitude when the number of feature…

---

## [Likelihood Ratio Test?](https://discourse.julialang.org/t/likelihood-ratio-test/38217)

<div class="topic-metadata">

**Author:** [@bjoerges](https://discourse.julialang.org/u/bjoerges)\
**Replies:** 10\
**Last updated:** [February 1, 2021, 3:23pm UTC](https://discourse.julialang.org/t/likelihood-ratio-test/38217 "2021-02-01T15:23:57Z")

</div>

Hey everyone, I am looking for a way to do Likelihood Ratio Tests in Julia. I have found a function in MixedModels.jl (which is also what I am using to fit my models), but it doesn’t seem to work for Generalized Linear …

---

## [How to subset dataframe by multiple rows with same value?](https://discourse.julialang.org/t/how-to-subset-dataframe-by-multiple-rows-with-same-value/54294)

<div class="topic-metadata">

**Author:** [@Shazman](https://discourse.julialang.org/u/Shazman)\
**Replies:** 6\
**Last updated:** [February 1, 2021, 11:48am UTC](https://discourse.julialang.org/t/how-to-subset-dataframe-by-multiple-rows-with-same-value/54294 "2021-02-01T11:48:01Z")

</div>

Hi I want to subset a dataframe(df) by multiple rows with value (10004). I tried the command, but it does not work, I am wondering why ? df\[in(\[10004\]).(df.HWG), :\] 28,960 rows × 11 columns (omitted printing of 2 co…

---

## [Auxiliary GLM Statistics](https://discourse.julialang.org/t/auxiliary-glm-statistics/14722)

<div class="topic-metadata">

**Author:** [@iwelch](https://discourse.julialang.org/u/iwelch)\
**Replies:** 9\
**Last updated:** [January 27, 2021, 6:35pm UTC](https://discourse.julialang.org/t/auxiliary-glm-statistics/14722 "2021-01-27T18:35:30Z")

</div>

Julia 1.0. GLM 1.0.0. My Basic Example julia\> using GLM, StatsBase julia\> y=\[1:10;\] ; X= hcat( fill(1,10), y.^2, y.^3 ); lm1= GLM.lm(X, y); Basic Questions is the sigma (i.e., sqrt(deviance(lm)/dof\_residual(lm)) al…

---

## [Generate a zero mean gaussian vector](https://discourse.julialang.org/t/generate-a-zero-mean-gaussian-vector/54056)

<div class="topic-metadata">

**Author:** [@cheng\_chen](https://discourse.julialang.org/u/cheng_chen)\
**Replies:** 1\
**Last updated:** [January 27, 2021, 2:12pm UTC](https://discourse.julialang.org/t/generate-a-zero-mean-gaussian-vector/54056 "2021-01-27T14:12:37Z")

</div>

How can I generate a zero mean 2 by 1 Gaussian vector with covariance Identity matrix?

---

## [How to make this AFT regression stable?](https://discourse.julialang.org/t/how-to-make-this-aft-regression-stable/53685)

<div class="topic-metadata">

**Author:** [@mkarikom](https://discourse.julialang.org/u/mkarikom)\
**Replies:** 5\
**Last updated:** [January 22, 2021, 1:01am UTC](https://discourse.julialang.org/t/how-to-make-this-aft-regression-stable/53685 "2021-01-22T01:01:00Z")

</div>

I have adapted this excellent post by @aaowens to fit a right-censored AFT regression model, but for some reason the likelihood is not compatible with forwarddiff. using Distributions,Optim function AFT(pvec,data) f…

---

## [Can anyone tell me what is this distribution?](https://discourse.julialang.org/t/can-anyone-tell-me-what-is-this-distribution/53726)

<div class="topic-metadata">

**Author:** [@abulak](https://discourse.julialang.org/u/abulak)\
**Replies:** 0\
**Last updated:** [January 21, 2021, 2:14pm UTC](https://discourse.julialang.org/t/can-anyone-tell-me-what-is-this-distribution/53726 "2021-01-21T14:14:15Z")

</div>

In my hobby project I arrived at the following problem: I’m given a family of random variables: Wₖ for which I need to estimate cdf: (the above is histogram approximation of pdf from 10^7 samples) From theoretical con…

---

## [Need advice when data almost (but not quite) fit log normal dist](https://discourse.julialang.org/t/need-advice-when-data-almost-but-not-quite-fit-log-normal-dist/53287)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 7\
**Last updated:** [January 13, 2021, 6:38pm UTC](https://discourse.julialang.org/t/need-advice-when-data-almost-but-not-quite-fit-log-normal-dist/53287 "2021-01-13T18:38:37Z")

</div>

This isn’t really a question about Julia so maybe this should be tagged as Offtopic, but I’m working with a data set that has several series of data that all closely follow a log normal distribution, but not quite. The p…

---

## [What shape should data be in to perform a linear fit with GLM?](https://discourse.julialang.org/t/what-shape-should-data-be-in-to-perform-a-linear-fit-with-glm/53120)

<div class="topic-metadata">

**Author:** [@alice](https://discourse.julialang.org/u/alice)\
**Replies:** 7\
**Last updated:** [January 11, 2021, 9:22am UTC](https://discourse.julialang.org/t/what-shape-should-data-be-in-to-perform-a-linear-fit-with-glm/53120 "2021-01-11T09:22:15Z")

</div>

I am trying to fit some data with a variable number of dimensions using a linear fit. I am currently trying to use GLM, but I am happy to use any other library. I’ve been running in to a lot of trouble getting the data …

---

## [How to fit a normal approximation to data in Julia](https://discourse.julialang.org/t/how-to-fit-a-normal-approximation-to-data-in-julia/35045)

<div class="topic-metadata">

**Author:** [@lesshaste](https://discourse.julialang.org/u/lesshaste)\
**Replies:** 21\
**Last updated:** [January 8, 2021, 7:32pm UTC](https://discourse.julialang.org/t/how-to-fit-a-normal-approximation-to-data-in-julia/35045 "2021-01-08T19:32:19Z")

</div>

I have an unnormalized array of counts for a probability mass function (derived from sampling from the true function). The array should be indexed from 0. 12 is the number of times a 0 occurred in the sample. 199 is th…

---

## [Why is statistics so difficult?](https://discourse.julialang.org/t/why-is-statistics-so-difficult/52416)

<div class="topic-metadata">

**Author:** [@rikh](https://discourse.julialang.org/u/rikh)\
**Replies:** 40\
**Last updated:** [January 6, 2021, 1:33pm UTC](https://discourse.julialang.org/t/why-is-statistics-so-difficult/52416 "2021-01-06T13:33:50Z")

</div>

Hi all! I see a few very interesting discussions going on on this forum, and would like to ask for your opinions on the difficulty of statistics. I see many people hanging around here with much experience with stats, and…

---

## [Meaning of Δ in different \`apply!\` methods is different in different Flux optimizers?](https://discourse.julialang.org/t/meaning-of-in-different-apply-methods-is-different-in-different-flux-optimizers/52811)

<div class="topic-metadata">

**Author:** [@gxyd](https://discourse.julialang.org/u/gxyd)\
**Replies:** 0\
**Last updated:** [January 4, 2021, 7:57am UTC](https://discourse.julialang.org/t/meaning-of-in-different-apply-methods-is-different-in-different-flux-optimizers/52811 "2021-01-04T07:57:08Z")

</div>

Hi, I want to understand the meaning of Δ used in different methods (i.e. for different optimizers) in Flux. For example: in Descent, the meaning of Δ seems to be: \[∂f/∂θᵢ\] i=1,2,3,… While looking at the meaning of Δ …

---

## [Compute integral of a beta distribution](https://discourse.julialang.org/t/compute-integral-of-a-beta-distribution/52607)

<div class="topic-metadata">

**Author:** [@TBu](https://discourse.julialang.org/u/TBu)\
**Replies:** 7\
**Last updated:** [December 30, 2020, 12:40pm UTC](https://discourse.julialang.org/t/compute-integral-of-a-beta-distribution/52607 "2020-12-30T12:40:54Z")

</div>

Hi everyone, I used the Distributions package in order to generate a beta distribution, which works quite nicely: using Distributions d = Beta(10,15) However - is it possible to integrate over an interval of the distr…

---

## [Density Distributions](https://discourse.julialang.org/t/density-distributions/52224)

<div class="topic-metadata">

**Author:** [@lorrp1](https://discourse.julialang.org/u/lorrp1)\
**Replies:** 6\
**Last updated:** [December 22, 2020, 12:28pm UTC](https://discourse.julialang.org/t/density-distributions/52224 "2020-12-22T12:28:25Z")

</div>

Hello, im trying to get the array of the discrete probability distribution of a normal distribution: i have tried: d = Normal() #standard a = rand(d,1000) x = pdf.(d, a) #even with just pdf(d,a) plot(x) but its not wo…

---

## [Package for time series?](https://discourse.julialang.org/t/package-for-time-series/29555)

<div class="topic-metadata">

**Author:** [@Donkeyfish87](https://discourse.julialang.org/u/Donkeyfish87)\
**Replies:** 10\
**Last updated:** [December 14, 2020, 5:59pm UTC](https://discourse.julialang.org/t/package-for-time-series/29555 "2020-12-14T17:59:42Z")

</div>

Are there any Julia packages for time series analysis? Like arma, arima, vector autoregressions etc. I’ve tried looking around a bit, but didn’t find anything. Thanks.

---

## [High Energy Physics (HEP) folks on Julia?](https://discourse.julialang.org/t/high-energy-physics-hep-folks-on-julia/15330)

<div class="topic-metadata">

**Author:** [@jling](https://discourse.julialang.org/u/jling)\
**Replies:** 100\
**Last updated:** [December 12, 2020, 1:07pm UTC](https://discourse.julialang.org/t/high-energy-physics-hep-folks-on-julia/15330 "2020-12-12T13:07:32Z")

</div>

Maybe it’s not the most suitable domain but imo close enough. For reference, ROOT refers to ROOT/CERN Really enjoyed using Julia for a while and trying to apply it to more study/research replated works. I noticed the…

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

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