Last updated: 2023-04-13

Checks: 6 1

Knit directory: Serreze-T1D_Workflow/

This reproducible R Markdown analysis was created with workflowr (version 1.7.0). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.


The R Markdown is untracked by Git. To know which version of the R Markdown file created these results, you’ll want to first commit it to the Git repo. If you’re still working on the analysis, you can ignore this warning. When you’re finished, you can run wflow_publish to commit the R Markdown file and build the HTML.

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

The command set.seed(20220210) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

Nice! There were no cached chunks for this analysis, so you can be confident that you successfully produced the results during this run.

Great job! Using relative paths to the files within your workflowr project makes it easier to run your code on other machines.

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version 51e3845. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:


Ignored files:
    Ignored:    .DS_Store
    Ignored:    analysis/.DS_Store

Untracked files:
    Untracked:  analysis/0.1.1_preparing.data_bqc_7.batches_myo.R
    Untracked:  analysis/0.1.1_preparing.data_bqc_7.batches_myo.Rmd
    Untracked:  analysis/0.1_samples_batch_20230124.Rmd
    Untracked:  analysis/0.1_samples_batch_20230225.Rmd
    Untracked:  analysis/0.1_samples_batch_20230313.Rmd
    Untracked:  analysis/0.2_haplotype_comparison_bqc_7.batches_myo_minprob.R
    Untracked:  analysis/0.2_haplotype_comparison_bqc_7.batches_myo_minprob.Rmd
    Untracked:  analysis/2.1_sample_bqc_7.batches_myo.R
    Untracked:  analysis/2.1_sample_bqc_7.batches_myo.Rmd
    Untracked:  analysis/2.1_sample_bqc_7.batches_myo.Rmd.R
    Untracked:  analysis/2.2.1_snp_qc_7.batches_myo.R
    Untracked:  analysis/2.2.1_snp_qc_7.batches_myo.Rmd
    Untracked:  analysis/2.2.1_snp_qc_7.batches_myo_mis.R
    Untracked:  analysis/2.2.1_snp_qc_7.batches_myo_mis.Rmd
    Untracked:  analysis/2.4_preparing.data_aqc_7.batches_myo.Rmd
    Untracked:  analysis/2.4_preparing.data_aqc_7.batches_myo_mis.Rmd
    Untracked:  analysis/2.4_preparing.data_aqc_7.batches_myo_misd.R
    Untracked:  analysis/3.1_phenotype.qc_corrected_7.batches_myo.Rmd
    Untracked:  analysis/3.1_phenotype.qc_corrected_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_het-ici-myo-yes.vs.het-ici-myo-no_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_het-ici-myo-yes.vs.het-ici-myo-no_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_het-ici-myo-yes.vs.het-ici-myo-no_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_het-ici.vs.het-pbs_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_het-ici.vs.het-pbs_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_het-ici.vs.het-pbs_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici-myo-yes.vs.ici-myo-no_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici-myo-yes.vs.ici-myo-no_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici-myo-yes.vs.ici-myo-no_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici-sick.vs.ici-eoi_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici-sick.vs.ici-eoi_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici-sick.vs.ici-eoi_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici.vs.pbs_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici.vs.pbs_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_ici.vs.pbs_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_myo-yes.vs.myo-no_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_myo-yes.vs.myo-no_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_myo-yes.vs.myo-no_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_pbs-myo-yes.vs.pbs-myo-no_snpsqc_dis_no-x_updated_7.batches_myo.R
    Untracked:  analysis/4.1.1_qtl.analysis_binary_pbs-myo-yes.vs.pbs-myo-no_snpsqc_dis_no-x_updated_7.batches_myo.Rmd
    Untracked:  analysis/4.1.1_qtl.analysis_binary_pbs-myo-yes.vs.pbs-myo-no_snpsqc_dis_no-x_updated_7.batches_myo_mis.Rmd
    Untracked:  analysis/debugging_comparing_old_Rmd/
    Untracked:  analysis/genotype.frequencies_het-ici-myo-yes.vs.het-ici-myo-no_7.batches_myo.Rmd
    Untracked:  analysis/genotype.frequencies_het-ici-myo-yes.vs.het-ici-myo-no_7.batches_myo_mis.Rmd
    Untracked:  analysis/genotype.frequencies_het-ici.vs.het-pbs_7.batches_myo.Rmd
    Untracked:  analysis/genotype.frequencies_het-ici.vs.het-pbs_7.batches_myo_mis.Rmd
    Untracked:  analysis/genotype.frequencies_ici-myo-yes.vs.ici-myo-no_7.batches_myo.Rmd
    Untracked:  analysis/genotype.frequencies_ici-myo-yes.vs.ici-myo-no_7.batches_myo_mis.Rmd
    Untracked:  analysis/genotype.frequencies_ici-sick.vs.ici-eoi_7.batches_myo.Rmd
    Untracked:  analysis/genotype.frequencies_ici-sick.vs.ici-eoi_7.batches_myo_mis.Rmd
    Untracked:  analysis/genotype.frequencies_ici.vs.pbs_7.batches_myo.Rmd
    Untracked:  analysis/genotype.frequencies_ici.vs.pbs_7.batches_myo_mis.Rmd
    Untracked:  analysis/genotype.frequencies_myo-yes.vs.myo-no_7.batches_myo.Rmd
    Untracked:  analysis/genotype.frequencies_pbs-myo-yes.vs.pbs-myo-no_7.batches_myo.Rmd
    Untracked:  analysis/index_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_het-ici-myo-yes.vs.het-ici-myo-no_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_het-ici-myo-yes.vs.het-ici-myo-no_7.batches_myo_mis.Rmd
    Untracked:  analysis/power.analysis_het-ici.vs.het-pbs_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_het-ici.vs.het-pbs_7.batches_myo_mis.Rmd
    Untracked:  analysis/power.analysis_ici-myo-yes.vs.ici-myo-no_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_ici-myo-yes.vs.ici-myo-no_7.batches_myo_mis.Rmd
    Untracked:  analysis/power.analysis_ici-sick.vs.ici-eoi_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_ici-sick.vs.ici-eoi_7.batches_myo_mis.Rmd
    Untracked:  analysis/power.analysis_ici.vs.pbs_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_ici.vs.pbs_7batches_myo_mis.Rmd
    Untracked:  analysis/power.analysis_myo-yes.vs.myo-no_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_myo-yes.vs.myo-no_7.batches_myo_mis.Rmd
    Untracked:  analysis/power.analysis_pbs-myo-yes.vs.pbs-myo-no_7.batches_myo.Rmd
    Untracked:  analysis/power.analysis_pbs-myo-yes.vs.pbs-myo-no_7.batches_myo_mis.Rmd
    Untracked:  data/GM_covar_4.batches_myo.csv
    Untracked:  data/GM_covar_7.batches_myo.csv
    Untracked:  data/bad_markers_all_4.batches_myo.RData
    Untracked:  data/bad_markers_all_7.batches_myo.RData
    Untracked:  data/e_4.batches_myo.RData
    Untracked:  data/e_7.batches_myo.RData
    Untracked:  data/e_snpg_samqc_4.batches_myo.RData
    Untracked:  data/e_snpg_samqc_7.batches_myo.RData
    Untracked:  data/errors_ind_4.batches_myo.RData
    Untracked:  data/errors_ind_7.batches_myo.RData
    Untracked:  data/genetic_map_4.batches_myo.csv
    Untracked:  data/genetic_map_7.batches_myo.csv
    Untracked:  data/genotype_errors_marker_4.batches_myo.RData
    Untracked:  data/genotype_errors_marker_7.batches_myo.RData
    Untracked:  data/genotype_freq_marker_4.batches_myo.RData
    Untracked:  data/genotype_freq_marker_7.batches_myo.RData
    Untracked:  data/gm_allqc_4.batches_myo.RData
    Untracked:  data/gm_allqc_4.batches_myo_mis.RData
    Untracked:  data/gm_allqc_7.batches_myo.RData
    Untracked:  data/gm_allqc_7.batches_myo_mis.RData
    Untracked:  data/gm_samqc_4.batches_myo.RData
    Untracked:  data/gm_samqc_7.batches_myo.RData
    Untracked:  data/gm_serreze.BC312.RData
    Untracked:  data/gm_serreze.BC366.RData
    Untracked:  data/gm_serreze_7.batches_myo.RData
    Untracked:  data/gm_serreze_bc_7.batches_myo.RData
    Untracked:  data/percent_missing_id_4.batches_myo.RData
    Untracked:  data/percent_missing_id_7.batches_myo.RData
    Untracked:  data/percent_missing_marker_4.batches_myo.RData
    Untracked:  data/percent_missing_marker_7.batches_myo.RData
    Untracked:  data/pheno_4.batches_myo.csv
    Untracked:  data/pheno_7.batches_myo.csv
    Untracked:  data/physical_map_4.batches_myo.csv
    Untracked:  data/physical_map_7.batches_myo.csv
    Untracked:  data/probs_36state_bc_7.batches_myo.rds
    Untracked:  data/probs_8state_bc_7.batches_myo.rds
    Untracked:  data/qc_info_bad_sample_4.batches_myo.RData
    Untracked:  data/qc_info_bad_sample_7.batches_myo.RData
    Untracked:  data/sample_geno_AHB_4.batches_myo.csv
    Untracked:  data/sample_geno_AHB_7.batches_myo.csv
    Untracked:  data/sample_geno_bc_4.batches_myo.csv
    Untracked:  data/sample_geno_bc_7.batches_myo.csv
    Untracked:  data/serreze_probs_4.batches_myo.rds
    Untracked:  data/serreze_probs_7.batches_myo.rds
    Untracked:  data/serreze_probs_allqc_4.batches_myo.rds
    Untracked:  data/serreze_probs_allqc_4.batches_myo_mis.rds
    Untracked:  data/serreze_probs_allqc_7.batches_myo.rds
    Untracked:  data/serreze_probs_allqc_7.batches_myo_mis.rds
    Untracked:  data/summary.cg_4.batches_myo.RData
    Untracked:  data/summary.cg_7.batches_myo.RData
    Untracked:  output/Percent_missing_genotype_data_4.batches_myo.pdf
    Untracked:  output/Percent_missing_genotype_data_7.batches_myo.pdf
    Untracked:  output/Percent_missing_genotype_data_per_marker_4.batches_myo.pdf
    Untracked:  output/Percent_missing_genotype_data_per_marker_7.batches_myo.pdf
    Untracked:  output/Proportion_matching_genotypes_before_removal_of_bad_samples_4.batches_myo.pdf
    Untracked:  output/Proportion_matching_genotypes_before_removal_of_bad_samples_7.batches_myo.pdf
    Untracked:  output/genotype_error_marker_4.batches_myo.pdf
    Untracked:  output/genotype_error_marker_7.batches_myo.pdf
    Untracked:  output/genotype_frequency_marker_4.batches_myo.pdf
    Untracked:  output/genotype_frequency_marker_7.batches_myo.pdf

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.


