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::deme_inbreeding_spcoef_vanilla(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 Dec. 22, 2024, 1:39 a.m.