CRAN Package Check Results for Maintainer ‘Woody Miao <miaolin at pku.edu.cn>’

Last updated on 2025-08-09 05:49:10 CEST.

Package ERROR NOTE OK
GRAB 4 2 7

Package GRAB

Current CRAN status: ERROR: 4, NOTE: 2, OK: 7

Version: 0.2.2
Check: examples
Result: ERROR Running examples in ‘GRAB-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: GRAB.POLMM > ### Title: POLMM method in GRAB package > ### Aliases: GRAB.POLMM > > ### ** Examples > > ### First, Read Data and Convert Phenotype to a Factor > library(dplyr) Attaching package: ‘dplyr’ The following objects are masked from ‘package:stats’: filter, lag The following objects are masked from ‘package:base’: intersect, setdiff, setequal, union > PhenoFile <- system.file("extdata", "simuPHENO.txt", package = "GRAB") > PhenoData <- data.table::fread(PhenoFile, header = TRUE) > PhenoData <- PhenoData %>% mutate(OrdinalPheno = factor(OrdinalPheno, + levels = c(0, 1, 2) + )) > > ### Step 1: Fit a null model > # If a sparse GRM is used in model fitting, SparseGRMFile is required. > # If SparseGRMFile isn't provided, GRAB.NullModel() will calculate dense GRM from GenoFile. > > SparseGRMFile <- system.file("SparseGRM", "SparseGRM.txt", package = "GRAB") > GenoFile <- system.file("extdata", "simuPLINK.bed", package = "GRAB") > obj.POLMM <- GRAB.NullModel( + formula = OrdinalPheno ~ AGE + GENDER, + data = PhenoData, + subjData = PhenoData$IID, + method = "POLMM", + traitType = "ordinal", + GenoFile = GenoFile, + SparseGRMFile = SparseGRMFile, + control = list( + showInfo = FALSE, + LOCO = FALSE, + tolTau = 0.2, + tolBeta = 0.1 + ) + ) Number of subjects in 'formula': 1000 Reading bim file: /home/hornik/tmp/R.check/r-devel-clang/Work/build/Packages/GRAB/extdata/simuPLINK.bim Reading fam file: /home/hornik/tmp/R.check/r-devel-clang/Work/build/Packages/GRAB/extdata/simuPLINK.fam setting PLINK object in CPP.... Reading bim file.... Reading fam file: /home/hornik/tmp/R.check/r-devel-clang/Work/build/Packages/GRAB/extdata/simuPLINK.fam m_N: 1000 Setting position of samples in PLINK files.... posSampleInPlink: 501 510 600 SampleInModel: Subj-1 Subj-10 Subj-100 SampleInPlink: f1_1 f1_2 f1_3 Number of samples: 1000 Based on the 'GenoFile' and 'GenoFileIndex',PLINKformat is used for genotype data. Sparse GRM is used when fitting a null model. The below are the list of control parameters used in null model fitting. 2-1301000.10.20.21001e-061001e-100.10.1200.0025auto10.010.1FALSEFALSEFALSE gIndex: 1016 gIndex: 678 gIndex: 128 gIndex: 929 gIndex: 470 gIndex: 298 gIndex: 269 gIndex: 596 gIndex: 329 gIndex: 36 gIndex: 728 gIndex: 877 gIndex: 484 gIndex: 873 gIndex: 381 gIndex: 974 gIndex: 800 gIndex: 851 gIndex: 930 gIndex: 325 gIndex: 1091 gIndex: 553 gIndex: 110 gIndex: 531 gIndex: 555 gIndex: 888 gIndex: 342 gIndex: 581 gIndex: 39 gIndex: 536 gIndex: 374 gIndex: 247 gIndex: 377 gIndex: 38 gIndex: 434 gIndex: 809 gIndex: 525 gIndex: 641 gIndex: 21 gIndex: 192 gIndex: 498 gIndex: 982 gIndex: 842 gIndex: 464 gIndex: 995 gIndex: 83 gIndex: 1052 gIndex: 251 gIndex: 732 gIndex: 619 gIndex: 556 gIndex: 286 gIndex: 613 gIndex: 328 gIndex: 486 gIndex: 854 gIndex: 850 gIndex: 857 gIndex: 618 gIndex: 575 gIndex: 989 gIndex: 1022 gIndex: 315 gIndex: 1051 gIndex: 649 gIndex: 1097 gIndex: 810 gIndex: 954 gIndex: 281 gIndex: 284 gIndex: 500 gIndex: 510 gIndex: 535 gIndex: 692 gIndex: 987 gIndex: 736 gIndex: 338 gIndex: 345 gIndex: 674 gIndex: 42 gIndex: 0 gIndex: 28 gIndex: 589 gIndex: 917 gIndex: 837 gIndex: 795 gIndex: 627 gIndex: 232 gIndex: 292 gIndex: 572 gIndex: 941 gIndex: 368 gIndex: 450 gIndex: 85 gIndex: 482 gIndex: 326 gIndex: 621 gIndex: 354 gIndex: 818 gIndex: 811 gIndex: 48 gIndex: 360 gIndex: 1037 gIndex: 241 gIndex: 439 gIndex: 757 gIndex: 816 gIndex: 817 gIndex: 246 gIndex: 750 gIndex: 218 gIndex: 134 It took 0.003075 seconds (0.003074 CPU seconds) to calculate 30 genKinbVec(). Start iteration ..... d2eps: -1.3071e+02 d1eps: 2.7770 deps: 0.0212 iter: 0 m_eps: 0 0.9821 eps0: 0 0.9609 m_tolEps: 1e-10 d2eps: -1.2536e+02 d1eps: 0.0574 deps: 0.0005 iter: 1 m_eps: 0 0.9826 eps0: 0 0.9821 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 2.5614e-05 deps: 2.0451e-07 iter: 2 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 5.0837e-12 deps: 4.0590e-14 iter: 3 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9826 beta: -37.5557 0.5975 0.3791 diffBeta: 0.00275036 Start updating tau... iter: 0 beta: -37.5557 0.5975 0.3791 tau: 0.0293943 d2eps: -1.2360e+02 d1eps: -2.2300 deps: -0.0180 iter: 0 m_eps: 0 0.9646 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2811e+02 d1eps: 0.0403 deps: 0.0003 iter: 1 m_eps: 0 0.9649 eps0: 0 0.9646 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2656e-05 deps: 9.8853e-08 iter: 2 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2440e-12 deps: 9.7169e-15 iter: 3 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9649 beta: -37.5479 0.5974 0.3811 diffBeta: 0.00231347 Start updating tau... iter: 1 beta: -37.5479 0.5974 0.3811 tau: 0.0136161 It took 0.070414 seconds (0.070407 CPU seconds) to fit the null POLMM.. Start estimating variance ratio.... nSNPs for CV: 20 CV: 5.84339e-06 Complete the null model fitting in package GRAB: Complete Time: 2025-08-07 17:15:58.789623 > > objPOLMMFile <- system.file("results", "objPOLMMFile.RData", package = "GRAB") > save(obj.POLMM, file = objPOLMMFile) Warning in gzfile(file, "wb") : cannot open compressed file '/home/hornik/tmp/R.check/r-devel-clang/Work/build/Packages/GRAB/results/objPOLMMFile.RData', probable reason 'Read-only file system' Error in gzfile(file, "wb") : cannot open the connection Calls: save -> gzfile Execution halted Flavor: r-devel-linux-x86_64-debian-clang

