diff --git a/R/PaulinoEx13.16.R b/R/PaulinoEx13.16.R
new file mode 100644
index 0000000000000000000000000000000000000000..d194a279f79b0d58226f8cee1ae6422f1b2d4552
--- /dev/null
+++ b/R/PaulinoEx13.16.R
@@ -0,0 +1,48 @@
+#' @name PaulinoEx13.16
+#' @title Frequências de pacientes com algum tipo de efeitos 
+#'     terapêuticos avaliados em três visitas.
+#' @description Um estudo feito por Molenberghs et al. (1997) trata 
+#'     sobre pacientes do foro psiquiátrico em tratamento com uma droga. 
+#'     Uma das respostas considerada em 3 ocasiões, (X1, X2, X3), indica 
+#'     os efeitos terapêuticos provocados pela droga (manifestados pelo 
+#'     desaparecimento de sintomas anteriormente presentes), registrada 
+#'     dicotomicamente nas categorias 1 (ausência de melhoria ou mesmo 
+#'     pioria ou ainda melhoria mínima sem alterações sintomáticas) e 2 
+#'     (melhoria ligeira ou pronunciada com desaparecimento parcial ou 
+#'     quase total dos sintomas)
+#' @format Um \code{data.frame} com 18 observações e 4 variáveis.
+#'
+#' \describe{
+#'
+#' \item{\code{vis1}}{Primeira visita feita.}
+#'
+#' \item{\code{vis2}}{Segunda visita feita.}
+#' 
+#' \item{\code{vis3}}{Terceira visita feita.}
+#' 
+#' \item{\code{freq}}{Frequência com os valores.}
+#'
+#' }
+#' @keywords TC
+#' @source PAULINO; SINGER (2006), Tabela 13.16, pág. 487.
+#' @references Molenberghs, G. e Goetghebeur, E. J. T. (1997). Simple
+#'     fitting algorithms for incomplete categorical data, Journal of 
+#'     the Royal Statistical Society, B 59:401-414.
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(PaulinoEx13.16)
+#' 
+#' dados <- xtabs(freq ~ vis1 + vis2 + vis3, data = PaulinoEx13.16)
+#' barchart(dados, horizontal = FALSE,
+#'    stack = FALSE,
+#'    scales = list(x = list(labels = c("1","2"))),
+#'    ylab = "Frequência",
+#'    xlab = "Visita 1",
+#'    main = "Frequências de pacientes com algum tipo de efeitos 
+#'    terapêuticos avaliados em três visitas.",
+#'    auto.key = list(corner = c(0.90, 0.95), 
+#'    text = c("1", "2", "omisso")))
+#'        
+NULL
diff --git a/data-raw/PaulinoEx13.16.txt b/data-raw/PaulinoEx13.16.txt
new file mode 100644
index 0000000000000000000000000000000000000000..6163a4627ee16b5af718db96525456cd4a5e0679
--- /dev/null
+++ b/data-raw/PaulinoEx13.16.txt
@@ -0,0 +1,19 @@
+vis1	vis2	vis3	freq
+1	1	1	73
+1	1	2	6
+1	1	omisso	13
+1	2	1	4
+1	2	2	3
+1	2	omisso	5
+1	omisso	1	
+1	omisso	2	
+1	omisso	omisso	10
+2	1	1	73
+2	1	2	3
+2	1	omisso	9
+2	2	1	46
+2	2	2	16
+2	2	omisso	17
+2	omisso	1	
+2	omisso	2	
+2	omisso	omisso	21
diff --git a/data/PaulinoEx13.16.rda b/data/PaulinoEx13.16.rda
new file mode 100644
index 0000000000000000000000000000000000000000..6308b0cae9376ddb787cbe72c71a9db5f2ec23ec
Binary files /dev/null and b/data/PaulinoEx13.16.rda differ
diff --git a/man/PaulinoEx13.16.Rd b/man/PaulinoEx13.16.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..bb210fc406af4433ff8611897891ecdfdc0eff36
--- /dev/null
+++ b/man/PaulinoEx13.16.Rd
@@ -0,0 +1,58 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PaulinoEx13.16.R
+\name{PaulinoEx13.16}
+\alias{PaulinoEx13.16}
+\title{Frequências de pacientes com algum tipo de efeitos 
+    terapêuticos avaliados em três visitas.}
+\format{Um \code{data.frame} com 18 observações e 4 variáveis.
+
+\describe{
+
+\item{\code{vis1}}{Primeira visita feita.}
+
+\item{\code{vis2}}{Segunda visita feita.}
+
+\item{\code{vis3}}{Terceira visita feita.}
+
+\item{\code{freq}}{Frequência com os valores.}
+
+}}
+\source{
+PAULINO; SINGER (2006), Tabela 13.16, pág. 487.
+}
+\description{
+Um estudo feito por Molenberghs et al. (1997) trata 
+    sobre pacientes do foro psiquiátrico em tratamento com uma droga. 
+    Uma das respostas considerada em 3 ocasiões, (X1, X2, X3), indica 
+    os efeitos terapêuticos provocados pela droga (manifestados pelo 
+    desaparecimento de sintomas anteriormente presentes), registrada 
+    dicotomicamente nas categorias 1 (ausência de melhoria ou mesmo 
+    pioria ou ainda melhoria mínima sem alterações sintomáticas) e 2 
+    (melhoria ligeira ou pronunciada com desaparecimento parcial ou 
+    quase total dos sintomas)
+}
+\examples{
+
+library(lattice)
+
+data(PaulinoEx13.16)
+
+dados <- xtabs(freq ~ vis1 + vis2 + vis3, data = PaulinoEx13.16)
+barchart(dados, horizontal = FALSE,
+   stack = FALSE,
+   scales = list(x = list(labels = c("1","2"))),
+   ylab = "Frequência",
+   xlab = "Visita 1",
+   main = "Frequências de pacientes com algum tipo de efeitos 
+   terapêuticos avaliados em três visitas.",
+   auto.key = list(corner = c(0.90, 0.95), 
+   text = c("1", "2", "omisso")))
+       
+}
+\references{
+Molenberghs, G. e Goetghebeur, E. J. T. (1997). Simple
+    fitting algorithms for incomplete categorical data, Journal of 
+    the Royal Statistical Society, B 59:401-414.
+}
+\keyword{TC}
+