Entering edit mode
11.2 years ago
SheelS
▴
40
Hi all,
I am doing the "GWAS" analysis and the method I used is "mmscore". I have parents and its offspring's trait data and GWA data.
Below is the model I used.
After creating a gwaa-class.....
>Model <- polygenic(weight ~ length, kin=gwasMydata, Mydata, trait.type = "gaussian")
>Modelan0 <- mmscore( Model,data=Mydata)
###Then I check my lambda inflation factor
>Modelan0
###Lambda
>37.94......(too big)
As I know the lambda should be around 1 (or a value very close to 1.00), I believe 39 is too big that should be something wrong,
###Does anyone know how to correct the inflation factor, or which means there are must something wrong in my input dataset?
By the way, the marker number in my dataset is around 10K and individual is around 650 including parents and offspring.
Really appreciate someone who can teach me to sort it out, thanks,