First of all I'd like to thank you for the tool, I find it really helpful as it addresses key challenge while using bulk data. I'm trying to reproduce your work of deconvolution analysis for RA patient data. I have downloaded all the required packages and have followed all the steps exactly in the same way as mentioned in the paper (all input file downloaded from the given location in the paper) but facing an error while performing concordance analysis at the FACS reference.
An Error is thrown after running " bMIND_B <- evaluate_concordance(genes,B,label_sc,deconv2$A[,1,],label) "
Error: Error in pre_exp[i, pre_label == "OA"] : incorrect number of dimensions
Calls: evaluate_concordance -> wilcox.test
Execution halted
System Configuration: CentOS 7
In addition to this, the step:
res_alg_all_ra_l2_norm_trace <- cell_deconvolve_trace(O=log2(exprs(bulk_ra_vs_oa)[rownames(tmp_r
a$type),]+1),
theta=Fra_ra$theta,
R=log2(tmp_ra$type+1),
alpha=0.5,beta=1,gamma=1,
max.iter=1000,solver = "admm",verbose=FALSE,No
rmalize=FALSE)
in the paper is converging in 21 steps but for me it is converging in 22 steps.
Hi,
First of all I'd like to thank you for the tool, I find it really helpful as it addresses key challenge while using bulk data. I'm trying to reproduce your work of deconvolution analysis for RA patient data. I have downloaded all the required packages and have followed all the steps exactly in the same way as mentioned in the paper (all input file downloaded from the given location in the paper) but facing an error while performing concordance analysis at the FACS reference. An Error is thrown after running " bMIND_B <- evaluate_concordance(genes,B,label_sc,deconv2$A[,1,],label) " Error: Error in pre_exp[i, pre_label == "OA"] : incorrect number of dimensions Calls: evaluate_concordance -> wilcox.test Execution halted
System Configuration: CentOS 7
In addition to this, the step: res_alg_all_ra_l2_norm_trace <- cell_deconvolve_trace(O=log2(exprs(bulk_ra_vs_oa)[rownames(tmp_r a$type),]+1), theta=Fra_ra$theta, R=log2(tmp_ra$type+1), alpha=0.5,beta=1,gamma=1, max.iter=1000,solver = "admm",verbose=FALSE,No rmalize=FALSE)
in the paper is converging in 21 steps but for me it is converging in 22 steps.
Could you please help me resolve this error?