copy number variation as a covariate in MatrixEQTL
1
0
Entering edit mode
6.9 years ago
mms140130 ▴ 60

I want to add the copy number variation as a covariate in the model by MatrixEQTL, but I have 20,000 rows and 1069 columns, this won't work as the number rows is greater than the number of samples

How can I fix that?

R gene snp • 1.6k views
ADD COMMENT
0
Entering edit mode
6.9 years ago
Fabio Marroni ★ 3.0k

I don't think that adding 20,000 covariates is a good idea. You should do a PCA analysis on the CNV data and then use the first X components (where X is a number between 2 and 10, let's say) as covariates. PCA can used to correct for population structure (doing PCA on the genotype data) or to incorporate unaccounted bias in genotypes (e.g. I also suggest that you use as a covariate the PCA of your expression data). There is lot of discussion on the use of PCA in GWAS. You will find useful references here.

ADD COMMENT

Login before adding your answer.

Traffic: 2048 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6