diff --git a/R/DemetrioEx2.12.5.R b/R/DemetrioEx2.12.5.R index f8beaa897b5abaadc78970d5dbe729625864be21..0b34ec77eb6f69f43e586c27411736f7f3567ed2 100644 --- a/R/DemetrioEx2.12.5.R +++ b/R/DemetrioEx2.12.5.R @@ -1,20 +1,21 @@ #' @name DemetrioEx2.12.5 -#' @title Estimativa de Quadrados Mínimos +#' @title Dados Genéricos para Regressão Segmentada #' -#' @description Os dados estão presentes no exercício com o objetivo de -#' estimar na prática o quadrado mínimo. +#' @description Dados para exercício analítico, com o objetivo de obter +#' as estimativas de mínimos quadrados dos parâmetros de um modelo +#' de regressão linear segmentada. #' #' @format Um \code{data.frame} de 5 linhas e 2 colunas. #' #' \describe{ #' -#' \item{\code{x}}{Sem interpretação.} +#' \item{\code{x}}{Variável independente, sem interpretação.} #' -#' \item{\code{y}}{Sem interpretação.} +#' \item{\code{y}}{Variável dependente, sem interpretação.} #' #' } #' -#' @keywords TODO +#' @keywords RegSeg #' #' @source Demétrio, C. G. B., & Zocchi, S. S. (2011). Modelos de #' Regressão. Piracicaba: ESALQ. (Exercício 2.12.5 pág. 60) @@ -26,9 +27,20 @@ #' library(lattice) #' #' xyplot(y ~ x, data = DemetrioEx2.12.5, -#' main = "x vs y", -#' xlab = "x", -#' ylab = "y", -#' type = c("p", "r"), col.line = 3) +#' cex = 1.1, pch = 19, +#' main = 'Regressão Segmentada', +#' grid = TRUE, +#' panel = function(x, y, ...) { +#' # Regressão Segmentada com ponto de corte conhecido +#' b <- DemetrioEx2.12.5$x[3] +#' m0 <- lm(y ~ x + I(pmax(x - b, 0)), +#' data = DemetrioEx2.12.5) +#' # Pontos que definem os dois segmentos +#' cx <- c(0, b, 8) +#' cy <- predict(m0, newdata = data.frame(x = cx)) +#' panel.xyplot(x, y, ...) +#' panel.segments(cx[1], cy[1], cx[2], cy[2]) +#' panel.segments(cx[2], cy[2], cx[3], cy[3]) +#' }) #' NULL diff --git a/man/DemetrioEx2.12.5.Rd b/man/DemetrioEx2.12.5.Rd index e49fedfb57aa4cb17d254d9582fa716c4b5492fb..f32c2adabb0fc5ebca808b75b3f0cd8cac51819c 100644 --- a/man/DemetrioEx2.12.5.Rd +++ b/man/DemetrioEx2.12.5.Rd @@ -2,14 +2,14 @@ % Please edit documentation in R/DemetrioEx2.12.5.R \name{DemetrioEx2.12.5} \alias{DemetrioEx2.12.5} -\title{Estimativa de Quadrados Mínimos} +\title{Dados Genéricos para Regressão Segmentada} \format{Um \code{data.frame} de 5 linhas e 2 colunas. \describe{ - \item{\code{x}}{Sem interpretação.} + \item{\code{x}}{Variável independente, sem interpretação.} - \item{\code{y}}{Sem interpretação.} + \item{\code{y}}{Variável dependente, sem interpretação.} }} \source{ @@ -17,8 +17,9 @@ Demétrio, C. G. B., & Zocchi, S. S. (2011). Modelos de Regressão. Piracicaba: ESALQ. (Exercício 2.12.5 pág. 60) } \description{ -Os dados estão presentes no exercício com o objetivo de - estimar na prática o quadrado mínimo. +Dados para exercício analítico, com o objetivo de obter + as estimativas de mínimos quadrados dos parâmetros de um modelo + de regressão linear segmentada. } \examples{ @@ -27,11 +28,22 @@ data(DemetrioEx2.12.5) library(lattice) xyplot(y ~ x, data = DemetrioEx2.12.5, - main = "x vs y", - xlab = "x", - ylab = "y", - type = c("p", "r"), col.line = 3) + cex = 1.1, pch = 19, + main = 'Regressão Segmentada', + grid = TRUE, + panel = function(x, y, ...) { + # Regressão Segmentada com ponto de corte conhecido + b <- DemetrioEx2.12.5$x[3] + m0 <- lm(y ~ x + I(pmax(x - b, 0)), + data = DemetrioEx2.12.5) + # Pontos que definem os dois segmentos + cx <- c(0, b, 8) + cy <- predict(m0, newdata = data.frame(x = cx)) + panel.xyplot(x, y, ...) + panel.segments(cx[1], cy[1], cx[2], cy[2]) + panel.segments(cx[2], cy[2], cx[3], cy[3]) + }) } -\keyword{TODO} +\keyword{RegSeg}