Version: 0.2.2
Check: examples
Result: ERROR Running examples in ‘GRAB-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: GRAB.POLMM > ### Title: POLMM method in GRAB package > ### Aliases: GRAB.POLMM > > ### ** Examples > > ### First, Read Data and Convert Phenotype to a Factor > library(dplyr) Attaching package: ‘dplyr’ The following objects are masked from ‘package:stats’: filter, lag The following objects are masked from ‘package:base’: intersect, setdiff, setequal, union > PhenoFile <- system.file("extdata", "simuPHENO.txt", package = "GRAB") > PhenoData <- data.table::fread(PhenoFile, header = TRUE) > PhenoData <- PhenoData %>% mutate(OrdinalPheno = factor(OrdinalPheno, + levels = c(0, 1, 2) + )) > > ### Step 1: Fit a null model > # If a sparse GRM is used in model fitting, SparseGRMFile is required. > # If SparseGRMFile isn't provided, GRAB.NullModel() will calculate dense GRM from GenoFile. > > SparseGRMFile <- system.file("SparseGRM", "SparseGRM.txt", package = "GRAB") > GenoFile <- system.file("extdata", "simuPLINK.bed", package = "GRAB") > obj.POLMM <- GRAB.NullModel( + formula = OrdinalPheno ~ AGE + GENDER, + data = PhenoData, + subjData = PhenoData$IID, + method = "POLMM", + traitType = "ordinal", + GenoFile = GenoFile, + SparseGRMFile = SparseGRMFile, + control = list( + showInfo = FALSE, + LOCO = FALSE, + tolTau = 0.2, + tolBeta = 0.1 + ) + ) Number of subjects in 'formula': 1000 Reading bim file: /home/hornik/tmp/R.check/r-devel-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.bim Reading fam file: /home/hornik/tmp/R.check/r-devel-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.fam setting PLINK object in CPP.... Reading bim file.... Reading fam file: /home/hornik/tmp/R.check/r-devel-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.fam m_N: 1000 Setting position of samples in PLINK files.... posSampleInPlink: 501 510 600 SampleInModel: Subj-1 Subj-10 Subj-100 SampleInPlink: f1_1 f1_2 f1_3 Number of samples: 1000 Based on the 'GenoFile' and 'GenoFileIndex',PLINKformat is used for genotype data. Sparse GRM is used when fitting a null model. The below are the list of control parameters used in null model fitting. 2-1301000.10.20.21001e-061001e-100.10.1200.0025auto10.010.1FALSEFALSEFALSE gIndex: 1016 gIndex: 678 gIndex: 128 gIndex: 929 gIndex: 470 gIndex: 298 gIndex: 269 gIndex: 596 gIndex: 329 gIndex: 36 gIndex: 728 gIndex: 877 gIndex: 484 gIndex: 873 gIndex: 381 gIndex: 974 gIndex: 800 gIndex: 851 gIndex: 930 gIndex: 325 gIndex: 1091 gIndex: 553 gIndex: 110 gIndex: 531 gIndex: 555 gIndex: 888 gIndex: 342 gIndex: 581 gIndex: 39 gIndex: 536 gIndex: 374 gIndex: 247 gIndex: 377 gIndex: 38 gIndex: 434 gIndex: 809 gIndex: 525 gIndex: 641 gIndex: 21 gIndex: 192 gIndex: 498 gIndex: 982 gIndex: 842 gIndex: 464 gIndex: 995 gIndex: 83 gIndex: 1052 gIndex: 251 gIndex: 732 gIndex: 619 gIndex: 556 gIndex: 286 gIndex: 613 gIndex: 328 gIndex: 486 gIndex: 854 gIndex: 850 gIndex: 857 gIndex: 618 gIndex: 575 gIndex: 989 gIndex: 1022 gIndex: 315 gIndex: 1051 gIndex: 649 gIndex: 1097 gIndex: 810 gIndex: 954 gIndex: 281 gIndex: 284 gIndex: 500 gIndex: 510 gIndex: 535 gIndex: 692 gIndex: 987 gIndex: 736 gIndex: 338 gIndex: 345 gIndex: 674 gIndex: 42 gIndex: 0 gIndex: 28 gIndex: 589 gIndex: 917 gIndex: 837 gIndex: 795 gIndex: 627 gIndex: 232 gIndex: 292 gIndex: 572 gIndex: 941 gIndex: 368 gIndex: 450 gIndex: 85 gIndex: 482 gIndex: 326 gIndex: 621 gIndex: 354 gIndex: 818 gIndex: 811 gIndex: 48 gIndex: 360 gIndex: 1037 gIndex: 241 gIndex: 439 gIndex: 757 gIndex: 816 gIndex: 817 gIndex: 246 gIndex: 750 gIndex: 218 gIndex: 134 It took 0.002471 seconds (0.002468 CPU seconds) to calculate 30 genKinbVec(). Start iteration ..... d2eps: -1.3071e+02 d1eps: 2.7770 deps: 0.0212 iter: 0 m_eps: 0 0.9821 eps0: 0 0.9609 m_tolEps: 1e-10 d2eps: -1.2536e+02 d1eps: 0.0574 deps: 0.0005 iter: 1 m_eps: 0 0.9826 eps0: 0 0.9821 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 2.5614e-05 deps: 2.0451e-07 iter: 2 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 5.0837e-12 deps: 4.0590e-14 iter: 3 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9826 beta: -37.5557 0.5975 0.3791 diffBeta: 0.00275036 Start updating tau... iter: 0 beta: -37.5557 0.5975 0.3791 tau: 0.0293943 d2eps: -1.2360e+02 d1eps: -2.2300 deps: -0.0180 iter: 0 m_eps: 0 0.9646 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2811e+02 d1eps: 0.0403 deps: 0.0003 iter: 1 m_eps: 0 0.9649 eps0: 0 0.9646 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2656e-05 deps: 9.8853e-08 iter: 2 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2440e-12 deps: 9.7169e-15 iter: 3 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9649 beta: -37.5479 0.5974 0.3811 diffBeta: 0.00231347 Start updating tau... iter: 1 beta: -37.5479 0.5974 0.3811 tau: 0.0136161 It took 0.042828 seconds (0.041628 CPU seconds) to fit the null POLMM.. Start estimating variance ratio.... nSNPs for CV: 20 CV: 5.84339e-06 Complete the null model fitting in package GRAB: Complete Time: 2025-08-08 13:19:11.500123 > > objPOLMMFile <- system.file("results", "objPOLMMFile.RData", package = "GRAB") > save(obj.POLMM, file = objPOLMMFile) Warning in gzfile(file, "wb") : cannot open compressed file '/home/hornik/tmp/R.check/r-devel-gcc/Work/build/Packages/GRAB/results/objPOLMMFile.RData', probable reason 'Read-only file system' Error in gzfile(file, "wb") : cannot open the connection Calls: save -> gzfile Execution halted Flavor: r-devel-linux-x86_64-debian-gcc

