R Under development (unstable) (2024-03-14 r86117 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

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.

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.

> ## BayesX REML testing
> library("BayesXsrc")
> reml <- run.bayesx("reml.prg", verbose = FALSE)
> fx1 <- read.table("reml_f_x1_pspline.res", header = TRUE)
> fx2 <- read.table("reml_f_x2_pspline.res", header = TRUE)
> print(round(head(fx1), digits = 3))
  intnr     x1  pmode ci95lower ci80lower   std ci80upper ci95upper pcat95
1     1 -2.999 -0.261    -0.464    -0.394 0.104    -0.128    -0.058     -1
2     2 -2.985 -0.268    -0.464    -0.396 0.100    -0.141    -0.073     -1
3     3 -2.981 -0.270    -0.464    -0.397 0.099    -0.144    -0.077     -1
4     4 -2.974 -0.274    -0.464    -0.398 0.097    -0.149    -0.084     -1
5     5 -2.963 -0.279    -0.464    -0.400 0.094    -0.159    -0.095     -1
6     6 -2.962 -0.280    -0.464    -0.400 0.094    -0.159    -0.096     -1
  pcat80
1     -1
2     -1
3     -1
4     -1
5     -1
6     -1
> print(round(head(fx2), digits = 3))
  intnr    x2  pmode ci95lower ci80lower   std ci80upper ci95upper pcat95
1     1 0.001 -0.381    -0.457    -0.431 0.039    -0.332    -0.305     -1
2     2 0.001 -0.381    -0.457    -0.431 0.039    -0.332    -0.305     -1
3     3 0.002 -0.381    -0.456    -0.430 0.039    -0.331    -0.305     -1
4     4 0.003 -0.380    -0.456    -0.430 0.039    -0.331    -0.305     -1
5     5 0.003 -0.380    -0.455    -0.429 0.039    -0.331    -0.304     -1
6     6 0.004 -0.379    -0.455    -0.429 0.038    -0.330    -0.304     -1
  pcat80
1     -1
2     -1
3     -1
4     -1
5     -1
6     -1
> 
> proc.time()
   user  system elapsed 
   0.15    0.03    0.51