error in seek_inflection.R

[florale]

Hi Salma and Manuel,

I’m exploring your package and have run into an error in a few cases. I’ve noticed these cases often include descents, but wanted to check whether that might be why a DLMO time wasn’t produced. I’ve included reproducible code below.

test_data <- data.table(
  datetime = as.POSIXct(c(
    "2023-11-24 09:00:00",
    "2023-11-24 10:00:00",
    "2023-11-24 11:00:00",
    "2023-11-24 12:00:00",
    "2023-11-23 13:00:00",
    "2023-11-23 14:00:00",
    "2023-11-23 15:00:00"
  ), tz = "UTC"),
  melatonin = c(1.0, 3.4, 7.3, 11.2, 13.5, 10.8, 12.0)
)
dlmo_test <- calculate_dlmo(
  data = test_data,
  threshold = 4,
  fine_flag = TRUE
)

Was wondering if you have any insights? Thank you!

[salmathalji]

Hi Flora @florale ,
Thank you for the reproducible example. I was able to reproduce the error and have pushed a bug fix to the repo

. The fix prevents non-finite residuals from causing that cryptic error. The updated version now handles this case properly. With the fix applied, your example runs successfully and returns:

• coarse DLMO: 09:21:00
• fine DLMO: 09:24:00

I’ve attached the resulting figures here for your reference. Please don't hesitate to get in touch if anything else comes up.

dlmo_result_plot_fine.pdf
dlmo_test_plot_coarse.pdf