I am trying to run SCDC with only one reference single cell rna seq data set, i.e. only one subject, e.g. sample H2 in Segerstolpe example data set.
SCDC_prop_ONE() runs until it encounters an error due to incompatible data in function SCDC_peval() line 258-260 on github, namely:
if (!is.data.frame(ptrue)){
ptrue <- as.data.frame.matrix(ptrue)
}
The following error message is displayed:
Error in seq_len(ncols) :
argument must be coercible to non-negative integer
In addition: Warning message:
In seq_len(ncols) : first element used of 'length.out' argument
It would be great if you could adapt the SCDC_peval() function. It ran smoothly when I used the full Segerstolpe example data set with function SCDC_prop().
Dear Meichen,
I am trying to run SCDC with only one reference single cell rna seq data set, i.e. only one subject, e.g. sample H2 in Segerstolpe example data set. SCDC_prop_ONE() runs until it encounters an error due to incompatible data in function SCDC_peval() line 258-260 on github, namely:
The following error message is displayed:
It would be great if you could adapt the SCDC_peval() function. It ran smoothly when I used the full Segerstolpe example data set with function SCDC_prop().
Best Beate
Here is my example workflow: