search

DoseResponse / drc

Fitting dose-response models in R
https://doseresponse.github.io/drc/
21 stars 16 forks source link

`predict` not returning confidence intervals for `twophase()` model #39

Closed Perugolate closed 4 months ago

Perugolate commented 4 months ago

similar to #23 , confidence intervals are not returned by predict, for models with two phase() functions:

spinach_model1_ll4 <- drm(SLOPE~DOSE, CURVE, data = spinach, fct = LL.4())
head(predict(spinach_model1_ll4, interval = "confidence"))

     Prediction    Lower    Upper
[1,]   1.879553 1.795308 1.963798
[2,]   1.879553 1.795308 1.963798
[3,]   1.879553 1.795308 1.963798
[4,]   1.186813 1.113500 1.260126
[5,]   1.186813 1.113500 1.260126
[6,]   1.186813 1.113500 1.260126
spinach_model1_tp <- drm(SLOPE~DOSE, CURVE, data = spinach, fct = twophase())
> head(predict(spinach_model1_tp, interval = "confidence"))

[1] 1.878608 1.878608 1.878608 1.202958 1.202958 1.202958
> sessionInfo()

R version 4.3.3 (2024-02-29)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Big Sur ... 10.16

Matrix products: default
BLAS/LAPACK: /Users/paul/mambaforge/envs/nls/lib/libopenblasp-r0.3.27.dylib;  LAPACK version 
3.12.0

locale:
[1] en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8

time zone: Europe/London
tzcode source: internal

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] dplyr_1.1.4   ggplot2_3.5.1 drc_3.0-1     MASS_7.3-60  

loaded via a namespace (and not attached):
 [1] bit_4.0.5        Matrix_1.6-5     gtable_0.3.5     crayon_1.5.2    
 [5] compiler_4.3.3   gtools_3.9.5     plotrix_3.8-4    tidyselect_1.2.1
 [9] parallel_4.3.3   splines_4.3.3    scales_1.3.0     lattice_0.22-6  
[13] TH.data_1.1-2    readr_2.1.5      R6_2.5.1         generics_0.1.3  
[17] tibble_3.2.1     car_3.1-2        munsell_0.5.1    tzdb_0.4.0      
[21] pillar_1.9.0     rlang_1.1.4      utf8_1.2.4       multcomp_1.4-25 
[25] bit64_4.0.5      cli_3.6.2        withr_3.0.0      magrittr_2.0.3  
[29] grid_4.3.3       vroom_1.6.5      mvtnorm_1.2-5    hms_1.1.3       
[33] sandwich_3.1-0   lifecycle_1.0.4  vctrs_0.6.5      glue_1.7.0      
[37] codetools_0.2-20 zoo_1.8-12       survival_3.7-0   abind_1.4-5     
[41] carData_3.0-5    fansi_1.0.6      colorspace_2.1-0 tools_4.3.3     
[45] pkgconfig_2.0.3
Perugolate commented 4 months ago

closing, as likely intended behaviour