# NOT RUN {
# Uses the Lucy data to draw a Bernoulli sample
data(Lucy)
attach(Lucy)
N <- dim(Lucy)[1]
n=400
prob=n/N
sam <- S.BE(N,prob)
# The information about the units in the sample is stored in an object called data
data <- Lucy[sam,]
attach(data)
names(data)
# The variables of interest are: Income, Employees and Taxes
# This information is stored in a data frame called estima
estima <- data.frame(Income, Employees, Taxes)
E.BE(estima,prob)
# }
Run the code above in your browser using DataLab