Last updated: 2022-04-25

Checks: 2 0

Knit directory: synovialscrnaseq/

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


Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

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 3443cc6. 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:    '/
    Ignored:    .Rhistory
    Ignored:    .Rproj.user/
    Ignored:    .empty/
    Ignored:    analysis/.Rhistory
    Ignored:    analysis/iSEE_interactive_document.html
    Ignored:    code/test_files/
    Ignored:    data/Culemann/
    Ignored:    data/E-MTAB-8322/
    Ignored:    data/Synovial scRNA-seq samples - Sheet1.csv
    Ignored:    data/Zhang_top20_singlecell_cluster_markers_fromGithub.csv
    Ignored:    data/findMarkers_results.rds
    Ignored:    data/findMarkers_results_v2.rds
    Ignored:    data/info/
    Ignored:    data/syn_sce_tidy_filtered.rds
    Ignored:    data/syn_sce_tidy_hvg.rds
    Ignored:    data/syn_sce_tidy_hvg_cms.rds
    Ignored:    docs/
    Ignored:    output/Figures_Paper/
    Ignored:    output/Sample_summaries_RA_comparisons.rds
    Ignored:    output/Sample_summaries_direct_dissociation.rds
    Ignored:    output/Sample_summaries_exvivo_treatment.rds
    Ignored:    output/count_matrix_unfiltered.mtx
    Ignored:    output/emptyDrops_result_v4.rds
    Ignored:    output/emptyDrops_result_v4_tmp.rds
    Ignored:    output/entropies_fstat_v4.rds
    Ignored:    output/entropies_fstat_v4_ec.rds
    Ignored:    output/entropies_fstat_v4_main.rds
    Ignored:    output/entropies_fstat_v4_mp.rds
    Ignored:    output/entropies_fstat_v4_sf.rds
    Ignored:    output/entropies_fstat_v4_tc.rds
    Ignored:    output/entropies_fstat_v5_ec.rds
    Ignored:    output/entropies_fstat_v5_main.rds
    Ignored:    output/entropies_fstat_v5_mp.rds
    Ignored:    output/entropies_fstat_v5_sf.rds
    Ignored:    output/entropies_fstat_v5_tc.rds
    Ignored:    output/findMarkers_results_v4.rds
    Ignored:    output/findMarkers_results_v4_ec.rds
    Ignored:    output/findMarkers_results_v4_main.rds
    Ignored:    output/findMarkers_results_v4_mp.rds
    Ignored:    output/findMarkers_results_v4_sf.rds
    Ignored:    output/findMarkers_results_v4_tc.rds
    Ignored:    output/findMarkers_results_v5_ec.rds
    Ignored:    output/findMarkers_results_v5_main.rds
    Ignored:    output/findMarkers_results_v5_mp.rds
    Ignored:    output/findMarkers_results_v5_sf.rds
    Ignored:    output/findMarkers_results_v5_tc.rds
    Ignored:    output/findMarkers_results_v6.rds
    Ignored:    output/findMarkers_results_v6_ec.rds
    Ignored:    output/findMarkers_results_v6_main.rds
    Ignored:    output/findMarkers_results_v6_mp.rds
    Ignored:    output/findMarkers_results_v6_sf.rds
    Ignored:    output/findMarkers_results_v6_tc.rds
    Ignored:    output/goana_results_v6_ec.rds
    Ignored:    output/goana_results_v6_mp.rds
    Ignored:    output/syn_v4/
    Ignored:    output/syn_v4_clustering_lookup.rds
    Ignored:    output/syn_v4_clustering_lookup_invivo.rds
    Ignored:    output/syn_v4_conos_obj.rds
    Ignored:    output/syn_v4_neiprop.rds
    Ignored:    output/syn_v4_neiprop_invivo.rds
    Ignored:    output/syn_v4_predref.rds
    Ignored:    output/syn_v4_predref_invivo.rds
    Ignored:    output/syn_v4_res_da_1.rds
    Ignored:    output/syn_v4_res_da_2.rds
    Ignored:    output/syn_v4_res_da_2_Accute_inflammation.rds
    Ignored:    output/syn_v4_res_da_2_Diagnosis.rds
    Ignored:    output/syn_v4_res_da_2_Diagnosis_main.rds
    Ignored:    output/syn_v4_res_da_2_Lymphoid_folicles.rds
    Ignored:    output/syn_v4_res_da_2_Pathotype.rds
    Ignored:    output/syn_v4_res_da_2_Sample_prep.rds
    Ignored:    output/syn_v4_res_da_2_Therapy.rds
    Ignored:    output/syn_v4_res_da_2_Treatment.rds
    Ignored:    output/syn_v4_res_da_2_Vascularisation_bin.rds
    Ignored:    output/syn_v4_res_da_Accute_inflammation_invivo.rds
    Ignored:    output/syn_v4_res_da_Diagnosis_invivo.rds
    Ignored:    output/syn_v4_res_da_Diagnosis_main_invivo.rds
    Ignored:    output/syn_v4_res_da_Lymphoid_folicles_invivo.rds
    Ignored:    output/syn_v4_res_da_Pathotype_invivo.rds
    Ignored:    output/syn_v4_res_da_Therapy_invivo.rds
    Ignored:    output/syn_v4_res_da_Vascularisation_bin_invivo.rds
    Ignored:    output/syn_v4_res_ds_2.rds
    Ignored:    output/syn_v4_res_ds_2_Accute_inflammation.rds
    Ignored:    output/syn_v4_res_ds_2_Diagnosis.rds
    Ignored:    output/syn_v4_res_ds_2_Diagnosis_main.rds
    Ignored:    output/syn_v4_res_ds_2_Lymphoid_folicles.rds
    Ignored:    output/syn_v4_res_ds_2_Pathotype.rds
    Ignored:    output/syn_v4_res_ds_2_Sample_prep.rds
    Ignored:    output/syn_v4_res_ds_2_Therapy.rds
    Ignored:    output/syn_v4_res_ds_2_Treatment.rds
    Ignored:    output/syn_v4_res_ds_2_Vascularisation_bin.rds
    Ignored:    output/syn_v4_res_ds_Accute_inflammation_invivo.rds
    Ignored:    output/syn_v4_res_ds_Diagnosis_invivo.rds
    Ignored:    output/syn_v4_res_ds_Diagnosis_main_invivo.rds
    Ignored:    output/syn_v4_res_ds_Lymphoid_folicles_invivo.rds
    Ignored:    output/syn_v4_res_ds_Pathotype_invivo.rds
    Ignored:    output/syn_v4_res_ds_Therapy_invivo.rds
    Ignored:    output/syn_v4_res_ds_Vascularisation_bin_invivo.rds
    Ignored:    output/syn_v4_sce.rds
    Ignored:    output/syn_v4_sce_annot_df_invivo.csv
    Ignored:    output/syn_v4_sce_annot_df_invivo.rds
    Ignored:    output/syn_v4_sce_conos.rds
    Ignored:    output/syn_v4_sce_ec.rds
    Ignored:    output/syn_v4_sce_ec_invivo.rds
    Ignored:    output/syn_v4_sce_emptyDrops_filt_invivo.rds
    Ignored:    output/syn_v4_sce_emptyDrops_invivo.rds
    Ignored:    output/syn_v4_sce_filtered.rds
    Ignored:    output/syn_v4_sce_filtered_invivo.rds
    Ignored:    output/syn_v4_sce_hdf5/
    Ignored:    output/syn_v4_sce_hvg.rds
    Ignored:    output/syn_v4_sce_hvg_cms.rds
    Ignored:    output/syn_v4_sce_hvg_cms_doublet.rds
    Ignored:    output/syn_v4_sce_hvg_cms_doublet_annot.rds
    Ignored:    output/syn_v4_sce_hvg_cms_doublet_annot_RA_only.rds
    Ignored:    output/syn_v4_sce_hvg_cms_doublet_annot_invivo.rds
    Ignored:    output/syn_v4_sce_hvg_cms_doublet_annot_manual_invivo.rds
    Ignored:    output/syn_v4_sce_hvg_cms_doublet_invivo.rds
    Ignored:    output/syn_v4_sce_hvg_invivo.rds
    Ignored:    output/syn_v4_sce_mp.rds
    Ignored:    output/syn_v4_sce_mp_invivo.rds
    Ignored:    output/syn_v4_sce_sf.rds
    Ignored:    output/syn_v4_sce_sf_invivo.rds
    Ignored:    output/syn_v4_sce_tc.rds
    Ignored:    output/syn_v4_sce_tc_invivo.rds
    Ignored:    output/syn_v4_swappedDrops_24300_after.rds
    Ignored:    output/syn_v4_swappedDrops_24300_before.rds
    Ignored:    output/syn_v4_swappedDrops_24555_before.rds
    Ignored:    output/syn_v4_swappedDrops_24793_after.rds
    Ignored:    output/syn_v4_swappedDrops_24793_before.rds
    Ignored:    output/syn_v4_swappedDrops_after.rds
    Ignored:    output/syn_v4_swappedDrops_before.rds
    Ignored:    output/syn_v5_annot_df_manual.rds
    Ignored:    output/syn_v5_cluster_cellid_match_invivo.rds
    Ignored:    output/syn_v5_clustering_lookup_invivo.rds
    Ignored:    output/syn_v5_clustering_lookup_multiple_invivo.rds
    Ignored:    output/syn_v5_res_da_Accute_inflammation_invivo.rds
    Ignored:    output/syn_v5_res_da_Diagnosis_invivo.rds
    Ignored:    output/syn_v5_res_da_Diagnosis_main_invivo.rds
    Ignored:    output/syn_v5_res_da_Lymphoid_folicles_invivo.rds
    Ignored:    output/syn_v5_res_da_Pathotype_invivo.rds
    Ignored:    output/syn_v5_res_da_Therapy_invivo.rds
    Ignored:    output/syn_v5_res_da_Vascularisation_bin_invivo.rds
    Ignored:    output/syn_v5_res_ds_Accute_inflammation_invivo.rds
    Ignored:    output/syn_v5_res_ds_Diagnosis_invivo.rds
    Ignored:    output/syn_v5_res_ds_Diagnosis_main_invivo.rds
    Ignored:    output/syn_v5_res_ds_Lymphoid_folicles_invivo.rds
    Ignored:    output/syn_v5_res_ds_Pathotype_invivo.rds
    Ignored:    output/syn_v5_res_ds_Therapy_invivo.rds
    Ignored:    output/syn_v5_res_ds_Vascularisation_bin_invivo.rds
    Ignored:    output/syn_v5_sce.rds
    Ignored:    output/syn_v5_sce_ec_invivo.rds
    Ignored:    output/syn_v5_sce_filtered_invivo.rds
    Ignored:    output/syn_v5_sce_hdf5/
    Ignored:    output/syn_v5_sce_hvg_cms_doublet_annot_manual_invivo.rds
    Ignored:    output/syn_v5_sce_hvg_cms_doublet_cmstest_invivo.rds
    Ignored:    output/syn_v5_sce_hvg_cms_doublet_invivo.rds
    Ignored:    output/syn_v5_sce_hvg_cms_doublet_subcluster_invivo.rds
    Ignored:    output/syn_v5_sce_hvg_invivo.rds
    Ignored:    output/syn_v5_sce_mp_invivo.rds
    Ignored:    output/syn_v5_sce_sf_invivo.rds
    Ignored:    output/syn_v5_sce_tc_invivo.rds
    Ignored:    output/syn_v5_vst_out_invivo.rds
    Ignored:    output/syn_v6_cluster_cellid_match_invivo.rds
    Ignored:    output/syn_v6_clustering_lookup_invivo.rds
    Ignored:    output/syn_v6_clustering_lookup_multiple_invivo.rds
    Ignored:    output/syn_v6_sce.rds
    Ignored:    output/syn_v6_sce_ec_invivo.rds
    Ignored:    output/syn_v6_sce_filtered_invivo.rds
    Ignored:    output/syn_v6_sce_hdf5/
    Ignored:    output/syn_v6_sce_hvg_cms_doublet_invivo.rds
    Ignored:    output/syn_v6_sce_hvg_cms_doublet_subcluster_invivo.rds
    Ignored:    output/syn_v6_sce_hvg_invivo.rds
    Ignored:    output/syn_v6_sce_hvg_marker_genes.rds
    Ignored:    output/syn_v6_sce_mp_invivo.rds
    Ignored:    output/syn_v6_sce_sf_invivo.rds
    Ignored:    output/syn_v6_sce_tc_invivo.rds
    Ignored:    output/syn_v6_sfig1.rds
    Ignored:    output/syn_v6_vst_out_invivo.rds

