GrahamDB / hzar

Hybrid Zone Analysis in R
1 stars 1 forks source link

Does model selection need to be conducted separately for each locus? #3

Open NACairns opened 5 years ago

NACairns commented 5 years ago

Dear Dr. Derryberry I have cline estimates for SNP hybrid index (overall and "diagnostic"), Cytb, morphology and for one transect call variables. I would like to look at a large number of my loci to get a feel for the range of estimates across loci as per Baldassarre et al. (2014)'s Figure 3 (https://onlinelibrary.wiley.com/doi/pdf/10.1111/evo.12457) and S3 (https://onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1111%2Fevo.12457&file=evo12457-sup-0003-TableS3.pdf). I have been trying to figure out a way to cycle through a list of input loci and run the whole process as a single function for each. Is there a built in way to do this in HZAR or does model selection for each locus need to be conducted separately. I apologies if I missed missed something obvious but I really appreciate any help you can provide me. Thank you, Best, Nick

NACairns commented 5 years ago

From Dr. Derryberry: ...I do have a series of helper scripts for repeated fitting of clines for multiple loci and character traits available at https://github.com/GrahamDB/MultiLocusHZAR. Hope that helps....