# Specific Domains

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

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

**Page:** 259

---

## [Recommended sampler for data with Categorical likelihood](https://discourse.julialang.org/t/recommended-sampler-for-data-with-categorical-likelihood/85820)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 1\
**Last updated:** [August 16, 2022, 9:36pm UTC](https://discourse.julialang.org/t/recommended-sampler-for-data-with-categorical-likelihood/85820 "2022-08-16T21:36:27Z")

</div>

The NUTS sampler takes a really long time in this MWE: using Distributions using Turing d = Categorical(\[0.7, 0.29, 0.01\]) N = 100 Y = rand(d, N) @model function test(Y) p ~ Dirichlet(\[2,2,2\]) for i in 1:N …

---

## [Visualising NetCDF using Raster.jl](https://discourse.julialang.org/t/visualising-netcdf-using-raster-jl/85761)

<div class="topic-metadata">

**Author:** [@Michaela\_Vollmer](https://discourse.julialang.org/u/Michaela_Vollmer)\
**Replies:** 10\
**Last updated:** [August 16, 2022, 9:29pm UTC](https://discourse.julialang.org/t/visualising-netcdf-using-raster-jl/85761 "2022-08-16T21:29:25Z")

</div>

Hello! I am trying to plot spatial data of a NetCDF file from a WFlow example. I have followed the steps outlined in the following question Raster.jl discourse post where a similar example is shown where the data can be …

---

## [Simulation example fail](https://discourse.julialang.org/t/simulation-example-fail/85823)

<div class="topic-metadata">

**Author:** [@tsugiharu](https://discourse.julialang.org/u/tsugiharu)\
**Replies:** 2\
**Last updated:** [August 16, 2022, 9:20pm UTC](https://discourse.julialang.org/t/simulation-example-fail/85823 "2022-08-16T21:20:21Z")

</div>

Hi all, I just fired up Julia for the first time. My goal is to use the JEMSS package for research purposes. First of all,I could install and run JEMSS. But,When I ran the Auckland City’s data with the file name unchan…

---

## [JuMPing with Unitful?](https://discourse.julialang.org/t/jumping-with-unitful/85839)

<div class="topic-metadata">

**Author:** [@bcon](https://discourse.julialang.org/u/bcon)\
**Replies:** 3\
**Last updated:** [August 16, 2022, 9:04pm UTC](https://discourse.julialang.org/t/jumping-with-unitful/85839 "2022-08-16T21:04:04Z")

</div>

The function I would like to optimize uses Unitful variables in its internal calculations, but when optimizing I receive MethodError: no method matching \_float(::Quantity{Int64, 𝐋, Unitful.FreeUnits{(mm,), 𝐋, nothing}}),…

---

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

---

## [Problem when intall MathTexEngine v0.5.0](https://discourse.julialang.org/t/problem-when-intall-mathtexengine-v0-5-0/85817)

<div class="topic-metadata">

**Author:** [@zhzhang](https://discourse.julialang.org/u/zhzhang)\
**Replies:** 4\
**Last updated:** [August 16, 2022, 1:30pm UTC](https://discourse.julialang.org/t/problem-when-intall-mathtexengine-v0-5-0/85817 "2022-08-16T13:30:48Z")

</div>

Hi all, I want to update the MathTexEngine from 0.4.3 to 0.5.0. However, the Pkg threw these errors. (@v1.7) pkg\> add MathTeXEngine@0.5 Resolving package versions... ERROR: Unsatisfiable requirements detected for pac…

---

## [ANN: upcoming refactoring of JuMP's nonlinear API](https://discourse.julialang.org/t/ann-upcoming-refactoring-of-jumps-nonlinear-api/83052)

<div class="topic-metadata">

**Author:** [@odow](https://discourse.julialang.org/u/odow)\
**Replies:** 9\
**Last updated:** [August 16, 2022, 8:24am UTC](https://discourse.julialang.org/t/ann-upcoming-refactoring-of-jumps-nonlinear-api/83052 "2022-08-16T08:24:23Z")

</div>

tl;dr: The next feature release of JuMP (probably v1.2.0) introduces a large refactoring of JuMP’s nonlinear API. If you have code that accessed private features of the nonlinear API, such as JuMP.\_Derivatives or mode…

---

## [Basins of attractros for 6D system DynamicalSystems.jl](https://discourse.julialang.org/t/basins-of-attractros-for-6d-system-dynamicalsystems-jl/85781)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 11\
**Last updated:** [August 16, 2022, 7:16am UTC](https://discourse.julialang.org/t/basins-of-attractros-for-6d-system-dynamicalsystems-jl/85781 "2022-08-16T07:16:24Z")

</div>

I trye found basins for attractors, but i have some problems. First problem is out of memmory when i use big length for xg1 = yg1 = zg1 = xg2 = yg2 = zg2 = range(-6.0, 6.0; length = 10) If i using for example length …

---

## [How to detect a container resize event in GTK](https://discourse.julialang.org/t/how-to-detect-a-container-resize-event-in-gtk/59624)

<div class="topic-metadata">

**Author:** [@johntrinder](https://discourse.julialang.org/u/johntrinder)\
**Replies:** 2\
**Last updated:** [August 16, 2022, 12:51am UTC](https://discourse.julialang.org/t/how-to-detect-a-container-resize-event-in-gtk/59624 "2022-08-16T00:51:28Z")

</div>

I’ve scoured google but came up with nothing. I’ve done the following: id = signal\_connect(alig, “check-resize”) do widget, event println(“The window was resized”) end where ‘alig’ is a GtkAlignment to which I add a…

---

## [User defined objective function from external pkl file](https://discourse.julialang.org/t/user-defined-objective-function-from-external-pkl-file/85759)

<div class="topic-metadata">

**Author:** [@chenr86](https://discourse.julialang.org/u/chenr86)\
**Replies:** 4\
**Last updated:** [August 15, 2022, 6:30pm UTC](https://discourse.julialang.org/t/user-defined-objective-function-from-external-pkl-file/85759 "2022-08-15T18:30:04Z")

</div>

Hello, I have a complex function wrapped in a pkl file and I want to use to define my objective function. Here is a small example that I want to realize and the problem I meet: using JuMP using PyCall py""" import pic…

---

## [DiffEqParamEstim.jl - Generalized Likelihood Inference](https://discourse.julialang.org/t/diffeqparamestim-jl-generalized-likelihood-inference/85709)

<div class="topic-metadata">

**Author:** [@xztt](https://discourse.julialang.org/u/xztt)\
**Replies:** 2\
**Last updated:** [August 15, 2022, 5:39pm UTC](https://discourse.julialang.org/t/diffeqparamestim-jl-generalized-likelihood-inference/85709 "2022-08-15T17:39:10Z")

</div>

Hi, I am looking to implement the method that is found here I can follow the method well but I can not for the life of me see a way in which I could extend this to a problem to a situation where we have got n\>2 paramete…

---

## [Implementing the tensor algebra functor](https://discourse.julialang.org/t/implementing-the-tensor-algebra-functor/64433)

<div class="topic-metadata">

**Author:** [@choquet](https://discourse.julialang.org/u/choquet)\
**Replies:** 0\
**Last updated:** [July 10, 2021, 5:32pm UTC](https://discourse.julialang.org/t/implementing-the-tensor-algebra-functor/64433 "2021-07-10T17:32:28Z")

</div>

I’m trying to implement the free tensor algebra functor T: given a vector space V it creates another vector space T(V). This vector space T(V) has many nice properties, for example, T(V) is not just a vector space but al…

---

## [Orbitdiagram DynamicalSystems.jl](https://discourse.julialang.org/t/orbitdiagram-dynamicalsystems-jl/85766)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 0\
**Last updated:** [August 15, 2022, 9:43am UTC](https://discourse.julialang.org/t/orbitdiagram-dynamicalsystems-jl/85766 "2022-08-15T09:43:04Z")

</div>

Parameter that changes at high values leads to synchronization, that is, a regular mode appears. Ideally, the diagram should show the transition from chaos to regularity. The diagram show this, but there are windows that…

---

## [Agents.elastic\_collision!() makes ageents move faster and faster](https://discourse.julialang.org/t/agents-elastic-collision-makes-ageents-move-faster-and-faster/84740)

<div class="topic-metadata">

**Author:** [@Wen-Wei\_Tseng](https://discourse.julialang.org/u/Wen-Wei_Tseng)\
**Replies:** 1\
**Last updated:** [August 14, 2022, 12:18pm UTC](https://discourse.julialang.org/t/agents-elastic-collision-makes-ageents-move-faster-and-faster/84740 "2022-08-14T12:18:41Z")

</div>

I tried replicating the results of social distancing example in Agents.jl Example Zoo. And the agents (balls) moved faster and faster in the output video. I was using Julia 1.8.0-rc3 (I’ve tried Julia 1.7.3, but the res…

---

## [Tullio without product](https://discourse.julialang.org/t/tullio-without-product/85708)

<div class="topic-metadata">

**Author:** [@qwerty](https://discourse.julialang.org/u/qwerty)\
**Replies:** 1\
**Last updated:** [August 13, 2022, 5:13pm UTC](https://discourse.julialang.org/t/tullio-without-product/85708 "2022-08-13T17:13:26Z")

</div>

How can I replace the \* with a function? @tullio c := a\[i\]\*b\[i\] a is an array of symbolic differential and b is an array of symbolic functions. I would like to get: a\[1\](b\[1\])+a\[2\](b\[2\]) Thank you

---

## [PyCall Holoviews output in Jupyter notebook](https://discourse.julialang.org/t/pycall-holoviews-output-in-jupyter-notebook/85707)

<div class="topic-metadata">

**Author:** [@ea42gh](https://discourse.julialang.org/u/ea42gh)\
**Replies:** 1\
**Last updated:** [August 13, 2022, 3:36pm UTC](https://discourse.julialang.org/t/pycall-holoviews-output-in-jupyter-notebook/85707 "2022-08-13T15:36:35Z")

</div>

Is it possible to get the pyviz stack working with PyCall in a Jupyter notebook? using PyCall hv = pyimport("holoviews") hv.extension("bokeh" ) t=range(-1,1,10) hv.Curve((t,t .^ 2)) currently yields PyObject :Curve \[…

---

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

---

## [Q: FiniteDifferences | gradient of a 3D scalar array](https://discourse.julialang.org/t/q-finitedifferences-gradient-of-a-3d-scalar-array/64065)

<div class="topic-metadata">

**Author:** [@Audrius-St](https://discourse.julialang.org/u/Audrius-St)\
**Replies:** 5\
**Last updated:** [August 13, 2022, 2:05am UTC](https://discourse.julialang.org/t/q-finitedifferences-gradient-of-a-3d-scalar-array/64065 "2022-08-13T02:05:53Z")

</div>

Hello, I have a 3D array consisting of numerical scalar values: imgVlm::Array{Float64, 3} The grid spacing along each axis is uniform, but the spacing values are different for each axis. FiniteDifferences seems to be…

---

## [How to take the 3D fft of a vector?](https://discourse.julialang.org/t/how-to-take-the-3d-fft-of-a-vector/85669)

<div class="topic-metadata">

**Author:** [@Nic\_Leblanc](https://discourse.julialang.org/u/Nic_Leblanc)\
**Replies:** 8\
**Last updated:** [August 12, 2022, 6:22pm UTC](https://discourse.julialang.org/t/how-to-take-the-3d-fft-of-a-vector/85669 "2022-08-12T18:22:44Z")

</div>

Hi all, I was wondering if there were a way to take the 3D fft of a vector. I am obviously using FFTW for my direct fft calculations using plans but I’m not sure for it to calculate the 3D fft and not just a 1D fft. The…

---

## [How to group/color lines per ID in Makie?](https://discourse.julialang.org/t/how-to-group-color-lines-per-id-in-makie/85622)

<div class="topic-metadata">

**Author:** [@aelmokadem](https://discourse.julialang.org/u/aelmokadem)\
**Replies:** 8\
**Last updated:** [August 12, 2022, 6:06pm UTC](https://discourse.julialang.org/t/how-to-group-color-lines-per-id-in-makie/85622 "2022-08-12T18:06:31Z")

</div>

Hi all. I am using CairoMakie.jl to plot a figure for longitudinal theophylline concentration data in 12 subjects simultaneously. The code I am using is Theoph = dataset("datasets", "Theoph") f = Figure() Axis(f\[1,1\]) C…

---

## [1000x1000 heatmap with pyplot() won't render in VS Code plots pane](https://discourse.julialang.org/t/1000x1000-heatmap-with-pyplot-wont-render-in-vs-code-plots-pane/85671)

<div class="topic-metadata">

**Author:** [@thompsonmj](https://discourse.julialang.org/u/thompsonmj)\
**Replies:** 0\
**Last updated:** [August 12, 2022, 5:25pm UTC](https://discourse.julialang.org/t/1000x1000-heatmap-with-pyplot-wont-render-in-vs-code-plots-pane/85671 "2022-08-12T17:25:02Z")

</div>

using Plots pyplot() @time heatmap(rand(1,1)) @time heatmap(rand(10,10)) @time heatmap(rand(100,100)) @time heatmap(rand(1000,1000)) Each takes up to a handful of milliseconds to run, but the 1k x 1k grid simply won’t r…

---

## [How to solve the problem of julia loading Moesk?](https://discourse.julialang.org/t/how-to-solve-the-problem-of-julia-loading-moesk/85467)

<div class="topic-metadata">

**Author:** [@sand\_feifsndali](https://discourse.julialang.org/u/sand_feifsndali)\
**Replies:** 2\
**Last updated:** [August 12, 2022, 3:13am UTC](https://discourse.julialang.org/t/how-to-solve-the-problem-of-julia-loading-moesk/85467 "2022-08-12T03:13:23Z")

</div>

Downloaded a MosekTools, but when using MosekTools, it will ask to build (“Mosek”), and there are various problems when building, as shown in the figure How to solve this problem, do I need to download a moesk softwa…

---

## [IloIfThen constraints for CPLEX in Julia](https://discourse.julialang.org/t/iloifthen-constraints-for-cplex-in-julia/85619)

<div class="topic-metadata">

**Author:** [@Deoch](https://discourse.julialang.org/u/Deoch)\
**Replies:** 3\
**Last updated:** [August 11, 2022, 11:58pm UTC](https://discourse.julialang.org/t/iloifthen-constraints-for-cplex-in-julia/85619 "2022-08-11T23:58:19Z")

</div>

Hello, I need to use the conditional constraints for a JuMP model with CPLEX, that would be IloIfThen in C++. Like if I have three variables x, y and z, a constraint like: If (x \>= 1) && (y \<= 2) then z == 1 I am usin…

---

## [Automatic update in Makie with Observable of vector](https://discourse.julialang.org/t/automatic-update-in-makie-with-observable-of-vector/85623)

<div class="topic-metadata">

**Author:** [@filchristou](https://discourse.julialang.org/u/filchristou)\
**Replies:** 5\
**Last updated:** [August 11, 2022, 8:28pm UTC](https://discourse.julialang.org/t/automatic-update-in-makie-with-observable-of-vector/85623 "2022-08-11T20:28:22Z")

</div>

I need to iteratively call a Makie recipe according to an observable value. Lets say I write a recipe MyPlot and I want to iteratively call Scatter inside. using Makie @recipe(MyPlot, x) do scene Theme( re…

---

## [X and y axis equal scaling](https://discourse.julialang.org/t/x-and-y-axis-equal-scaling/30366)

<div class="topic-metadata">

**Author:** [@squirrel](https://discourse.julialang.org/u/squirrel)\
**Replies:** 2\
**Last updated:** [August 11, 2022, 7:40pm UTC](https://discourse.julialang.org/t/x-and-y-axis-equal-scaling/30366 "2022-08-11T19:40:08Z")

</div>

How do I make the x and y axis have equal scaling? If can cheat it with size = (800,800), for example, but then if I add labels or a title, or change the axis limits, this requires a lot of adjusting. Is there any way I…

---

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

---

## [Lyapunov spectrum, HyperChaos, DynamicalSystems.jl](https://discourse.julialang.org/t/lyapunov-spectrum-hyperchaos-dynamicalsystems-jl/85539)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 8\
**Last updated:** [August 11, 2022, 1:19pm UTC](https://discourse.julialang.org/t/lyapunov-spectrum-hyperchaos-dynamicalsystems-jl/85539 "2022-08-11T13:19:53Z")

</div>

I calculate lyapunov spectrum in hyperchaos. Why i have such spectrum: \[0.015581182326147337, 0.0028325903213687266, 1.0762649529273773e-6, -0.007004756972893843, -5.7866575928839605, -11.883439707053556\] For hyp…

---

## [How do I do a fast cumulative integration](https://discourse.julialang.org/t/how-do-i-do-a-fast-cumulative-integration/85574)

<div class="topic-metadata">

**Author:** [@owiecc](https://discourse.julialang.org/u/owiecc)\
**Replies:** 8\
**Last updated:** [August 11, 2022, 1:14pm UTC](https://discourse.julialang.org/t/how-do-i-do-a-fast-cumulative-integration/85574 "2022-08-11T13:14:02Z")

</div>

What are my options for a fast cumulative integration? I can use one of many quadrature packages but this does not scale. using QuadGK f = cos x₁ = 0.0 int\_f(x₂) = quadgk(f, x₁, x₂)\[1\] using Plots plot(\[f, int\_f\], xlim…

---

## [Looking for co-authors for animal movement and epidemiology ABM](https://discourse.julialang.org/t/looking-for-co-authors-for-animal-movement-and-epidemiology-abm/85614)

<div class="topic-metadata">

**Author:** [@franfram](https://discourse.julialang.org/u/franfram)\
**Replies:** 1\
**Last updated:** [August 11, 2022, 1:13pm UTC](https://discourse.julialang.org/t/looking-for-co-authors-for-animal-movement-and-epidemiology-abm/85614 "2022-08-11T13:13:56Z")

</div>

Hello everyone, I just started building an ABM for animal movement and epidemiology (PhD Project) Is anyone working on something similar and would like to team up? Kindly, Francis

---

## [How to run many copies of a random function in parallel on GPU?](https://discourse.julialang.org/t/how-to-run-many-copies-of-a-random-function-in-parallel-on-gpu/85602)

<div class="topic-metadata">

**Author:** [@GeekyPeas](https://discourse.julialang.org/u/GeekyPeas)\
**Replies:** 5\
**Last updated:** [August 11, 2022, 7:02am UTC](https://discourse.julialang.org/t/how-to-run-many-copies-of-a-random-function-in-parallel-on-gpu/85602 "2022-08-11T07:02:25Z")

</div>

Consider this “embarrassingly parallel” and fairly general problem: You have a random function f (such as a random number generator, or a sampling procedure, or a run of a Markov Chain, or a Jump Process) that is imple…

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

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