Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Drop {mpcbs} tests and examples, because that package is broken in R (>= 4.2.0) #29

Closed
HenrikBengtsson opened this issue Nov 10, 2022 · 1 comment
Milestone

Comments

@HenrikBengtsson
Copy link
Owner

From R-hub checks:

About to run xvfb-run R CMD check aroma.core_3.3.0.tar.gz
'getOption("repos")' replaces Bioconductor standard repositories, see
'?repositories' for details
replacement repositories:
CRAN: https://cloud.r-project.org
* using log directory ‘/home/docker/aroma.core.Rcheck’
* using R Under development (unstable) (2022-11-09 r83320)
* using platform: x86_64-pc-linux-gnu (64-bit)
* using session charset: UTF-8
* checking for file ‘aroma.core/DESCRIPTION’ ... OK
* this is package ‘aroma.core’ version ‘3.3.0’
* checking package namespace information ... OK
* checking package dependencies ... NOTE
Packages suggested but not available for checking: 'EBImage', 'GLAD'
* 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 ‘aroma.core’ can be installed ... OK
* checking installed package size ... OK
* checking package 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 R 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 loading without being on the library search path ... 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
Package unavailable to check Rd xrefs: ‘GLAD’
* 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 examples ... ERROR
Running examples in ‘aroma.core-Ex.R’ failed
The error most likely occurred in:
> ### Name: segmentByMPCBS.RawGenomicSignals
> ### Title: Segment copy numbers using the multi-platform CBS (mpCBS) method
> ### Aliases: segmentByMPCBS.RawGenomicSignals
> ### RawGenomicSignals.segmentByMPCBS
> ### segmentByMPCBS,RawGenomicSignals-method
> ### Keywords: internal methods IO
>
> ### ** Examples
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Simulating copy-number data from multiple platforms
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Piecewise-constant copy-number state function
> cnState <- function(x, ...) {
+ n <- length(x)
+ mu <- double(n)
+ mu[20e6 <= x & x <= 30e6] <- +1
+ mu[65e6 <= x & x <= 80e6] <- -1
+ mu
+ } # cnState()
>
> xMax <- 100e6
>
> Js <- c(200, 400, 100)
> bs <- c(1, 1.4, 0.5)
> as <- c(0, +0.5, -0.5)
> sds <- c(0.5, 0.3, 0.8)
>
> cnList <- list()
> for (kk in seq_along(Js)) {
+ J <- Js[kk]
+ a <- as[kk]
+ b <- bs[kk]
+ sd <- sds[kk]
+ x <- sort(runif(J, max=xMax))
+ mu <- a + b * cnState(x)
+ eps <- rnorm(J, sd=sd)
+ y <- mu + eps
+ cn <- RawCopyNumbers(y, x)
+ cnList[[kk]] <- cn
+ }
>
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Merge platform data (record their sources in 'id')
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> cn <- Reduce(append, cnList)
> plot(cn, ylim=c(-3,3), col=cn$id)
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Segment
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> legend <- c()
>
> if (require("DNAcopy")) {
+ fit <- segmentByCBS(cn)
+ cnr <- extractCopyNumberRegions(fit)
+ print(cnr)
+ drawLevels(cnr, col="white", lwd=6)
+ drawLevels(cnr, col="red", lwd=3)
+ legend <- c(legend, red="CBS")
+ }
Loading required package: DNAcopy
CopyNumberRegions:
Number of regions: 5
>
>
> if (require("mpcbs")) {
+ fit <- local({
+ ## WORKAROUND: There's a _R_CHECK_LENGTH_1_LOGIC2_ bug in
+ ## mpcbs::mpcbs.mbic(). Until fixed, we need to disable this check
+ ## while calling segmentByMPCBS().
+ ovalue <- Sys.getenv("_R_CHECK_LENGTH_1_LOGIC2_")
+ on.exit(Sys.setenv("_R_CHECK_LENGTH_1_LOGIC2_"=ovalue))
+ ovalue <- Sys.setenv("_R_CHECK_LENGTH_1_LOGIC2_"="warn,verbose")
+ segmentByMPCBS(cn)
+ })
+ cnr <- extractCopyNumberRegions(fit)
+ print(cnr)
+ drawLevels(cnr, col="white", lwd=6)
+ drawLevels(cnr, col="blue", lwd=3)
+ legend <- c(legend, blue="MPCBS")
+ }
Loading required package: mpcbs
Loading required package: fields
Loading required package: spam
Spam version 2.9-1 (2022-08-07) is loaded.
Type 'help( Spam)' or 'demo( spam)' for a short introduction
and overview of this package.
Help for individual functions is also obtained by adding the
suffix '.spam' to the function name, e.g. 'help( chol.spam)'.
Attaching package: ‘spam’
The following object is masked from ‘package:aroma.core’:
display
The following object is masked from ‘package:R.utils’:
cleanup
The following objects are masked from ‘package:base’:
backsolve, forwardsolve
Loading required package: viridis
Loading required package: viridisLite
Try help(fields) to get started.
Error in is.na(rratio) || (length(rratio) != K) :
'length = 3' in coercion to 'logical(1)'
Calls: local ... system.time -> do.call -> mpcbs.mbic -> fcompute.max.ProjectedZ
Timing stopped at: 0 0 0.001
Execution halted
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
Running ‘PairedPSCNData,SEG.R’
Running ‘PairedPSCNData.R’
Running ‘RawCopyNumbers,states.R’
Running ‘RawCopyNumbers.R’
Running ‘RawGenomicSignals.R’
Running ‘RawGenomicSignals.SEG,MP.R’
ERROR
Running the tests in ‘tests/RawGenomicSignals.SEG,MP.R’ failed.
Last 13 lines of output:
cleanup

