Open YichaoOU opened 5 years ago
Currently, I'm using transcript_tpm = sleuth_to_matrix(so,"obs_norm","tpm") to get a normalized gene expression matrix and then use log2((transcript_tpm[,"treatment_mean"]+1e-100)/(transcript_tpm[,"control_mean"]+1e-100)) as the logFC.
transcript_tpm = sleuth_to_matrix(so,"obs_norm","tpm")
log2((transcript_tpm[,"treatment_mean"]+1e-100)/(transcript_tpm[,"control_mean"]+1e-100))
Is it correct?
Thanks, Yichao
Currently, I'm using
transcript_tpm = sleuth_to_matrix(so,"obs_norm","tpm")
to get a normalized gene expression matrix and then uselog2((transcript_tpm[,"treatment_mean"]+1e-100)/(transcript_tpm[,"control_mean"]+1e-100))
as the logFC.Is it correct?
Thanks, Yichao