无法在 CRAN solaris 机器中复制错误

问题描述

我在 CRAN 上有一个名为“metapack”的包，并收到 Ripley 教授的消息，必须解决错误才能将该包安全地保留在 CRAN 中。我在 GitHub 上的 R-CMD-check 上遇到了这个错误，但无法在我的任何本地机器上复制它来挽救我的生命。可重现的代码片段如下（我使用 Rcpp/RcppArmadillo）：

library(metapack)
data("cholesterol")
Outcome <- model.matrix(~ pldlc + phdlc + ptg,data = cholesterol)
SD <- model.matrix(~ sdldl + sdhdl + sdtg,data = cholesterol)
Trial <- cholesterol$Trial
Treat <- cholesterol$trt
Npt <- cholesterol$Npt
XCovariate <- model.matrix(~ 0 + bldlc + bhdlc + btg + age + durat + white + male + dm,data = cholesterol)
WCovariate <- model.matrix(~ trt,data = cholesterol)

fmodel <- 1
set.seed(2797542)
fit <- bayes.parobs(Outcome,SD,XCovariate,WCovariate,Treat,Trial,Npt,fmodel,mcmc = list(ndiscard = 1,nskip = 1,nkeep = 1),scale_x = TRUE,group = cholesterol$onstat,verbose = FALSE)

这给了我以下错误：

warning: chol(): given matrix is not symmetric  
error: chol(): decomposition failed

这不可能是真的，因为 (1) 它在其他机器上运行没有问题，(2) 每个 Cholesky 分解的矩阵都确保通过 Sig = 0.5 * (Sig + Sig.t()); 对称。即使使用 Sig = arma::symmatu(Sig) 似乎也不能保证对称。

有没有人遇到过类似的情况？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）