

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

.:[[ScripfuncaoLP]]

.:[[HelpLP]]