Last updated: 2021-07-01

Checks: 7 0

Knit directory: ebpmf_data_analysis/

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 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(20200511) 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 a87c087. 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:    .Rproj.user/
    Ignored:    analysis/fastTopics_on_sla_cache/
    Ignored:    analysis/investigate_F_proj_cache/
    Ignored:    analysis/prepare_sla_sim_cache/

Untracked files:
    Untracked:  Rplot.pdf
    Untracked:  analysis/droplet_simulated_0.10_compare_mle_eb.Rmd
    Untracked:  analysis/fastTopics_on_droplet_sim.Rmd
    Untracked:  analysis/fastTopics_on_sla_sim.Rmd
    Untracked:  analysis/investigate_F_proj_cache_old/
    Untracked:  analysis/prepare_droplet.Rmd
    Untracked:  analysis/prepare_droplet_sim.Rmd
    Untracked:  analysis/prepare_sim_data.Rmd
    Untracked:  code/hpmf.R
    Untracked:  data/GSE103354_Trachea_droplet_UMIcounts.txt.gz
    Untracked:  data/SLA/SLA_preprocessed_minzhe.RData
    Untracked:  data/SLA/docword.sla_minzhe.txt
    Untracked:  data/SLA/vocab.sla_minzhe.txt
    Untracked:  data/pbmc_68k/
    Untracked:  data/sim/droplet_simulated.Rds
    Untracked:  data/sim/droplet_simulated_0.10.Rds
    Untracked:  data/sim/droplet_simulated_0.20.Rds
    Untracked:  data/sim/pbmc_68k_simulated.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter1000.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter1500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter2500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter3000.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter3500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter4000.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter4500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_ebpmf_wbg_initLF50_K20_maxiter500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_init_nmf_K20_iter50.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter1000.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter1500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter2500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter3000.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter3500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter4000.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter4500.Rds
    Untracked:  output/NIPS/v0.4.5/nips_pmf_initLF50_K20_maxiter500.Rds
    Untracked:  output/News/v0.4.5/news_ebpmf_wbg_initLF50_K20_maxiter1000.Rds
    Untracked:  output/News/v0.4.5/news_ebpmf_wbg_initLF50_K20_maxiter1500.Rds
    Untracked:  output/News/v0.4.5/news_ebpmf_wbg_initLF50_K20_maxiter500.Rds
    Untracked:  output/News/v0.4.5/news_init_nmf_K20_iter50.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter1000.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter1500.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter2500.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter3000.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter3500.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter4000.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter4500.Rds
    Untracked:  output/News/v0.4.5/news_pmf_initLF50_K20_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/Rplots.pdf
    Untracked:  output/SLA/v0.4.5/correct_kl.R
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K100_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K100_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K100_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K100_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K10_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K50_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_ebpmf_wbg_initLF50_K5_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_init_nmf_K100_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_init_nmf_K10_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_init_nmf_K50_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_init_nmf_K5_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K100_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K10_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K50_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_ebpmf_wbg_initLF50_K5_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_init_nmf_K100_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_init_nmf_K10_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_init_nmf_K3_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_init_nmf_K50_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_init_nmf_K5_iter50.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K100_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K10_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K3_maxiter30.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K50_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_minzhe_pmf_initLF50_K5_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K100_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K10_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K50_maxiter500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter1000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter1500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter2500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter3000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter3500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter4000.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter4500.Rds
    Untracked:  output/SLA/v0.4.5/sla_pmf_initLF50_K5_maxiter500.Rds
    Untracked:  output/fastTopics_fit/fit__fastTopics_k.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_data_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k10.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k11.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k12.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k13.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k14.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k15.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k16.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k17.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k18.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k19.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k2.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k20.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k3.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k4.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k5.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k6.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k7.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k8.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_0.10_ebpmf_wbg_init_truth_k9.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k10.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k11.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k12.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k13.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k14.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k15.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k16.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k17.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k18.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k19.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k2.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k20.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k3.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k4.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k5.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k6.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k7.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k8.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.10_fastTopics_init_truth_k9.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k10.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k11.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k12.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k13.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k14.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k15.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k16.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k17.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k18.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k19.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k2.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k20.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k3.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k4.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k5.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k6.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k7.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k8.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted__0.20_fastTopics_init_truth_k9.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_ebpmf_wbg_k10_init_truth.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_evals.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k10.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k10_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k11.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k11_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k12.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k12_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k13.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k13_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k14.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k14_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k15.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k15_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k16.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k16_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k17.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k17_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k18.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k18_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k19.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k19_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k2.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k20.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k20_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k2_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k3.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k3_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k4.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k4_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k5.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k5_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k6.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k6_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k7.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k7_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k8.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k8_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k9.Rds.out
    Untracked:  output/fastTopics_fit/fit_droplet_simualted_fastTopics_init_truth_k9_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.10_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_0.20_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_droplet_simulated_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_pbmc_68k_simulated_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_sla_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_evals.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k10.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k10.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k11.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k11.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k12.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k12.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k13.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k13.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k13_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k14.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k14.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k15.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k15.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k16.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k16.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k17.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k17.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k18.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k18.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k19.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k19.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k2.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k2.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k20.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k20.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k3.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k3.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k4.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k4.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k5.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k5.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k6.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k6.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k7.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k7.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k8.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k8.Rout
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k9.Rds
    Untracked:  output/fastTopics_fit/fit_sla_sim_fastTopics_k9.Rout
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k10.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k10.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k11.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k11.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k12.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k12.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k13.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k13.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k14.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k14.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k15.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k15.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k16.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k16.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k17.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k17.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k18.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k18.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k19.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k19.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k2.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k2.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k20.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k20.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k3.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k3.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k4.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k4.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k5.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k5.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k6.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k6.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k7.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k7.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k8.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k8.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k9.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k9.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_evals.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k10.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k10.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k10_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k11.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k11.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k11_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k12.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k12.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k12_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k13.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k13.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k13_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k14.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k14.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k14_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k15.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k15.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k15_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k16.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k16.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k16_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k17.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k17.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k17_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k18.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k18.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k18_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k19.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k19.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k19_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k2.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k2.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k20.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k20.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k20_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k2_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k3.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k3.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k3_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k4.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k4.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k4_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k5.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k5.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k5_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k6.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k6.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k6_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k7.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k7.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k7_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k8.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k8.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k8_comparison.pdf
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k9.Rds
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k9.Rds.out
    Untracked:  output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k9_comparison.pdf
    Untracked:  output/fit_kos_coneNMF_F.Rout
    Untracked:  output/sim/v0.4.5/exper1/fit_sim_bg_block_n1100_p2100_K50_pmf_maxiter_5000_from_truth.Rout
    Untracked:  output/sim/v0.4.5/exper1/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper1/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter10_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper1/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter2000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper1/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter3000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper1/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter4000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper1/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter5_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000_from_truth.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000_pmf_bg_K50_maxiter1000_from_truth.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000_pmf_bg_K50_maxiter10_from_truth_scaled0.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_pmf_bg_maxiter_5000_from_truth.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_pmf_maxiter_5000.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_pmf_maxiter_5000_from_truth.Rout
    Untracked:  output/sim/v0.4.5/exper2/fit_sim_bg_block_n1100_p2100_K50_pmf_maxiter_5000_from_truth.Rout_old
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1000.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1000_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1000_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1000_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter100_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter10_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1500.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1500_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1500_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1500_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter1_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2000_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2000_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2000_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2500.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2500_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2500_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter2500_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3000.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3000_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3000_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3000_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter300_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3500.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3500_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3500_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3500_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4000.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4000_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4000_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4000_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter400_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4500.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4500_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4500_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter4500_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter500.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter5000_pmf_bg_K50_maxiter10_from_truth_scaled0.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter5000_pmf_bg_K50_maxiter10_from_truth_scaled1.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter5_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter2000.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter2000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter3000.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter3000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter3_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter4000.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_K50_maxiter4000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter1500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter2000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter2500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter3000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter3500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter4000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter4500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper2/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter5000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000.Rout
    Untracked:  output/sim/v0.4.5/exper3/fit_sim_bg_block_n1100_p2100_K50_pmf_bg_maxiter_5000_from_truth.Rout
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter10_init_random2.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter20_init_random2.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter1000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter1500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter2000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter2500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter3000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter3500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter4000_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter4500_from_truth.Rds
    Untracked:  output/sim/v0.4.5/exper3/sim_bg_block_n1100_p2100_K50_pmf_bg_K50_maxiter500_from_truth.Rds
    Untracked:  script/ConeFactor_kos_K20_2.R
    Untracked:  script/Rplots.pdf
    Untracked:  script/droplet.RData
    Untracked:  script/eval_fit_simulated.R
    Untracked:  script/eval_fit_simulated.sh
    Untracked:  script/eval_fit_simulated_plot.R
    Untracked:  script/eval_fit_simulated_plot.sh
    Untracked:  script/fit_droplet_fastTopics.R
    Untracked:  script/fit_droplet_fastTopics.sbatch
    Untracked:  script/fit_droplet_fastTopics.sh
    Untracked:  script/fit_kos_coneNMF_F.R
    Untracked:  script/fit_kos_coneNMF_F.sbatch
    Untracked:  script/fit_news_ebpmf_wbg.sbatch
    Untracked:  script/fit_news_ebpmf_wbg_initLF.R
    Untracked:  script/fit_news_pmf.sbatch
    Untracked:  script/fit_news_pmf_initLF.R
    Untracked:  script/fit_news_pmf_initialization.R
    Untracked:  script/fit_nips_ebpmf_wbg.sbatch
    Untracked:  script/fit_nips_ebpmf_wbg_initLF.R
    Untracked:  script/fit_nips_pmf.sbatch
    Untracked:  script/fit_nips_pmf_initLF.R
    Untracked:  script/fit_nips_pmf_initialization.R
    Untracked:  script/fit_pbmc_68k_fastTopics.R
    Untracked:  script/fit_pbmc_68k_fastTopics.sbatch
    Untracked:  script/fit_pbmc_68k_fastTopics.sh
    Untracked:  script/fit_sim_ebpmf_wbg2.sbatch
    Untracked:  script/fit_sim_pmf_bg_from_truth.R
    Untracked:  script/fit_sim_pmf_bg_from_truth.sbatch
    Untracked:  script/fit_sim_pmf_bg_from_truth_exper3.sbatch
    Untracked:  script/fit_simulated_ebpmf_wbg_init.R
    Untracked:  script/fit_simulated_ebpmf_wbg_init.sbatch
    Untracked:  script/fit_simulated_ebpmf_wbg_init.sh
    Untracked:  script/fit_simulated_fastTopics.R
    Untracked:  script/fit_simulated_fastTopics.sbatch
    Untracked:  script/fit_simulated_fastTopics.sh
    Untracked:  script/fit_simulated_fastTopics2.sh
    Untracked:  script/fit_simulated_fastTopics3.sh
    Untracked:  script/fit_simulated_fastTopics4.sh
    Untracked:  script/fit_simulated_fastTopics_init.R
    Untracked:  script/fit_simulated_fastTopics_init.sbatch
    Untracked:  script/fit_simulated_fastTopics_init.sh
    Untracked:  script/fit_sla_ebpmf_wbg2.sbatch
    Untracked:  script/fit_sla_ebpmf_wbg3.sbatch
    Untracked:  script/fit_sla_fastTopics.R
    Untracked:  script/fit_sla_fastTopics.sbatch
    Untracked:  script/fit_sla_pmf2.sbatch
    Untracked:  script/fit_sla_pmf3.sbatch
    Untracked:  script/fit_sla_sim_fastTopics.R
    Untracked:  script/fit_sla_sim_fastTopics.sbatch
    Untracked:  script/fit_sla_sim_fastTopics.sh
    Untracked:  script/pbmc_68k.RData
    Untracked:  script/prepare_68k_pbmc.R
    Untracked:  script/prepare_droplet.R
    Untracked:  script/simulate_droplet.R
    Untracked:  script/simulate_pbmc_68k.R
    Untracked:  script/slurm-11984946.out
    Untracked:  script/slurm-11984947.out
    Untracked:  script/slurm-11984948.out
    Untracked:  script/slurm-11984950.out
    Untracked:  script/slurm-11984951.out
    Untracked:  script/slurm-11984952.out
    Untracked:  script/slurm-11984953.out
    Untracked:  script/slurm-11984954.out
    Untracked:  script/slurm-11984955.out
    Untracked:  script/slurm-11984956.out
    Untracked:  script/slurm-11984957.out
    Untracked:  script/slurm-11984958.out
    Untracked:  script/slurm-11984959.out
    Untracked:  script/slurm-11984960.out
    Untracked:  script/slurm-11984961.out
    Untracked:  script/slurm-11984962.out
    Untracked:  script/slurm-11984963.out
    Untracked:  script/slurm-11984964.out
    Untracked:  script/slurm-11984965.out
    Untracked:  script/slurm-11984966.out
    Untracked:  script/slurm-11984967.out
    Untracked:  script/slurm-11984968.out
    Untracked:  script/slurm-11984969.out
    Untracked:  script/slurm-11984970.out
    Untracked:  script/slurm-11984971.out
    Untracked:  script/slurm-11984972.out
    Untracked:  script/slurm-11984973.out
    Untracked:  script/slurm-11984974.out
    Untracked:  script/slurm-11984975.out
    Untracked:  script/slurm-11984976.out
    Untracked:  script/slurm-11984977.out
    Untracked:  script/slurm-11984978.out
    Untracked:  script/slurm-11984979.out
    Untracked:  script/slurm-11984980.out
    Untracked:  script/slurm-11984981.out
    Untracked:  script/slurm-11984982.out
    Untracked:  script/slurm-11984983.out
    Untracked:  script/slurm-11984984.out
    Untracked:  script/slurm-11987069.out
    Untracked:  script/slurm-11987070.out
    Untracked:  script/slurm-11987071.out
    Untracked:  script/slurm-11987072.out
    Untracked:  script/slurm-11987073.out
    Untracked:  script/slurm-11987074.out
    Untracked:  script/slurm-11987075.out
    Untracked:  script/slurm-11987077.out
    Untracked:  script/slurm-11987079.out
    Untracked:  script/slurm-11987080.out
    Untracked:  script/slurm-11987081.out
    Untracked:  script/slurm-11987082.out
    Untracked:  script/slurm-11987083.out
    Untracked:  script/slurm-11987084.out
    Untracked:  script/slurm-11987085.out
    Untracked:  script/slurm-11987086.out
    Untracked:  script/slurm-11987087.out
    Untracked:  script/slurm-11987088.out
    Untracked:  script/slurm-11987089.out
    Untracked:  script/slurm-11987090.out
    Untracked:  script/slurm-11987091.out
    Untracked:  script/slurm-11987092.out
    Untracked:  script/slurm-11987093.out
    Untracked:  script/slurm-11987094.out
    Untracked:  script/slurm-11987095.out
    Untracked:  script/slurm-11987096.out
    Untracked:  script/slurm-11987097.out
    Untracked:  script/slurm-11987098.out
    Untracked:  script/slurm-11987099.out
    Untracked:  script/slurm-11987100.out
    Untracked:  script/slurm-11987101.out
    Untracked:  script/slurm-11987102.out
    Untracked:  script/slurm-11987103.out
    Untracked:  script/slurm-11987104.out
    Untracked:  script/slurm-11987105.out
    Untracked:  script/slurm-11987106.out
    Untracked:  script/slurm-11987107.out
    Untracked:  script/slurm-11987108.out
    Untracked:  script/slurm-11987401.out
    Untracked:  script/slurm-11987402.out
    Untracked:  script/slurm-11987403.out
    Untracked:  script/slurm-11987404.out
    Untracked:  script/slurm-11987405.out
    Untracked:  script/slurm-11987406.out
    Untracked:  script/slurm-11987407.out
    Untracked:  script/slurm-11987408.out
    Untracked:  script/slurm-11987409.out
    Untracked:  script/slurm-11987410.out
    Untracked:  script/slurm-11987411.out
    Untracked:  script/slurm-11987412.out
    Untracked:  script/slurm-11987413.out
    Untracked:  script/slurm-11987414.out
    Untracked:  script/slurm-11987415.out
    Untracked:  script/slurm-11987416.out
    Untracked:  script/slurm-11987417.out
    Untracked:  script/slurm-11987418.out
    Untracked:  script/slurm-11987419.out
    Untracked:  script/slurm-11987420.out
    Untracked:  script/slurm-11987421.out
    Untracked:  script/slurm-11987422.out
    Untracked:  script/slurm-11987423.out
    Untracked:  script/slurm-11987424.out
    Untracked:  script/slurm-11987425.out
    Untracked:  script/slurm-11987426.out
    Untracked:  script/slurm-11987427.out
    Untracked:  script/slurm-11987428.out
    Untracked:  script/slurm-11987429.out
    Untracked:  script/slurm-11987430.out
    Untracked:  script/slurm-11987431.out
    Untracked:  script/slurm-11987432.out
    Untracked:  script/slurm-11987433.out
    Untracked:  script/slurm-11987434.out
    Untracked:  script/slurm-11987435.out
    Untracked:  script/slurm-11987436.out
    Untracked:  script/slurm-11987437.out
    Untracked:  script/slurm-11987438.out
    Untracked:  script/smallsim2.R
    Untracked:  test.R