Untracked files:
    Untracked:  analysis/test.Rmd
    Untracked:  code/rebuild_ezRun.R
    Untracked:  nonhosted_public/
    Untracked:  singRstudio.sh.bak

Unstaged changes:
    Modified:   README.md
    Modified:   analysis/index.Rmd
    Modified:   analysis/scRNAseq_complete_04-2_celltype_markers.Rmd
    Modified:   analysis/scRNAseq_complete_04-2_celltype_markers_subcelltypes.Rmd
    Modified:   analysis/scRNAseq_complete_Figures.Rmd
    Modified:   code/create_hdf5.R
    Modified:   data/syn_sce_tidy_hvg_cms_doublet.rds

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.


These are the previous versions of the repository in which changes were made to the R Markdown (analysis/about.Rmd) and HTML (public/about.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File Version Author Date Message
html 3443cc6 Reto Gerber 2022-04-25 Update
Rmd e88c23e Reto Gerber 2021-07-12 add ambient RNA analysis
html e88c23e Reto Gerber 2021-07-12 add ambient RNA analysis
html a28e4c4 retogerber 2021-01-29 update website
html 04c48dc retogerber 2021-01-29 add main website files
Rmd 0e7eabd retogerber 2021-01-05 Start workflowr project.

Scripts

The complete analysis was run using singularity. Install instructions can be found here. To recreate the analysis first download the container:

singularity pull library://retogerber/default/bioc_python_synovial_v2.sif

and then run the following (adjustments needed, read below)

R_LIBS_USER='/usr/local/lib/R/site-library,/usr/local/lib/R/library'
export R_LIBS_USER
SINGULARITY_BIND='PROJECT_DIR:PROJECT_DIR,RAWDATA_DIR/data:PROJECT_DIR/data_server,TMP_DIR/tmp:TMP_DIR/tmp,TMP_DIR/run:/run' 
export SINGULARITY_BIND
singularity exec bioc_python_synovial_v2.sif R -e "workflowr::wflow_build()"

where R_LIBS_USER makes sure that only R libraries within the container are used and SINGULARITY_BIND makes all the binds between the host and the container. Replace PROJECT_DIR with the directory of the project (e.g. ~/protocol_synovial) and RAWDATA_DIR with the directory of the raw data, which contains a directory called data containing the raw data. TMP_DIR is a temporary directory to allow R to function properly and needs two directories within, tmp and run. Those can be created with mkdir -p TMP_DIR/{tmp,run}.

To run individual scripts run

singularity exec bioc_python_synovial_v2.sif R -e "workflowr::wflow_build(SCRIPT_NAME)"

replacing SCRIPT_NAME with the name of the script to run.

To speed the analysis up the number of cores was set between 10-20, if you want to change that number adjust the variable n_workers at the top within each script keeping in mind that this might influence the results slightly (random seed).

The Singularity build file (for the image) can be found here.

Data