R Under development (unstable) (2024-04-27 r86487 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.

> library(NewmanOmics)
> 
> ### pairedStat
> data(LungPair)
> lung <- as.matrix(log2(1 + LungPair))
> 
> set.seed(12345)
> normal <- lung[, 1, drop=FALSE]
> tumor  <- lung[, 2, drop=FALSE]
> ps <- pairedStat(normal, tumor)
> summary(ps@p.values)
 TCGA.38.4625.01A.01R.1758.07
 Min.   :0.0000              
 1st Qu.:0.2723              
 Median :0.5709              
 Mean   :0.5564              
 3rd Qu.:0.8499              
 Max.   :1.0000              
> 
> fm <- fitMix3(as.vector(ps@p.values), print.level = 0)
Warning message:
In nlm(logLikelihood, mle, q = datavec, z = Z, stepmax = 10000,  :
  NA/Inf replaced by maximum positive value
> ### reduce displayed digits to avoid differences in NLM
> ### on different platforms
> ### FAILS anyway on M1Max and NoLD tests
> ## formatC(fm@mle, digits=3, format="e")
> ## round(fm@psi, 5)
> ## But the plots are indistinguishable.
> plot(fm)
> 
> 
> proc.time()
   user  system elapsed 
   8.12    0.31    8.37