Abnormal results of interaction model

[yrjia1015]

Abnormal results of interaction model

[yrjia1015]

Thank you for your gracious and patient response.

Best regards,
Yuran

[francois-a]

Yes that's expected — to avoid artifacts due to outlier effects, a MAF filter that depends on the interaction term is applied (i.e., variants with MAF ≥ 0.05 in the samples belonging to each of the top and bottom halves of the interaction term). This can be adjusted with --maf_threshold_interaction. Please see section "Identification of cell type interaction eQTLs and sQTLs" in the supplement of https://www.science.org/doi/10.1126/science.aaz8528 for additional details.

[yrjia1015]

Thank you for your reply. Under the same phenotype and genotype background, changing only the cell score covariate of the interaction term can also lead to differences in the overall SNP gene pair of the test.

Can I understand that when testing each gene, only variations belonging to the upper and lower parts of the cell score distribution were considered, and the frequency of secondary alleles for these variations was not less than 5%?

This phenomenon is not caused by operational errors.

[francois-a]

That's almost right. It's the samples that were split into those belonging to the upper and lower parts of the cell score distribution, and MAF was computed separately for each of those two groups. Only variants that passed the MAF filter in both groups were retained. If you want to avoid this, you can set --maf_threshold_interaction 0,