The following objects are masked from 'package:base':

backsolve, forwardsolve

Loading required package: viridis
Loading required package: viridisLite

Try help(fields) to get started.
Error in is.na(rratio) || (length(rratio) != K) :
'length = 3' in coercion to 'logical(1)'
Calls: local ... system.time -> do.call -> mpcbs.mbic -> fcompute.max.ProjectedZ
Timing stopped at: 0 0 0.001
Execution halted
* checking PDF version of manual ... OK
* DONE
Status: 2 ERRORs, 2 NOTEs
See ‘/home/docker/aroma.core.Rcheck/00check.log’ for details.
@HenrikBengtsson HenrikBengtsson added this to the Next release milestone Nov 10, 2022
@HenrikBengtsson HenrikBengtsson changed the title TESTS: Drop {mpcps} tests, because that package is broken in R (>= 4.2.0) Drop {mpcps} tests and examples, because that package is broken in R (>= 4.2.0) Nov 10, 2022
@HenrikBengtsson
Copy link
Owner Author

There is (at least) one _R_CHECK_LENGTH_1_CONDITION_ bug and (at least) one _R_CHECK_LENGTH_1_LOGIC2_ bug:

$ wget http://download.r-forge.r-project.org/src/contrib/mpcbs_1.1.1.tar.gz
$ R CMD check --as-cran mpcbs_1.1.1.tar.gz 
* using log directory ‘/home/hb/repositories/aroma.core/.local/mpcbs.Rcheck’
* using R version 4.2.2 (2022-10-31)
* using platform: x86_64-pc-linux-gnu (64-bit)
* using session charset: UTF-8
* using option ‘--as-cran’
* checking for file ‘mpcbs/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘mpcbs’ version ‘1.1.1’
* checking CRAN incoming feasibility ... NOTE
Maintainer: ‘Nancy R. Zhang <nzhang@stanford.edu>’

New submission

Non-FOSS package license (What license is it under?)

The Title field should be in title case. Current version is:
‘Multi-sample and Multiplatform CNV.’
In title case that is:
‘Multi-Sample and Multiplatform CNV.’

The Date field is over a month old.

This build time stamp is over a month old.
* 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 serialization versions ... OK
* checking whether package ‘mpcbs’ can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking for future file timestamps ... OK
* checking DESCRIPTION meta-information ... NOTE
Malformed Title field: should not end in a period.
Non-standard license specification:
  What license is it under?
Standardizable: FALSE
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking R 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 loading without being on the library search path ... OK
* checking use of S3 registration ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... WARNING
merge:
  function(x, y, ...)
merge.pos:
  function(pos, anchors)

plot:
  function(x, ...)
plot.crossplatform:
  function(pos, y, anchor, chpts, ranking, yhat, platform.names,
           num.panels, xlim, ...)

