| Back to Multiple platform build/check report for BioC 3.22: simplified long |
|
This page was generated on 2025-08-15 12:06 -0400 (Fri, 15 Aug 2025).
| Hostname | OS | Arch (*) | R version | Installed pkgs |
|---|---|---|---|---|
| nebbiolo2 | Linux (Ubuntu 24.04.3 LTS) | x86_64 | 4.5.1 (2025-06-13) -- "Great Square Root" | 4818 |
| palomino8 | Windows Server 2022 Datacenter | x64 | 4.5.1 (2025-06-13 ucrt) -- "Great Square Root" | 4554 |
| lconway | macOS 12.7.1 Monterey | x86_64 | 4.5.1 (2025-06-13) -- "Great Square Root" | 4595 |
| kjohnson3 | macOS 13.7.7 Ventura | arm64 | 4.5.1 Patched (2025-06-14 r88325) -- "Great Square Root" | 4537 |
| taishan | Linux (openEuler 24.03 LTS) | aarch64 | 4.5.0 (2025-04-11) -- "How About a Twenty-Six" | 4535 |
| Click on any hostname to see more info about the system (e.g. compilers) (*) as reported by 'uname -p', except on Windows and Mac OS X | ||||
| Package 575/2317 | Hostname | OS / Arch | INSTALL | BUILD | CHECK | BUILD BIN | ||||||||
| diffHic 1.41.0 (landing page) Aaron Lun
| nebbiolo2 | Linux (Ubuntu 24.04.3 LTS) / x86_64 | OK | OK | OK | |||||||||
| palomino8 | Windows Server 2022 Datacenter / x64 | OK | OK | OK | OK | |||||||||
| lconway | macOS 12.7.1 Monterey / x86_64 | OK | OK | OK | OK | |||||||||
| kjohnson3 | macOS 13.7.7 Ventura / arm64 | OK | OK | OK | OK | |||||||||
| taishan | Linux (openEuler 24.03 LTS) / aarch64 | OK | OK | OK | ||||||||||
|
To the developers/maintainers of the diffHic package: - Allow up to 24 hours (and sometimes 48 hours) for your latest push to git@git.bioconductor.org:packages/diffHic.git to reflect on this report. See Troubleshooting Build Report for more information. - Use the following Renviron settings to reproduce errors and warnings. - If 'R CMD check' started to fail recently on the Linux builder(s) over a missing dependency, add the missing dependency to 'Suggests:' in your DESCRIPTION file. See Renviron.bioc for more information. |
| Package: diffHic |
| Version: 1.41.0 |
| Command: F:\biocbuild\bbs-3.22-bioc\R\bin\R.exe CMD check --no-multiarch --install=check:diffHic.install-out.txt --library=F:\biocbuild\bbs-3.22-bioc\R\library --no-vignettes --timings diffHic_1.41.0.tar.gz |
| StartedAt: 2025-08-15 02:37:14 -0400 (Fri, 15 Aug 2025) |
| EndedAt: 2025-08-15 02:46:06 -0400 (Fri, 15 Aug 2025) |
| EllapsedTime: 531.2 seconds |
| RetCode: 0 |
| Status: OK |
| CheckDir: diffHic.Rcheck |
| Warnings: 0 |
##############################################################################
##############################################################################
###
### Running command:
###
### F:\biocbuild\bbs-3.22-bioc\R\bin\R.exe CMD check --no-multiarch --install=check:diffHic.install-out.txt --library=F:\biocbuild\bbs-3.22-bioc\R\library --no-vignettes --timings diffHic_1.41.0.tar.gz
###
##############################################################################
##############################################################################
* using log directory 'F:/biocbuild/bbs-3.22-bioc/meat/diffHic.Rcheck'
* using R version 4.5.1 (2025-06-13 ucrt)
* using platform: x86_64-w64-mingw32
* R was compiled by
gcc.exe (GCC) 14.2.0
GNU Fortran (GCC) 14.2.0
* running under: Windows Server 2022 x64 (build 20348)
* using session charset: UTF-8
* using option '--no-vignettes'
* checking for file 'diffHic/DESCRIPTION' ... OK
* this is package 'diffHic' version '1.41.0'
* 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 hidden files and directories ... OK
* checking for portable file names ... OK
* checking whether package 'diffHic' can be installed ... OK
* used C++ compiler: 'G__~1.EXE (GCC) 14.2.0'
* checking installed package size ... INFO
installed size is 14.9Mb
sub-directories of 1Mb or more:
doc 1.3Mb
libs 12.4Mb
* 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 ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... 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 ... OK
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... NOTE
Found the following Rd file(s) with Rd \link{} targets missing package
anchors:
annotatePairs.Rd: findOverlaps
connectCounts.Rd: findOverlaps, linkOverlaps
consolidatePairs.Rd: combineTests
correctedContact.Rd: mglmOneGroup
cutGenome.Rd: BSgenome-class, GRanges-class, matchPattern
diClusters.Rd: controlClusterFDR
filterDiag.Rd: pairdist
filterPeaks.Rd: aveLogCPM
filters.Rd: loessFit, scaledAverage
normalizeCNV.Rd: locfit, lp
readMTX2IntSet.Rd: InteractionSet-class, GInteractions
Please provide package anchors for all Rd \link{} targets not in the
package itself and the base packages.
* 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 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 ... INFO
GNU make is a SystemRequirements.
* checking for portable use of $(BLAS_LIBS) and $(LAPACK_LIBS) ... OK
* checking use of PKG_*FLAGS in Makefiles ... OK
* checking compiled code ... NOTE
Note: information on .o files for x64 is not available
File 'F:/biocbuild/bbs-3.22-bioc/R/library/diffHic/libs/x64/diffHic.dll':
Found '_assert', possibly from 'assert' (C)
Found '_exit', possibly from '_exit' (C)
Found 'abort', possibly from 'abort' (C), 'runtime' (Fortran)
Found 'exit', possibly from 'exit' (C), 'stop' (Fortran)
Compiled code should not call entry points which might terminate R nor
write to stdout/stderr instead of to the console, nor use Fortran I/O
nor system RNGs nor [v]sprintf. The detected symbols are linked into
the code but might come from libraries and not actually be called.
See 'Writing portable packages' in the 'Writing R Extensions' manual.
* checking sizes of PDF files under 'inst/doc' ... OK
* checking installed files from 'inst/doc' ... OK
* checking files in 'vignettes' ... OK
* checking examples ... OK
Examples with CPU (user + system) or elapsed time > 5s
user system elapsed
normalizeCNV 9.14 0.14 9.28
cutGenome 7.23 0.17 7.78
* checking for unstated dependencies in 'tests' ... OK
* checking tests ...
Running 'test-basic.R'
Comparing 'test-basic.Rout' to 'test-basic.Rout.save' ... OK
Running 'testthat.R'
OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE
Status: 2 NOTEs
See
'F:/biocbuild/bbs-3.22-bioc/meat/diffHic.Rcheck/00check.log'
for details.
diffHic.Rcheck/00install.out
##############################################################################
##############################################################################
###
### Running command:
###
### F:\biocbuild\bbs-3.22-bioc\R\bin\R.exe CMD INSTALL diffHic
###
##############################################################################
##############################################################################
* installing to library 'F:/biocbuild/bbs-3.22-bioc/R/library'
* installing *source* package 'diffHic' ...
** this is package 'diffHic' version '1.41.0'
** using staged installation
** libs
using C++ compiler: 'G__~1.EXE (GCC) 14.2.0'
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c binner.cpp -o binner.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c check_input.cpp -o check_input.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c cluster_2d.cpp -o cluster_2d.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c count_background.cpp -o count_background.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c count_connect.cpp -o count_connect.o
count_connect.cpp: In function 'SEXPREC* count_connect(SEXP, SEXP)':
count_connect.cpp:57:30: warning: comparison of integer expressions of different signedness: 'size_t' {aka 'long long unsigned int'} and 'const int' [-Wsign-compare]
57 | for (size_t odex=0; odex < ncombos; ++odex) {
| ~~~~~^~~~~~~~~
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c count_patch.cpp -o count_patch.o
count_patch.cpp: In function 'SEXPREC* count_patch(SEXP, SEXP, SEXP, SEXP, SEXP)':
count_patch.cpp:45:33: warning: comparison of integer expressions of different signedness: 'size_t' {aka 'long long unsigned int'} and 'const int' [-Wsign-compare]
45 | for (size_t vecdex=0; vecdex<ncombos; ++vecdex) {
| ~~~~~~^~~~~~~~
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c directionality.cpp -o directionality.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c init.cpp -o init.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c iterative_correction.cpp -o iterative_correction.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c neighbors.cpp -o neighbors.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c pair_stats.cpp -o pair_stats.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c quadrant_bg.cpp -o quadrant_bg.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c report_hic_pairs.cpp -o report_hic_pairs.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c trended_filter.cpp -o trended_filter.o
g++ -std=gnu++17 -I"F:/biocbuild/bbs-3.22-bioc/R/include" -DNDEBUG -D_FILE_OFFSET_BITS=64 -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/include' -I'F:/biocbuild/bbs-3.22-bioc/R/library/Rcpp/include' -I"C:/rtools45/x86_64-w64-mingw32.static.posix/include" -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c utils.cpp -o utils.o
g++ -std=gnu++17 -shared -s -static-libgcc -o diffHic.dll tmp.def binner.o check_input.o cluster_2d.o count_background.o count_connect.o count_patch.o directionality.o init.o iterative_correction.o neighbors.o pair_stats.o quadrant_bg.o report_hic_pairs.o trended_filter.o utils.o F:/biocbuild/bbs-3.22-bioc/R/library/Rhtslib/usrlib/x64/libhts.a -lm -lbz2 -llzma -lcurl -lnghttp2 -lpsl -lbrotlidec -lbrotlicommon -lbcrypt -lidn2 -lunistring -liconv -lssl -lcrypto -lz -lcrypt32 -lwsock32 -lwldap32 -lssh2 -lgcrypt -lgpg-error -lws2_32 -lzstd -lregex -LC:/rtools45/x86_64-w64-mingw32.static.posix/lib/x64 -LC:/rtools45/x86_64-w64-mingw32.static.posix/lib -LF:/biocbuild/bbs-3.22-bioc/R/bin/x64 -lR
installing to F:/biocbuild/bbs-3.22-bioc/R/library/00LOCK-diffHic/00new/diffHic/libs/x64
** R
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** installing vignettes
** testing if installed package can be loaded from temporary location
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (diffHic)
diffHic.Rcheck/tests/test-basic.Rout
R version 4.5.1 (2025-06-13 ucrt) -- "Great Square Root"
Copyright (C) 2025 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.
> # These are just placeholders for the real things in inst/tests.
>
> suppressWarnings(suppressPackageStartupMessages(require(diffHic)))
>
> hic.file <- system.file("exdata", "hic_sort.bam", package="diffHic")
> cuts <- readRDS(system.file("exdata", "cuts.rds", package="diffHic"))
> param <- pairParam(fragments=cuts)
>
> # Setting up the parameters
> fout <- "output.h5"
> preparePairs(hic.file, param, file=fout)
$pairs
total marked filtered mapped
32 7 3 22
$same.id
dangling self.circle
4 1
$singles
[1] 2
$chimeras
total mapped multi invalid
12 8 7 5
> head(getPairData(fout, param))
length orientation insert
1 40 1 40
2 80 1 80
3 60 2 54
4 40 3 160
5 20 1 66
6 80 0 114
>
> loadChromos(fout)
anchor1 anchor2
1 chrA chrA
2 chrB chrA
3 chrB chrB
> head(loadData(fout, "chrA", "chrA"))
anchor1.id anchor2.id anchor1.pos anchor2.pos anchor1.len anchor2.len
1 2 1 49 14 -5 10
2 2 1 79 9 -10 10
3 2 1 65 21 10 -10
4 4 1 156 6 -10 -10
5 4 2 141 80 -5 10
6 4 2 154 50 10 10
> head(loadData(fout, "chrA", "chrB"))
anchor1.id anchor2.id anchor1.pos anchor2.pos anchor1.len anchor2.len
1 5 2 19 44 10 -10
2 5 3 14 105 10 -10
3 6 1 65 19 5 10
4 6 2 24 90 -5 5
5 6 3 24 100 -5 -10
6 6 3 24 95 -5 -10
Warning message:
In value[[3L]](cond) : anchor definitions are reversed
>
> # Loading the counts.
> data <- squareCounts(fout, param, width=50, filter=1)
> data
class: InteractionSet
dim: 10 1
metadata(2): param width
assays(1): counts
rownames: NULL
rowData names(0):
colnames: NULL
colData names(1): totals
type: ReverseStrictGInteractions
regions: 6
>
> margins <- marginCounts(fout, param, width=50)
> margins
class: RangedSummarizedExperiment
dim: 6 1
metadata(1): param
assays(1): counts
rownames: NULL
rowData names(1): nfrags
colnames: NULL
colData names(1): totals
> totalCounts(fout, param)
[1] 17
>
> regions <- GRanges("chrA", IRanges(c(1, 100, 150), c(20, 140, 160)))
> connectCounts(fout, param, regions=regions, filter=1L)
class: InteractionSet
dim: 2 1
metadata(1): param
assays(1): counts
rownames: NULL
rowData names(0):
colnames: NULL
colData names(1): totals
type: ReverseStrictGInteractions
regions: 3
>
> # Checking some values.
> head(getArea(data))
[1] 2208 2304 2208 2208 3312 3174
> head(pairdist(data))
[1] 47 140 93 47 NA NA
>
> anchors(data, type="first")
GRanges object with 10 ranges and 1 metadata column:
seqnames ranges strand | nfrags
<Rle> <IRanges> <Rle> | <integer>
[1] chrA 49-94 * | 1
[2] chrA 141-188 * | 1
[3] chrA 141-188 * | 1
[4] chrA 141-188 * | 1
[5] chrB 1-69 * | 2
[6] chrB 1-69 * | 2
[7] chrB 1-69 * | 2
[8] chrB 1-69 * | 2
[9] chrB 70-92 * | 1
[10] chrB 70-92 * | 1
-------
seqinfo: 2 sequences from an unspecified genome
> anchors(data, type="second")
GRanges object with 10 ranges and 1 metadata column:
seqnames ranges strand | nfrags
<Rle> <IRanges> <Rle> | <integer>
[1] chrA 1-48 * | 1
[2] chrA 1-48 * | 1
[3] chrA 49-94 * | 1
[4] chrA 95-140 * | 1
[5] chrA 1-48 * | 1
[6] chrA 49-94 * | 1
[7] chrA 95-140 * | 1
[8] chrB 1-69 * | 2
[9] chrB 1-69 * | 2
[10] chrB 70-92 * | 1
-------
seqinfo: 2 sequences from an unspecified genome
> assay(data)
[,1]
[1,] 3
[2,] 1
[3,] 2
[4,] 1
[5,] 1
[6,] 2
[7,] 3
[8,] 1
[9,] 2
[10,] 1
> regions(data)
GRanges object with 6 ranges and 1 metadata column:
seqnames ranges strand | nfrags
<Rle> <IRanges> <Rle> | <integer>
[1] chrA 1-48 * | 1
[2] chrA 49-94 * | 1
[3] chrA 95-140 * | 1
[4] chrA 141-188 * | 1
[5] chrB 1-69 * | 2
[6] chrB 70-92 * | 1
-------
seqinfo: 2 sequences from an unspecified genome
>
> data$totals
[1] 17
> colData(data)
DataFrame with 1 row and 1 column
totals
<integer>
1 17
> metadata(data)
$param
Genome contains 7 restriction fragments across 2 chromosomes
No discard regions are specified
No limits on chromosomes for read extraction
No cap on the read pairs per pair of restriction fragments
$width
[1] 50
>
> asDGEList(data)
An object of class "DGEList"
$counts
Sample1
1 3
2 1
3 2
4 1
5 1
6 2
7 3
8 1
9 2
10 1
$samples
group lib.size norm.factors
Sample1 1 17 1
> asDGEList(data, lib.size=20)$samples
group lib.size norm.factors
Sample1 1 20 1
> asDGEList(data, norm.factors=2, group="a")$samples
group lib.size norm.factors
Sample1 a 17 2
Warning message:
In DGEList.default(group = "a", lib.size = 17L, norm.factors = 2, :
norm factors don't multiply to 1
>
> # Playing around with some bin counts.
> stuff <- correctedContact(data)
> head(stuff$truth)
[1] 1.517032e+04 6.123523e-01 6.123523e-01 8.162749e-01 2.041633e-01
[6] 2.041633e-01
>
> data.large <- squareCounts(fout, param, width=100, filter=1)
> boxed <- boxPairs(larger=data.large, smaller=data)
> head(boxed$indices$larger)
[1] 1 2 3 4 5 6
> head(boxed$indices$smaller)
[1] 1 2 2 3 4 4
>
> head(enrichedPairs(data))
class: InteractionSet
dim: 6 1
metadata(2): param width
assays(5): counts quadrant vertical horizontal surrounding
rownames: NULL
rowData names(4): N.quadrant N.vertical N.horizontal N.surrounding
colnames: NULL
colData names(1): totals
type: ReverseStrictGInteractions
regions: 6
> head(clusterPairs(data, tol=10)$indices[[1]])
[1] 1 2 2 2 3 3
>
> # End.
>
> unlink(fout)
>
> proc.time()
user system elapsed
15.34 1.07 16.42
diffHic.Rcheck/tests/test-basic.Rout.save
R version 4.4.0 (2024-04-24 ucrt) -- "Puppy Cup"
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.
> # These are just placeholders for the real things in inst/tests.
>
> suppressWarnings(suppressPackageStartupMessages(require(diffHic)))
>
> hic.file <- system.file("exdata", "hic_sort.bam", package="diffHic")
> cuts <- readRDS(system.file("exdata", "cuts.rds", package="diffHic"))
> param <- pairParam(fragments=cuts)
>
> # Setting up the parameters
> fout <- "output.h5"
> preparePairs(hic.file, param, file=fout)
$pairs
total marked filtered mapped
32 7 3 22
$same.id
dangling self.circle
4 1
$singles
[1] 2
$chimeras
total mapped multi invalid
12 8 7 5
> head(getPairData(fout, param))
length orientation insert
1 40 1 40
2 80 1 80
3 60 2 54
4 40 3 160
5 20 1 66
6 80 0 114
>
> loadChromos(fout)
anchor1 anchor2
1 chrA chrA
2 chrB chrA
3 chrB chrB
> head(loadData(fout, "chrA", "chrA"))
anchor1.id anchor2.id anchor1.pos anchor2.pos anchor1.len anchor2.len
1 2 1 49 14 -5 10
2 2 1 79 9 -10 10
3 2 1 65 21 10 -10
4 4 1 156 6 -10 -10
5 4 2 141 80 -5 10
6 4 2 154 50 10 10
> head(loadData(fout, "chrA", "chrB"))
anchor1.id anchor2.id anchor1.pos anchor2.pos anchor1.len anchor2.len
1 5 2 19 44 10 -10
2 5 3 14 105 10 -10
3 6 1 65 19 5 10
4 6 2 24 90 -5 5
5 6 3 24 100 -5 -10
6 6 3 24 95 -5 -10
Warning message:
In value[[3L]](cond) : anchor definitions are reversed
>
> # Loading the counts.
> data <- squareCounts(fout, param, width=50, filter=1)
> data
class: InteractionSet
dim: 10 1
metadata(2): param width
assays(1): counts
rownames: NULL
rowData names(0):
colnames: NULL
colData names(1): totals
type: ReverseStrictGInteractions
regions: 6
>
> margins <- marginCounts(fout, param, width=50)
> margins
class: RangedSummarizedExperiment
dim: 6 1
metadata(1): param
assays(1): counts
rownames: NULL
rowData names(1): nfrags
colnames: NULL
colData names(1): totals
> totalCounts(fout, param)
[1] 17
>
> regions <- GRanges("chrA", IRanges(c(1, 100, 150), c(20, 140, 160)))
> connectCounts(fout, param, regions=regions, filter=1L)
class: InteractionSet
dim: 2 1
metadata(1): param
assays(1): counts
rownames: NULL
rowData names(0):
colnames: NULL
colData names(1): totals
type: ReverseStrictGInteractions
regions: 3
>
> # Checking some values.
> head(getArea(data))
[1] 2208 2304 2208 2208 3312 3174
> head(pairdist(data))
[1] 47 140 93 47 NA NA
>
> anchors(data, type="first")
GRanges object with 10 ranges and 1 metadata column:
seqnames ranges strand | nfrags
<Rle> <IRanges> <Rle> | <integer>
[1] chrA 49-94 * | 1
[2] chrA 141-188 * | 1
[3] chrA 141-188 * | 1
[4] chrA 141-188 * | 1
[5] chrB 1-69 * | 2
[6] chrB 1-69 * | 2
[7] chrB 1-69 * | 2
[8] chrB 1-69 * | 2
[9] chrB 70-92 * | 1
[10] chrB 70-92 * | 1
-------
seqinfo: 2 sequences from an unspecified genome
> anchors(data, type="second")
GRanges object with 10 ranges and 1 metadata column:
seqnames ranges strand | nfrags
<Rle> <IRanges> <Rle> | <integer>
[1] chrA 1-48 * | 1
[2] chrA 1-48 * | 1
[3] chrA 49-94 * | 1
[4] chrA 95-140 * | 1
[5] chrA 1-48 * | 1
[6] chrA 49-94 * | 1
[7] chrA 95-140 * | 1
[8] chrB 1-69 * | 2
[9] chrB 1-69 * | 2
[10] chrB 70-92 * | 1
-------
seqinfo: 2 sequences from an unspecified genome
> assay(data)
[,1]
[1,] 3
[2,] 1
[3,] 2
[4,] 1
[5,] 1
[6,] 2
[7,] 3
[8,] 1
[9,] 2
[10,] 1
> regions(data)
GRanges object with 6 ranges and 1 metadata column:
seqnames ranges strand | nfrags
<Rle> <IRanges> <Rle> | <integer>
[1] chrA 1-48 * | 1
[2] chrA 49-94 * | 1
[3] chrA 95-140 * | 1
[4] chrA 141-188 * | 1
[5] chrB 1-69 * | 2
[6] chrB 70-92 * | 1
-------
seqinfo: 2 sequences from an unspecified genome
>
> data$totals
[1] 17
> colData(data)
DataFrame with 1 row and 1 column
totals
<integer>
1 17
> metadata(data)
$param
Genome contains 7 restriction fragments across 2 chromosomes
No discard regions are specified
No limits on chromosomes for read extraction
No cap on the read pairs per pair of restriction fragments
$width
[1] 50
>
> asDGEList(data)
An object of class "DGEList"
$counts
Sample1
1 3
2 1
3 2
4 1
5 1
6 2
7 3
8 1
9 2
10 1
$samples
group lib.size norm.factors
Sample1 1 17 1
> asDGEList(data, lib.size=20)$samples
group lib.size norm.factors
Sample1 1 20 1
> asDGEList(data, norm.factors=2, group="a")$samples
group lib.size norm.factors
Sample1 a 17 2
Warning message:
In DGEList.default(group = "a", lib.size = 17L, norm.factors = 2, :
norm factors don't multiply to 1
>
> # Playing around with some bin counts.
> stuff <- correctedContact(data)
> head(stuff$truth)
[1] 1.517032e+04 6.123523e-01 6.123523e-01 8.162749e-01 2.041633e-01
[6] 2.041633e-01
>
> data.large <- squareCounts(fout, param, width=100, filter=1)
> boxed <- boxPairs(larger=data.large, smaller=data)
> head(boxed$indices$larger)
[1] 1 2 3 4 5 6
> head(boxed$indices$smaller)
[1] 1 2 2 3 4 4
>
> head(enrichedPairs(data))
class: InteractionSet
dim: 6 1
metadata(2): param width
assays(5): counts quadrant vertical horizontal surrounding
rownames: NULL
rowData names(4): N.quadrant N.vertical N.horizontal N.surrounding
colnames: NULL
colData names(1): totals
type: ReverseStrictGInteractions
regions: 6
> head(clusterPairs(data, tol=10)$indices[[1]])
[1] 1 2 2 2 3 3
>
> # End.
>
> unlink(fout)
>
> proc.time()
user system elapsed
4.00 0.42 6.53
diffHic.Rcheck/tests/testthat.Rout
R version 4.5.1 (2025-06-13 ucrt) -- "Great Square Root"
Copyright (C) 2025 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.
> library(testthat)
> library(diffHic)
Loading required package: GenomicRanges
Loading required package: stats4
Loading required package: BiocGenerics
Loading required package: generics
Attaching package: 'generics'
The following objects are masked from 'package:base':
as.difftime, as.factor, as.ordered, intersect, is.element, setdiff,
setequal, union
Attaching package: 'BiocGenerics'
The following objects are masked from 'package:stats':
IQR, mad, sd, var, xtabs
The following objects are masked from 'package:base':
Filter, Find, Map, Position, Reduce, anyDuplicated, aperm, append,
as.data.frame, basename, cbind, colnames, dirname, do.call,
duplicated, eval, evalq, get, grep, grepl, is.unsorted, lapply,
mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int,
rank, rbind, rownames, sapply, saveRDS, table, tapply, unique,
unsplit, which.max, which.min
Loading required package: S4Vectors
Attaching package: 'S4Vectors'
The following object is masked from 'package:utils':
findMatches
The following objects are masked from 'package:base':
I, expand.grid, unname
Loading required package: IRanges
Attaching package: 'IRanges'
The following object is masked from 'package:grDevices':
windows
Loading required package: Seqinfo
Loading required package: InteractionSet
Loading required package: SummarizedExperiment
Loading required package: MatrixGenerics
Loading required package: matrixStats
Attaching package: 'MatrixGenerics'
The following objects are masked from 'package:matrixStats':
colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse,
colCounts, colCummaxs, colCummins, colCumprods, colCumsums,
colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs,
colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats,
colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds,
colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads,
colWeightedMeans, colWeightedMedians, colWeightedSds,
colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet,
rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods,
rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps,
rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins,
rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks,
rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs, rowVars,
rowWeightedMads, rowWeightedMeans, rowWeightedMedians,
rowWeightedSds, rowWeightedVars
Loading required package: Biobase
Welcome to Bioconductor
Vignettes contain introductory material; view with
'browseVignettes()'. To cite Bioconductor, see
'citation("Biobase")', and for packages 'citation("pkgname")'.
Attaching package: 'Biobase'
The following object is masked from 'package:MatrixGenerics':
rowMedians
The following objects are masked from 'package:matrixStats':
anyMissing, rowMedians
>
> test_check("diffHic")
[ FAIL 0 | WARN 0 | SKIP 0 | PASS 427 ]
>
> proc.time()
user system elapsed
57.53 1.79 69.59
diffHic.Rcheck/diffHic-Ex.timings
| name | user | system | elapsed | |
| DNaseHiC | 0.16 | 0.05 | 0.27 | |
| annotatePairs | 0.53 | 0.01 | 0.55 | |
| boxPairs | 0.76 | 0.02 | 0.78 | |
| clusterPairs | 1.06 | 0.01 | 1.08 | |
| compartmentalize | 1.75 | 0.02 | 1.76 | |
| connectCounts | 1.83 | 0.05 | 1.97 | |
| consolidatePairs | 0.45 | 0.01 | 0.50 | |
| correctedContact | 0.55 | 0.00 | 0.55 | |
| cutGenome | 7.23 | 0.17 | 7.78 | |
| diClusters | 1.97 | 0.00 | 1.97 | |
| diffHicUsersGuide | 0 | 0 | 0 | |
| domainDirections | 0.19 | 0.00 | 0.22 | |
| enrichedPairs | 0.48 | 0.00 | 0.48 | |
| extractPatch | 0.44 | 0.00 | 0.49 | |
| filterDiag | 0.11 | 0.00 | 0.14 | |
| filterPeaks | 0.19 | 0.00 | 0.19 | |
| filters | 0.29 | 0.00 | 0.29 | |
| getArea | 0.21 | 0.00 | 0.21 | |
| getPairData | 0.06 | 0.02 | 0.10 | |
| loadData | 0.08 | 0.00 | 0.10 | |
| marginCounts | 0.68 | 0.00 | 0.72 | |
| mergeCMs | 0.49 | 0.00 | 0.56 | |
| mergePairs | 0.22 | 0.03 | 0.36 | |
| neighborCounts | 0.11 | 0.00 | 0.14 | |
| normalizeCNV | 9.14 | 0.14 | 9.28 | |
| pairParam | 0.12 | 0.00 | 0.13 | |
| plotDI | 0.61 | 0.02 | 0.62 | |
| plotPlaid | 1.55 | 0.03 | 1.89 | |
| preparePairs | 0.25 | 0.12 | 0.63 | |
| prunePairs | 0.19 | 0.11 | 0.50 | |
| readMTX2IntSet | 0.48 | 0.03 | 0.51 | |
| savePairs | 0.06 | 0.00 | 0.30 | |
| squareCounts | 0.74 | 0.03 | 1.37 | |
| totalCounts | 0.47 | 0.02 | 0.58 | |