# Statistics

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

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

**Page:** 25

---

## [LoadError: MethodError: no method matching quantile(::Distributions.JohnsonSU{Float64}, ::Float64)](https://discourse.julialang.org/t/loaderror-methoderror-no-method-matching-quantile-distributions-johnsonsu-float64-float64/3371)

<div class="topic-metadata">

**Author:** [@Wael\_Sboui](https://discourse.julialang.org/u/Wael_Sboui)\
**Replies:** 2\
**Last updated:** [April 26, 2017, 8:03pm UTC](https://discourse.julialang.org/t/loaderror-methoderror-no-method-matching-quantile-distributions-johnsonsu-float64-float64/3371 "2017-04-26T20:03:06Z")

</div>

Hi smart community, I faced this problem after having added a univariate continuous to Distributions.jl and univariates.jl Always working on Julia 0.5.1. However, I didn’t face any similar problem when I was on Julia 0.…

---

## [LoadError: error in method definition: function Base.min must be explicitly imported to be extended](https://discourse.julialang.org/t/loaderror-error-in-method-definition-function-base-min-must-be-explicitly-imported-to-be-extended/3362)

<div class="topic-metadata">

**Author:** [@Wael\_Sboui](https://discourse.julialang.org/u/Wael_Sboui)\
**Replies:** 3\
**Last updated:** [April 25, 2017, 2:03pm UTC](https://discourse.julialang.org/t/loaderror-error-in-method-definition-function-base-min-must-be-explicitly-imported-to-be-extended/3362 "2017-04-25T14:03:55Z")

</div>

Hi everybody, I’m working under Julia 0.5.1 and I face this problem below: Error evaluating AirplaneRealStats5.9.jl LoadError: LoadError: error in method definition: function Base.min must be explicitly imported to be …

---

## [Binomial distribution without Distributions.jl?](https://discourse.julialang.org/t/binomial-distribution-without-distributions-jl/3176)

<div class="topic-metadata">

**Author:** [@anon94023334](https://discourse.julialang.org/u/anon94023334)\
**Replies:** 11\
**Last updated:** [April 11, 2017, 11:57pm UTC](https://discourse.julialang.org/t/binomial-distribution-without-distributions-jl/3176 "2017-04-11T23:57:19Z")

</div>

Hi all, This is not intended to be a knock on the excellent Distributions.jl package - we’ve used it happily for several years in LightGraphs and it’s working wonderfully. We currently have an initiative to make LightG…

---

## [Projection Pursuit Regression](https://discourse.julialang.org/t/projection-pursuit-regression/3147)

<div class="topic-metadata">

**Author:** [@babaq](https://discourse.julialang.org/u/babaq)\
**Replies:** 1\
**Last updated:** [April 11, 2017, 7:49am UTC](https://discourse.julialang.org/t/projection-pursuit-regression/3147 "2017-04-11T07:49:34Z")

</div>

Hi, I am considering using Projection Pursuit Regression in a project to analyze data. After some search in the internet, it seems only R provides ppr method in stats package. I am wandering is there some Julia package …

---

## [(approximate) PDFs/CDFs suitable for automatic differentiation?](https://discourse.julialang.org/t/approximate-pdfs-cdfs-suitable-for-automatic-differentiation/3027)

<div class="topic-metadata">

**Author:** [@evanfields](https://discourse.julialang.org/u/evanfields)\
**Replies:** 6\
**Last updated:** [April 3, 2017, 9:56pm UTC](https://discourse.julialang.org/t/approximate-pdfs-cdfs-suitable-for-automatic-differentiation/3027 "2017-04-03T21:56:42Z")

</div>

I’m working on a project where, given some observations of transformed random variables, we would like to estimate the parameters of the distribution those observations came from. As a simple example, suppose Y = min(X\_…

---

## [Get columns of DataFrames.jl ModelMatrix object](https://discourse.julialang.org/t/get-columns-of-dataframes-jl-modelmatrix-object/2780)

<div class="topic-metadata">

**Author:** [@sbohor3](https://discourse.julialang.org/u/sbohor3)\
**Replies:** 12\
**Last updated:** [March 23, 2017, 5:36pm UTC](https://discourse.julialang.org/t/get-columns-of-dataframes-jl-modelmatrix-object/2780 "2017-03-23T17:36:32Z")

</div>

Is there a way to extract column names of an object created from ModelMatrix() so that I can use those columns as row names for coefficients? Thanks.

---

## [Equivalent Julia package to R package \`prob\`](https://discourse.julialang.org/t/equivalent-julia-package-to-r-package-prob/2829)

<div class="topic-metadata">

**Author:** [@AjayTalati](https://discourse.julialang.org/u/AjayTalati)\
**Replies:** 9\
**Last updated:** [March 23, 2017, 2:56pm UTC](https://discourse.julialang.org/t/equivalent-julia-package-to-r-package-prob/2829 "2017-03-23T14:56:13Z")

</div>

Hello :slight\_smile:, I’m new to Julia, and have been using R previously. Could I ask if there’s an equivalent package in Julia to the package, prob in R. https://cran.r-project.org/web/packages/prob/vignettes/prob.pdf …

---

## [Is there a dropna() for DataFrame?](https://discourse.julialang.org/t/is-there-a-dropna-for-dataframe/1777)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 23\
**Last updated:** [March 22, 2017, 2:11pm UTC](https://discourse.julialang.org/t/is-there-a-dropna-for-dataframe/1777 "2017-03-22T14:11:47Z")

</div>

How to drop the rows of a DataFrame that contain NA in any column? dropna() apparently only works with DataArray?

---

## [Evaluating explicitly the kernel for a Gaussian/Normal distribution](https://discourse.julialang.org/t/evaluating-explicitly-the-kernel-for-a-gaussian-normal-distribution/2632)

<div class="topic-metadata">

**Author:** [@Bernal](https://discourse.julialang.org/u/Bernal)\
**Replies:** 1\
**Last updated:** [March 13, 2017, 9:02am UTC](https://discourse.julialang.org/t/evaluating-explicitly-the-kernel-for-a-gaussian-normal-distribution/2632 "2017-03-13T09:02:46Z")

</div>

Hi folks, probably I’m attempting to discover “how to boil water” however I have faced some weird problems building a vector of values matching a Gaussian bell (I do NOT want to create a random variables vector that is…

---

## [How does Loess decide how many sub-intervals to use for local regression?](https://discourse.julialang.org/t/how-does-loess-decide-how-many-sub-intervals-to-use-for-local-regression/2615)

<div class="topic-metadata">

**Author:** [@KrishnaKanhaiya](https://discourse.julialang.org/u/KrishnaKanhaiya)\
**Replies:** 0\
**Last updated:** [March 12, 2017, 5:34am UTC](https://discourse.julialang.org/t/how-does-loess-decide-how-many-sub-intervals-to-use-for-local-regression/2615 "2017-03-12T05:34:05Z")

</div>

Hi, I want to do local regression at each and every input point. For example: using Loess xs = sort(rand(100)) #100 input points p = 1 ys = exp(p\*tpoints) + 0.1 \* rand(100) model = loess(xs,ys) when I do model.bs …

---

## [My Julia videos on YouTube](https://discourse.julialang.org/t/my-julia-videos-on-youtube/2437)

<div class="topic-metadata">

**Author:** [@parthasarathy\_gangul](https://discourse.julialang.org/u/parthasarathy_gangul)\
**Replies:** 0\
**Last updated:** [March 4, 2017, 12:17pm UTC](https://discourse.julialang.org/t/my-julia-videos-on-youtube/2437 "2017-03-04T12:17:46Z")

</div>

See my videos on Julia on YouTube. Feel free to contribute. Here is the link: Parthasarathy Ganguly

---

## [Plot histogram with Plots](https://discourse.julialang.org/t/plot-histogram-with-plots/2304)

<div class="topic-metadata">

**Author:** [@Andrea\_Tersigni](https://discourse.julialang.org/u/Andrea_Tersigni)\
**Replies:** 3\
**Last updated:** [February 25, 2017, 11:52pm UTC](https://discourse.julialang.org/t/plot-histogram-with-plots/2304 "2017-02-25T23:52:50Z")

</div>

“In convertToAnyVector, could not handle the argument types: (Float64,)”

---

## [Creation of R formula with interpolated Julia Symbol (RCall)](https://discourse.julialang.org/t/creation-of-r-formula-with-interpolated-julia-symbol-rcall/2240)

<div class="topic-metadata">

**Author:** [@Diego\_Javier\_Zea](https://discourse.julialang.org/u/Diego_Javier_Zea)\
**Replies:** 1\
**Last updated:** [February 23, 2017, 9:20am UTC](https://discourse.julialang.org/t/creation-of-r-formula-with-interpolated-julia-symbol-rcall/2240 "2017-02-23T09:20:49Z")

</div>

I have Julia Symbols with the name of the columns and I want to interpolate those to R (using RCall) to run a R’s lm(). I found and alternative using R’s as.formula to construct the formula from a string. But I was wait…

---

## [Construct a generic discrete distribution in Distributions.jl?](https://discourse.julialang.org/t/construct-a-generic-discrete-distribution-in-distributions-jl/2238)

<div class="topic-metadata">

**Author:** [@xiaoweiz](https://discourse.julialang.org/u/xiaoweiz)\
**Replies:** 11\
**Last updated:** [February 23, 2017, 6:17am UTC](https://discourse.julialang.org/t/construct-a-generic-discrete-distribution-in-distributions-jl/2238 "2017-02-23T06:17:50Z")

</div>

This may be a silly question. But what is the easies/most “julian” way to construct a generic discrete distribution type given a vector of outcomes and a vector of probabilities of the same length? It seems a workaround…

---

## [Multiple testing correction (e.g. FDR)?](https://discourse.julialang.org/t/multiple-testing-correction-e-g-fdr/2213)

<div class="topic-metadata">

**Author:** [@Ali\_Vahdati](https://discourse.julialang.org/u/Ali_Vahdati)\
**Replies:** 2\
**Last updated:** [February 21, 2017, 1:27pm UTC](https://discourse.julialang.org/t/multiple-testing-correction-e-g-fdr/2213 "2017-02-21T13:27:16Z")

</div>

Hi, I may have missed it, but is there already any function for multiple testing correction, such as FDR and Bonferroni? Thanks in advance, Ali

---

## [Parallel random number generation](https://discourse.julialang.org/t/parallel-random-number-generation/1950)

<div class="topic-metadata">

**Author:** [@kkmann](https://discourse.julialang.org/u/kkmann)\
**Replies:** 11\
**Last updated:** [February 10, 2017, 1:39pm UTC](https://discourse.julialang.org/t/parallel-random-number-generation/1950 "2017-02-10T13:39:43Z")

</div>

Hey, is the default random generation of the Distributions package safe for parallel use? I.e., are multiple streams of random numbers generated on different workers guaranteed to be independent? If not, what would be t…

---

## [Kolmogorov–Smirnov test, retrieve output](https://discourse.julialang.org/t/kolmogorov-smirnov-test-retrieve-output/1480)

<div class="topic-metadata">

**Author:** [@evolgenomology](https://discourse.julialang.org/u/evolgenomology)\
**Replies:** 2\
**Last updated:** [January 14, 2017, 1:42pm UTC](https://discourse.julialang.org/t/kolmogorov-smirnov-test-retrieve-output/1480 "2017-01-14T13:42:33Z")

</div>

Hi, I would like to use the KS test on 2000 datasets and each time I run it just retrieve the output values into an array. Being quite new to Julia, I wonder how (and if) this is possible. Currently, I am doing Approxim…

---

## [MixedModels fit error](https://discourse.julialang.org/t/mixedmodels-fit-error/1416)

<div class="topic-metadata">

**Author:** [@JordanCuevas](https://discourse.julialang.org/u/JordanCuevas)\
**Replies:** 6\
**Last updated:** [January 13, 2017, 5:19pm UTC](https://discourse.julialang.org/t/mixedmodels-fit-error/1416 "2017-01-13T17:19:06Z")

</div>

I am trying to fit a fairly complicated mixed model using Julia and the MixedModels package, on a dataset that has ~3.5 million rows (see screen shot for model specification and errors). The model is nested, with the r…

---

## [Formula creation without direct name specification](https://discourse.julialang.org/t/formula-creation-without-direct-name-specification/1448)

<div class="topic-metadata">

**Author:** [@evoclue](https://discourse.julialang.org/u/evoclue)\
**Replies:** 2\
**Last updated:** [January 12, 2017, 8:06pm UTC](https://discourse.julialang.org/t/formula-creation-without-direct-name-specification/1448 "2017-01-12T20:06:21Z")

</div>

Say I have a DataFrame called mydata where names(mydata) is Symbol\[5001\] consisting of :y, :x1, :x2, :x3, :x4, … (the other 4995 symbols) … , :x5000. I want to fit models with the formula ModelFrame(y ~ x1, mydata), Mod…

---

## [Which MCMC library for continuous time stochastic process inference?](https://discourse.julialang.org/t/which-mcmc-library-for-continuous-time-stochastic-process-inference/1422)

<div class="topic-metadata">

**Author:** [@nbren12](https://discourse.julialang.org/u/nbren12)\
**Replies:** 12\
**Last updated:** [January 12, 2017, 8:14am UTC](https://discourse.julialang.org/t/which-mcmc-library-for-continuous-time-stochastic-process-inference/1422 "2017-01-12T08:14:42Z")

</div>

I am just getting started in julia, and there seem to be many options for doing bayesian inference. I have been looking at Mamba and Klara, but I am not sure which is better for what I want to do: inference using continu…

---

## [Query on histc and histogram counts](https://discourse.julialang.org/t/query-on-histc-and-histogram-counts/1371)

<div class="topic-metadata">

**Author:** [@Donald\_Lacombe](https://discourse.julialang.org/u/Donald_Lacombe)\
**Replies:** 5\
**Last updated:** [January 9, 2017, 7:52pm UTC](https://discourse.julialang.org/t/query-on-histc-and-histogram-counts/1371 "2017-01-09T19:52:17Z")

</div>

Dear Julia Users, I want to count the number of times a zip code appears in a vector and found the following topic: I ran the following code: using StatsBase using DataFrames df1 = readtable("Test2.csv") Out\[4\]: ir…

---

## [Fitting a Distribution to existing data](https://discourse.julialang.org/t/fitting-a-distribution-to-existing-data/928)

<div class="topic-metadata">

**Author:** [@bluesmoon](https://discourse.julialang.org/u/bluesmoon)\
**Replies:** 16\
**Last updated:** [January 5, 2017, 7:26pm UTC](https://discourse.julialang.org/t/fitting-a-distribution-to-existing-data/928 "2017-01-05T19:26:18Z")

</div>

Good $(period\_of\_day) to all, Let me preface this question by saying that there is a good chance (p \> 0.8) that I don’t know what I’m talking about, so feel free to redirect me. I’ll explain what I have and what I need…

---

## [Hist() and hits() in Julia?](https://discourse.julialang.org/t/hist-and-hits-in-julia/1224)

<div class="topic-metadata">

**Author:** [@mgiugliano](https://discourse.julialang.org/u/mgiugliano)\
**Replies:** 2\
**Last updated:** [December 30, 2016, 11:45pm UTC](https://discourse.julialang.org/t/hist-and-hits-in-julia/1224 "2016-12-30T23:45:26Z")

</div>

Given a set of elements, I would love to efficiently partition their numerical values into bins, and then return the count in each bin, as well as the bin edges. MATLAB used to have super simple functions for this (i.e.…

---

## [New package: ForecastEval](https://discourse.julialang.org/t/new-package-forecasteval/774)

<div class="topic-metadata">

**Author:** [@colintbowers](https://discourse.julialang.org/u/colintbowers)\
**Replies:** 0\
**Last updated:** [December 7, 2016, 5:03am UTC](https://discourse.julialang.org/t/new-package-forecasteval/774 "2016-12-07T05:03:16Z")

</div>

I’ve just published a new package, ForecastEval, to the official METADATA. The package implements several of the statistical tests from the forecast evaluation literature including: Diebold, Mariano (1995) “Comparing P…

---

## [ERROR: type CoefTable has no field mat](https://discourse.julialang.org/t/error-type-coeftable-has-no-field-mat/688)

<div class="topic-metadata">

**Author:** [@Harish\_Kumar](https://discourse.julialang.org/u/Harish_Kumar)\
**Replies:** 1\
**Last updated:** [December 2, 2016, 2:05pm UTC](https://discourse.julialang.org/t/error-type-coeftable-has-no-field-mat/688 "2016-12-02T14:05:57Z")

</div>

After getting estimate for lmm, I am using coeftable to covert in to DF. But in 0.5 Julia i am getting below error ERROR: type CoefTable has no field mat. Is it due to any library mismatch? LMM output… Fixed-effects …

---

## [How to get p value from Chi Square table](https://discourse.julialang.org/t/how-to-get-p-value-from-chi-square-table/649)

<div class="topic-metadata">

**Author:** [@Hongwei](https://discourse.julialang.org/u/Hongwei)\
**Replies:** 6\
**Last updated:** [November 30, 2016, 5:09pm UTC](https://discourse.julialang.org/t/how-to-get-p-value-from-chi-square-table/649 "2016-11-30T17:09:27Z")

</div>

Hi guys, Given cquantile(Chisq(1), .005), we can get its Chi square table value 7.8 Is there anyway to get the corresponding p value when given Chi square table value and freedom? Thanks, Hongwei

---

## [Step-by-step hidden Markov chain log-likelihood calculation](https://discourse.julialang.org/t/step-by-step-hidden-markov-chain-log-likelihood-calculation/527)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 3\
**Last updated:** [November 24, 2016, 2:14pm UTC](https://discourse.julialang.org/t/step-by-step-hidden-markov-chain-log-likelihood-calculation/527 "2016-11-24T14:14:47Z")

</div>

I have a HMC for which the transition and observation matrices are time varying (and have some 0 elements). I would like to calculate the log-likelihood (for Bayesian inference). Which library should I use? I looked at …

---

## [PCA in MultivariateStats](https://discourse.julialang.org/t/pca-in-multivariatestats/307)

<div class="topic-metadata">

**Author:** [@jianghaizhu](https://discourse.julialang.org/u/jianghaizhu)\
**Replies:** 0\
**Last updated:** [November 14, 2016, 10:05pm UTC](https://discourse.julialang.org/t/pca-in-multivariatestats/307 "2016-11-14T22:05:12Z")

</div>

I am using PCA in MultivariateStats and trying to covert an old Matlab script to Julia. I am dealing with a serials of images. First I want to make sure I got the matrix right. I reshaped each image to a vector and put n…

---

## [Defining a factor on DataFrames](https://discourse.julialang.org/t/defining-a-factor-on-dataframes/268)

<div class="topic-metadata">

**Author:** [@hpoit](https://discourse.julialang.org/u/hpoit)\
**Replies:** 1\
**Last updated:** [November 13, 2016, 10:07am UTC](https://discourse.julialang.org/t/defining-a-factor-on-dataframes/268 "2016-11-13T10:07:53Z")

</div>

Hello! Would anyone know how to define a factor on DataFrames according to

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