ema.r

dim(popB[[7]])

sampleA<-sample(popA$id,75, replace=FALSE) # mean(popA)=0.7
sampleB<-sample(popB[[7]]$id, 75, replace=FALSE) # mean(popB[[7]])=0.7

popB7id <- setdiff(popB[[7]]$id, sampleB)
popB7 <- popB[[7]][popB[[7]]$id %in% popB7id, ]

sampleB_pure <- list()

while (length(sampleB_pure) < 75) {
  sampleB2<-sample(popB7$id,1)
  if (!(sampleB2 %in% sampleA)) {
    sampleB_pure <- c(sampleB_pure, sampleB2)
  }
  popB7id <- setdiff(popB7$id,sampleB2)
  popB7 <- popB7[popB7$id %in% popB7id, ]
}

stopifnot(length(intersect(sampleB_pure, sampleA)) == 0)