set.seed(1234)

library(Matching)

# population
data(lalonde)
lalonde$PNR <- 1:(length(lalonde$age))

# sample of population
population_ID <- sample(1:length(lalonde$PNR), 350)
sample <- lalonde[population,]