update:
  function(object, ...)
update.bic.baseline:
  function(this.y, segbag, carriers, this.sum.deltasq,
           this.sum.deltasqwin, newchpt, samp.ord, win, y.var,
           PRIOR.CARRIER.PROB, PRIOR.SEGCOUNT, do.plots)

See section ‘Generic functions and methods’ in the ‘Writing R
Extensions’ manual.

Found the following apparent S3 methods exported but not registered:
  merge.pos plot.crossplatform update.bic.baseline
See section ‘Registering S3 methods’ in the ‘Writing R Extensions’
manual.
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... [11s/11s] NOTE
NB: .First.lib is obsolete and will not be used in R >= 3.0.0

Chisq.Contrib.test: no visible global function definition for ‘runif’
Chisq.Contrib.test: no visible global function definition for ‘par’
Chisq.Contrib.test: no visible global function definition for ‘image’
ComputeZ.R: no visible binding for global variable ‘y’
Remove.Overlap.test: no visible global function definition for ‘par’
Remove.Overlap.test: no visible global function definition for ‘runif’
Remove.Overlap.test: no visible global function definition for ‘image’
compute.var: no visible binding for global variable ‘median’
compute.var: no visible global function definition for ‘median’
computeBeta: no visible global function definition for ‘integrate’
computeMoments: no visible global function definition for ‘integrate’
computeZ.onechange: no visible global function definition for ‘pchisq’
computeZ.onechange.sample: no visible global function definition for
  ‘pchisq’
computeZ.squarewave.sample: no visible global function definition for
  ‘pchisq’
computeZ.test: no visible global function definition for ‘par’
computeZ.test: no visible global function definition for ‘rnorm’
computeZ.test: no visible binding for global variable ‘z’
computeZ.test: no visible global function definition for ‘image’
fcompute.max.Z: no visible binding for global variable ‘var’
fscan.max: no visible global function definition for ‘lines’
fscan.max: no visible global function definition for ‘points’
fscan.max: no visible binding for global variable ‘MIN.SNPS’
getCutoffEpidemicChangeLocationKnown: no visible global function
  definition for ‘qchisq’
mpcbs: no visible global function definition for ‘pdf’
mpcbs: no visible global function definition for ‘par’
mpcbs: no visible global function definition for ‘segments’
mpcbs: no visible global function definition for ‘dev.off’
mpcbs.mbic: no visible global function definition for ‘pdf’
mpcbs.mbic: no visible global function definition for ‘par’
mpcbs.mbic: no visible global function definition for ‘segments’
mpcbs.mbic: no visible global function definition for ‘dev.off’
mpscan: no visible global function definition for ‘par’
mpscan: no visible global function definition for ‘lines’
mpscan: no visible global function definition for ‘segments’
plot.crossplatform: no visible global function definition for ‘par’
plot.crossplatform: no visible global function definition for ‘lines’
plot.crossplatform: no visible global function definition for
  ‘segments’
plot.crossplatform: no visible global function definition for ‘text’
pvalueMultisampleWeightedChisq: no visible global function definition
  for ‘integrate’
pvalueMultisampleWeightedChisqold: no visible global function
  definition for ‘integrate’
pvalueSumChisq: no visible global function definition for ‘integrate’
scatter.ci: no visible global function definition for ‘segments’
scatter.ci: no visible global function definition for ‘points’
scatter.cpfit: no visible global function definition for ‘par’
scatter.cpfit: no visible binding for global variable ‘sd’
scatter.cpfit: no visible global function definition for ‘abline’
scatter.cpfit: no visible global function definition for ‘points’
update.bic.baseline: no visible global function definition for ‘par’
update.bic.baseline: no visible global function definition for ‘points’
update.bic.baseline: no visible global function definition for
  ‘segments’
vu: no visible global function definition for ‘pnorm’
wpca: no visible global function definition for ‘lm’
wpca: no visible global function definition for ‘par’
wpca: no visible global function definition for ‘abline’
wpca: no visible global function definition for ‘points’
wpca: no visible global function definition for ‘lines’
wpca: no visible global function definition for ‘legend’
Undefined global functions or variables:
  MIN.SNPS abline dev.off image integrate legend lines lm median par
  pchisq pdf pnorm points qchisq rnorm runif sd segments text var y z
