R Under development (unstable) (2026-01-07 r89288) -- "Unsuffered Consequences" Copyright (C) 2026 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > pkgname <- "spBFA" > source(file.path(R.home("share"), "R", "examples-header.R")) > options(warn = 1) > library('spBFA') > > base::assign(".oldSearch", base::search(), pos = 'CheckExEnv') > base::assign(".old_wd", base::getwd(), pos = 'CheckExEnv') > cleanEx() > nameEx("bfa_sp") > ### * bfa_sp > > flush(stderr()); flush(stdout()) > > ### Name: bfa_sp > ### Title: Spatial factor analysis using a Bayesian hierarchical model. > ### Aliases: bfa_sp > > ### ** Examples > > > > > > cleanEx() > nameEx("diagnostics") > ### * diagnostics > > flush(stderr()); flush(stdout()) > > ### Name: diagnostics > ### Title: diagnostics > ### Aliases: diagnostics > > ### ** Examples > > ###Load pre-computed regression results > data(reg.bfa_sp) > > ###Compute and print diagnostics > diags <- diagnostics(reg.bfa_sp) 2 Calculating Log-Lik: 0%.. /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/subview_meat.hpp:1413:54: runtime error: reference binding to null pointer of type 'const double' #0 0x7b5c5f51972f in arma::subview::colptr(unsigned int) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/subview_meat.hpp:1413 #1 0x7b5c5f51972f in arma::subview_row::subview_row(arma::Mat const&, unsigned int) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/subview_meat.hpp:4271 #2 0x7b5c5f528e7b in arma::Mat::row(unsigned int) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/Mat_meat.hpp:3350 #3 0x7b5c5f528e7b in GetLogLik(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>, int, bool) /data/gannet/ripley/R/packages/tests-gcc-SAN/spBFA/src/DIAG_GetLogLik.cpp:63 #4 0x7b5c5f71414f in _spBFA_GetLogLik /data/gannet/ripley/R/packages/tests-gcc-SAN/spBFA/src/RcppExports.cpp:63 #5 0x00000074f306 in R_doDotCall /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:763 #6 0x000000753b8a in do_dotcall /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:1437 #7 0x0000008ceba2 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8132 #8 0x0000008bbdeb in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7515 #9 0x00000085c642 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167 #10 0x00000087238a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2389 #11 0x00000087644a in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2302 #12 0x00000085cce3 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2324 #13 0x00000085cce3 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1280 #14 0x00000088abfe in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3581 #15 0x00000085d106 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1232 #16 0x000000a02488 in Rf_ReplIteration /data/gannet/ripley/R/svn/R-devel/src/main/main.c:264 #17 0x000000a02488 in R_ReplConsole /data/gannet/ripley/R/svn/R-devel/src/main/main.c:317 #18 0x000000a1056a in run_Rmainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1235 #19 0x000000a10602 in Rf_mainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1242 #20 0x0000004131bf in main /data/gannet/ripley/R/svn/R-devel/src/main/Rmain.c:29 #21 0x7f5c79e115f4 in __libc_start_call_main (/lib64/libc.so.6+0x35f4) (BuildId: a1dda014206b55b07f58fe8db80121b752dc3d03) #22 0x7f5c79e116a7 in __libc_start_main@@GLIBC_2.34 (/lib64/libc.so.6+0x36a7) (BuildId: a1dda014206b55b07f58fe8db80121b752dc3d03) #23 0x000000413ba4 in _start (/data/gannet/ripley/R/gcc-SAN3/bin/exec/R+0x413ba4) (BuildId: 09d1e93a0fe813e632cb484c467b0eaf790fc71f) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/access.hpp:26:100: runtime error: reference binding to null pointer of type 'double' #0 0x7b5c5f5195e8 in double& arma::access::rw(double const&) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/access.hpp:26 #1 0x7b5c5f5195e8 in arma::subview::colptr(unsigned int) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/subview_meat.hpp:1413 #2 0x7b5c5f5195e8 in arma::subview_row::subview_row(arma::Mat const&, unsigned int) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/subview_meat.hpp:4271 #3 0x7b5c5f528e7b in arma::Mat::row(unsigned int) /data/gannet/ripley/R/test-dev/RcppArmadillo/include/current/armadillo_bits/Mat_meat.hpp:3350 #4 0x7b5c5f528e7b in GetLogLik(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>, int, bool) /data/gannet/ripley/R/packages/tests-gcc-SAN/spBFA/src/DIAG_GetLogLik.cpp:63 #5 0x7b5c5f71414f in _spBFA_GetLogLik /data/gannet/ripley/R/packages/tests-gcc-SAN/spBFA/src/RcppExports.cpp:63 #6 0x00000074f306 in R_doDotCall /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:763 #7 0x000000753b8a in do_dotcall /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:1437 #8 0x0000008ceba2 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8132 #9 0x0000008bbdeb in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7515 #10 0x00000085c642 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167 #11 0x00000087238a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2389 #12 0x00000087644a in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2302 #13 0x00000085cce3 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2324 #14 0x00000085cce3 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1280 #15 0x00000088abfe in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3581 #16 0x00000085d106 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1232 #17 0x000000a02488 in Rf_ReplIteration /data/gannet/ripley/R/svn/R-devel/src/main/main.c:264 #18 0x000000a02488 in R_ReplConsole /data/gannet/ripley/R/svn/R-devel/src/main/main.c:317 #19 0x000000a1056a in run_Rmainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1235 #20 0x000000a10602 in Rf_mainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1242 #21 0x0000004131bf in main /data/gannet/ripley/R/svn/R-devel/src/main/Rmain.c:29 #22 0x7f5c79e115f4 in __libc_start_call_main (/lib64/libc.so.6+0x35f4) (BuildId: a1dda014206b55b07f58fe8db80121b752dc3d03) #23 0x7f5c79e116a7 in __libc_start_main@@GLIBC_2.34 (/lib64/libc.so.6+0x36a7) (BuildId: a1dda014206b55b07f58fe8db80121b752dc3d03) #24 0x000000413ba4 in _start (/data/gannet/ripley/R/gcc-SAN3/bin/exec/R+0x413ba4) (BuildId: 09d1e93a0fe813e632cb484c467b0eaf790fc71f) 25%.. 50%.. 75%.. 100%.. Done! Calculating PPD: 0%.. 25%.. 50%.. 75%.. 100%.. Done! > print(unlist(diags)) dic.dic dic.pd dinf.p dinf.g dinf.dinf -100.57865 124.26469 214.96291 66.53240 281.49531 waic.waic waic.p_waic waic.lppd waic.p_waic_1 -92.93843 84.67741 131.14662 37.44991 > > > > > cleanEx() > nameEx("is.spBFA") > ### * is.spBFA > > flush(stderr()); flush(stdout()) > > ### Name: is.spBFA > ### Title: is.spBFA > ### Aliases: is.spBFA > > ### ** Examples > > ###Load pre-computed results > data(reg.bfa_sp) > > ###Test function > is.spBFA(reg.bfa_sp) [1] TRUE > > > > > cleanEx() > nameEx("predict.spBFA") > ### * predict.spBFA > > flush(stderr()); flush(stdout()) > > ### Name: predict.spBFA > ### Title: predict.spBFA > ### Aliases: predict.spBFA > > ### ** Examples > > ###Load pre-computed regression results > data(reg.bfa_sp) > > ###Compute predictions > pred <- predict(reg.bfa_sp, NewTimes = 3) Krigging Eta: 0%.. 25%.. 50%.. 75%.. 100%.. Done! Krigging Y: 0%.. 25%.. 50%.. 75%.. 100%.. Done! > pred.observations <- pred$Y$Y10 # observed data predictions > pred.krig <- pred$Eta$Eta10 # krigged parameters > > > > > ### *