sunduanchen / Scissor

Scissor package
GNU General Public License v3.0
179 stars 31 forks source link

cutoff parameter does not limit number of cells chosen #40

Open kpruss opened 2 years ago

kpruss commented 2 years ago

I have not been able to successfully implement a cut-off threshold while calling Scissor.

For example, if alpha = 0.05 leads to 40% of cells being selected by Scissor and I want to decrease the proportion of Scissor-selected cells

cutoff_threshold = 0.2 alpha_threshold = 0.05

info <- Scissor(bulk_dataset = duod_counts, sc_dataset = seur_duod, phenotype = phenotype, alpha = alpha_threshold, cutoff = cutoff_threshold, tag = tag, family = "gaussian", Save_file = paste0('sbatch_outputs/grid_search/', date, alphathreshold, "", cutoff_threshold, protein, '.RData'))

The number and identity of cells selected is exactly the same as that chosen without the inclusion of the cutoff parameter.