devmotion

using Turing
using StatsPlots

using Random
using Statistics

function demo(N::Int; n::Int = 10)
    # observation noise
    σ² = 0.3

devmotion

using Turing
using StatsPlots

using Random
using Statistics

function demo(N::Int; n::Int = 10)
    # observation noise
    σ² = 0.3

devmotion

using DelayDiffEq, DiffEqParamEstim, BlackBoxOptim, DataFrames, LsqFit
using Plots
gr()

include("importData.jl")
include("plot.jl")

# import data from the path, in which:
# pop: population data
# g1, g2: g1 and g2 data

devmotion

using MuladdMacro, EllipsisNotation, BenchmarkTools

# In place
@generated function f_loop!(out, dt, y, k, b, idxs::Nothing, ::Type{Val{N}}) where N
    sumexpr = :(+($((:(k[$i][idx]*b[$i]) for i in 1:N)...)))

    quote
        @muladd @inbounds for idx in eachindex(out)
            out[idx] = y[idx] + dt*$sumexpr
        end

devmotion

using MacroTools: namify
using Base.Meta: isexpr

## Utilities

# Split struct definition (neglecting constructors)

function splitstruct(structdef)
    # Split type definition
    isexpr(structdef, :struct, 3) || error("Not a type definition:", structdef)