In nickbrazeau/discent: Estimation of Deme Inbreeding Spatial Coefficients with Gradient Descent

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(discent) 

***

data("IBD_simulation_data")

start_params <- rep(0.3, 3)
names(start_params) <- 1:3
start_params <- c(start_params, "m"= 250)

#run
ret <- discent::disc(discdat = IBD_simulation_data,
                     start_params = start_params,
                     learningrate = 1e-3,
                     lambda = 0.01,
                     m_lowerbound = 0,
                     m_upperbound = Inf,
                     b1 = 0.9,
                     b2 = 0.999,
                     e = 1e-8,
                     steps = 1e4,
                     thin = 1,
                     normalize_geodist = TRUE,
                     report_progress = TRUE,
                     return_verbose = TRUE)

nickbrazeau/discent documentation built on June 13, 2025, 5:13 a.m.