Consider adding
  importFrom("grDevices", "dev.off", "pdf")
  importFrom("graphics", "abline", "image", "legend", "lines", "par",
             "points", "segments", "text")
  importFrom("stats", "integrate", "lm", "median", "pchisq", "pnorm",
             "qchisq", "rnorm", "runif", "sd", "var")
to your NAMESPACE file.
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd line widths ... NOTE
Rd file 'cross.platform.consensus.Rd':
  \examples lines wider than 100 characters:
     consensus.cn = res$consensus.cn[[1]]    # consensus.cn[t] is consensus copy number estimate at anchor$merged.pos[t]
     plot.crossplatform(pos,y, yhat=seg$yhat, anchor=anchor,platform.names=platform.names, num.panels=K+1, col="darkgray")
     plot.crossplatform(pos,y, yhat=seg$yhat, anchor=anchor,platform.names=platform.names, num.panels=K+1, xlim=c(xmin,xmax), col="darkgray" ... [TRUNCATED]

Rd file 'mpcbs-package.Rd':
  \examples lines wider than 100 characters:
     consensus.cn = res$consensus.cn[[1]]    # consensus.cn[t] is consensus copy number estimate at anchor$merged.pos[t]
     plot.crossplatform(pos,y, yhat=seg$yhat, anchor=anchor,platform.names=platform.names, num.panels=K+1, col="darkgray")
     plot.crossplatform(pos,y, yhat=seg$yhat, anchor=anchor,platform.names=platform.names, num.panels=K+1, xlim=c(xmin,xmax), col="darkgray" ... [TRUNCATED]
             xlab="Genome Position", ylab="Log Ratio", main="Consensus Copy Number Estimate", cex.lab=2, xlim=c(xmin,xmax))

Rd file 'mpcbs.Rd':
  \examples lines wider than 100 characters:
     #plot.crossplatform(pos,y, yhat=seg$yhat, chpts=seg$chpts, anchor=anchor, platform.names=platform.names, col="darkgray")

Rd file 'mpcbs.mbic.Rd':
  \examples lines wider than 100 characters:
     plot.crossplatform(pos,y, yhat=seg$yhat, chpts=seg$chpts, anchor=anchor, platform.names=platform.names, col="darkgray")

Rd file 'plot.crossplatform.Rd':
  \usage lines wider than 90 characters:
     plot.crossplatform(pos, y, anchor, chpts = NULL, ranking = NULL, yhat = NULL, platform.names, num.panels = NA, xlim = NULL, ...)
  \examples lines wider than 100 characters:
     plot.crossplatform(pos,y, yhat=seg$yhat, anchor=anchor, platform.names=platform.names, col="darkgray")
     plot.crossplatform(pos,y, yhat=seg$yhat, chpts=seg$chpts, anchor=anchor, platform.names=platform.names, col="darkgray")
     plot.crossplatform(pos,y, yhat=seg$yhat, chpts=seg$chpts, anchor=anchor, platform.names=platform.names, col="darkgray", xlim=c(5e7,1e8) ... [TRUNCATED]
     plot.crossplatform(pos,y, yhat=seg$yhat, chpts=seg$chpts, anchor=anchor, platform.names=platform.names, col="darkgray", xlim=c(5e7,6e7) ... [TRUNCATED]

These lines will be truncated in the PDF manual.
* checking Rd cross-references ... OK
* checking for missing documentation entries ... WARNING
Undocumented code objects:
  ‘ComputeProjectedZ.fromS.R.segments’
Undocumented data sets:
  ‘mpcbs.example’
All user-level objects in a package should have documentation entries.
See chapter ‘Writing R documentation files’ in the ‘Writing R
Extensions’ manual.
* checking for code/documentation mismatches ... WARNING
Variables with usage in documentation object 'mpcbsexample' but not in code:
  ‘mpcbsexample’

Codoc mismatches from documentation object 'compute.var':
compute.var
  Code: function(y, use.mean = TRUE)
  Docs: function(y)
  Argument names in code not in docs:
    use.mean

Codoc mismatches from documentation object 'merge.pos':
merge.pos
  Code: function(pos, anchors = NULL)
  Docs: function(pos, anchors = 1)
  Mismatches in argument default values:
    Name: 'anchors' Code: NULL Docs: 1

Codoc mismatches from documentation object 'mpcbs':
mpcbs
  Code: function(y, pos, anchor, win, Z = NULL, MIN.SNPs = 2,
                 MIN.BP.LEN = 10000, ALPHA = 0, rratio = NULL, f =
                 NULL, MAX.CHPTS = 50, WCHISQ.CUTOFF = NA,
                 platform.names = NULL, plots = TRUE, plotspdf = NULL,
                 use.filtered.scan = TRUE, useProjectedChisquare =
                 TRUE)
  Docs: function(y, pos, anchors)
  Argument names in code not in docs:
    anchor win Z MIN.SNPs MIN.BP.LEN ALPHA rratio f MAX.CHPTS
    WCHISQ.CUTOFF platform.names plots plotspdf use.filtered.scan
    useProjectedChisquare
  Argument names in docs not in code:
    anchors
  Mismatches in argument names:
    Position: 3 Code: anchor Docs: anchors

Codoc mismatches from documentation object 'mpcbs.mbic':
mpcbs.mbic
  Code: function(y, pos, anchor, MIN.SNPs = 2, MIN.BP.LEN = 1000,
                 rratio = NULL, f = NULL, MAX.CHPTS = 30,
                 platform.names = NULL, plots = TRUE, plotspdf = NULL,
                 use.filtered.scan = TRUE)
  Docs: function(y, pos, anchor, MIN.SNPs = 2, MAX.CHPTS = 30,
                 platform.names, plots = TRUE)
  Argument names in code not in docs:
    MIN.BP.LEN rratio f plotspdf use.filtered.scan
  Mismatches in argument names:
    Position: 5 Code: MIN.BP.LEN Docs: MAX.CHPTS
    Position: 6 Code: rratio Docs: platform.names
    Position: 7 Code: f Docs: plots
  Mismatches in argument default values:
    Name: 'platform.names' Code: NULL Docs: 

Codoc mismatches from documentation object 'mpscan':
mpscan
  Code: function(y, pos, anchors, win, Z = NULL, MIN.SNPs = 1, ALPHA =
                 0, WCHISQ.CUTOFF = NA)
  Docs: function(y, pos, anchors, win)
  Argument names in code not in docs:
    Z MIN.SNPs ALPHA WCHISQ.CUTOFF

* checking Rd \usage sections ... WARNING
Documented arguments not in \usage in documentation object 'compute.var':
  ‘use.mean’

Undocumented arguments in documentation object 'mpcbs'
  ‘anchors’
Documented arguments not in \usage in documentation object 'mpcbs':
  ‘anchor’ ‘MIN.SNPs’ ‘MIN.BP.LEN’ ‘WCHISQ.CUTOFF’ ‘rratio’
  ‘platform.names’ ‘plots’ ‘plotspdf’ ‘use.filtered.scan’

Documented arguments not in \usage in documentation object 'mpcbs.mbic':
  ‘MIN.BP.LEN’ ‘rratio’ ‘plotspdf’ ‘use.filtered.scan’

Undocumented arguments in documentation object 'mpscan'
  ‘anchors’
Documented arguments not in \usage in documentation object 'mpscan':
  ‘anchor’ ‘MIN.SNPs’ ‘WCHISQ.CUTOFF’

Functions with \usage entries need to have the appropriate \alias
entries, and all their arguments documented.
The \usage entries must correspond to syntactically valid R code.
See chapter ‘Writing R documentation files’ in the ‘Writing R
Extensions’ manual.
S3 methods shown with full name in documentation object 'merge.pos':
  ‘merge.pos’

S3 methods shown with full name in documentation object 'plot.crossplatform':
  ‘plot.crossplatform’

The \usage entries for S3 methods should use the \method markup and not
their full name.
See chapter ‘Writing R documentation files’ in the ‘Writing R
Extensions’ manual.
* 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 data for ASCII and uncompressed saves ... OK
* checking line endings in C/C++/Fortran sources/headers ... OK
* checking pragmas in C/C++ headers and code ... OK
* checking compilation flags used ... OK
* checking compiled code ... NOTE
File ‘mpcbs/libs/mpcbs.so’:
  Found no calls to: ‘R_registerRoutines’, ‘R_useDynamicSymbols’

It is good practice to register native routines and to disable symbol
search.

See ‘Writing portable packages’ in the ‘Writing R Extensions’ manual.
* checking examples ... ERROR
Running examples in ‘mpcbs-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: cross.platform.consensus
> ### Title: Compute the cross platform consensus copy number given a
> ###   segmentation.
> ### Aliases: cross.platform.consensus
> ### Keywords: weighted least squares
> 
> ### ** Examples
> 
> 
> data(mpcbs.example)
> 
> # there are 3 platforms represented in this data example: Illumina, Affymetrix, and Agilent.
> names(mpcbs.example)
[1] "illu" "affy" "agil"
> 
> # K is the number of platforms.
> K=3
> 
> # Store the chromosome positions in vector pos,
> # the intensities in vector y:
> 
> pos=vector("list",K)
> pos[[1]] = mpcbs.example$illu[,1]
> pos[[2]] = mpcbs.example$affy[,1]
> pos[[3]] = mpcbs.example$agil[,1]
> 
> y = vector("list",K)
> y[[1]] = mpcbs.example$illu[,2]
> y[[2]] = mpcbs.example$affy[,2]
> y[[3]] = mpcbs.example$agil[,2]
> 
> # Names of the platforms:
> platform.names=c("Illumina","Affymetrix","Agilent")
> 
> # Get the anchor set.
> anchor = merge.pos(pos)
> 
> # Perform the segmentation.
> seg<-mpcbs.mbic(y,pos,anchor, MAX.CHPTS=10, platform.names=platform.names,plots=TRUE)
Using anchors containing 142404 positions.
Use filtered scan:  TRUE , maximum change-points:  10 
 ----------- FAILURE REPORT -------------- 
 --- failure: length > 1 in coercion to logical ---
 --- srcref --- 
: 
 --- package (from environment) --- 
mpcbs
 --- call from context --- 
fcompute.max.ProjectedZ(S, SST, imap, win, rratio, MIN.SNPs, 
    f = f)
 --- call from argument --- 
is.na(rratio) || (length(rratio) != K)
 --- R stacktrace ---
where 1: fcompute.max.ProjectedZ(S, SST, imap, win, rratio, MIN.SNPs, 
    f = f)
where 2: mpcbs.mbic(y, pos, anchor, MAX.CHPTS = 10, platform.names = platform.names, 
    plots = TRUE)

 --- value of length: 3 type: logical ---
[1] FALSE FALSE FALSE
 --- function from context --- 
function (this.S, this.SST, this.imap, win, rratio, MIN.SNPs, 
    f = NULL) 
{
    K = ncol(this.S)
    this.T = nrow(this.S)
    if (is.na(rratio) || (length(rratio) != K)) {
        rratio = rep(1, K)
    }
    win = min(win, this.T - 1)
    temp = fscan.max(this.S, this.SST, this.imap = this.imap, 
        MIN.SNPs = MIN.SNPs, rratio = rratio, f = f, use.Project.statistic = TRUE, 
        verbose = FALSE)
    bestchpt = temp$seg
    bestZ = temp$maxZ
    list(bestchpt = bestchpt, bestZ = bestZ)
}
<bytecode: 0x55ad735e2ba0>
<environment: namespace:mpcbs>
 --- function search by body ---
Function fcompute.max.ProjectedZ in namespace mpcbs has this body.
 ----------- END OF FAILURE REPORT -------------- 
Error in is.na(rratio) || (length(rratio) != K) : 
  'length(x) = 3 > 1' in coercion to 'logical(1)'
Calls: mpcbs.mbic -> fcompute.max.ProjectedZ
Execution halted
* checking PDF version of manual ... OK
* checking HTML version of manual ... OK
* checking for non-standard things in the check directory ... OK
* checking for detritus in the temp directory ... OK
* DONE

Status: 1 ERROR, 4 WARNINGs, 5 NOTEs
See
  ‘/home/hb/repositories/aroma.core/.local/mpcbs.Rcheck/00check.log’
for details.

@HenrikBengtsson HenrikBengtsson changed the title Drop {mpcps} tests and examples, because that package is broken in R (>= 4.2.0) Drop {mpcbs} tests and examples, because that package is broken in R (>= 4.2.0) Nov 10, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant