Assuming Michael is correct, you can use setdiff(): > set.seed(42) > current <- sample.int(4500, 4495) # All but 5 numbers used > setdiff(1:4500, current) # Find which numbers are left [1] 905 1252 2508 3192 4484
------------------------------------- David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -----Original Message----- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Michael Dewey Sent: Thursday, January 21, 2016 7:34 AM To: maryam firoozi; r-help@r-project.org Subject: Re: [R] Fwd: Hello Maryam See below On 20/01/2016 21:26, maryam firoozi via R-help wrote: > Hello, > i made a population about 4500 individual. So you mean nearly 4500, not exactly 4500? this has two sex(female and male).they had pedigree. > i wanted to enter new indiviual but their ID of indiviual mustnot be same > perivous and their ID number mustnot be bigger than 4500. > first population's ID number is 1:4500. If you tabulate ID you will get a table where all the entries should be 1 except for the ID which you have not yet used which will be zero. ?tabulate _NOT_ ?table All you need to do then is find them ?which should work > how can i handel it? > > Sent from my iPhone > > Begin forwarded message: > >> From: maryam firoozi <firoozi_maryam6...@yahoo.com> >> Date: January 18, 2016 at 1:14:56 AM GMT+3:30 >> To: r-help@r-project.org >> >> hello, >> we want to do genomic blup in r.i know that use pedigree package. >> the formule is >> gblup( P~1,data=ped[,c('ID','P')],M=M,lambda=1/h2-1) >> P:phenotype variance >> ped:pedigree >> M: matrix marker or genotype >> my ped has 4500 ID.but my M has 9000 individual.becasue i have two row for >> each ID in M matrix becasue each ID has two haplotype.how can i solve it.the >> formula didnt solve. >> sincerely >> >> >> Sent from my iPhone > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Michael http://www.dewey.myzen.co.uk/home.html ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
