# NOT RUN {
data('BigCity')
Households <- BigCity %>% group_by(HHID) %>%
summarise(Stratum = unique(Stratum),
PSU = unique(PSU),
Persons = n(),
Income = sum(Income),
Expenditure = sum(Expenditure))
attach(Households)
UI <- levels(as.factor(Households$PSU))
NI <- length(UI)
nI <- 100
samI <- S.SI(NI, nI)
sampleI <- UI[samI]
CityI <- Households[which(Households$PSU %in% sampleI), ]
attach(CityI)
area <- as.factor(CityI$PSU)
estima <- data.frame(CityI$Persons, CityI$Income, CityI$Expenditure)
E.1SI(NI, nI, estima, area)
# }
Run the code above in your browser using DataLab