Unstaged changes:
    Modified:   code/misc.R
    Deleted:    data/sim/docword.sim_bg_block_n1100_p2100_K50.txt
    Deleted:    data/sim/init.sim_bg_block_n1100_p2100_K50.Rds
    Deleted:    data/sim/truth.sim_bg_block_n1100_p2100_K50.Rds
    Modified:   output/fastTopics_fit/fit_droplet_simulated_fastTopics_k19_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k10_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k11_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k12_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k2_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k3_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k4_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k5_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k6_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k7_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k8_comparison.pdf
    Modified:   output/fastTopics_fit/fit_sla_sim_fastTopics_k9_comparison.pdf
    Deleted:    output/sim/v0.4.5/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000.Rout
    Deleted:    output/sim/v0.4.5/fit_sim_bg_block_n1100_p2100_K50_ebpmf_wbg_maxiter_5000_from_truth.Rout
    Deleted:    output/sim/v0.4.5/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter3.Rds
    Deleted:    output/sim/v0.4.5/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter5000.Rds
    Deleted:    output/sim/v0.4.5/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter5000_from_truth.Rds
    Deleted:    output/sim/v0.4.5/sim_bg_block_n1100_p2100_K50_ebpmf_wbg_K50_maxiter50_from_truth2.Rds
    Modified:   script/fit_kos_NMF_F.R
    Modified:   script/fit_sim_ebpmf_wbg.sbatch
    Modified:   script/fit_sim_ebpmf_wbg_from_pmf_bg.sbatch
    Modified:   script/fit_sim_ebpmf_wbg_from_pmf_bg2.sbatch

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/sla_simulated_compare_mle_eb.Rmd) and HTML (docs/sla_simulated_compare_mle_eb.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
Rmd a87c087 Zihao 2021-07-01 update sla_simulated_compare_mle_eb.Rmd
html 87da1d8 Zihao 2021-07-01 Build site.
Rmd 872d6c2 Zihao 2021-07-01 update sla_simulated_compare_mle_eb.Rmd
html 0e449f1 Zihao 2021-07-01 Build site.
Rmd e0bcfa6 Zihao 2021-07-01 update sla_simulated_compare_mle_eb.Rmd
html 37a48ab Zihao 2021-06-30 Build site.
Rmd 938181c Zihao 2021-06-30 update sla_simulated_compare_mle_eb.Rmd
html a616a65 Zihao 2021-06-30 Build site.
Rmd 983a327 Zihao 2021-06-30 sla_simulated_compare_mle_eb.Rmd

Introduction

  • I thin sla dataset by half to get training and test dataset (see prepare_sla_sim). I first fit MLE (fastTopics) on original dataset, and call it truth . (Probably not a good name but it’s convenient…). Then I fit MLE and EB on thinned dataset. A good model should have good likelihood on test data, and should be close to truth.
  • I initialize MLE (fastTopics) and EB (ebpmf-wbg) methods using the truth.
  • EB approach is much better than MLE in the two cases shown below: \(K = 6\) & \(K = 20\).
  • Note that the mean matrix is slightly different bewteen the original data and train/test, as I remove some columns with 0 counts on training.
rm(list = ls())
source("code/misc.R")
library(fastTopics)
library(ebpmf.alpha)

datafile = "data/sim/sla_simulated.Rds"
data = readRDS(datafile)
data$Y_train = as(data$Y_train, "sparseMatrix")
data$Y_test = as(as.matrix(data$Y_test), "sparseMatrix") ## the way i generate simulated data is a bit stupid

dim(data$Y_train)
[1] 3207 8030
dim(data$Y_test)
[1] 3207 8030
sum(data$Y_train != 0)/(nrow(data$Y_train) * ncol(data$Y_train))
[1] 0.004549318
sum(data$Y_test != 0)/(nrow(data$Y_test) * ncol(data$Y_test))
[1] 0.004498682

\(K = 6\)

k = 6
truth.file = sprintf("output/fastTopics_fit/fit_sla_fastTopics_k%d.Rds", k)
fitted.mle.file = sprintf("output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k%d.Rds", k)
fitted.eb.file = sprintf("output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k%d.Rds", k)


truth = readRDS(truth.file)
fitted.mle = readRDS(fitted.mle.file)
fitted.eb = readRDS(fitted.eb.file)

## get "truth"
truth = fastTopics::multinom2poisson(truth)
truth = list(L = truth$L, F = truth$F[data$w_idx,])
class(truth) <- c("poisson_nmf_fit","list")
truth = fastTopics::poisson2multinom(truth) ## s is not right, but this does not affect computation of log-likelihood


fitted.mle.m = fastTopics::poisson2multinom(fitted.mle)
fitted.eb.p = ebpmf2fastTopics_pmf(fitted.eb)
fitted.eb.m = fastTopics::poisson2multinom(fitted.eb.p)

idx.mle = match_topics(truth$F, fitted.mle.m$F)
idx.eb = match_topics(truth$F, fitted.eb.m$F)

## see if there is 1-1 correspondence bewteen topics
idx.mle
[1] 1 2 3 4 5 6
idx.eb
[1] 1 2 3 4 5 6

compare log-lik on training/test

loglik.train.truth = fastTopics::loglik_multinom_topic_model(X = data$Y_train, truth)
loglik.train.mle = fastTopics::loglik_multinom_topic_model(X = data$Y_train, fitted.mle.m)
loglik.train.eb = fastTopics::loglik_multinom_topic_model(X = data$Y_train, fitted.eb.m)

loglik.test.truth = fastTopics::loglik_multinom_topic_model(X = data$Y_test, truth)
loglik.test.mle = fastTopics::loglik_multinom_topic_model(X = data$Y_test, fitted.mle.m)
loglik.test.eb = fastTopics::loglik_multinom_topic_model(X = data$Y_test, fitted.eb.m)


## training error
###  truth          mle             eb
c(mean(loglik.train.truth), mean(loglik.train.mle), mean(loglik.train.eb))
[1] -187.0736 -184.3535 -188.8222
## test error
###  truth          mle             eb
c(mean(loglik.test.truth), mean(loglik.test.mle), mean(loglik.test.eb))
[1] -183.5079 -204.8046 -191.4603

MLE vs truth

model = fitted.mle.m
idx = match_topics(truth$F, model$F)
l_mse <- c()
f_mse <- c()

par(mfrow = c(k, 2))
for(i in 1:k){
  j = idx[i]
  plot(truth$L[,i], model$L[,j], xlab = "loading (truth)", ylab = "loading (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  plot(truth$F[,i], model$F[,j], xlab = "factor (truth)", ylab = "factor (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  l_mse = c(l_mse, mean((truth$L[,i] - model$L[,j])^2))
  f_mse = c(f_mse, mean((truth$F[,i] - model$F[,j])^2))
}

Version Author Date
a616a65 Zihao 2021-06-30
l_mse.mle = l_mse
f_mse.mle = f_mse

EB vs truth

model = fitted.eb.m
idx = match_topics(truth$F, model$F)
l_mse <- c()
f_mse <- c()

par(mfrow = c(k, 2))
for(i in 1:k){
  j = idx[i]
  plot(truth$L[,i], model$L[,j], xlab = "loading (truth)", ylab = "loading (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  plot(truth$F[,i], model$F[,j], xlab = "factor (truth)", ylab = "factor (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  l_mse = c(l_mse, mean((truth$L[,i] - model$L[,j])^2))
  f_mse = c(f_mse, mean((truth$F[,i] - model$F[,j])^2))
}

Version Author Date
a616a65 Zihao 2021-06-30
l_mse.eb = l_mse
f_mse.eb = f_mse

compare MSE from truth

round(l_mse.eb/l_mse.mle, digits = 2)
[1] 0.55 0.64 0.74 0.56 0.40 0.35
round(f_mse.eb/f_mse.mle, digits = 2)
[1] 0.71 1.21 1.29 1.02 0.91 0.65

look at EB closely

\(g(L)\)

get_prior_summary(gs = fitted.eb$qg$gls, log10 = FALSE)

Version Author Date
0e449f1 Zihao 2021-07-01

\(g(F)\)

get_prior_summary(gs = fitted.eb$qg$gfs, log10 = FALSE)

Version Author Date
0e449f1 Zihao 2021-07-01

\(l0, f0, q(l), q(f)\)

par(mfrow = c(k, 2))
hist(log(fitted.eb$l0), breaks = 30, xlab = "log(l0)", main = "log(l0)")
hist(log(fitted.eb$f0), breaks = 30, xlab = "log(f0)", main = "log(f0)")
for(i in 1:k){
  hist(log(fitted.eb$qg$qls_mean[,i]), breaks = 30, xlab = "log(l)", main = sprintf("topic %d", i))
  hist(log(fitted.eb$qg$qfs_mean[,i]), breaks = 30, xlab = "log(f)", main = sprintf("topic %d", i))
}

Version Author Date
0e449f1 Zihao 2021-07-01
a616a65 Zihao 2021-06-30

Version Author Date
0e449f1 Zihao 2021-07-01

\(w\)

fitted.eb$w
[1] 3.864248 2.497604 2.850447 3.505004 5.351958 3.435125

\(K = 20\)

k = 20
truth.file = sprintf("output/fastTopics_fit/fit_sla_fastTopics_k%d.Rds", k)
fitted.mle.file = sprintf("output/fastTopics_fit/fit_sla_simualted_fastTopics_init_truth_k%d.Rds", k)
fitted.eb.file = sprintf("output/fastTopics_fit/fit_sla_simualted_ebpmf_wbg_init_truth_k%d.Rds", k)

truth = readRDS(truth.file)
fitted.mle = readRDS(fitted.mle.file)
fitted.eb = readRDS(fitted.eb.file)

## get "truth"
truth = fastTopics::multinom2poisson(truth)
truth = list(L = truth$L, F = truth$F[data$w_idx,])
class(truth) <- c("poisson_nmf_fit","list")
truth = fastTopics::poisson2multinom(truth) ## s is not right, but this does not affect computation of log-likelihood


fitted.mle.m = fastTopics::poisson2multinom(fitted.mle)
fitted.eb.p = ebpmf2fastTopics_pmf(fitted.eb)
fitted.eb.m = fastTopics::poisson2multinom(fitted.eb.p)

idx.mle = match_topics(truth$F, fitted.mle.m$F)
idx.eb = match_topics(truth$F, fitted.eb.m$F)

## see if there is 1-1 correspondence bewteen topics
idx.mle
 [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20
idx.eb
 [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20

compare log-lik on training/test

loglik.train.truth = fastTopics::loglik_multinom_topic_model(X = data$Y_train, truth)
loglik.train.mle = fastTopics::loglik_multinom_topic_model(X = data$Y_train, fitted.mle.m)
loglik.train.eb = fastTopics::loglik_multinom_topic_model(X = data$Y_train, fitted.eb.m)

loglik.test.truth = fastTopics::loglik_multinom_topic_model(X = data$Y_test, truth)
loglik.test.mle = fastTopics::loglik_multinom_topic_model(X = data$Y_test, fitted.mle.m)
loglik.test.eb = fastTopics::loglik_multinom_topic_model(X = data$Y_test, fitted.eb.m)


## training error
###  truth          mle             eb
c(mean(loglik.train.truth), mean(loglik.train.mle), mean(loglik.train.eb))
[1] -172.6338 -167.1411 -174.7843
## test error
###  truth          mle             eb
c(mean(loglik.test.truth), mean(loglik.test.mle), mean(loglik.test.eb))
[1] -169.4004 -218.8784 -188.2133

MLE vs truth

model = fitted.mle.m
idx = match_topics(truth$F, model$F)
l_mse <- c()
f_mse <- c()

par(mfrow = c(k, 2))
for(i in 1:k){
  j = idx[i]
  plot(truth$L[,i], model$L[,j], xlab = "loading (truth)", ylab = "loading (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  plot(truth$F[,i], model$F[,j], xlab = "factor (truth)", ylab = "factor (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  l_mse = c(l_mse, mean((truth$L[,i] - model$L[,j])^2))
  f_mse = c(f_mse, mean((truth$F[,i] - model$F[,j])^2))
}

Version Author Date
0e449f1 Zihao 2021-07-01
l_mse.mle = l_mse
f_mse.mle = f_mse

EB vs truth

model = fitted.eb.m
idx = match_topics(truth$F, model$F)
l_mse <- c()
f_mse <- c()

par(mfrow = c(k, 2))
for(i in 1:k){
  j = idx[i]
  plot(truth$L[,i], model$L[,j], xlab = "loading (truth)", ylab = "loading (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  plot(truth$F[,i], model$F[,j], xlab = "factor (truth)", ylab = "factor (fitted)", main = sprintf("topic %d", i))
  abline(a = 0, b = 1, col = "blue")
  l_mse = c(l_mse, mean((truth$L[,i] - model$L[,j])^2))
  f_mse = c(f_mse, mean((truth$F[,i] - model$F[,j])^2))
}

Version Author Date
0e449f1 Zihao 2021-07-01
l_mse.eb = l_mse
f_mse.eb = f_mse

compare MSE from truth

round(l_mse.eb/l_mse.mle, digits = 2)
 [1] 0.37 0.26 0.37 0.41 0.39 0.28 0.63 0.28 0.32 0.46 0.58 0.40 0.41 0.44
[15] 0.40 0.91 0.45 0.56 0.50 0.51
round(f_mse.eb/f_mse.mle, digits = 2)
 [1] 0.67 0.54 0.61 0.54 0.47 0.72 1.09 0.69 0.64 0.39 0.62 1.02 0.56 0.69
[15] 0.33 0.35 0.91 0.75 0.49 0.82

look at EB closely

\(g(L)\)

get_prior_summary(gs = fitted.eb$qg$gls, log10 = FALSE)

Version Author Date
0e449f1 Zihao 2021-07-01

\(g(F)\)

get_prior_summary(gs = fitted.eb$qg$gfs, log10 = FALSE)

Version Author Date
0e449f1 Zihao 2021-07-01

\(l0, f0, q(l), q(f)\)

par(mfrow = c(k, 2))
hist(log(fitted.eb$l0), breaks = 30, xlab = "log(l0)", main = "log(l0)")
hist(log(fitted.eb$f0), breaks = 30, xlab = "log(f0)", main = "log(f0)")
for(i in 1:k){
  hist(log(fitted.eb$qg$qls_mean[,i]), breaks = 30, xlab = "log(l)", main = sprintf("topic %d", i))
  hist(log(fitted.eb$qg$qfs_mean[,i]), breaks = 30, xlab = "log(f)", main = sprintf("topic %d", i))
}

Version Author Date
0e449f1 Zihao 2021-07-01

Version Author Date
0e449f1 Zihao 2021-07-01

\(w\)

fitted.eb$w
 [1]  402.3158  485.2620  480.9474  420.1032  362.1982 1342.4510  343.3042
 [8]  394.2655  449.9989  429.4477  521.7581  291.1336  451.3496  383.4733
[15]  731.5467  735.3171  475.3372  258.9959  316.9660  321.8124

sessionInfo()
R version 3.6.1 (2019-07-05)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Scientific Linux 7.4 (Nitrogen)

Matrix products: default
BLAS/LAPACK: /software/openblas-0.2.19-el7-x86_64/lib/libopenblas_haswellp-r0.2.19.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

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

other attached packages:
[1] ebpmf.alpha_0.4.5 fastTopics_0.5-53 pheatmap_1.0.12   workflowr_1.6.2  

loaded via a namespace (and not attached):
 [1] httr_1.4.1         tidyr_1.1.0        jsonlite_1.6      
 [4] viridisLite_0.3.0  gtools_3.8.1       RcppParallel_4.4.3
 [7] assertthat_0.2.1   mixsqp_0.3-43      yaml_2.2.0        
[10] progress_1.2.2     ggrepel_0.9.1      pillar_1.4.2      
[13] backports_1.1.4    lattice_0.20-38    quantreg_5.41     
[16] glue_1.3.1         quadprog_1.5-7     digest_0.6.20     
[19] RColorBrewer_1.1-2 promises_1.0.1     colorspace_1.4-1  
[22] cowplot_1.0.0      htmltools_0.3.6    httpuv_1.5.1      
[25] Matrix_1.2-18      pkgconfig_2.0.2    invgamma_1.1      
[28] SparseM_1.77       purrr_0.3.4        scales_1.1.0      
[31] whisker_0.3-2      later_0.8.0        Rtsne_0.15        
[34] MatrixModels_0.4-1 git2r_0.26.1       tibble_2.1.3      
[37] ggplot2_3.3.3      ashr_2.2-47        lazyeval_0.2.2    
[40] magrittr_1.5       crayon_1.3.4       mcmc_0.9-6        
[43] evaluate_0.14      ebpm_0.0.1.1       fs_1.3.1          
[46] MASS_7.3-51.4      truncnorm_1.0-8    tools_3.6.1       
[49] data.table_1.13.2  prettyunits_1.0.2  hms_0.5.3         
[52] lifecycle_0.1.0    stringr_1.4.0      MCMCpack_1.4-5    
[55] plotly_4.9.0       munsell_0.5.0      irlba_2.3.3       
[58] compiler_3.6.1     rlang_0.4.10       grid_3.6.1        
[61] htmlwidgets_1.3    rmarkdown_1.13     gtable_0.3.0      
[64] R6_2.4.0           knitr_1.23         dplyr_0.8.3       
[67] uwot_0.1.5         rprojroot_1.3-2    stringi_1.4.3     
[70] SQUAREM_2017.10-1  Rcpp_1.0.5         vctrs_0.3.1       
[73] tidyselect_1.1.0   xfun_0.8           coda_0.19-3