cran / OptM

:exclamation: This is a read-only mirror of the CRAN R package repository. OptM — Estimating the Optimal Number of Migration Edges from 'Treemix'
2 stars 0 forks source link

optM error in R Studio #1

Closed Asrix closed 2 years ago

Asrix commented 2 years ago

Hi, I am trying to run optM on my data. folder2 has all of the output files from your suggested run: for m in {1..10}; do for i in {1..10}; do s=$RANDOM do treemix -i AC.filt.recode.treemix.frq.gz -o treemix/AC.${i}.${m} -global -m ${m} -k 1 -seed ${s}; done; done When I try to run test.optM = optM(folder2) I get the following:

Finished reading .llik, .modelcov.goz, and .cov.gz files. No output file will be saved. To save an output file, run with 'tsv = "file.tsv"' All migration edges will be included. m ranges between 0 and 10. The average number of iterations per m (not including m = 0) was 10. Make sure these values are correct... Analyzing the treemix results using the Evanno method. Error in optM(folder2) : Error: It is recommended to run more than 2 iterations for each m.

I do have 10 iterations for each m in folder2 and nothing else in that folder. R info: RStudio 2022.02.0+443 "Prairie Trillium" Release (9f7969398b90468440a501cf065295d9050bb776, 2022-02-16) for macOS Mozilla/5.0 (Macintosh; Intel Mac OS X 11_1_0) AppleWebKit/537.36 (KHTML, like Gecko) QtWebEngine/5.12.10 Chrome/69.0.3497.128 Safari/537.36 R version 4.1.2 (2021-11-01) Platform: x86_64-apple-darwin17.0 (64-bit) Running under: macOS Big Sur 11.1

Matrix products: default LAPACK: /Library/Frameworks/R.framework/Versions/4.1/Resources/lib/libRlapack.dylib

Random number generation: RNG: Mersenne-Twister Normal: Inversion Sample: Rounding

locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

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

other attached packages: [1] OptM_0.1.6 R.utils_2.11.0
[3] R.oo_1.24.0 R.methodsS3_1.8.1
[5] SeqArray_1.34.0 pegas_1.1
[7] ape_5.6-2 poppr_2.9.3
[9] adegenet_2.1.5 ade4_1.7-18
[11] vcfR_1.12.0 RColorBrewer_1.1-2
[13] VariantAnnotation_1.38.0 Rsamtools_2.8.0
[15] Biostrings_2.62.0 XVector_0.34.0
[17] SummarizedExperiment_1.22.0 Biobase_2.52.0
[19] GenomicRanges_1.46.1 GenomeInfoDb_1.30.1
[21] IRanges_2.28.0 S4Vectors_0.32.3
[23] MatrixGenerics_1.4.3 matrixStats_0.61.0
[25] BiocGenerics_0.40.0 SNPRelate_1.28.0
[27] gdsfmt_1.28.1 forcats_0.5.1
[29] stringr_1.4.0 dplyr_1.0.8
[31] purrr_0.3.4 readr_2.1.2
[33] tidyr_1.2.0 tibble_3.1.6
[35] ggplot2_3.3.5 tidyverse_1.3.1
[37] rgl_0.108.3

gaborcsardi commented 2 years ago

Hi, this is a read only mirror of CRAN, please see the package authors in the DESCRIPTION file. Look for Maintainer, BugReports and URL. Thanks!

mxt2022202059 commented 3 months ago

Hi, I'm having the same problem, what was your previous solution?