Version: 0.2.2
Check: examples
Result: ERROR Running examples in ‘GRAB-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: GRAB.POLMM > ### Title: POLMM method in GRAB package > ### Aliases: GRAB.POLMM > > ### ** Examples > > ### First, Read Data and Convert Phenotype to a Factor > library(dplyr) Attaching package: ‘dplyr’ The following objects are masked from ‘package:stats’: filter, lag The following objects are masked from ‘package:base’: intersect, setdiff, setequal, union > PhenoFile <- system.file("extdata", "simuPHENO.txt", package = "GRAB") > PhenoData <- data.table::fread(PhenoFile, header = TRUE) > PhenoData <- PhenoData %>% mutate(OrdinalPheno = factor(OrdinalPheno, + levels = c(0, 1, 2) + )) > > ### Step 1: Fit a null model > # If a sparse GRM is used in model fitting, SparseGRMFile is required. > # If SparseGRMFile isn't provided, GRAB.NullModel() will calculate dense GRM from GenoFile. > > SparseGRMFile <- system.file("SparseGRM", "SparseGRM.txt", package = "GRAB") > GenoFile <- system.file("extdata", "simuPLINK.bed", package = "GRAB") > obj.POLMM <- GRAB.NullModel( + formula = OrdinalPheno ~ AGE + GENDER, + data = PhenoData, + subjData = PhenoData$IID, + method = "POLMM", + traitType = "ordinal", + GenoFile = GenoFile, + SparseGRMFile = SparseGRMFile, + control = list( + showInfo = FALSE, + LOCO = FALSE, + tolTau = 0.2, + tolBeta = 0.1 + ) + ) Number of subjects in 'formula': 1000 Reading bim file: /home/hornik/tmp/R.check/r-patched-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.bim Reading fam file: /home/hornik/tmp/R.check/r-patched-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.fam setting PLINK object in CPP.... Reading bim file.... Reading fam file: /home/hornik/tmp/R.check/r-patched-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.fam m_N: 1000 Setting position of samples in PLINK files.... posSampleInPlink: 501 510 600 SampleInModel: Subj-1 Subj-10 Subj-100 SampleInPlink: f1_1 f1_2 f1_3 Number of samples: 1000 Based on the 'GenoFile' and 'GenoFileIndex',PLINKformat is used for genotype data. Sparse GRM is used when fitting a null model. The below are the list of control parameters used in null model fitting. 2-1301000.10.20.21001e-061001e-100.10.1200.0025auto10.010.1FALSEFALSEFALSE gIndex: 1016 gIndex: 678 gIndex: 128 gIndex: 929 gIndex: 470 gIndex: 298 gIndex: 269 gIndex: 596 gIndex: 329 gIndex: 36 gIndex: 728 gIndex: 877 gIndex: 484 gIndex: 873 gIndex: 381 gIndex: 974 gIndex: 800 gIndex: 851 gIndex: 930 gIndex: 325 gIndex: 1091 gIndex: 553 gIndex: 110 gIndex: 531 gIndex: 555 gIndex: 888 gIndex: 342 gIndex: 581 gIndex: 39 gIndex: 536 gIndex: 374 gIndex: 247 gIndex: 377 gIndex: 38 gIndex: 434 gIndex: 809 gIndex: 525 gIndex: 641 gIndex: 21 gIndex: 192 gIndex: 498 gIndex: 982 gIndex: 842 gIndex: 464 gIndex: 995 gIndex: 83 gIndex: 1052 gIndex: 251 gIndex: 732 gIndex: 619 gIndex: 556 gIndex: 286 gIndex: 613 gIndex: 328 gIndex: 486 gIndex: 854 gIndex: 850 gIndex: 857 gIndex: 618 gIndex: 575 gIndex: 989 gIndex: 1022 gIndex: 315 gIndex: 1051 gIndex: 649 gIndex: 1097 gIndex: 810 gIndex: 954 gIndex: 281 gIndex: 284 gIndex: 500 gIndex: 510 gIndex: 535 gIndex: 692 gIndex: 987 gIndex: 736 gIndex: 338 gIndex: 345 gIndex: 674 gIndex: 42 gIndex: 0 gIndex: 28 gIndex: 589 gIndex: 917 gIndex: 837 gIndex: 795 gIndex: 627 gIndex: 232 gIndex: 292 gIndex: 572 gIndex: 941 gIndex: 368 gIndex: 450 gIndex: 85 gIndex: 482 gIndex: 326 gIndex: 621 gIndex: 354 gIndex: 818 gIndex: 811 gIndex: 48 gIndex: 360 gIndex: 1037 gIndex: 241 gIndex: 439 gIndex: 757 gIndex: 816 gIndex: 817 gIndex: 246 gIndex: 750 gIndex: 218 gIndex: 134 It took 0.003503 seconds (0.003502 CPU seconds) to calculate 30 genKinbVec(). Start iteration ..... d2eps: -1.3071e+02 d1eps: 2.7770 deps: 0.0212 iter: 0 m_eps: 0 0.9821 eps0: 0 0.9609 m_tolEps: 1e-10 d2eps: -1.2536e+02 d1eps: 0.0574 deps: 0.0005 iter: 1 m_eps: 0 0.9826 eps0: 0 0.9821 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 2.5614e-05 deps: 2.0451e-07 iter: 2 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 5.0837e-12 deps: 4.0590e-14 iter: 3 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9826 beta: -37.5557 0.5975 0.3791 diffBeta: 0.00275036 Start updating tau... iter: 0 beta: -37.5557 0.5975 0.3791 tau: 0.0293943 d2eps: -1.2360e+02 d1eps: -2.2300 deps: -0.0180 iter: 0 m_eps: 0 0.9646 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2811e+02 d1eps: 0.0403 deps: 0.0003 iter: 1 m_eps: 0 0.9649 eps0: 0 0.9646 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2656e-05 deps: 9.8853e-08 iter: 2 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2440e-12 deps: 9.7169e-15 iter: 3 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9649 beta: -37.5479 0.5974 0.3811 diffBeta: 0.00231347 Start updating tau... iter: 1 beta: -37.5479 0.5974 0.3811 tau: 0.0136161 It took 0.106994 seconds (0.067012 CPU seconds) to fit the null POLMM.. Start estimating variance ratio.... nSNPs for CV: 20 CV: 5.84339e-06 Complete the null model fitting in package GRAB: Complete Time: 2025-08-06 05:31:40.248063 > > objPOLMMFile <- system.file("results", "objPOLMMFile.RData", package = "GRAB") > save(obj.POLMM, file = objPOLMMFile) Warning in gzfile(file, "wb") : cannot open compressed file '/home/hornik/tmp/R.check/r-patched-gcc/Work/build/Packages/GRAB/results/objPOLMMFile.RData', probable reason 'Read-only file system' Error in gzfile(file, "wb") : cannot open the connection Calls: save -> gzfile Execution halted Flavor: r-patched-linux-x86_64

