cozygene/TCA

cell-type-specific methylation for associated sites

antgomo opened this issue · 1 comments

Hi,

i am using tcasub and tensor to get the cell-type-specific methylation for a set of associated sites for plot purposes

I am following this

##subset the estimates of the model parameters
mdl.tca.sub <- tcasub(tca.mdl, features = rownames(betas.test), log_file = paste0(imid,".tcasub.log"))


# estimate cell-type-specific methylation for the associated site

Z_hat <- tensor(X = betas.test[my.sel.sites,], mdl.tca.sub, log_file=paste0(imid,".tensor.log"),parallel=T,num_cores=5)

My question is, if I ran previously tca for tca.mdlwith sex and age as a covariates, the cell-type-specific-methylation for this sites will be adjusted for sex and age, or do i have to adjust later? i.e, the signal contains still sex and age or not?

Thanks in advance

E-R commented

The 'tensor' function uses the covariates in your tca model (tca.mdl) when estimating the cell-type-specific levels, so the values are not adjusted for these covariates. It means that if you want to use Z_hat in a down-stream analysis then you should adjust for these covariates if appropriate.