<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)

}

.:Scrip Funcao LP

.:Help funcao LP