<div dir="ltr">Hi there, <div><br></div><div>it is a bit hard to figure things out as we don't see the outputs of several commands and there is no reproducible example, but your import of the data looks good to me. It could be too many missing data in one of your groups, or a fixed locus? </div><div><br></div><div>Maybe worth checking allele counts per populations: </div><div>tab(genind2genpop(<span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">pinus.genind</span>))</div><div><br></div><div>Best</div><div>Thibaut</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><br>--<br>Dr Thibaut Jombart<br>Lecturer, Department of Infectious Disease Epidemiology, Imperial College London<br>Head of RECON: <a href="http://repidemicsconsortium.org" target="_blank">repidemicsconsortium.org</a><br>WHO Consultant - outbreak analysis</div><div><a href="https://thibautjombart.netlify.com" target="_blank">https://thibautjombart.netlify.com</a><br>Twitter: @TeebzR<br>+44(0)20 7594 3658</div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On 24 June 2018 at 02:31, Marcelo Laia <span dir="ltr"><<a href="mailto:marcelolaia@gmail.com" target="_blank">marcelolaia@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear all,<br>
<br>
I try to use adegenet for the first time.<br>
<br>
library("ape")<br>
library("pegas")<br>
library("seqinr")<br>
library("ggplot2")<br>
library("adegenet")<br>
<br>
pinus.data <- read.table("Analises_pinus_<wbr>individual.csv", header=TRUE,<br>
                         sep="\t", row.names = 1)<br>
head(pinus.data)<br>
> head(pinus.data)<br>
   PTTX3081 RIPT0031 PTTX3011 PTTX2037 PTTX311 Populacao<br>
1    186243   235284   131155   190190  191191    POP001<br>
3    177235   235235   131155   190190  191195    POP003<br>
15   177239   235235   131155   190190  191191    POP015<br>
16   177254   235235   131155   190190  191195    POP016<br>
17   177235   235235       NA   190190  191195    POP017<br>
19   177235   235235   131155   190190  191195    POP019<br>
> <br>
<br>
pinus.genind <- df2genind(pinus.data[, 1:5], ploidy=2, ncode=3, NA.char = "NA",<br>
                          pop = pinus.data$Populacao, type = "codom")<br>
pinus.genind<br>
tab(pinus.genind)<br>
popNames(pinus.genind)<br>
indNames(pinus.genind)<br>
locNames(pinus.genind)<br>
<br>
> is.genind(pinus.genind)<br>
[1] TRUE<br>
> <br>
<br>
> Fst(as.loci(pinus.genind))<br>
         Fit Fst Fis<br>
PTTX3081 NaN NaN NaN<br>
RIPT0031 NaN NaN NaN<br>
PTTX3011 NaN NaN NaN<br>
PTTX2037 NaN NaN NaN<br>
PTTX311  NaN NaN NaN<br>
<br>
The col "Populacao" in file Analises_pinus_individual.csv is the pop slot. I<br>
consider each individual from a different population.<br>
<br>
Please, could you help me?<br>
<br>
Thank you very much!<br>
<span class="HOEnZb"><font color="#888888"><br>
-- <br>
Laia, ML<br>
______________________________<wbr>_________________<br>
adegenet-forum mailing list<br>
<a href="mailto:adegenet-forum@lists.r-forge.r-project.org">adegenet-forum@lists.r-forge.<wbr>r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum" rel="noreferrer" target="_blank">https://lists.r-forge.r-<wbr>project.org/cgi-bin/mailman/<wbr>listinfo/adegenet-forum</a><br>
</font></span></blockquote></div><br></div>