<code>
##funcao
LP ← function(data) {
N <- c()
for (i in 1:length(unique(data$inq))) {
tmp <- subset(data, inq == i)
N[i] <- with(tmp, ( (sum(n1) + 1)*(sum(n2) + 1) / (sum(m2) + 1) ) - 1)
}
x11() # abre um dispositivo grafico
par(mfrow = c(3,1))
plot((N),main=("Plot de Lincoln Peterson"), xlab="inquerito")
barplot((N),main=("Barplot de Lincoln Peterson"))
boxplot((N),main=("Boxplot de Lincoln Peterson"))
return(N)
}