Hello everyone!
I got 90 individuals genotyped in 3 replicates (genotyping by sequencing), so as a result, my vcf file contains information on 270 samples. The phenotype data were collected for 90 individuals. Then I try to run the mixed linear model in TASSEL and each of 3 genotypes replicates got the same phenotype data. As I understand this is not ok since the phenotype variance getting smaller artificially.
And my problem solution now looks as follows:
- I make the permutations, changing the phenotype data for each of 3 genotype replicates randomly for 500 times and running the MLM to find significant associations on random combinations.
- Then I calculate p-value as follows: I divide the number of iterations with the number of significant SNPs equal to or greater than the number of significant SNPs obtained with using original phenotype-genotype combination by the number of iterations.
I obtain a pretty small p-value, around 0.02 and then also make a binomial test in order to ensure that it is unexpectable to obtain that number of significant SNPs by random chance.
If you ever had a similar situation please share some advice. I also would be grateful if you share your opinion.
Thank you.
0 answers
No answers yet.
Log in to answer this question.
-
Significance testing of top vs. random SNPs
written by am29 •
-
random forests regression permuation importance
written by bsp017 •
-
Mixed Linear Model and GWAS study
written by Kumar
-
How to formally test the replicate consistency having genotype data?
written by rimgubaev
-
How to interpret Centred IBS kinship matrix produced by TASSEL?
written by rimgubaev
-
How to make association studies using haplotype (LD) blocks instead of SNPs?
written by rimgubaev
-
Has anyone dealt with "Specified element is out of bounds" error produced by TASSEL?
written by rimgubaev
-
Which test is applied in order to generate p-values in TASSEL when using Mixed Linear Model?
written by rimgubaev
-
How to build a phylogenetic tree using fasta obtained from VCF file?
written by rimgubaev
-
SNP genotyping and analysis of association with phenotype
written by orzech_mag
Please format your post better so it is easier to read - split into paragraphs with breaks where there is a slight change in context. When the content is a "wall of text", it discourages people from reading and by extension, getting to the core of your problem.
I've tried to format your post as best as I can, but please feel free to edit it and make any changes you think would be appopriate.
That sounds not appropriate, and I think you should merge those replicates to keep one record per individual rather than 3. Are the genotypes between replicates mostly concordant?
Yes, most of them (>75%) cluster together but others do not (I find out it by building an NJ tree). It s a plant sample so each replicate corresponds to a plant from a certain line, so it is biological replicate not technical. I tried to merge genotypes but in that case, I lose a huge amount of meaningful SNPs.