diff --git a/R/music2.R b/R/music2.R index c393006..cf75a2a 100644 --- a/R/music2.R +++ b/R/music2.R @@ -337,7 +337,7 @@ music2_prop_t_statistics = function(bulk.control.mtx, bulk.case.mtx, sc.sce, clu MOD0=MOD0+mod0 id_d = sample(colnames(bulk.case),round(ncol(bulk.case)*sample_prop)) - case_s = exprs(bulk.case[exp_genel, colnames(bulk.case) %in% id_d]) + case_s = bulk.case[exp_genel, colnames(bulk.case) %in% id_d] prop_d = prop_CASE[colnames(case_s),] mod1 = apply(case_s, 1, function(x){ @@ -627,4 +627,4 @@ music2_prop_toast = function(bulk.control.mtx, bulk.case.mtx, sc.sce, clusters, return(list('Est.prop' = prop_all,'convergence'=TRUE,'n.iter'=iter,'DE.genes'=lr))} else{ return(list('Est.prop' = prop_all,'convergence'=FALSE,'id.not.converge'=rownames(prop_case_ini)))} -} \ No newline at end of file +}