GWAS of ordinal phenotypes with MCMCglmm model

Question

[rthapa]

Hi,

I want to do genome-wide association studies of ordinal phenotypes (1,2,3,4,5). Since my phenotypic data is not continuous, I cannot use a generalized linear model. The few R packages I found are MASS:polr, ordinal package with clm function and MCMCglmm. It looks like polr and ordinal package cannot take random effects. Does anyone have experience in using any ordinal model for GWAS of ordinal traits? Did anyone use MCMCglmm model for association studies? Thanks