Hi, <div><br></div><div>Thank you for your reply Thibaut. I was able to use your suggestion to produce the 'epid' vector from H3N2.</div><div><br></div><div>I have a follow-up question that will likely highlight my novice abilities with R. When you say that this vector can by "passed" to the 'pop' part of the genind object, I think this is the step I am now having trouble doing successfully. I have created a vector with the desired population names listed for each individual (with then number of entries equaling the number of individuals, as you stated). However, I am not able to append this to the 'pop' or 'other' component of the genind object. </div>
<div><br></div><div>Calling "pop('my genind object')" displays the name of the last unique individual in each separate population group (i.e. separated by POP in the Genepop file). There are 326 such entries.</div>
<div><br></div><div>I created a vector with 326 rows/entries with my preferred names for the scatterplot labels, as described in the previous email ("new labels vector").</div><div><br></div><div>Attempting "pop(my genind object) <- new labels vector", I receive an error:</div>
<div><br></div><div>"Error in `pop<-`(`*tmp*`, value = list(V1 = c(12L, 12L, 12L, 12L, 12L,  : </div><div>  wrong length for population factor"</div><div><br></div><div>Again, any insight you can provide on passing the new population name vector to the 'pop' part of the genind object would be greatly appreciated. Ultimately, I am hoping to replicate the amendment of the scatterplot that you did with the H3N2 epidemic year to change the labels within each ellipse (i.e. page 18 in the DAPC tutorial). Thank you again. </div>
<div><br></div><div>Jonathan</div><div><br></div><div>    <br>
<br><div class="gmail_quote">On Tue, Nov 13, 2012 at 5:28 AM, Jombart, Thibaut <span dir="ltr"><<a href="mailto:t.jombart@imperial.ac.uk" target="_blank">t.jombart@imperial.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Hello,<br>
<br>
Just to clarify, R uses no 'file' for storing objects during a session, everything is stored in the RAM.<br>
<br>
The rule in R is to type the name of an object to print its content. You can use '@' or '$' to access the content of genind objects.<br>
In your case, you're looking for:<br>
##<br>
H3N2@other$epid<br>
##<br>
<br>
which is simply a vector of numeric. You can pass any vector of character, numeric, or factor to 'pop' to change the population information of a genind. The only requirement is for this vector to have a length which is the number of individuals in the study.<br>


<br>
Cheers<br>
<br>
Thibaut<br>
<br>
<br>
________________________________________<br>
From: <a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org" target="_blank">adegenet-forum-bounces@lists.r-forge.r-project.org</a> [<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 J. Richardson [<a href="mailto:jrichardson4@gmail.com" target="_blank">jrichardson4@gmail.com</a>]<br>


Sent: 12 November 2012 21:05<br>
To: <a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.r-project.org</a><br>
Subject: [adegenet-forum] customize labels within centroid<br>
<div><div><br>
Hi Thibaut and Adegenet users,<br>
<br>
I have what I think is a simple/quick question, but have not been able to figure out on my own or via the archives. I have a dataset (Genepop formatted) I have been working with to conduct DAPC and create a scatterplot. I have a legend with each population/sampling site name listed (there are 13), which roughly corresponds to the population/site names within the Genepop file (between each "POP" line). I would like to include these same labels within each of the 13 centroids on the scatterplot, but currently when I set "clab=1" the labels are displayed "P01", "P02"... "P13" rather than the names of the populations/sites in the data file or the object created to put this names in the legend.<br>


<br>
Do you know how I can change these labels, either manually with a new object for the "clab=" command to use, or for "scatter" to refer to in the Genepop input file? I see in the H3N2 example, that the labels were amended to the epidemic years - I think using the "pop(H3N2) <- H3N2$other$epid" line. But I can't figure out how to display the H3N2 data to see how the "epid" item is incorporated into the file.<br>


<br>
Thank you in advance for any advice or insights.<br>
<br>
Jonathan<br>
</div></div></blockquote></div><br></div>