Version: 0.2.2
Check: examples
Result: ERROR Running examples in ‘GRAB-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: GRAB.POLMM > ### Title: POLMM method in GRAB package > ### Aliases: GRAB.POLMM > > ### ** Examples > > ### First, Read Data and Convert Phenotype to a Factor > library(dplyr) Attaching package: ‘dplyr’ The following objects are masked from ‘package:stats’: filter, lag The following objects are masked from ‘package:base’: intersect, setdiff, setequal, union > PhenoFile <- system.file("extdata", "simuPHENO.txt", package = "GRAB") > PhenoData <- data.table::fread(PhenoFile, header = TRUE) > PhenoData <- PhenoData %>% mutate(OrdinalPheno = factor(OrdinalPheno, + levels = c(0, 1, 2) + )) > > ### Step 1: Fit a null model > # If a sparse GRM is used in model fitting, SparseGRMFile is required. > # If SparseGRMFile isn't provided, GRAB.NullModel() will calculate dense GRM from GenoFile. > > SparseGRMFile <- system.file("SparseGRM", "SparseGRM.txt", package = "GRAB") > GenoFile <- system.file("extdata", "simuPLINK.bed", package = "GRAB") > obj.POLMM <- GRAB.NullModel( + formula = OrdinalPheno ~ AGE + GENDER, + data = PhenoData, + subjData = PhenoData$IID, + method = "POLMM", + traitType = "ordinal", + GenoFile = GenoFile, + SparseGRMFile = SparseGRMFile, + control = list( + showInfo = FALSE, + LOCO = FALSE, + tolTau = 0.2, + tolBeta = 0.1 + ) + ) Number of subjects in 'formula': 1000 Reading bim file: /home/hornik/tmp/R.check/r-release-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.bim Reading fam file: /home/hornik/tmp/R.check/r-release-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.fam setting PLINK object in CPP.... Reading bim file.... Reading fam file: /home/hornik/tmp/R.check/r-release-gcc/Work/build/Packages/GRAB/extdata/simuPLINK.fam m_N: 1000 Setting position of samples in PLINK files.... posSampleInPlink: 501 510 600 SampleInModel: Subj-1 Subj-10 Subj-100 SampleInPlink: f1_1 f1_2 f1_3 Number of samples: 1000 Based on the 'GenoFile' and 'GenoFileIndex',PLINKformat is used for genotype data. Sparse GRM is used when fitting a null model. The below are the list of control parameters used in null model fitting. 2-1301000.10.20.21001e-061001e-100.10.1200.0025auto10.010.1FALSEFALSEFALSE gIndex: 1016 gIndex: 678 gIndex: 128 gIndex: 929 gIndex: 470 gIndex: 298 gIndex: 269 gIndex: 596 gIndex: 329 gIndex: 36 gIndex: 728 gIndex: 877 gIndex: 484 gIndex: 873 gIndex: 381 gIndex: 974 gIndex: 800 gIndex: 851 gIndex: 930 gIndex: 325 gIndex: 1091 gIndex: 553 gIndex: 110 gIndex: 531 gIndex: 555 gIndex: 888 gIndex: 342 gIndex: 581 gIndex: 39 gIndex: 536 gIndex: 374 gIndex: 247 gIndex: 377 gIndex: 38 gIndex: 434 gIndex: 809 gIndex: 525 gIndex: 641 gIndex: 21 gIndex: 192 gIndex: 498 gIndex: 982 gIndex: 842 gIndex: 464 gIndex: 995 gIndex: 83 gIndex: 1052 gIndex: 251 gIndex: 732 gIndex: 619 gIndex: 556 gIndex: 286 gIndex: 613 gIndex: 328 gIndex: 486 gIndex: 854 gIndex: 850 gIndex: 857 gIndex: 618 gIndex: 575 gIndex: 989 gIndex: 1022 gIndex: 315 gIndex: 1051 gIndex: 649 gIndex: 1097 gIndex: 810 gIndex: 954 gIndex: 281 gIndex: 284 gIndex: 500 gIndex: 510 gIndex: 535 gIndex: 692 gIndex: 987 gIndex: 736 gIndex: 338 gIndex: 345 gIndex: 674 gIndex: 42 gIndex: 0 gIndex: 28 gIndex: 589 gIndex: 917 gIndex: 837 gIndex: 795 gIndex: 627 gIndex: 232 gIndex: 292 gIndex: 572 gIndex: 941 gIndex: 368 gIndex: 450 gIndex: 85 gIndex: 482 gIndex: 326 gIndex: 621 gIndex: 354 gIndex: 818 gIndex: 811 gIndex: 48 gIndex: 360 gIndex: 1037 gIndex: 241 gIndex: 439 gIndex: 757 gIndex: 816 gIndex: 817 gIndex: 246 gIndex: 750 gIndex: 218 gIndex: 134 It took 0.005931 seconds (0.003393 CPU seconds) to calculate 30 genKinbVec(). Start iteration ..... d2eps: -1.3071e+02 d1eps: 2.7770 deps: 0.0212 iter: 0 m_eps: 0 0.9821 eps0: 0 0.9609 m_tolEps: 1e-10 d2eps: -1.2536e+02 d1eps: 0.0574 deps: 0.0005 iter: 1 m_eps: 0 0.9826 eps0: 0 0.9821 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 2.5614e-05 deps: 2.0451e-07 iter: 2 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2525e+02 d1eps: 5.0837e-12 deps: 4.0590e-14 iter: 3 m_eps: 0 0.9826 eps0: 0 0.9826 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9826 beta: -37.5557 0.5975 0.3791 diffBeta: 0.00275036 Start updating tau... iter: 0 beta: -37.5557 0.5975 0.3791 tau: 0.0293943 d2eps: -1.2360e+02 d1eps: -2.2300 deps: -0.0180 iter: 0 m_eps: 0 0.9646 eps0: 0 0.9826 m_tolEps: 1e-10 d2eps: -1.2811e+02 d1eps: 0.0403 deps: 0.0003 iter: 1 m_eps: 0 0.9649 eps0: 0 0.9646 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2656e-05 deps: 9.8853e-08 iter: 2 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 d2eps: -1.2803e+02 d1eps: 1.2440e-12 deps: 9.7169e-15 iter: 3 m_eps: 0 0.9649 eps0: 0 0.9649 m_tolEps: 1e-10 UpdateEps iter: 3 eps: 0 0.9649 beta: -37.5479 0.5974 0.3811 diffBeta: 0.00231347 Start updating tau... iter: 1 beta: -37.5479 0.5974 0.3811 tau: 0.0136161 It took 0.073767 seconds (0.073541 CPU seconds) to fit the null POLMM.. Start estimating variance ratio.... nSNPs for CV: 20 CV: 5.84339e-06 Complete the null model fitting in package GRAB: Complete Time: 2025-08-02 05:17:44.613291 > > objPOLMMFile <- system.file("results", "objPOLMMFile.RData", package = "GRAB") > save(obj.POLMM, file = objPOLMMFile) Warning in gzfile(file, "wb") : cannot open compressed file '/home/hornik/tmp/R.check/r-release-gcc/Work/build/Packages/GRAB/results/objPOLMMFile.RData', probable reason 'Read-only file system' Error in gzfile(file, "wb") : cannot open the connection Calls: save -> gzfile Execution halted Flavor: r-release-linux-x86_64

Version: 0.2.2
Check: installed package size
Result: NOTE installed size is 17.7Mb sub-directories of 1Mb or more: extdata 2.1Mb libs 14.9Mb Flavors: r-oldrel-macos-arm64, r-oldrel-macos-x86_64

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.