Error running spans_procedure

[lukedurell]

I run the following code:

library(pmartR)
library(pmartRdata)

mypep <- edata_transform(omicsData = pep_object, data_scale = "log2")

mypep <- group_designation(omicsData = mypep, main_effects = "Phenotype")

spans_procedure(mypep)

I see the following error:

> spans_procedure(mypep)
Error: Error in function boost::math::cdf(const chi_squared_distribution<double>&, double): Chi Square parameter was -1.1657341758564144e-15, but must be > 0 !

sessionInfo() returns:

R version 4.3.1 (2023-06-16)
Platform: aarch64-apple-darwin20 (64-bit)
Running under: macOS Ventura 13.5

Matrix products: default
BLAS:   /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib 
LAPACK: /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/lib/libRlapack.dylib;  LAPACK version 3.11.0

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

time zone: America/Chicago
tzcode source: internal

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

other attached packages:
[1] pmartRdata_2.0.0 pmartR_2.4.0    

loaded via a namespace (and not attached):
 [1] utf8_1.2.3        R6_2.5.1         
 [3] codetools_0.2-19  tidyselect_1.2.0 
 [5] magrittr_2.0.3    iterators_1.0.14 
 [7] glue_1.6.2        tibble_3.2.1     
 [9] pkgconfig_2.0.3   dplyr_1.1.2      
[11] generics_0.1.3    lifecycle_1.0.3  
[13] cli_3.6.1         foreach_1.5.2    
[15] fansi_1.0.4       vctrs_0.6.3      
[17] withr_2.5.0       compiler_4.3.1   
[19] rstudioapi_0.15.0 tools_4.3.1      
[21] pillar_1.9.0      Rcpp_1.0.11      
[23] rlang_1.1.1

[clabornd]

Thanks Luke, I can't reproduce this on my machine so there's gotta be some weird thing with a new package version or R version. I'll try reproducing with R 4.3.

[evanglass]

I cannot reproduce this either (with R version 4.3.2). Perhaps it's an Apple M chip thing?

[clabornd]

Yea this is an arm64 thing...it is showing up on some tests I'm running on my fork: https://github.com/clabornd/pmartR/actions/runs/9132936942

This doesn't show up on the CRAN checks since we are skipping SPANS tests to save time there, but is a real issue.