There are no past versions. Publish this analysis with wflow_publish() to start tracking its development.


Genotype Summary by Sample

Neogen_Sample_ID Missing AA AC AG AT CC CG GC GG TA TC TG TT Total_Genotyped Percent_Genotyped MURGIGV01_SNP.Count cr
D928-ICI-Myo 5355 29098 1897 7859 11 30153 8 6 30454 6 7741 1792 28879 137904 96.26% 143259 0.962620149519402
D266-PBS-EOI 5327 28883 2078 8242 6 29861 5 4 30129 1 8064 1939 28720 137932 96.28% 143259 0.962815599717993
D947-PBS-EOI 5312 29644 1725 6939 6 30660 6 6 30960 1 6856 1759 29385 137947 96.29% 143259 0.962920305181524
D268-PBS-EOI 5321 30037 1520 6316 3 31105 9 2 31409 6 6238 1504 29789 137938 96.29% 143259 0.962857481903406
D948-PBS-EOI 5293 30376 1443 5815 3 31312 3 5 31702 2 5744 1397 30164 137966 96.31% 143259 0.963052932101997
D271-PBS-EOI 5256 28019 2301 9672 7 29157 8 5 29314 4 9325 2248 27943 138003 96.33% 143259 0.963311205578707
D1706-ICI-Myo 5290 30287 1417 5824 6 31406 5 4 31720 2 5721 1447 30130 137969 96.31% 143259 0.963073873194703
D307-ICI-EOI 5282 29700 1702 6816 11 30733 9 8 31091 3 6694 1599 29611 137977 96.31% 143259 0.963129716108587
D1708-ICI-Myo 5245 30442 1400 5737 8 31452 6 4 31734 0 5566 1361 30304 138014 96.34% 143259 0.963387989585297
D308-ICI-EOI 5263 30105 1568 6195 6 31030 4 5 31430 3 6218 1490 29942 137996 96.33% 143259 0.963262343029059
D1737-ICI-Myo 5233 30084 1552 6305 7 31050 9 9 31393 6 6262 1462 29887 138026 96.35% 143259 0.963471753956121
D322-PBS-EOI 5217 29406 1750 7424 9 30458 4 4 30733 2 7286 1720 29246 138042 96.36% 143259 0.963583439783888
D1731-ICI-Myo 5338 30217 1505 6000 5 31143 6 7 31549 2 5968 1491 30028 137921 96.27% 143259 0.962738815711404
D1739-ICI-Myo 5376 31119 1123 4567 5 32071 4 4 32416 1 4533 1071 30969 137883 96.25% 143259 0.962473561870458
D1771-ICI-Myo 5329 29744 1668 6679 6 30844 5 4 31221 4 6581 1536 29638 137930 96.28% 143259 0.962801638989522
D1774-ICI-Myo 5330 30321 1403 5918 7 31382 5 3 31596 1 5774 1434 30085 137929 96.28% 143259 0.962794658625287
7451-ICI-EOI 5299 30207 1587 6126 3 31139 4 5 31450 3 5934 1453 30049 137960 96.30% 143259 0.963011049916585
D249-PBS-EOI 5264 29060 1924 7934 9 30101 11 4 30367 4 7782 1862 28937 137995 96.33% 143259 0.963031991009291
D249-PBS-EOI 5264 29060 1924 7934 9 30101 11 4 30367 4 7782 1862 28937 137995 96.33% 143259 0.963255362664824
D249-PBS-EOI 5264 29060 1924 7934 9 30101 11 4 30367 4 7782 1862 28937 137995 96.33% 143259 0.963255362664824
D249-PBS-EOI 5264 29060 1924 7934 9 30101 11 4 30367 4 7782 1862 28937 137995 96.33% 143259 0.963031991009291

