SAUCIE icon indicating copy to clipboard operation
SAUCIE copied to clipboard

About 'sinh'

Open yuhuang-cst opened this issue 5 years ago • 1 comments

  1. The cell-gene matrix from the same batch should be stored in the same '.csv' located in 'input_dir'. Right?
  2. The values of cell-gene matrix stored in '.csv' represent the raw counts, because the 'get_data' function will transform them with 'math.asinh'. Right?
  3. Why 'output_batch_correction' function doesn't take 'sinh' to get imputated raw count matrix? https://github.com/KrishnaswamyLab/SAUCIE/blob/c2e59683ddf401f07d4c226a420b367181934715/SAUCIE.py#L159
  4. If '3.' holds, the data used for clustering will be transformed by 'math.asinh' twice: one is in the batch correction stage, the other is in cluster stage.

yuhuang-cst avatar Nov 07 '19 03:11 yuhuang-cst

It seems that for clustering, SAUCIE uses batch-corrected data (asinh values) if SAUCIE was also directed to do batch-correction. Raw data will be used if batch-correction is not specified. Thus, asinh transformation will occur only once, unless raw data are asinh values, which it seems will be so if SAUCIE is used separately to first batch-correct and then cluster (see issue #22). So having an option to output non-asinh data will be good.

(SAUCIE.py)
# CLUSTERING
if args.cluster:
    if args.batch_correct:
        input_files = sorted(glob.glob(os.path.join(args.output_dir, 'batch_corrected', '*.{}'.format(args.format))))
    else:
        input_files = rawfiles

ashipde avatar Apr 23 '20 01:04 ashipde