[adegenet-forum] repooling random rows from genind objects

Spencer Bruce goatsrunfaster at gmail.com
Thu Oct 23 17:03:54 CEST 2014


they both say Null, if I just type them into R.

Just to be clear these genind objects contains microsat data for 11 loci
for thousands of individuals.

I'm rather new to R, so I apologize if I'm missing something obvious here...

On Thu, Oct 23, 2014 at 10:54 AM, Jombart, Thibaut <t.jombart at imperial.ac.uk
> wrote:

>
>
>  What does nrow(F1) and other nrow(...)'s say?
>
>
>
>
>   ------------------------------
> *From:* Spencer Bruce [goatsrunfaster at gmail.com]
> *Sent:* 23 October 2014 15:52
> *To:* Jombart, Thibaut
> *Cc:* adegenet-forum at lists.r-forge.r-project.org
> *Subject:* Re: [adegenet-forum] repooling random rows from genind objects
>
>   Error message:
>
>  Error in sample.int(length(x), size, replace, prob) :
>   invalid first argument
>
> On Thu, Oct 23, 2014 at 10:50 AM, Jombart, Thibaut <
> t.jombart at imperial.ac.uk> wrote:
>
>>
>> Hello,
>> hard to figure out what is wrong without the error message..
>> Cheers
>> Thibaut
>>  ------------------------------
>> *From:* adegenet-forum-bounces at lists.r-forge.r-project.org [
>> adegenet-forum-bounces at lists.r-forge.r-project.org] on behalf of Spencer
>> Bruce [goatsrunfaster at gmail.com]
>> *Sent:* 23 October 2014 15:29
>> *To:* adegenet-forum at lists.r-forge.r-project.org
>> *Subject:* [adegenet-forum] repooling random rows from genind objects
>>
>>    Hello All!
>>
>>  I have three seperate populations as genind objects. What I would like
>> to do is pull a certain number of random individuals from each, to form a
>> new single genind population.
>>
>>  I would then like individuals from this new genind population to mate
>> randomly, producing another genind object which would contain their
>> offspring.
>>
>>  Below is the code I came up with (which does not work):
>>
>>  Year1 <- repool(F1[sample(nrow(F1), 500), ], pop1[sample(nrow(pop1),
>> 750), ], pop2[sample(nrow(pop2), 750), ], n=2000)
>>
>>  Year2 <- hybridize(Year1[sample(nrow(Year1), 1000), ],
>> Year1[sample(nrow(Year1), 1000), ], n=2000)
>>
>>
>>  any help would be greatly appreciated!
>>
>>  Best,
>> Spencer
>>
>>  --
>> Spencer A Bruce
>> 200 Washington St.
>> Troy, NY 12180
>> 518 225 0787
>>
>
>
>
>  --
> Spencer A Bruce
> 200 Washington St.
> Troy, NY 12180
> 518 225 0787
>



-- 
Spencer A Bruce
200 Washington St.
Troy, NY 12180
518 225 0787
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20141023/c51ccf09/attachment-0001.html>


More information about the adegenet-forum mailing list