R Under development (unstable) (2025-09-01 r88761 ucrt) -- "Unsuffered Consequences" 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(albatross) > > x <- feem(matrix(1:42, 6), 1:6, 1:7) > feemscale(x) excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scale") [1] 12.26784 attr(,"class") [1] "feem" > > y <- replicate(2, x, FALSE) > feemscale(y) [[1]] excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scale") [1] 12.26784 attr(,"class") [1] "feem" [[2]] excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scale") [1] 12.26784 attr(,"class") [1] "feem" > > z <- feemcube(y, TRUE) > feemscale(z) , , 1 excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 , , 2 excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scales") [1] 12.26784 12.26784 attr(,"class") [1] "feemcube" > > z[1, 1, ] <- NA > tools::assertError(feemscale(z), verbose = TRUE) Asserted error: While processing "1": Scale function returned NA: does it need na.rm = TRUE? > feemscale(z, multiway::sumsq, na.rm = TRUE) , , 1 excitation emission 1 2 3 4 5 1 NA 0.0002736085 0.0005081301 0.0007426517 0.0009771732 2 7.817386e-05 0.0003126954 0.0005472170 0.0007817386 0.0010162602 3 1.172608e-04 0.0003517824 0.0005863039 0.0008208255 0.0010553471 4 1.563477e-04 0.0003908693 0.0006253909 0.0008599124 0.0010944340 5 1.954346e-04 0.0004299562 0.0006644778 0.0008989994 0.0011335210 6 2.345216e-04 0.0004690432 0.0007035647 0.0009380863 0.0011726079 excitation emission 6 7 1 0.001211695 0.001446216 2 0.001250782 0.001485303 3 0.001289869 0.001524390 4 0.001328956 0.001563477 5 0.001368043 0.001602564 6 0.001407129 0.001641651 , , 2 excitation emission 1 2 3 4 5 1 NA 0.0002736085 0.0005081301 0.0007426517 0.0009771732 2 7.817386e-05 0.0003126954 0.0005472170 0.0007817386 0.0010162602 3 1.172608e-04 0.0003517824 0.0005863039 0.0008208255 0.0010553471 4 1.563477e-04 0.0003908693 0.0006253909 0.0008599124 0.0010944340 5 1.954346e-04 0.0004299562 0.0006644778 0.0008989994 0.0011335210 6 2.345216e-04 0.0004690432 0.0007035647 0.0009380863 0.0011726079 excitation emission 6 7 1 0.001211695 0.001446216 2 0.001250782 0.001485303 3 0.001289869 0.001524390 4 0.001328956 0.001563477 5 0.001368043 0.001602564 6 0.001407129 0.001641651 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scales") [1] 25584 25584 attr(,"class") [1] "feemcube" > > feemscale(y, progress = TRUE) ================================================================================ [[1]] excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scale") [1] 12.26784 attr(,"class") [1] "feem" [[2]] excitation emission 1 2 3 4 5 6 7 1 0.08151391 0.5705974 1.059681 1.548764 2.037848 2.526931 3.016015 2 0.16302783 0.6521113 1.141195 1.630278 2.119362 2.608445 3.097529 3 0.24454174 0.7336252 1.222709 1.711792 2.200876 2.689959 3.179043 4 0.32605566 0.8151391 1.304223 1.793306 2.282390 2.771473 3.260557 5 0.40756957 0.8966531 1.385737 1.874820 2.363904 2.852987 3.342070 6 0.48908349 0.9781670 1.467250 1.956334 2.445417 2.934501 3.423584 attr(,"emission") [1] 1 2 3 4 5 6 attr(,"excitation") [1] 1 2 3 4 5 6 7 attr(,"scale") [1] 12.26784 attr(,"class") [1] "feem" > > proc.time() user system elapsed 0.35 0.17 0.45