* using log directory ‘/data/gannet/ripley/R/packages/tests-OpenBLAS/robsurvey.Rcheck’ * using R Under development (unstable) (2026-01-28 r89344) * using platform: x86_64-pc-linux-gnu * R was compiled by gcc (GCC) 15.1.1 20250521 (Red Hat 15.1.1-2) GNU Fortran (GCC) 15.1.1 20250521 (Red Hat 15.1.1-2) * running under: Fedora Linux 42 (Workstation Edition) * using session charset: UTF-8 * using option ‘--no-stop-on-test-error’ * checking for file ‘robsurvey/DESCRIPTION’ ... OK * checking extension type ... Package * this is package ‘robsurvey’ version ‘0.7-1’ * package encoding: UTF-8 * checking package namespace information ... OK * checking package dependencies ... OK * checking if this is a source package ... OK * checking if there is a namespace ... OK * checking for executable files ... OK * checking for hidden files and directories ... OK * checking for portable file names ... OK * checking for sufficient/correct file permissions ... OK * checking whether package ‘robsurvey’ can be installed ... [23s/77s] OK * used C compiler: ‘gcc (GCC) 15.1.1 20250521 (Red Hat 15.1.1-2)’ * checking package directory ... OK * checking ‘build’ directory ... OK * checking DESCRIPTION meta-information ... OK * checking top-level files ... OK * checking for left-over files ... OK * checking index information ... OK * checking package subdirectories ... OK * checking code files for non-ASCII characters ... OK * checking R files for syntax errors ... OK * checking whether the package can be loaded ... [3s/13s] OK * checking whether the package can be loaded with stated dependencies ... [4s/15s] OK * checking whether the package can be unloaded cleanly ... [4s/12s] OK * checking whether the namespace can be loaded with stated dependencies ... OK * checking whether the namespace can be unloaded cleanly ... [4s/16s] OK * checking loading without being on the library search path ... [4s/12s] OK * checking whether startup messages can be suppressed ... OK * checking use of S3 registration ... OK * checking dependencies in R code ... OK * checking S3 generic/method consistency ... OK * checking replacement functions ... OK * checking foreign function calls ... OK * checking R code for possible problems ... [29s/113s] OK * checking Rd files ... OK * checking Rd metadata ... OK * checking Rd line widths ... OK * checking Rd cross-references ... OK * checking for missing documentation entries ... OK * checking for code/documentation mismatches ... OK * checking Rd \usage sections ... OK * checking Rd contents ... OK * checking for unstated dependencies in examples ... OK * checking contents of ‘data’ directory ... OK * checking data for non-ASCII characters ... OK * checking LazyData ... OK * checking data for ASCII and uncompressed saves ... OK * checking line endings in C/C++/Fortran sources/headers ... OK * checking line endings in Makefiles ... OK * checking compilation flags in Makevars ... OK * checking for GNU extensions in Makefiles ... OK * checking for portable use of $(BLAS_LIBS) and $(LAPACK_LIBS) ... OK * checking use of PKG_*FLAGS in Makefiles ... OK * checking use of SHLIB_OPENMP_*FLAGS in Makefiles ... OK * checking pragmas in C/C++ headers and code ... OK * checking compilation flags used ... OK * checking compiled code ... OK * checking sizes of PDF files under ‘inst/doc’ ... OK * checking installed files from ‘inst/doc’ ... OK * checking files in ‘vignettes’ ... OK * checking examples ... [6s/19s] OK * checking for unstated dependencies in ‘tests’ ... OK * checking tests ... Running ‘basic.R’ [4s/20s] Running ‘check_functions.R’ Running ‘domain.R’ [4s/19s] Running ‘greg.R’ [4s/15s] Running ‘huber2.R’ [4s/13s] Running ‘ratio.R’ [4s/13s] Running ‘reg.R’ [4s/10s] [24s/93s] ERROR Running the tests in ‘tests/reg.R’ failed. Complete output: > suppressPackageStartupMessages(library("survey")) > library("robsurvey", quietly = TRUE) > attach(workplace) > > # design object > dn <- svydesign(ids = ~ID, strata = ~strat, fpc = ~fpc, weights = ~weight, + data = workplace) > > source("check_functions.R") > > #=============================================================================== > # 1 Test against MASS::rlm > #=============================================================================== > if (requireNamespace("MASS", quietly = TRUE)) { + library(MASS) + # make a copy of MASS::rlm + rlm_mod <- MASS:::rlm.default + # replace wmad-function with our weighted mad + body(rlm_mod)[[4]] <- substitute(wmad <- function(x, w) weighted_mad(x, w)) + + #--------------------------------------------------------------------------- + # Huber regression M-estimator + est <- svyreg_huberM(payroll ~ employment, dn, k = 1.345, tol = 1e-9) + ref <- rlm_mod(x = est$model$x, y = est$model$y, weights = est$model$w, + k = 1.345, scale.est = "MAD", method = "M", wt.method = "case", + acc = 1e-9, maxit = 50, test.vec = "coef") + all_equal(coef(est), coef(ref), tolerance = 1e-7, + label = "Huber regression regression: coefficients") + # model-based covariance matrix + all_equal(vcov(est, "model"), vcov(ref), tolerance = 1e-7, + label = "Huber regression M-est: model-based cov") + # design-based covariance matrix (test against fixed values) + ref <- matrix(c(61704317.084028, -3367311.1407500, -3367311.140750, + 870905.6027803), ncol = 2) + colnames(ref) <- rownames(ref) <- c("(Intercept)", "employment") + all_equal(vcov(est), ref, + label = "Huber regression M-est: design-based cov") + + #--------------------------------------------------------------------------- + # Tukey regression M-estimator + est <- svyreg_tukeyM(payroll ~ employment, dn, k = 4.6, tol = 1e-9) + ref <- rlm_mod(x = est$model$x, y = est$model$y, weights = est$model$w, + c = 4.6, psi = psi.bisquare, scale.est = "MAD", method = "M", + wt.method = "case", acc = 1e-9, maxit = 50, test.vec = "coef") + all_equal(coef(est), coef(ref), tolerance = 1e-7, + label = "Tukey regression regression: coefficients") + # model-based covariance matrix + all_equal(vcov(est, "model"), vcov(ref), tolerance = 1e-7, + label = "Tukey regression M-est: model-based cov") + # design-based covariance matrix (test against fixed values) + ref <- matrix(c(43513457.824988, -1421272.808406, -1421272.808406, + 589121.583227), ncol = 2) + colnames(ref) <- rownames(ref) <- c("(Intercept)", "employment") + all_equal(vcov(est), ref, + label = "Tukey regression M-est: design-based cov") + } Error: Huber regression M-est: model-based cov: Mean relative difference: 1 Execution halted * checking for unstated dependencies in vignettes ... OK * checking package vignettes ... OK * checking re-building of vignette outputs ... [48s/194s] OK * checking PDF version of manual ... [13s/40s] OK * checking HTML version of manual ... [14s/57s] OK * checking for non-standard things in the check directory ... OK * checking for detritus in the temp directory ... OK * DONE Status: 1 ERROR See ‘/data/gannet/ripley/R/packages/tests-OpenBLAS/robsurvey.Rcheck/00check.log’ for details. Command exited with non-zero status 1 Time 15:32.11, 224.35 + 25.03