<div dir="ltr">Hi Thibault,<div><br></div><div>     Thank you very much for your help. I just downloaded the dev version but I get the same results using summary(x) than with the previous version. summary gives me expected and observed heterozygosity but per loci:</div><div><br></div><div>summary(x)$Hobs</div><div><br></div><div>     What I would like to have is Hobs <u>per population</u> instead. Something similar to the function Hs but for Hobs. Is there anything like that in adegenet? I tried basic.stats from Hierfstat as Maria suggested, but I got an error message, probably due to something wrong in the way I formatted the data.</div><div><br></div><div>Simon</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 23, 2015 at 6:58 AM, Maria Guerrina <span dir="ltr"><<a href="mailto:maria.guerrina@edu.unige.it" target="_blank">maria.guerrina@edu.unige.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi Thibaut!<div>I have a problem similar to Simon.</div><div>I have tried to use the function 'summary', but I get 0 both for the observed heterozigosity both for the expected one...</div><div>I have tried to use the function 'basic.stats' in hierfstat, but it doesn't give a result for Ho, in my case with my data.</div><div><br></div><div>any idea?</div><div>Thanks!</div><div>Maria<br><div>
<span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-size:medium;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><div style="word-wrap:break-word"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><div style="word-wrap:break-word"><div>--</div><div>Maria Guerrina PhD</div><div><div>Università di Genova</div><div>DISTAV</div><div>Corso Dogali 1M</div><div>I - 16136 GENOVA (Italy)</div></div><div><a href="mailto:maria.guerrina@edu.unige.it" target="_blank">maria.guerrina@edu.unige.it</a></div><div><br></div></div></span><br></div></span><br></span><br>
</div>
<br><div><div><div class="h5"><div>On 23/nov/2015, at 13.11, Jombart, Thibaut wrote:</div><br></div></div><blockquote type="cite"><span style="border-collapse:separate;font-family:Helvetica;font-size:medium;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div><div><div class="h5"><div style="direction:ltr;font-family:Tahoma;color:rgb(0,0,0);font-size:10pt">Hello,<span> </span><br><br>have you tried 'summary'? Note, the current devel version is better (previous summary was returning invisibly a list).<span> </span><br><br>Cheers<br>Thibaut<br><div><br><div style="font-family:Tahoma;font-size:13px"><div><font size="2"><span style="font-size:10pt"><div><br></div></span></font></div></div></div><div style="font-family:'Times New Roman';color:rgb(0,0,0);font-size:16px"><hr><div style="direction:ltr"><font color="#000000" face="Tahoma" size="2"><b>From:</b><span> </span><a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org" target="_blank">adegenet-forum-bounces@lists.r-forge.r-project.org</a><span> </span>[<a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org" target="_blank">adegenet-forum-bounces@lists.r-forge.r-project.org</a>] on behalf of Simon Nadeau [<a href="mailto:simon.nadeau.ubc@gmail.com" target="_blank">simon.nadeau.ubc@gmail.com</a>]<br><b>Sent:</b><span> </span>19 November 2015 23:00<br><b>To:</b><span> </span><a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.r-project.org</a><br><b>Subject:</b><span> </span>[adegenet-forum] How to get observed heterozygosity per population<br></font><br></div><div></div><div><div dir="ltr">Hi,<div><br></div><div>      I have recently been using the various features of adegenet unsing genind objects and I works great so far. However, I can't seem to be able to find easy summary stats:</div><div><br></div><div>Observed heterozygosity per population: this would be very handy but I only found the function Hs (expected heterozygosity).</div><div><br></div><div>Global allele frequencies per loci: there is a nice histogram of minor allele frequencies on page 30 of the adegenet genomic tutorial, but I can't find how to generate the same plot using genind objects.</div><div><br></div><div>Thank you very much for your help,</div><div><br></div><div>Simon<br clear="all"><div><br></div>--<span> </span><br><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><font color="#3d85c6" size="4">Simon Nadeau, M. Sc.</font></div><div><font color="#3d85c6"><font size="3">Biologiste / Biologist</font></font></div><div><strong><font color="#007f00" size="3"></font></strong><br></div><div><div>Ressources naturelles Canada, Service canadien des forêts, Centre de foresterie des Laurentides / Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre</div><div>1055 rue du P.E.P.S., C.P. 10380 succ. Sainte-Foy, Québec (Québec), Canada G1V 4C7</div></div><div><a href="mailto:Simon.Nadeau.UBC@gmail.com" style="font-size:12.8px" target="_blank">Simon.Nadeau@a</a><span style="font-size:12.8px;color:rgb(0,0,255)"><a href="http://lumni.ubc.ca" target="_blank">lumni.UBC.ca</a>, </span><u style="color:rgb(0,0,255);font-size:12.8px"><span> </span><div style="display:inline!important"><a href="mailto:Simon.Nadeau@RNCan-NRCan.gc.ca" target="_blank">Simon.Nadeau@canada.ca</a></div></u></div><div><span style="font-size:12.8px">Tel: <a href="tel:%28604%29%20349-5196" value="+16043495196" target="_blank">(604) 349-5196</a></span><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>_______________________________________________<br>adegenet-forum mailing list<br><a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.r-project.org</a><br><a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum</a></div></span></blockquote></div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>
<font size="4" color="#3d85c6">Simon Nadeau, M. Sc.</font></div><div><font color="#3d85c6"><font size="3">Biologiste / Biologist</font></font></div><div><strong><font color="#007f00" size="3"></font></strong><br></div><div><div>Ressources naturelles Canada, Service canadien des forêts, Centre de foresterie des Laurentides / Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre</div><div>1055 rue du P.E.P.S., C.P. 10380 succ. Sainte-Foy, Québec (Québec), Canada G1V 4C7</div></div><div><a href="mailto:Simon.Nadeau.UBC@gmail.com" style="font-size:12.8000001907349px" target="_blank">Simon.Nadeau@a</a><span style="font-size:12.8000001907349px;color:rgb(0,0,255)"><a href="http://lumni.ubc.ca" target="_blank">lumni.UBC.ca</a>, </span><u style="color:rgb(0,0,255);font-size:12.8000001907349px"><div style="display:inline!important"><a href="mailto:Simon.Nadeau@RNCan-NRCan.gc.ca" target="_blank">Simon.Nadeau@canada.ca</a></div></u></div><div><span style="font-size:12.8000001907349px">Tel: (604) 349-5196</span><br></div></div></div></div></div></div></div></div></div></div></div></div>
</div>