Question: QTCAT Quantitative Trait Cluster Association Test problem with reading SNP data
0
gravatar for tauseef.kiani
18 months ago by
tauseef.kiani0 wrote:

Hello everyone, I'm running a GWAS analysis. First I used GAPIT and wanted to use PCs as population control covariates. To determine the appropriate number of PCs to include, I used BIC values which recommended 0 PCs (highest BIC value) which was consistent with PC plots as they also show no clear population stratification.

As my population seems unstructured, i wanted to used another association tool QTCAT:Quantitative Trait Cluster Association Test which is for unstructured populations and removes the need for adjustment for population structure in association analysis. It is also an R package just like GAPIT.

the genotype file requirement for this package is diploid Hapmap. now i have the diploid hapmap from TASSEL, but the problem is the file contains so many NN for missing values, which is not accepted my QTCAT, hence, all SNPs with NN are removed/not considered which removes 2/3rd of genotype file and considers only 1/3rd, for which it is highly unlikely to get any association.

Question is, How can i use whole Hapmap file without removing NN? is there any way to Impute it first? another requirement is the snp data should be of snp.matrix class which only has homozygous and heterozygous form of allels. NN tells it there's another form, hence not considered and removed.

qtcat gapit gwas hapmap • 481 views
ADD COMMENTlink modified 10 months ago by BioRyder180 • written 18 months ago by tauseef.kiani0
0
gravatar for BioRyder
10 months ago by
BioRyder180
India
BioRyder180 wrote:

You should mention the missing genotype format ('NN' f) while reading the genotype file in QTCAT. Below are sample code for considering NN while reading file.

gfile <- system.file("genotype_file.csv", package = "qtcat")

snp <- read.snpData(gfile, sep = ",",na.string = "NN")

ADD COMMENTlink modified 10 months ago • written 10 months ago by BioRyder180
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2266 users visited in the last hour