ID error in cross-sectional analysis

[CynFP]

Hi! I want to apply ALASCA for a design that does not include repeated measurements similar to the first model in example 3. To test it, I first downloaded the files of example 3 (https://figshare.com/articles/software/ALASCA_An_R_package_for_longitudinal_and_cross-sectional_analysis_of_multivariate_data_by_ASCA-based_methods/21362979/1) and when running the model (02.ex3.part1) it generated an error with the IDs, as follows:

INFO [2023-01-15 11:17:47] Initializing ALASCA (v1.0.7, 2022-12-10) WARN [2023-01-15 11:17:47] Guessing effects: disease INFO [2023-01-15 11:17:47] Will use linear models! INFO [2023-01-15 11:17:47] Will use Rfast! WARN [2023-01-15 11:17:47] Converting IDs to integer values WARN [2023-01-15 11:17:47] The disease column is used for stratification fstcore package v0.9.14 (OpenMP detected, using 20 threads) INFO [2023-01-15 11:17:48] Scaling data with sdref ... WARN [2023-01-15 11:17:48] The scaling sdref has been replaced by sdt1 as there is only one effect term. This corresponds to the column disease INFO [2023-01-15 11:17:50] Calculating LM coefficients INFO [2023-01-15 11:17:50] Reducing the number of dimensions with PCA INFO [2023-01-15 11:17:54] Keeping 111 components from initial PCA, explaining 95.11 % of variation. The limit can be changed with reduce_dimensions.limit INFO [2023-01-15 11:17:54] -> Finished the reduction of dimensions! Error in bmerge(i, x, leftcols, rightcols, roll, rollends, nomatch, mult, : Incompatible join types: x.ID (integer) and i.V1 (character)

The same error appears when trying to apply it to my data. However, when I run the script 02.ex3.part2 , example 1 and example 2, which are repeated measurements, they were processed without any problem.

How could I solve this problem? Thank you so much for your time and help.

Sincerely, Cynthia

sessionInfo() R version 4.2.2 (2022-10-31 ucrt) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 10 x64 (build 22621)

Matrix products: default

locale: [1] LC_COLLATE=Spanish_Mexico.utf8 LC_CTYPE=Spanish_Mexico.utf8
[3] LC_MONETARY=Spanish_Mexico.utf8 LC_NUMERIC=C
[5] LC_TIME=Spanish_Mexico.utf8

attached base packages: [1] stats graphics grDevices utils datasets methods base

other attached packages: [1] fstcore_0.9.14 ggrepel_0.9.2 ggplot2_3.4.0 ALASCA_1.0.7
[5] data.table_1.14.6 here_1.0.1

loaded via a namespace (and not attached): [1] Rcpp_1.0.9 log4r_0.4.3 pillar_1.8.1 compiler_4.2.2
[5] tools_4.2.2 bit_4.0.5 lifecycle_1.0.3 tibble_3.1.8
[9] gtable_0.3.1 pkgconfig_2.0.3 rlang_1.0.6 cli_3.4.1
[13] DBI_1.1.3 rstudioapi_0.14 parallel_4.2.2 duckdb_0.6.1
[17] withr_2.5.0 dplyr_1.0.10 generics_0.1.3 vctrs_0.5.1
[21] hms_1.1.2 RcppZiggurat_0.1.6 Rfast_2.0.6 rprojroot_2.0.3
[25] bit64_4.0.5 grid_4.2.2 tidyselect_1.2.0 glue_1.6.2
[29] R6_2.5.1 fansi_1.0.3 vroom_1.6.0 tzdb_0.3.0
[33] readr_2.1.3 magrittr_2.0.3 scales_1.2.1 ellipsis_0.3.2
[37] fst_0.9.8 assertthat_0.2.1 colorspace_2.0-3 utf8_1.2.2
[41] munsell_0.5.0 crayon_1.5.2

[andjar]

Hi Cynthia,

Thank you for testing ALASCA and for making me aware of this issue! Interestingly, I got another error when running example 3 part 1 (Error in crossprod(x, y) : requires numeric/complex matrix/vector arguments). I have corrected that bug in version 1.0.8, and the example now runs correctly for me. Can you please re-try with version 1.0.8 and see if the error persists?

Best, Anders

[CynFP]

Hello Anders!

I have run the example again and it works perfect with version 1.0.8. Thank you so much for your time and help.

Thanks Cynthia

[andjar]

Perfect! Thanks for the update, and please let me know if you experience more problems!

Best, Anders