openanalytics / BIGL

GNU General Public License v3.0
7 stars 8 forks source link

Vingette produces Error in numericDeriv #8

Closed sickman2000 closed 1 year ago

sickman2000 commented 2 years ago

Hi,

I started to go through you vignette to test your package. Within the analysis workflow at:

Parameter estimates will now satisfy equality:

constraints$matrix %*% pars == constraints$vector

fitMarginals(data, transforms = transforms,constraints = constraints)

I get the error:

Error in numericDeriv(form[[3L]], names(ind), env) : Missing value or an infinity produced when evaluating the model

any ideas ? Best Sigurd

Session Info

R version 4.0.5 (2021-03-31) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1

Matrix products: default

locale: [1] LC_COLLATE=German_Austria.1252 LC_CTYPE=German_Austria.1252 LC_MONETARY=German_Austria.1252 LC_NUMERIC=C
[5] LC_TIME=German_Austria.1252

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

other attached packages: [1] minpack.lm_1.2-1 ggplot2_3.3.5 rgl_0.108.3 knitr_1.37 BIGL_1.6.5

loaded via a namespace (and not attached): [1] fs_1.5.0 bit64_4.0.5 progress_1.2.2 httr_1.4.2 numDeriv_2016.8-1.1 tools_4.0.5
[7] R6_2.5.0 affyio_1.60.0 DBI_1.1.1 BiocGenerics_0.36.0 colorspace_2.0-0 ade4_1.7-16
[13] withr_2.4.1 tidyselect_1.1.0 prettyunits_1.1.1 bit_4.0.4 curl_4.3 compiler_4.0.5
[19] preprocessCore_1.52.1 Biobase_2.50.0 airway_1.10.0 sandwich_3.0-0 labeling_0.4.2 scales_1.1.1
[25] DEoptimR_1.0-10 mvtnorm_1.1-1 robustbase_0.93-9 affy_1.68.0 multcompView_0.1-8 pkgdown_1.6.1
[31] digest_0.6.27 foreign_0.8-81 rmarkdown_2.6 rio_0.5.16 pkgconfig_2.0.3 htmltools_0.5.1
[37] plotrix_3.8-1 fastmap_1.1.0 htmlwidgets_1.5.3 rlang_0.4.10 readxl_1.3.1 RSQLite_2.2.2
[43] farver_2.0.3 generics_0.1.0 zoo_1.8-8 jsonlite_1.7.2 gtools_3.8.2 dplyr_1.0.3
[49] zip_2.1.1 car_3.0-10 magrittr_2.0.1 Matrix_1.3-2 Rcpp_1.0.6 munsell_0.5.0
[55] S4Vectors_0.28.1 abind_1.4-5 lifecycle_1.0.0 stringi_1.5.3 multcomp_1.4-16 yaml_2.2.1
[61] nleqslv_3.3.2 carData_3.0-4 MASS_7.3-53.1 zlibbioc_1.36.0 plyr_1.8.6 grid_4.0.5
[67] affxparser_1.62.0 blob_1.2.1 parallel_4.0.5 forcats_0.5.1 crayon_1.4.1 lattice_0.20-41
[73] splines_4.0.5 haven_2.3.1 annotate_1.68.0 hms_1.0.0 pillar_1.4.7 codetools_0.2-18
[79] stats4_4.0.5 XML_3.99-0.6 glue_1.4.2 drc_3.0-1 evaluate_0.14 data.table_1.14.0
[85] BiocManager_1.30.10 vctrs_0.3.6 cellranger_1.1.0 gtable_0.3.0 purrr_0.3.4 assertthat_0.2.1
[91] cachem_1.0.3 xfun_0.29 openxlsx_4.2.3 ALL_1.32.0 xtable_1.8-4 survival_3.2-10
[97] tibble_3.0.5 AnnotationDbi_1.52.0 memoise_2.0.0 IRanges_2.24.1 TH.data_1.0-10 ellipsis_0.3.1

mnazarov commented 1 year ago

The error is due to 'bad' fixed value for b, it is changed in the latest version (1.6.7). Note that this code is not actually evaluated when vignette is run (it has eval = FALSE).