Gender Summary by Sample

Neogen_Sample_ID Provided.Sex Inferred.Sex
D928-ICI-Myo F F
D266-PBS-EOI F F
D947-PBS-EOI F F
D268-PBS-EOI F F
D948-PBS-EOI F F
D271-PBS-EOI F F
D1706-ICI-Myo F F
D307-ICI-EOI F F
D1708-ICI-Myo F F
D308-ICI-EOI F F
D1737-ICI-Myo F F
D322-PBS-EOI F F
D1731-ICI-Myo F F
D1739-ICI-Myo F F
D1771-ICI-Myo F F
D1774-ICI-Myo F F
7451-ICI-EOI F F
D249-PBS-EOI F F
D249-PBS-EOI F F
D249-PBS-EOI F F
D249-PBS-EOI F F

Phenotype Summary by Sample

Neogen_Sample_ID sex_pheno Inferred.Sex_pheno Myocarditis.Status age.of.clinical.phenotype..weeks.
D928-ICI-Myo F F YES
D266-PBS-EOI F F NO
D947-PBS-EOI F F YES
D268-PBS-EOI F F NO
D948-PBS-EOI F F YES
D271-PBS-EOI F F NO
D1706-ICI-Myo F F YES 13
D307-ICI-EOI F F NO
D1708-ICI-Myo F F YES 13
D308-ICI-EOI F F NO
D1737-ICI-Myo F F YES 13
D322-PBS-EOI F F NO
D1731-ICI-Myo F F YES 13
D1739-ICI-Myo F F YES 13
D1771-ICI-Myo F F YES 13
D1774-ICI-Myo F F YES 13
7451-ICI-EOI F F NO
D249-PBS-EOI F F NO
D249-PBS-EOI F F NO
D249-PBS-EOI F F NO
D249-PBS-EOI F F NO

