I have #genomic variants called using #FreeBayes. I want to identify sex-linked regions of the genome. As a complement to visualization I’m considering doing some type of #statistical analysis, in part driven by wanting to learn more #stats. I’m less interested in allele differences (GWAS) and more in #heterozygosity, coverage depth, etc. Are there standard models to identify region differences? Could a #MixedModel be used for this?

Hey #mixedmodel #rstats peeps - anyone know of a package or function that does for #glmmTMB what merTools::predictInterval() does for #lme4? #wantingToMoveEverythingOver