* using log directory ‘/data/gannet/ripley/R/packages/tests-LENGTH1/ALS.Rcheck’ * using R Under development (unstable) (2022-04-26 r82260) * using platform: x86_64-pc-linux-gnu (64-bit) * using session charset: UTF-8 * using option ‘--no-stop-on-test-error’ * checking for file ‘ALS/DESCRIPTION’ ... OK * checking extension type ... Package * this is package ‘ALS’ version ‘0.0.6’ * 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 whether package ‘ALS’ can be installed ... 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 use of S3 registration ... 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 ... [6s/18s] OK * checking Rd files ... OK * checking Rd metadata ... OK * checking Rd line widths ... OK * checking Rd cross-references ... OK * 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 contents of ‘data’ directory ... OK * checking data for non-ASCII characters ... OK * checking data for ASCII and uncompressed saves ... OK * checking installed files from ‘inst/doc’ ... OK * checking examples ... ERROR Running examples in ‘ALS-Ex.R’ failed The error most likely occurred in: > ### Name: als > ### Title: alternating least squares multivariate curve resolution > ### (MCR-ALS) > ### Aliases: als > ### Keywords: optimize > > ### ** Examples > > ## load 2 matrix datasets into variables d1 and d2 > ## load starting values for elution profiles > ## into variables Cstart1 and Cstart2 > ## load time labels as x, m/z values as x2 > data(multiex) > > ## starting values for elution profiles > matplot(x,Cstart1,type="l") > matplot(x,Cstart2,type="l",add=TRUE) > > ## using MCR-ALS, improve estimates for mass spectra S and the two > ## matrices of elution profiles > ## apply unimodality constraints to the elution profile estimates > ## note that the starting estimates for S just contain a dummy matrix > > test0 <- als(CList=list(Cstart1,Cstart2),S=matrix(1,nrow=400,ncol=2), + PsiList=list(d1,d2), x=x, x2=x2, uniC=TRUE, normS=0) Initial RSS 3.039967e+13 Iteration (opt. S): 1, RSS: 1.330703e+12, RD: 0.9562264 ----------- FAILURE REPORT -------------- --- failure: length > 1 in coercion to logical --- --- srcref --- : --- package (from environment) --- ALS --- call from context --- als(CList = list(Cstart1, Cstart2), S = matrix(1, nrow = 400, ncol = 2), PsiList = list(d1, d2), x = x, x2 = x2, uniC = TRUE, normS = 0) --- call from argument --- normS > S && (iter%%2 != b) --- R stacktrace --- where 1: als(CList = list(Cstart1, Cstart2), S = matrix(1, nrow = 400, ncol = 2), PsiList = list(d1, d2), x = x, x2 = x2, uniC = TRUE, normS = 0) --- value of length: 800 type: logical --- [,1] [,2] [1,] FALSE FALSE [2,] FALSE FALSE [3,] FALSE FALSE [4,] FALSE FALSE [5,] FALSE FALSE [6,] FALSE FALSE [7,] FALSE FALSE [8,] FALSE FALSE [9,] FALSE FALSE [10,] FALSE FALSE [11,] FALSE FALSE [12,] FALSE FALSE [13,] FALSE FALSE [14,] FALSE FALSE [15,] FALSE FALSE [16,] FALSE FALSE [17,] FALSE FALSE [18,] FALSE FALSE [19,] FALSE FALSE [20,] FALSE FALSE [21,] FALSE FALSE [22,] FALSE FALSE [23,] FALSE FALSE [24,] FALSE FALSE [25,] FALSE FALSE [26,] FALSE FALSE [27,] FALSE FALSE [28,] FALSE FALSE [29,] FALSE FALSE [30,] FALSE FALSE [31,] FALSE FALSE [32,] FALSE FALSE [33,] FALSE FALSE [34,] FALSE FALSE [35,] FALSE FALSE [36,] FALSE FALSE [37,] FALSE FALSE [38,] FALSE FALSE [39,] FALSE FALSE [40,] FALSE FALSE [41,] FALSE FALSE [42,] FALSE FALSE [43,] FALSE FALSE [44,] FALSE FALSE [45,] FALSE FALSE [46,] FALSE FALSE [47,] FALSE FALSE [48,] FALSE FALSE [49,] FALSE FALSE [50,] FALSE FALSE [51,] FALSE FALSE [52,] FALSE FALSE [53,] FALSE FALSE [54,] FALSE FALSE [55,] FALSE FALSE [56,] FALSE FALSE [57,] FALSE FALSE [58,] FALSE FALSE [59,] FALSE FALSE [60,] FALSE FALSE [61,] FALSE FALSE [62,] FALSE FALSE [63,] FALSE FALSE [64,] FALSE FALSE [65,] FALSE FALSE [66,] FALSE FALSE [67,] FALSE FALSE [68,] FALSE FALSE [69,] FALSE FALSE [70,] FALSE FALSE [71,] FALSE FALSE [72,] FALSE FALSE [73,] FALSE FALSE [74,] FALSE FALSE [75,] FALSE FALSE [76,] FALSE FALSE [77,] FALSE FALSE [78,] FALSE FALSE [79,] FALSE FALSE [80,] FALSE FALSE [81,] FALSE FALSE [82,] FALSE FALSE [83,] FALSE FALSE [84,] FALSE FALSE [85,] FALSE FALSE [86,] FALSE FALSE [87,] FALSE FALSE [88,] FALSE FALSE [89,] FALSE FALSE [90,] FALSE FALSE [91,] FALSE FALSE [92,] FALSE FALSE [93,] FALSE FALSE [94,] FALSE FALSE [95,] FALSE FALSE [96,] FALSE FALSE [97,] FALSE FALSE [98,] FALSE FALSE [99,] FALSE FALSE [100,] FALSE FALSE [101,] FALSE FALSE [102,] FALSE FALSE [103,] FALSE FALSE [104,] FALSE FALSE [105,] FALSE FALSE [106,] FALSE FALSE [107,] FALSE FALSE [108,] FALSE FALSE [109,] FALSE FALSE [110,] FALSE FALSE [111,] FALSE FALSE [112,] FALSE FALSE [113,] FALSE FALSE [114,] FALSE FALSE [115,] FALSE FALSE [116,] FALSE FALSE [117,] FALSE FALSE [118,] FALSE FALSE [119,] FALSE FALSE [120,] FALSE FALSE [121,] FALSE FALSE [122,] FALSE FALSE [123,] FALSE FALSE [124,] FALSE FALSE [125,] FALSE FALSE [126,] FALSE FALSE [127,] FALSE FALSE [128,] FALSE FALSE [129,] FALSE FALSE [130,] FALSE FALSE [131,] FALSE FALSE [132,] FALSE FALSE [133,] FALSE FALSE [134,] FALSE FALSE [135,] FALSE FALSE [136,] FALSE FALSE [137,] FALSE FALSE [138,] FALSE FALSE [139,] FALSE FALSE [140,] FALSE FALSE [141,] FALSE FALSE [142,] FALSE FALSE [143,] FALSE FALSE [144,] FALSE FALSE [145,] FALSE FALSE [146,] FALSE FALSE [147,] FALSE FALSE [148,] FALSE FALSE [149,] FALSE FALSE [150,] FALSE FALSE [151,] FALSE FALSE [152,] FALSE FALSE [153,] FALSE FALSE [154,] FALSE FALSE [155,] FALSE FALSE [156,] FALSE FALSE [157,] FALSE FALSE [158,] FALSE FALSE [159,] FALSE FALSE [160,] FALSE FALSE [161,] FALSE FALSE [162,] FALSE FALSE [163,] FALSE FALSE [164,] FALSE FALSE [165,] FALSE FALSE [166,] FALSE FALSE [167,] FALSE FALSE [168,] FALSE FALSE [169,] FALSE FALSE [170,] FALSE FALSE [171,] FALSE FALSE [172,] FALSE FALSE [173,] FALSE FALSE [174,] FALSE FALSE [175,] FALSE FALSE [176,] FALSE FALSE [177,] FALSE FALSE [178,] FALSE FALSE [179,] FALSE FALSE [180,] FALSE FALSE [181,] FALSE FALSE [182,] FALSE FALSE [183,] FALSE FALSE [184,] FALSE FALSE [185,] FALSE FALSE [186,] FALSE FALSE [187,] FALSE FALSE [188,] FALSE FALSE [189,] FALSE FALSE [190,] FALSE FALSE [191,] FALSE FALSE [192,] FALSE FALSE [193,] FALSE FALSE [194,] FALSE FALSE [195,] FALSE FALSE [196,] FALSE FALSE [197,] FALSE FALSE [198,] FALSE FALSE [199,] FALSE FALSE [200,] FALSE FALSE [201,] FALSE FALSE [202,] FALSE FALSE [203,] FALSE FALSE [204,] FALSE FALSE [205,] FALSE FALSE [206,] FALSE FALSE [207,] FALSE FALSE [208,] FALSE FALSE [209,] FALSE FALSE [210,] FALSE FALSE [211,] FALSE FALSE [212,] FALSE FALSE [213,] FALSE FALSE [214,] FALSE FALSE [215,] FALSE FALSE [216,] FALSE FALSE [217,] FALSE FALSE [218,] FALSE FALSE [219,] FALSE FALSE [220,] FALSE FALSE [221,] FALSE FALSE [222,] FALSE FALSE [223,] FALSE FALSE [224,] FALSE FALSE [225,] FALSE FALSE [226,] FALSE FALSE [227,] FALSE FALSE [228,] FALSE FALSE [229,] FALSE FALSE [230,] FALSE FALSE [231,] FALSE FALSE [232,] FALSE FALSE [233,] FALSE FALSE [234,] FALSE FALSE [235,] FALSE FALSE [236,] FALSE FALSE [237,] FALSE FALSE [238,] FALSE FALSE [239,] FALSE FALSE [240,] FALSE FALSE [241,] FALSE FALSE [242,] FALSE FALSE [243,] FALSE FALSE [244,] FALSE FALSE [245,] FALSE FALSE [246,] FALSE FALSE [247,] FALSE FALSE [248,] FALSE FALSE [249,] FALSE FALSE [250,] FALSE FALSE [251,] FALSE FALSE [252,] FALSE FALSE [253,] FALSE FALSE [254,] FALSE FALSE [255,] FALSE FALSE [256,] FALSE FALSE [257,] FALSE FALSE [258,] FALSE FALSE [259,] FALSE FALSE [260,] FALSE FALSE [261,] FALSE FALSE [262,] FALSE FALSE [263,] FALSE FALSE [264,] FALSE FALSE [265,] FALSE FALSE [266,] FALSE FALSE [267,] FALSE FALSE [268,] FALSE FALSE [269,] FALSE FALSE [270,] FALSE FALSE [271,] FALSE FALSE [272,] FALSE FALSE [273,] FALSE FALSE [274,] FALSE FALSE [275,] FALSE FALSE [276,] FALSE FALSE [277,] FALSE FALSE [278,] FALSE FALSE [279,] FALSE FALSE [280,] FALSE FALSE [281,] FALSE FALSE [282,] FALSE FALSE [283,] FALSE FALSE [284,] FALSE FALSE [285,] FALSE FALSE [286,] FALSE FALSE [287,] FALSE FALSE [288,] FALSE FALSE [289,] FALSE FALSE [290,] FALSE FALSE [291,] FALSE FALSE [292,] FALSE FALSE [293,] FALSE FALSE [294,] FALSE FALSE [295,] FALSE FALSE [296,] FALSE FALSE [297,] FALSE FALSE [298,] FALSE FALSE [299,] FALSE FALSE [300,] FALSE FALSE [301,] FALSE FALSE [302,] FALSE FALSE [303,] FALSE FALSE [304,] FALSE FALSE [305,] FALSE FALSE [306,] FALSE FALSE [307,] FALSE FALSE [308,] FALSE FALSE [309,] FALSE FALSE [310,] FALSE FALSE [311,] FALSE FALSE [312,] FALSE FALSE [313,] FALSE FALSE [314,] FALSE FALSE [315,] FALSE FALSE [316,] FALSE FALSE [317,] FALSE FALSE [318,] FALSE FALSE [319,] FALSE FALSE [320,] FALSE FALSE [321,] FALSE FALSE [322,] FALSE FALSE [323,] FALSE FALSE [324,] FALSE FALSE [325,] FALSE FALSE [326,] FALSE FALSE [327,] FALSE FALSE [328,] FALSE FALSE [329,] FALSE FALSE [330,] FALSE FALSE [331,] FALSE FALSE [332,] FALSE FALSE [333,] FALSE FALSE [334,] FALSE FALSE [335,] FALSE FALSE [336,] FALSE FALSE [337,] FALSE FALSE [338,] FALSE FALSE [339,] FALSE FALSE [340,] FALSE FALSE [341,] FALSE FALSE [342,] FALSE FALSE [343,] FALSE FALSE [344,] FALSE FALSE [345,] FALSE FALSE [346,] FALSE FALSE [347,] FALSE FALSE [348,] FALSE FALSE [349,] FALSE FALSE [350,] FALSE FALSE [351,] FALSE FALSE [352,] FALSE FALSE [353,] FALSE FALSE [354,] FALSE FALSE [355,] FALSE FALSE [356,] FALSE FALSE [357,] FALSE FALSE [358,] FALSE FALSE [359,] FALSE FALSE [360,] FALSE FALSE [361,] FALSE FALSE [362,] FALSE FALSE [363,] FALSE FALSE [364,] FALSE FALSE [365,] FALSE FALSE [366,] FALSE FALSE [367,] FALSE FALSE [368,] FALSE FALSE [369,] FALSE FALSE [370,] FALSE FALSE [371,] FALSE FALSE [372,] FALSE FALSE [373,] FALSE FALSE [374,] FALSE FALSE [375,] FALSE FALSE [376,] FALSE FALSE [377,] FALSE FALSE [378,] FALSE FALSE [379,] FALSE FALSE [380,] FALSE FALSE [381,] FALSE FALSE [382,] FALSE FALSE [383,] FALSE FALSE [384,] FALSE FALSE [385,] FALSE FALSE [386,] FALSE FALSE [387,] FALSE FALSE [388,] FALSE FALSE [389,] FALSE FALSE [390,] FALSE FALSE [391,] FALSE FALSE [392,] FALSE FALSE [393,] FALSE FALSE [394,] FALSE FALSE [395,] FALSE FALSE [396,] FALSE FALSE [397,] FALSE FALSE [398,] FALSE FALSE [399,] FALSE FALSE [400,] FALSE FALSE --- function from context --- function (CList, PsiList, S = matrix(), WList = list(), thresh = 0.001, maxiter = 100, forcemaxiter = FALSE, optS1st = TRUE, x = 1:nrow(CList[[1]]), x2 = 1:nrow(S), baseline = FALSE, fixed = vector("list", length(PsiList)), uniC = FALSE, uniS = FALSE, nonnegC = TRUE, nonnegS = TRUE, normS = 0, closureC = list()) { RD <- 10^20 PsiAll <- do.call("rbind", PsiList) resid <- vector("list", length(PsiList)) if (length(WList) == 0) { WList <- vector("list", length(PsiList)) for (i in 1:length(PsiList)) WList[[i]] <- matrix(1, nrow(PsiList[[1]]), ncol(PsiList[[1]])) } W <- do.call("rbind", WList) for (i in 1:length(PsiList)) resid[[i]] <- matrix(0, nrow(PsiList[[i]]), ncol(PsiList[[i]])) for (j in 1:length(PsiList)) { for (i in 1:nrow(PsiList[[j]])) { resid[[j]][i, ] <- PsiList[[j]][i, ] - CList[[j]][i, ] %*% t(S * WList[[j]][i, ]) } } initialrss <- oldrss <- sum(unlist(resid)^2) cat("Initial RSS", initialrss, "\n") iter <- 1 b <- if (optS1st) 1 else 0 oneMore <- FALSE while (((RD > thresh || forcemaxiter) && maxiter >= iter) || oneMore) { if (iter%%2 == b) S <- getS(CList, PsiAll, S, W, baseline, uniS, nonnegS, normS, x2) else CList <- getCList(S, PsiList, CList, WList, resid, x, baseline, fixed, uniC, nonnegC, closureC) for (j in 1:length(PsiList)) { for (i in 1:nrow(PsiList[[j]])) { resid[[j]][i, ] <- PsiList[[j]][i, ] - CList[[j]][i, ] %*% t(S * WList[[j]][i, ]) } } rss <- sum(unlist(resid)^2) RD <- ((oldrss - rss)/oldrss) oldrss <- rss typ <- if (iter%%2 == b) "S" else "C" cat("Iteration (opt. ", typ, "): ", iter, ", RSS: ", rss, ", RD: ", RD, "\n", sep = "") iter <- iter + 1 oneMore <- (normS > S && (iter%%2 != b) && maxiter != 1) || (length(closureC) > 0 && (iter%%2 == b)) } cat("Initial RSS / Final RSS =", initialrss, "/", rss, "=", initialrss/rss, "\n") return(list(CList = CList, S = S, rss = rss, resid = resid, iter = iter)) } --- function search by body --- Function als in namespace ALS has this body. ----------- END OF FAILURE REPORT -------------- Fatal error: length > 1 in coercion to logical * checking PDF 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 See ‘/data/gannet/ripley/R/packages/tests-LENGTH1/ALS.Rcheck/00check.log’ for details. Command exited with non-zero status 1 Time 1:44.44, 32.92 + 6.35