All Problematic Samples

Neogen_Sample_ID no_pheno low_call.rate different_sex

R version 4.2.2 (2022-10-31)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Big Sur ... 10.16

Matrix products: default
BLAS:   /Library/Frameworks/R.framework/Versions/4.2/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/4.2/Resources/lib/libRlapack.dylib

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] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] tibble_3.1.8      readxl_1.4.1      cluster_2.1.4     dplyr_1.0.10     
 [5] optparse_1.7.3    rhdf5_2.40.0      mclust_6.0.0      tidyr_1.2.1      
 [9] data.table_1.14.6 knitr_1.41        kableExtra_1.3.4  workflowr_1.7.0  

loaded via a namespace (and not attached):
 [1] tidyselect_1.2.0   xfun_0.35          bslib_0.4.1        purrr_0.3.5       
 [5] colorspace_2.0-3   vctrs_0.5.1        generics_0.1.3     getopt_1.20.3     
 [9] htmltools_0.5.3    viridisLite_0.4.1  yaml_2.3.6         utf8_1.2.2        
[13] rlang_1.0.6        jquerylib_0.1.4    later_1.3.0        pillar_1.8.1      
[17] withr_2.5.0        DBI_1.1.3          glue_1.6.2         lifecycle_1.0.3   
[21] stringr_1.5.0      cellranger_1.1.0   munsell_0.5.0      rvest_1.0.3       
[25] evaluate_0.18      callr_3.7.3        fastmap_1.1.0      httpuv_1.6.6      
[29] ps_1.7.2           fansi_1.0.3        Rcpp_1.0.9         promises_1.2.0.1  
[33] scales_1.2.1       cachem_1.0.6       webshot_0.5.4      jsonlite_1.8.4    
[37] fs_1.5.2           systemfonts_1.0.4  digest_0.6.30      stringi_1.7.8     
[41] processx_3.8.0     getPass_0.2-2      rprojroot_2.0.3    rhdf5filters_1.8.0
[45] cli_3.4.1          tools_4.2.2        magrittr_2.0.3     sass_0.4.4        
[49] whisker_0.4.1      pkgconfig_2.0.3    ellipsis_0.3.2     xml2_1.3.3        
[53] assertthat_0.2.1   rmarkdown_2.18     svglite_2.1.0      httr_1.4.4        
[57] rstudioapi_0.14    Rhdf5lib_1.18.2    R6_2.5.1           git2r_0.30.1      
[61] compiler_4.2.2