[R] Sampling procedure
Bert Gunter
gunter.berton at gene.com
Thu Oct 15 17:22:09 CEST 2009
If I understand what is wanted correctly, this can be a one-liner! -- think
whole objects:
splitup <- function(x,n.groups)
#split x into n.groups mutually exclusive sets
{
lx <- length(x)
if(n.groups >= lx) stop("Number of groups greater than vector length")
x <- x[sample(lx,lx)]
split(x,seq_len(n.groups))
}
## testit
> splitup(1:71,9)
$`1`
[1] 22 26 38 50 65 60 9 27
$`2`
[1] 24 2 69 28 71 31 41 13
$`3`
[1] 16 47 63 45 23 1 8 32
$`4`
[1] 34 39 64 35 7 19 4 55
$`5`
[1] 54 10 37 68 6 17 70 18
$`6`
[1] 61 11 5 46 33 43 14 56
$`7`
[1] 42 44 12 62 66 48 57 58
$`8`
[1] 21 40 30 29 20 49 52 67
$`9`
[1] 59 15 25 51 3 36 53
Cheers,
Bert Gunter
Genentech Nonclinical Statistics
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of David Winsemius
Sent: Thursday, October 15, 2009 7:55 AM
To: Marcio Resende
Cc: r-help at r-project.org
Subject: Re: [R] Sampling procedure
On Oct 15, 2009, at 10:19 AM, Marcio Resende wrote:
>
> I would like to divide a vector in 9 groups in a way that each
> number is
> present in only one group.
> In a vector of 783 I would like to divide in 9 different groups of 87
>
> Example <- matrix(c(1:783),ncol = 1)
> Example <- matrix(c(1:783),ncol = 1)
> Grp1 <- sample(Example, 87, replace=FALSE)
> Grp2 <- sample(Example[-Grp1], 87, replace=FALSE)
> Grp3 <- sample(Example[-c(Grp1, Grp2)], 87, replace=FALSE)
# lather, rinse , repeat
> s1 <- as.matrix(sample(Example,87, re = FALSE))
> Example <- Example[-s1]
> s2 <- as.matrix(sample(Example,87, re = FALSE))
> #however I don´t know how to remove the second group from the
> "Example" to
> continue sampling.
>
#Don't mess up the original
> There is probably an easy and faster way to do this.
> Could anybody help me?
> Thanks
--
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
______________________________________________
R-help at r-project.org mailing list
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.
More information about the R-help
mailing list