<html><body><div style="font-family: trebuchet ms,sans-serif; font-size: 12pt; color: #000000"><div>Perhaps something like this?</div><div style="padding-left: 30px;"><br data-mce-bogus="1"></div><div style="padding-left: 30px;"><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">library(adegenet)</span><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">library(pegas)</span><br><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">data(nancycats)</span><br><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">lvls <- levels(pop(nancycats))</span><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">out <- vector("list", length(lvls)) # prepare list for output</span><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">names(out) <- lvls</span><br><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">for (i in lvls) {</span><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;"> out[[i]] <- hw.test(nancycats[pop(nancycats) == i])</span><br><span style="font-family: "courier new", courier, monaco, monospace, sans-serif;" data-mce-style="font-family: 'courier new', courier, monaco, monospace, sans-serif;">}</span><br></div><div><br data-mce-bogus="1"></div><div>Cheers,</div><div>Roman</div><div><br></div><div data-marker="__SIG_PRE__">----<br>In god we trust, all others bring data.</div><br><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>From: </b>"Felipe Hernández" <fhernandeu@uc.cl><br><b>To: </b>adegenet-forum@lists.r-forge.r-project.org<br><b>Sent: </b>Thursday, March 2, 2017 6:04:42 PM<br><b>Subject: </b>[adegenet-forum] HWE calculation<br></div><br><div data-marker="__QUOTED_TEXT__"><div dir="ltr"><div><span style="font-size:13px">Hi everyone, </span><div style="font-size:13px"><br></div><div style="font-size:13px">I have used the adegenet/pegas packages to estimate HWE using a genind object. I got the results from chi-squared and exact test based on Monte Carlo for my whole set of loci (52 loci) across all my 29 populations. I wonder if there is any option to calculate HWE and get exact test Monte Carlo results, but for each population separately. Sorry if the question is so basic, but I would appreciate any helpful advice!</div></div><div style="font-size:13px"><br></div><div style="font-size:13px">Best,</div><div style="font-size:13px">Felipe</div><div style="font-size:13px"><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Felipe Hernández<div>Médico Veterinario (DVM), MSc.<div>PhD. Candidate<div>Interdisciplinary Ecology Program</div><div>School of Natural Resources and Environment</div><div>Wildlife Ecology and Conservation Department</div><div>University of Florida</div></div></div></div></div></div></div></div></div>
</div>
<br>_______________________________________________<br>adegenet-forum mailing list<br>adegenet-forum@lists.r-forge.r-project.org<br>https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum<br></div></div></body></html>