# Specific Domains

**URL:** https://discourse.julialang.org/c/domain/10.md?no_subcategories=false&page=387

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

**Page:** 388

---

## [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…

---

## [Mixture models with binary indicator variables](https://discourse.julialang.org/t/mixture-models-with-binary-indicator-variables/55219)

<div class="topic-metadata">

**Author:** [@drbenvincent](https://discourse.julialang.org/u/drbenvincent)\
**Replies:** 11\
**Last updated:** [February 18, 2021, 12:43pm UTC](https://discourse.julialang.org/t/mixture-models-with-binary-indicator-variables/55219 "2021-02-18T12:43:08Z")

</div>

I’m having trouble with a number of models which all involve a mixture model with binary indicator variables. I’m reasonably sure they are all suffering from the same underlying issue, and I’m semi-convinced the problem …

---

## [Makie - Is there a \`show\_axis=false\` equivalent for the \`Axis\` object?](https://discourse.julialang.org/t/makie-is-there-a-show-axis-false-equivalent-for-the-axis-object/55535)

<div class="topic-metadata">

**Author:** [@jonalm](https://discourse.julialang.org/u/jonalm)\
**Replies:** 1\
**Last updated:** [February 18, 2021, 11:04am UTC](https://discourse.julialang.org/t/makie-is-there-a-show-axis-false-equivalent-for-the-axis-object/55535 "2021-02-18T11:04:39Z")

</div>

Is it possible to mimic the old show\_axis=false behavior for Axis objects? potentially relevant issue: show\_axis parameter is ignored · Issue #88 · jkrumbiegel/MakieLayout.jl · GitHub

---

## [Understanding the performance of Zygote](https://discourse.julialang.org/t/understanding-the-performance-of-zygote/55447)

<div class="topic-metadata">

**Author:** [@Himanshu\_Chaudhary](https://discourse.julialang.org/u/Himanshu_Chaudhary)\
**Replies:** 1\
**Last updated:** [February 18, 2021, 10:37am UTC](https://discourse.julialang.org/t/understanding-the-performance-of-zygote/55447 "2021-02-18T10:37:27Z")

</div>

I remember reading somewhere that the cost of automatic differentiation is usually less than 2 times the cost of evaluating the function itself. I tried to benchmark a small code and the results that I got were a little…

---

## [Querying WRDS Data using Julia](https://discourse.julialang.org/t/querying-wrds-data-using-julia/31835)

<div class="topic-metadata">

**Author:** [@Yifan\_Liu](https://discourse.julialang.org/u/Yifan_Liu)\
**Replies:** 13\
**Last updated:** [February 18, 2021, 3:44am UTC](https://discourse.julialang.org/t/querying-wrds-data-using-julia/31835 "2021-02-18T03:44:39Z")

</div>

I have written a short tutorial about querying WRDS data using Julia. A user only needs to input her username and password to work with WRDS remotely in Julia. Basically, I translated the existing tutorial about querying…

---

## [Makie - How to pad data inside an axis?](https://discourse.julialang.org/t/makie-how-to-pad-data-inside-an-axis/55488)

<div class="topic-metadata">

**Author:** [@Kolaru](https://discourse.julialang.org/u/Kolaru)\
**Replies:** 2\
**Last updated:** [February 18, 2021, 12:19am UTC](https://discourse.julialang.org/t/makie-how-to-pad-data-inside-an-axis/55488 "2021-02-18T00:19:23Z")

</div>

I currently have the following and I would like to put a bit more space above the curves (generated with density for the second and third). How can I achieve that? If I could access the current xlims or ylims I could…

---

## [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…

---

## [Query derivative of DAE solution](https://discourse.julialang.org/t/query-derivative-of-dae-solution/55150)

<div class="topic-metadata">

**Author:** [@ohmsweetohm1](https://discourse.julialang.org/u/ohmsweetohm1)\
**Replies:** 3\
**Last updated:** [February 17, 2021, 7:08pm UTC](https://discourse.julialang.org/t/query-derivative-of-dae-solution/55150 "2021-02-17T19:08:03Z")

</div>

Is there a way to query the derivative of the solution when solving a DAE? function f(out, du, u, p, t) out\[1\] = - 0.04u\[1\] + 1e4 \* u\[2\] \* u\[3\] - du\[1\] out\[2\] = + 0.04u\[1\] - 3e7 \* u\[2\]^2 - 1e4 \* u\[2\] \* u\[3\] - du…

---

## [DSP Filter application](https://discourse.julialang.org/t/dsp-filter-application/55156)

<div class="topic-metadata">

**Author:** [@roble](https://discourse.julialang.org/u/roble)\
**Replies:** 13\
**Last updated:** [February 17, 2021, 4:28pm UTC](https://discourse.julialang.org/t/dsp-filter-application/55156 "2021-02-17T16:28:45Z")

</div>

Hi, I am trying to use the filt(f, x\[, si\]) function from the DSP.jl module. My raw data looks like this: On this data I apply a lowpass filter: df.data\_raw # my raw data w0 = 1/120 # cutoff frequency df.data\_lp= …

---

## [Logistic Regression Mixture Model](https://discourse.julialang.org/t/logistic-regression-mixture-model/55406)

<div class="topic-metadata">

**Author:** [@Carter](https://discourse.julialang.org/u/Carter)\
**Replies:** 2\
**Last updated:** [February 17, 2021, 2:20pm UTC](https://discourse.julialang.org/t/logistic-regression-mixture-model/55406 "2021-02-17T14:20:23Z")

</div>

I’m trying to implement robust logistic regression fromchapter 10.5 of PML. The model is p(y|\\mathbf{x}) = \\pi Ber(y|0.5) + (1-\\pi) Ber(y|\\sigma(\\mathbf{w}^\\top\\mathbf{x})) where \\pi are the mixture weights drawn from a …

---

## [Trying to use textbook with ChemEquations.jl package](https://discourse.julialang.org/t/trying-to-use-textbook-with-chemequations-jl-package/55356)

<div class="topic-metadata">

**Author:** [@22natarajank](https://discourse.julialang.org/u/22natarajank)\
**Replies:** 2\
**Last updated:** [February 17, 2021, 2:07pm UTC](https://discourse.julialang.org/t/trying-to-use-textbook-with-chemequations-jl-package/55356 "2021-02-17T14:07:19Z")

</div>

I have been learning Julia along side my AP chemistry class and was trying the ChemEquations.jl. It uses keyword ce to define a chemical Equations like this. balance(ce"Al{3} + Cl{-1} = AlCl3") What is ce above as a la…

---

## [Plot function with multiple arguments](https://discourse.julialang.org/t/plot-function-with-multiple-arguments/55411)

<div class="topic-metadata">

**Author:** [@Evan](https://discourse.julialang.org/u/Evan)\
**Replies:** 3\
**Last updated:** [February 17, 2021, 9:12am UTC](https://discourse.julialang.org/t/plot-function-with-multiple-arguments/55411 "2021-02-17T09:12:13Z")

</div>

I would like to plot a function that takes multiple inputs or arguments. An example for illustrative purpose: using Plots function dgp\_DD\_t(x, treatment=1) y = (2+ treatment)\*x end plot(dgp\_DD\_t,0,50) But what i wa…

---

## [Lag with WGLMakie](https://discourse.julialang.org/t/lag-with-wglmakie/55399)

<div class="topic-metadata">

**Author:** [@yakir12](https://discourse.julialang.org/u/yakir12)\
**Replies:** 9\
**Last updated:** [February 17, 2021, 10:51am UTC](https://discourse.julialang.org/t/lag-with-wglmakie/55399 "2021-02-17T10:51:20Z")

</div>

In the following MWE (requires a connected webcam to work), I experience a substantial lag and a surprisingly low frame-rate: using VideoIO, WGLMakie, JSServe, Observables cam = VideoIO.opencamera() # open the webcam i…

---

## [JuMP/Simplex - basic question: derivatives](https://discourse.julialang.org/t/jump-simplex-basic-question-derivatives/55416)

<div class="topic-metadata">

**Author:** [@TheMEGuy](https://discourse.julialang.org/u/TheMEGuy)\
**Replies:** 3\
**Last updated:** [February 16, 2021, 10:27pm UTC](https://discourse.julialang.org/t/jump-simplex-basic-question-derivatives/55416 "2021-02-16T22:27:47Z")

</div>

Hi, I am new to optimization and I have a question for you: these might be naive questions but when I register a nonlinear function I have two options: to provide JuMP with the first two derivations to use the auto d…

---

## [Using GPU via PyCall causes non-reusable memory allocation](https://discourse.julialang.org/t/using-gpu-via-pycall-causes-non-reusable-memory-allocation/55140)

<div class="topic-metadata">

**Author:** [@JosePereiraUA](https://discourse.julialang.org/u/JosePereiraUA)\
**Replies:** 4\
**Last updated:** [February 16, 2021, 6:33pm UTC](https://discourse.julialang.org/t/using-gpu-via-pycall-causes-non-reusable-memory-allocation/55140 "2021-02-16T18:33:14Z")

</div>

My name is José Pereira, I’m a portuguese PhD student currently using Julia to develop a light-weight protein design package. I’m trying to implement TorchANI, a known ML model for molecular energy calculation, developed…

---

## [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.

---

## [Please help with my first research paper (on distributed computing, implemented in Julia)](https://discourse.julialang.org/t/please-help-with-my-first-research-paper-on-distributed-computing-implemented-in-julia/55395)

<div class="topic-metadata">

**Author:** [@tisztamo](https://discourse.julialang.org/u/tisztamo)\
**Replies:** 3\
**Last updated:** [February 16, 2021, 3:15pm UTC](https://discourse.julialang.org/t/please-help-with-my-first-research-paper-on-distributed-computing-implemented-in-julia/55395 "2021-02-16T15:15:22Z")

</div>

I am finishing my very first research paper at the age of 42 :slight\_smile: , with the help of a friend. It is about automatically restructuring complex, long-running multi-threaded/distributed computations on the fly, …

---

## [Cannot capture log with JuMP](https://discourse.julialang.org/t/cannot-capture-log-with-jump/55365)

<div class="topic-metadata">

**Author:** [@Henrique\_Becker](https://discourse.julialang.org/u/Henrique_Becker)\
**Replies:** 5\
**Last updated:** [February 16, 2021, 2:43pm UTC](https://discourse.julialang.org/t/cannot-capture-log-with-jump/55365 "2021-02-16T14:43:48Z")

</div>

I wanna capture CPLEX log in a file when it executes. However, when solving MIP instances the last bit of the log is never captured by redirect\_stdout and redirect\_stderr. I do not know the reason, I only know it happens…

---

## [Gradient of sum](https://discourse.julialang.org/t/gradient-of-sum/55389)

<div class="topic-metadata">

**Author:** [@e3c6](https://discourse.julialang.org/u/e3c6)\
**Replies:** 4\
**Last updated:** [February 16, 2021, 12:06pm UTC](https://discourse.julialang.org/t/gradient-of-sum/55389 "2021-02-16T12:06:39Z")

</div>

T = randn(10,10) ps = Flux.params(T) opt=ADAM() for iter = 1:100 gs = gradient(ps) do sum(T) end Flux.update!(opt, ps, gs) end throws an error: ArgumentError: Cannot setindex! to 0.00099999999000000…

---

## [DiffEqCallbacks problem with DDE](https://discourse.julialang.org/t/diffeqcallbacks-problem-with-dde/54921)

<div class="topic-metadata">

**Author:** [@zamk](https://discourse.julialang.org/u/zamk)\
**Replies:** 10\
**Last updated:** [February 16, 2021, 10:50am UTC](https://discourse.julialang.org/t/diffeqcallbacks-problem-with-dde/54921 "2021-02-16T10:50:48Z")

</div>

I have created a SavingCallback which prints out and saves the du, here is my code: prob = DDEProblem(model,u0,h,(0.0,120.0),p) alg = MethodOfSteps(Tsit5()) saved\_values = SavedValues(Float64,Array{Float64}) saving = Sa…

---

## [TextModels without GPU](https://discourse.julialang.org/t/textmodels-without-gpu/55374)

<div class="topic-metadata">

**Author:** [@gpucce](https://discourse.julialang.org/u/gpucce)\
**Replies:** 0\
**Last updated:** [February 16, 2021, 9:01am UTC](https://discourse.julialang.org/t/textmodels-without-gpu/55374 "2021-02-16T09:01:46Z")

</div>

I can’t add TextModels on a machine without GPU. More specifically I am trying to use in windows subsystem linux, whereas on a linux system with GPU it works fine. Does anybody know if this is expected and if there is a …

---

## [Strange behavior of JuMP 0.21.5 while using CPLEX 20.1](https://discourse.julialang.org/t/strange-behavior-of-jump-0-21-5-while-using-cplex-20-1/55329)

<div class="topic-metadata">

**Author:** [@ahmadreza-marandi](https://discourse.julialang.org/u/ahmadreza-marandi)\
**Replies:** 4\
**Last updated:** [February 16, 2021, 8:02am UTC](https://discourse.julialang.org/t/strange-behavior-of-jump-0-21-5-while-using-cplex-20-1/55329 "2021-02-16T08:02:20Z")

</div>

Dear all, I have faced a strange behavior of CPLEX 20.1 while solving a nonconvex quadratic problem. As I couldn’t replicate the issue with small instances, I am attaching the data here: data information data (in .ma…

---

## [Wrapping a simple R package](https://discourse.julialang.org/t/wrapping-a-simple-r-package/55349)

<div class="topic-metadata">

**Author:** [@Albert\_Zevelev](https://discourse.julialang.org/u/Albert_Zevelev)\
**Replies:** 2\
**Last updated:** [February 16, 2021, 4:47am UTC](https://discourse.julialang.org/t/wrapping-a-simple-r-package/55349 "2021-02-16T04:47:37Z")

</div>

I want to wrap a simple ML package: sprintr. In R the code is straight-forward: library(sprintr); m = cv.sprinter(x = X, y = y) pred = predict(m, newdata = XH) In Julia it works in 1-step: X = randn(100,5) y = 7.0 .…

---

## [This file wont precompile](https://discourse.julialang.org/t/this-file-wont-precompile/55362)

<div class="topic-metadata">

**Author:** [@DevJac](https://discourse.julialang.org/u/DevJac)\
**Replies:** 1\
**Last updated:** [February 15, 2021, 11:38pm UTC](https://discourse.julialang.org/t/this-file-wont-precompile/55362 "2021-02-15T23:38:35Z")

</div>

using Flux using Statistics using Networks export PolicyNetwork, advantage ... I have this source file in the src directory of my project. src is in the JULIA\_LOAD\_PATH, but when I attempt to use Networks Julia pre-c…

---

## [Choice of the integrator for mixed first-second order system of ODE](https://discourse.julialang.org/t/choice-of-the-integrator-for-mixed-first-second-order-system-of-ode/55330)

<div class="topic-metadata">

**Author:** [@Gregstrq](https://discourse.julialang.org/u/Gregstrq)\
**Replies:** 1\
**Last updated:** [February 15, 2021, 10:40pm UTC](https://discourse.julialang.org/t/choice-of-the-integrator-for-mixed-first-second-order-system-of-ode/55330 "2021-02-15T22:40:04Z")

</div>

I have a system of ODE with the following peculiar functional structure. It is second order for one variable, but first order for another one (actually, A and B are some vectors): A'' = f(A, B) B' = g(A, B) I need the…

---

## [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…

---

## [BoundsError](https://discourse.julialang.org/t/boundserror/53961)

<div class="topic-metadata">

**Author:** [@ashwanimalviya](https://discourse.julialang.org/u/ashwanimalviya)\
**Replies:** 16\
**Last updated:** [February 15, 2021, 4:11pm UTC](https://discourse.julialang.org/t/boundserror/53961 "2021-02-15T16:11:25Z")

</div>

I have created a function as function ShortestPath(DAG, K2Score) for i = 1:length(K2Score) ind = findall(K2Score .== 0) #K2Score\[1,ind\] ~= root, K2Score\[1,ind\] = 0.01 end Finding the root. G = DAG M,N = size\[DAG…

---

## [Apply Flux model to a matrix](https://discourse.julialang.org/t/apply-flux-model-to-a-matrix/55324)

<div class="topic-metadata">

**Author:** [@anon74562486](https://discourse.julialang.org/u/anon74562486)\
**Replies:** 0\
**Last updated:** [February 15, 2021, 3:18pm UTC](https://discourse.julialang.org/t/apply-flux-model-to-a-matrix/55324 "2021-02-15T15:18:26Z")

</div>

I have a program that outputs a matrix and I have also a NN model that takes as input one element array. I want to apply the NN to every element of the matrix. I have tried in this way my\_model(x) = reshape(reshape(x,…

---

## [Nice way to do parallel matrix multiplication](https://discourse.julialang.org/t/nice-way-to-do-parallel-matrix-multiplication/55310)

<div class="topic-metadata">

**Author:** [@Ribeiro](https://discourse.julialang.org/u/Ribeiro)\
**Replies:** 0\
**Last updated:** [February 15, 2021, 12:14pm UTC](https://discourse.julialang.org/t/nice-way-to-do-parallel-matrix-multiplication/55310 "2021-02-15T12:14:18Z")

</div>

Hello, I’m trying to do parallel matrix multiplication in a for loop. I’m using Threads.@threads, but this also applies to a GPU version of the code and maybe a distributed (MPI-like) version in the future. Within my c…

---

## [Strange issues fitting GMMs in Turing.jl by extending a fairly simple example](https://discourse.julialang.org/t/strange-issues-fitting-gmms-in-turing-jl-by-extending-a-fairly-simple-example/55232)

<div class="topic-metadata">

**Author:** [@HarrisonWilde](https://discourse.julialang.org/u/HarrisonWilde)\
**Replies:** 10\
**Last updated:** [February 15, 2021, 9:54am UTC](https://discourse.julialang.org/t/strange-issues-fitting-gmms-in-turing-jl-by-extending-a-fairly-simple-example/55232 "2021-02-15T09:54:29Z")

</div>

I have been extending the Turing example on Bayesian GMMs today in a number of ways, trying it out with non-symmetric 2D Gaussian locations and more clusters than the two in the example I have linked, and namely with som…

[Previous page](https://discourse.julialang.org/c/domain/10.md?no_subcategories=false&page=386)

[Next page](https://discourse.julialang.org/c/domain/10.md?no_subcategories=false&page=388)