diff --git a/R/SchwaabTb2.2.1.R b/R/SchwaabTb2.2.1.R index 83f5b476dfa91b5b36cf2698b54311817bfb82fb..4fb918a300647cbb51d7baaa6578d70c9a5187ad 100644 --- a/R/SchwaabTb2.2.1.R +++ b/R/SchwaabTb2.2.1.R @@ -1,37 +1,42 @@ #' @name SchwaabTb2.2.1 -#' @title Efeito da Concentração e Temperatura Sobre Uma Taxa de Reação -#' @description São escolhidos três níveis de temperatura e cinco níveis -#' de concentração distintos. São feitas trẽs réplicas para cada -#' condição experimental. -#' @format Um \code{data.frame} com 45 observações e 3 variáveis, em que +#' @title Efeito da Concentração e Temperatura Sobre uma Taxa de Reação +#' @description Foram escolhidos três níveis de temperatura e cinco níveis +#' de concentração distintos. São feitas três réplicas para cada +#' condição experimental. +#' @format Um \code{data.frame} com 45 observações e 4 variáveis, em que #' #' \describe{ #' -#' \item{\code{temp(°C)}}{Fator numérico que representa a temperatura +#' \item{\code{temp(°C)}}{Fator numérico que representa a temperatura #' em que a reação ocorreu.} #' #' \item{\code{conc(mol/l)}}{Fator numérico que representa o valor da #' concentração.} -#' -#' \item{\code{z1}}{Fator numérico que representa o valor da variável -#' normalizada para a temperatura.} -#' -#' \item{\code{z2}}{Fator numérico que representa o valor da variável -#' normalizada para a concentração.} -#' -#' \item{\code{ordem}}{Fator numérico inteiro que representa o número do -#' ordenamento aleatório.} -#' -#' \item{\code{e}}{Fator numérico que representa um sinal aleatório com -#' distribuição normal, média igual a zero e desvio padrão igual a -#' 0.1.} -#' +#' +#' \item{\code{z1}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A +#' TEMP, APENAS COM RÓTULOS DIFERENTES] Fator numérico que +#' representa o valor da variável normalizada para a temperatura.} +#' +#' \item{\code{z2}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A +#' CONC, APENAS COM RÓTULOS DIFERENTES] Fator numérico que +#' representa o valor da variável normalizada para a concentração.} +#' +#' \item{\code{ordem}}{[OU É FATOR OU NUMÉRICO OU INTEIRO. EVITE INICIAR +#' A DESCRIÇÃO DE UMA VARIÁVEL COM SUA REPRESENTAÇÃO NO R (INTEIRO +#' QUE, NUMÉRICO QUE, FATOR QUE, ETC.)] Fator numérico inteiro que +#' representa o número do ordenamento aleatório.} +#' +#' \item{\code{e}}{[ESSE NÃO É UM RESULTADO? POIS DIFICILMENTE SE +#' OBSERVA UM SINAL ALEATÓRIO] Fator numérico que representa um +#' sinal aleatório com distribuição normal, média igual a zero e +#' desvio padrão igual a 0.1.} +#' #' \item{\code{rc(mol/l.h)}}{Fator numérico que representa o valor da -#' taxa de reação em função da concentração.} -#' +#' taxa de reação em função da concentração.} +#' #' \item{\code{rexp(mol/l.h)}}{Fator numérico que representa o valor da -#' taxa de reação em função da temperatura.} -#' +#' taxa de reação em função da temperatura.} +#' #' } #' @keywords DBC #' @source Schwaab, M., Pinto, J. C.(2007). Análise de Dados @@ -40,16 +45,17 @@ #' #' data(SchwaabTb2.2.1) #' str(SchwaabTb2.2.1) +#' #' library(lattice) -#' -#' xyplot(rc ~ rexp, -#' groups = temp , -#' auto.key = list(title = "Temperaturas", -#' cex.title = 1, -#' columns = 3), -#' type = c("a", "p"), -#' data = SchwaabTb2.2.1, -#' ylab = "Reação em função da temperatura ", -#' xlab = "Reação em função da concentração") -#' -NULL \ No newline at end of file +#' xyplot(rc ~ rexp, +#' groups = temp , +#' auto.key = list( +#' title = "Temperaturas", +#' cex.title = 1, +#' columns = 3), +#' type = c("a", "p"), +#' data = SchwaabTb2.2.1, +#' ylab = "Reação em função da temperatura ", +#' xlab = "Reação em função da concentração") +#' +NULL diff --git a/R/SchwaabTb2.2.15.R b/R/SchwaabTb2.2.15.R index 0c3bc2be1f5da9a0a08557b89bf6c93b4f34ca0f..7bc578ea77012c3fa609abbcda2beebf31066234 100644 --- a/R/SchwaabTb2.2.15.R +++ b/R/SchwaabTb2.2.15.R @@ -6,7 +6,7 @@ #' temperatura e potencial hidrogeniônico (pH) que resultam em uma #' maior atividade enzimática, expressas como gramas de produto #' produzido por um miligrama de enzima em uma hora. -#' @format Um \code{data.frame} com 27 observações e 6 variáveis, em que +#' @format Um \code{data.frame} com 27 observações e 4 variáveis, em que #' #' \describe{ #' diff --git a/R/SchwaabTb2.2.7.R b/R/SchwaabTb2.2.7.R index 4a0920e63efb461d8ba00bde974a2c850a854b63..4b9f0bccf73f3471ed276bdfbcde7c6ef669b6ba 100644 --- a/R/SchwaabTb2.2.7.R +++ b/R/SchwaabTb2.2.7.R @@ -3,9 +3,9 @@ #' Reação #' @description Resultados de um experimento em arranjo fatorial #' \eqn{2^2}, realizado sob delineamento inteiramente casualizado -#' com três réplicadas de cada condição experimental. Os fatores +#' com três réplicas de cada condição experimental. Os fatores #' experimentados foram temperatura e concentração de uma solução -#' não informada, com interesse de avaliar qual o efeito casuado +#' não informada, com interesse de avaliar qual o efeito causado #' sobre a taxa de reação. #' @format Um \code{data.frame} com 12 observações e 4 variáveis, em que #' diff --git a/R/SchwaabTb2.3.5.R b/R/SchwaabTb2.3.5.R index c7c98988d86dfbe7448ca7f5e9f0ac3c439570af..daad3f6b34113786cc208127119acaf8259c1c39 100644 --- a/R/SchwaabTb2.3.5.R +++ b/R/SchwaabTb2.3.5.R @@ -7,51 +7,58 @@ #' #' \describe{ #' -#' \item{\code{temp(°C)}}{Fator numérico que representa a temperatura +#' \item{\code{temp(°C)}}{Fator numérico que representa a temperatura #' em que a reação ocorreu.} -#' -#' \item{\code{z1}}{Fator numérico que representa o valor da variável -#' normalizada para a temperatura.} -#' -#' \item{\code{z2}}{Fator numérico que representa o valor da variável -#' normalizada para a atividade enzimática.} -#' -#' \item{\code{z3}}{Fator numérico que representa o valor da variável -#' normalizada para a velocidade de agitação.} -#' +#' +#' \item{\code{z1}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A +#' TEMP, APENAS COM RÓTULOS DIFERENTES] Fator numérico que +#' representa o valor da variável normalizada para a temperatura.} +#' +#' \item{\code{z2}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A +#' CONC, APENAS COM RÓTULOS DIFERENTES] Fator numérico que +#' representa o valor da variável normalizada para a atividade +#' enzimática.} +#' +#' \item{\code{z3}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A RPM, +#' APENAS COM RÓTULOS DIFERENTES] Fator numérico que representa o +#' valor da variável normalizada para a velocidade de agitação.} +#' #' \item{\code{ordem}}{Fator numérico inteiro que representa o número do #' ordenamento aleatório.} -#' -#' \item{\code{e}}{Fator numérico inteiro que representa o erro -#' experimental.} -#' -#' \item{\code{rc}}{Fator numérico inteiro que representa o valor da -#' reação em função da concentração.} -#' -#' \item{\code{rexp}}{Fator numérico inteiro que representa o valor da -#' reação em função temperatura.} -#' +#' +#' \item{\code{e}}{[ESSE NÃO É UM RESULTADO? POIS DIFICILMENTE SE +#' OBSERVA UM SINAL ALEATÓRIO] Fator numérico inteiro que representa +#' o erro experimental.} +#' +#' \item{\code{rc}}{[PROVAVELMENTE É RESULTADO DE ALGUMA FORMULA FÍSICA +#' TEORICA] Fator numérico inteiro que representa o valor da reação +#' em função da concentração.} +#' +#' \item{\code{rexp}}{[VERIFICAR SE HÁ INDICAÇÃO DE UNIDADE DE MEDIDA, +#' POIS ESTÁ DIFERENTE DE SCHWAABTB2.2.7] Fator numérico inteiro +#' que representa o valor da reação em função temperatura.} +#' #' \item{\code{rpm}}{Fator numérico inteiro que representa o valor da -#' velocidade de agitação.} +#' velocidade de agitação.} #' } -#' @keywords FAT +#' @keywords FAT3K #' @source Schwaab, M., Pinto, J. C.(2007). Análise de Dados #' Experimentais II. Rio de Janeiro, RJ. COPPE UFRJ e-papers. #' @examples #' #' data(SchwaabTb2.3.5) #' str(SchwaabTb2.3.5) -#' +#' #' library(lattice) -#' -#' xyplot(rc ~ rexp, -#' groups = rpm , -#' auto.key = list(title = "rpm", -#' cex.title = 1, -#' columns = 3), -#' type = c("a", "p"), -#' data = SchwaabTb2.3.5, -#' ylab = "Função temperatura ", -#' xlab = "Função da concentração") -#' -NULL \ No newline at end of file +#' +#' xyplot(rc ~ rexp, +#' groups = rpm , +#' auto.key = list(title = "rpm", +#' cex.title = 1, +#' columns = 3), +#' type = c("a", "p"), +#' data = SchwaabTb2.3.5, +#' ylab = "Função temperatura ", +#' xlab = "Função da concentração") +#' +NULL diff --git a/R/SchwaabTb5.5.3.R b/R/SchwaabTb5.5.3.R index 1f714959b0835048522b7b874b9a3a5c42f5ada8..45a3d4a9831ae4085c2abbb5b7709057801d81de 100644 --- a/R/SchwaabTb5.5.3.R +++ b/R/SchwaabTb5.5.3.R @@ -1,37 +1,57 @@ #' @name SchwaabTb5.5.3 -#' @title Reação de Isomeração -#' @description Modelo cinético que representa uma reação de isomeração -#' por uma cinética de primeira ordem. +#' @title Modelo Cinético de uma Reação de Isomerização +#' @description Dados resultantes de um estudo sobre uma reação de +#' isomerização. Foram registrados no estudo o tempo da reação +#' (\code{tempo}), a temperatura sob a qual ela ocorreu +#' (\code{temp})e a fração que restou do reagente após o tempo +#' decorrido (\code{reag}). O interesse no estudo é relacionar a +#' fração que restou do reagente com o tempo e temperatura da +#' reação. +#' +#' Para expressar tal relação existe um modelo cinético, descrito +#' por uma cinética de primeira ordem, que representa uma reação de +#' isomerização conforme descrito na equação abaixo +#' \deqn{\texttt{reag}_i = \exp \left ( -(\alpha 10^{17})\, +#' \texttt{tempo}_i \exp \left ( -\frac{\beta}{\texttt{temp}_i} +#' \right )\right )} +#' em que o índice \eqn{i} representa cada realização do experimento +#' e \eqn{\alpha} e \eqn{\beta} são parâmetros a serem estimados. +#' #' @format Um \code{data.frame} com 41 observações e 3 variáveis, em que #' #' \describe{ #' -#' \item{\code{tempo(s)}}{Fator numérico que representa o tempo em que a -#' reação ocorreu.} +#' \item{\code{tempo}}{Tempo decorrido da reação, mensurado em segundos.} #' -#' \item{\code{temp(K)}}{Fator numérico que representa o valor da -#' temperatura.} -#' -#' \item{\code{reag}}{Fator numérico que representa o valor da -#' fração do reagente.} +#' \item{\code{temp}}{Temperatura sob a qual a reação ocorreu, mensurada +#' na escala kelvin.} +#' +#' \item{\code{reag}}{Fração que restou do reagente após a reação.} #' } -#' @keywords DBC +#' +#' @keywords RNL RM #' @source Schwaab, M., Pinto, J. C.(2007). Análise de Dados #' Experimentais I. Rio de Janeiro, RJ. COPPE UFRJ e-papers. #' @examples #' #' data(SchwaabTb5.5.3) #' str(SchwaabTb5.5.3) +#' #' library(lattice) +#' splom(~SchwaabTb5.5.3) +#' +#' # Ajuste do modelo cinético +#' model <- nls(reag ~ exp(-(alpha * 1e+17) * tempo * exp(-beta/temp)), +#' start = list(alpha = 0.8, beta = 27500), +#' data = SchwaabTb5.5.3) +#' summary(model) +#' +#' # Observados vs preditos +#' xyplot(SchwaabTb5.5.3$reag ~ fitted(model), +#' type = c("g", "p", "r"), +#' panel = function(x, y, ...) { +#' panel.xyplot(x, y, ...) +#' panel.abline(0, 1) +#' }) #' -#' xyplot(reag ~ temp, -#' groups = tempo , -#' type = c("a", "p"), -#' auto.key = list(title = "Tempo(s)", -#' cex.title = 1, -#' columns = 6), -#' data = SchwaabTb5.5.3, -#' ylab = "Reação", -#' xlab = "Temperatura(k)") -#' -NULL \ No newline at end of file +NULL diff --git a/man/SchwaabEx3.3.Rd b/man/SchwaabEx3.3.Rd index c0bd951c95817ded52dd4be60811d05b40394f6b..042ab11df0628b6aeefeb95c697196b5d3c15cfe 100644 --- a/man/SchwaabEx3.3.Rd +++ b/man/SchwaabEx3.3.Rd @@ -2,24 +2,24 @@ % Please edit documentation in R/SchwaabEx3.3.R \name{SchwaabEx3.3} \alias{SchwaabEx3.3} -\title{Avaliação de Operadores numa Empresa Quimica} +\title{Avaliação de Operadores numa Empresa Química} \format{Um \code{data.frame} com 80 observações e 2 variáveis, em que \describe{ -\item{\code{turma}}{Fator de 4 níveis qualitativos, que são as - turmas avaliadas.} +\item{\code{turma}}{Fator de 4 níveis que representa a turma do + operador avaliado.} -\item{\code{dados}}{Fator numérico que representa os valores dos - dados de conversão do reator.} +\item{\code{desemp}}{Escore de desempenho do avaliador na operação do + reator.} }} \source{ -Schwaab, M., Pinto, J. C.(2007). Análise de Dados +Schwaab, M., Pinto, J. C.(2007) . Análise de Dados Experimentais I. Rio de Janeiro, RJ. COPPE UFRJ e-papers. } \description{ -Quatro turmas de operadores devem ser avaliados em +Estudo que avaliou operadores de quatro turmas em relação ao desempenho na operação de um reator. } \examples{ @@ -27,12 +27,15 @@ Quatro turmas de operadores devem ser avaliados em data(SchwaabEx3.3) str(SchwaabEx3.3) +meds <- with(SchwaabEx3.3, tapply(desemp, turma, summary)) +do.call(cbind, meds) + library(lattice) -bwplot(dados ~ turma, - data =SchwaabEx3.3 , - ylab = "Dados de conversão no reator", - xlab = "Turmas") +bwplot(desemp ~ turma, + data = SchwaabEx3.3, + ylab = "Desempenho do operador", + xlab = "Turma") } -\keyword{DBC} +\keyword{AASI} diff --git a/man/SchwaabTb2.2.1.Rd b/man/SchwaabTb2.2.1.Rd index 16e522d03a6454e80cc9c05cc3d030f98e0957fe..b9216367afc5cd6c362befde5edd9e02e915e172 100644 --- a/man/SchwaabTb2.2.1.Rd +++ b/man/SchwaabTb2.2.1.Rd @@ -2,61 +2,67 @@ % Please edit documentation in R/SchwaabTb2.2.1.R \name{SchwaabTb2.2.1} \alias{SchwaabTb2.2.1} -\title{Efeito da Concentração e Temperatura Sobre Uma Taxa de Reação} -\format{Um \code{data.frame} com 45 observações e 3 variáveis, em que +\title{Efeito da Concentração e Temperatura Sobre uma Taxa de Reação} +\format{Um \code{data.frame} com 45 observações e 4 variáveis, em que \describe{ -\item{\code{temp(°C)}}{Fator numérico que representa a temperatura +\item{\code{temp(°C)}}{Fator numérico que representa a temperatura em que a reação ocorreu.} \item{\code{conc(mol/l)}}{Fator numérico que representa o valor da concentração.} - -\item{\code{z1}}{Fator numérico que representa o valor da variável - normalizada para a temperatura.} - -\item{\code{z2}}{Fator numérico que representa o valor da variável - normalizada para a concentração.} - -\item{\code{ordem}}{Fator numérico inteiro que representa o número do - ordenamento aleatório.} - -\item{\code{e}}{Fator numérico que representa um sinal aleatório com - distribuição normal, média igual a zero e desvio padrão igual a - 0.1.} - + +\item{\code{z1}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A + TEMP, APENAS COM RÓTULOS DIFERENTES] Fator numérico que + representa o valor da variável normalizada para a temperatura.} + +\item{\code{z2}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A + CONC, APENAS COM RÓTULOS DIFERENTES] Fator numérico que + representa o valor da variável normalizada para a concentração.} + +\item{\code{ordem}}{[OU É FATOR OU NUMÉRICO OU INTEIRO. EVITE INICIAR + A DESCRIÇÃO DE UMA VARIÁVEL COM SUA REPRESENTAÇÃO NO R (INTEIRO + QUE, NUMÉRICO QUE, FATOR QUE, ETC.)] Fator numérico inteiro que + representa o número do ordenamento aleatório.} + +\item{\code{e}}{[ESSE NÃO É UM RESULTADO? POIS DIFICILMENTE SE + OBSERVA UM SINAL ALEATÓRIO] Fator numérico que representa um + sinal aleatório com distribuição normal, média igual a zero e + desvio padrão igual a 0.1.} + \item{\code{rc(mol/l.h)}}{Fator numérico que representa o valor da - taxa de reação em função da concentração.} - + taxa de reação em função da concentração.} + \item{\code{rexp(mol/l.h)}}{Fator numérico que representa o valor da - taxa de reação em função da temperatura.} - + taxa de reação em função da temperatura.} + }} \source{ Schwaab, M., Pinto, J. C.(2007). Análise de Dados Experimentais I. Rio de Janeiro, RJ. COPPE UFRJ e-papers. } \description{ -São escolhidos três níveis de temperatura e cinco níveis - de concentração distintos. São feitas trẽs réplicas para cada +Foram escolhidos três níveis de temperatura e cinco níveis + de concentração distintos. São feitas três réplicas para cada condição experimental. } \examples{ data(SchwaabTb2.2.1) str(SchwaabTb2.2.1) -library(lattice) - xyplot(rc ~ rexp, - groups = temp , - auto.key = list(title = "Temperaturas", - cex.title = 1, - columns = 3), - type = c("a", "p"), - data = SchwaabTb2.2.1, - ylab = "Reação em função da temperatura ", - xlab = "Reação em função da concentração") +library(lattice) +xyplot(rc ~ rexp, + groups = temp , + auto.key = list( + title = "Temperaturas", + cex.title = 1, + columns = 3), + type = c("a", "p"), + data = SchwaabTb2.2.1, + ylab = "Reação em função da temperatura ", + xlab = "Reação em função da concentração") } \keyword{DBC} diff --git a/man/SchwaabTb2.2.15.Rd b/man/SchwaabTb2.2.15.Rd index f10d32f7373f1419f7f47ecb34b702247f4a4afb..fb90f73a5708318d5d39bb5fa0001422da38f99b 100644 --- a/man/SchwaabTb2.2.15.Rd +++ b/man/SchwaabTb2.2.15.Rd @@ -2,55 +2,62 @@ % Please edit documentation in R/SchwaabTb2.2.15.R \name{SchwaabTb2.2.15} \alias{SchwaabTb2.2.15} -\title{Atividade Enzimática} -\format{Um \code{data.frame} com 27 observações e 6 variáveis, em que +\title{Atividade Enzimática em Problema Biotecnológico} +\format{Um \code{data.frame} com 27 observações e 4 variáveis, em que \describe{ -\item{\code{temp(°C)}}{Fator numérico que representa a temperatura - em que a reação ocorreu.} - -\item{\code{ph}}{Fator numérico que representa o grau de atividade da - enzima.} - -\item{\code{z1}}{Fator numérico que representa o valor da variável - normalizada para a temperatura.} - -\item{\code{z2}}{Fator numérico que representa o valor da variável - normalizada para a atividade enzimática.} - -\item{\code{ordem}}{Fator numérico inteiro que representa o número do - ordenamento aleatório.} - -\item{\code{ativ(g/mgh)}}{Fator numérico que representa o valor da - atividade enzimática.} - +\item{\code{ordem}}{Identifica a ordem de coleta das observações.} + +\item{\code{temp}}{Níveis de temperatura experimentados, registrados + em graus Celsius.} + +\item{\code{ph}}{Níveis de potencial hidrogeniônico (pH) + experimentados.} + +\item{\code{ativ}}{Atividade enzimática do solo, mensurada pelas + gramas de produto produzido por um miligrama de enzima em uma + hora \eqn{g\,(mg\,h)^{-1}}.} + }} \source{ Schwaab, M., Pinto, J. C.(2007). Análise de Dados Experimentais II. Rio de Janeiro, RJ. COPPE UFRJ e-papers. } \description{ -Plano experimental a três níveis para identificar pontos - de máxima atividade enzimática em um problema biotecnológico. +Resultados de um experimento fatorial \eqn{3^2}, + realizado sob delineamento inteiramente casualizado com três + repetições, com o objetivo de identificar as combinações de + temperatura e potencial hidrogeniônico (pH) que resultam em uma + maior atividade enzimática, expressas como gramas de produto + produzido por um miligrama de enzima em uma hora. +} +\note{ +Embora a variável \code{ativ} represente produção, valores + negativos foram registrados. Não há justificativas sobre tais + valores na descrição do experimento. } \examples{ data(SchwaabTb2.2.15) str(SchwaabTb2.2.15) -library(lattice) +means <- aggregate(ativ ~ temp + ph, data = SchwaabTb2.2.15, mean) +xtabs(ativ ~ ., data = means) - xyplot(temp ~ ativ, - groups = ph , - auto.key = list(title = "pH", - cex.title = 1, - columns = 3), - type = c("a", "p"), - data = SchwaabTb2.2.15, - ylab = "Temperatura ", - xlab = "Atividade") +library(lattice) +xyplot(ativ ~ temp, groups = ph, + data = SchwaabTb2.2.15, + type = c("g", "p", "a"), + xlab = expression("Temperatura em"~degree*"C"), + ylab = expression("Atividade enzimátiva em"~g~(mg~h)^{-1}), + auto.key = list( + title = "pH", + cex.title = 1.1, + column = 3, + lines = TRUE + )) } -\keyword{FAT} +\keyword{FAT2K} diff --git a/man/SchwaabTb2.2.7.Rd b/man/SchwaabTb2.2.7.Rd index 000cad571ba5a6db138a6b9c171178d5ff77f922..fac89be660bae85c1f21c0e25e5d08ece3550e09 100644 --- a/man/SchwaabTb2.2.7.Rd +++ b/man/SchwaabTb2.2.7.Rd @@ -2,63 +2,69 @@ % Please edit documentation in R/SchwaabTb2.2.7.R \name{SchwaabTb2.2.7} \alias{SchwaabTb2.2.7} -\title{Efeito da Concentração e Temperatura Sobre Uma Taxa de Reação} -\format{Um \code{data.frame} com 12 observações e 8 variáveis, em que +\title{Efeito da Concentração Molar e Temperatura Sobre Taxa de + Reação} +\format{Um \code{data.frame} com 12 observações e 4 variáveis, em que \describe{ -\item{\code{temp(°C)}}{Fator numérico que representa a temperatura - em que a reação ocorreu.} + \item{\code{ordem}}{Identifica a ordem de coleta das observações.} + +\item{\code{temp}}{Níveis de temperatura em que a reação ocorreu, + registrados em graus Celsius.} + +\item{\code{conc}}{Níveis de concentração em \eqn{mol/L} da solução.} + +\item{\code{reac}}{Taxa de reação, mensurada em \eqn{mol/L.h}.} -\item{\code{conc(mol/l)}}{Fator numérico que representa o valor da - concentração.} - -\item{\code{z1}}{Fator numérico que representa o valor da variável - normalizada para a temperatura.} - -\item{\code{z2}}{Fator numérico que representa o valor da variável - normalizada para a concentração.} - -\item{\code{ordem}}{Fator numérico inteiro que representa o número do - ordenamento aleatório.} - -\item{\code{e}}{Fator numérico que representa um sinal aleatório com - distribuição normal, média igual a zero e desvio padrão igual a - 0.1.} - -\item{\code{rc(mol/l.h)}}{Fator numérico que representa o valor da - taxa de reação em função da concentração.} - -\item{\code{rexp(mol/l.h)}}{Fator numérico que representa o valor da - taxa de reação em função da temperatura.} - }} \source{ Schwaab, M., Pinto, J. C.(2007). Análise de Dados Experimentais II. Rio de Janeiro, RJ. COPPE UFRJ e-papers. } \description{ -São escolhidos três níveis de temperatura e cinco níveis - de concentração distintos. São feitas trẽs réplicas para cada - condição experimental fatorial. +Resultados de um experimento em arranjo fatorial + \eqn{2^2}, realizado sob delineamento inteiramente casualizado + com três réplicas de cada condição experimental. Os fatores + experimentados foram temperatura e concentração de uma solução + não informada, com interesse de avaliar qual o efeito causado + sobre a taxa de reação. +} +\note{ +Na tabela 2.7 do Capítulo 2 de Schwaab e Pinto (2007) são + exibidos também quantidades denominadas como \eqn{\varepsilon}, + erro aleatório e \eqn{R^c}. Essas são quantidades decorrentes da + equação de Arrhenius, que nesse experimento, descreve + teoricamente a taxa de reação como + \deqn{R^c = 1000\exp \left ( -\frac{2500}{\texttt{temp}+273,15} + \right ) \texttt{conc}} + Assim o valor observado da taxa de reação, \code{reac} pode ser + obtido com \eqn{R^c + \varepsilon}. Nesse \emph{dataset} optou-se + somente pela inclusão da taxa de reação observada uma vez que + essa é a quantidade efetivamente utilizada nas análises e as + demais podem ser obtidas a partir dela. } \examples{ data(SchwaabTb2.2.7) str(SchwaabTb2.2.7) -library(lattice) +means <- aggregate(reac ~ temp + conc, data = SchwaabTb2.2.7, mean) +xtabs(reac ~ ., data = means) - xyplot(rc ~ rexp, - groups = temp , - auto.key = list(title = "Temperaturas", - cex.title = 1, - columns = 2), - type = c("a", "p"), - data = SchwaabTb2.2.7, - ylab = "Reação em função da temperatura ", - xlab = "Reação em função da concentração") +library(lattice) +xyplot(reac ~ temp, groups = conc, + data = SchwaabTb2.2.7, + type = c("g", "p", "a"), + xlab = expression("Temperatura em"~degree*"C"), + ylab = "Taxa de reação", + auto.key = list( + corner = c(0.1, 0.95), + title = "Concentração em mol/L", + cex.title = 1, + lines = TRUE + )) } -\keyword{FAT} +\keyword{FAT2K} diff --git a/man/SchwaabTb2.3.5.Rd b/man/SchwaabTb2.3.5.Rd index e6ea492726775b9ed50fed754445ec36c49485c9..d578774fc3b94cd99eec68a86eaaeea79a53a6f6 100644 --- a/man/SchwaabTb2.3.5.Rd +++ b/man/SchwaabTb2.3.5.Rd @@ -7,32 +7,39 @@ \describe{ -\item{\code{temp(°C)}}{Fator numérico que representa a temperatura +\item{\code{temp(°C)}}{Fator numérico que representa a temperatura em que a reação ocorreu.} - -\item{\code{z1}}{Fator numérico que representa o valor da variável - normalizada para a temperatura.} - -\item{\code{z2}}{Fator numérico que representa o valor da variável - normalizada para a atividade enzimática.} - -\item{\code{z3}}{Fator numérico que representa o valor da variável - normalizada para a velocidade de agitação.} - + +\item{\code{z1}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A + TEMP, APENAS COM RÓTULOS DIFERENTES] Fator numérico que + representa o valor da variável normalizada para a temperatura.} + +\item{\code{z2}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A + CONC, APENAS COM RÓTULOS DIFERENTES] Fator numérico que + representa o valor da variável normalizada para a atividade + enzimática.} + +\item{\code{z3}}{[DEVE SER REMOVIDA POIS ESSA VARIÁVEL É IGUAL A RPM, + APENAS COM RÓTULOS DIFERENTES] Fator numérico que representa o + valor da variável normalizada para a velocidade de agitação.} + \item{\code{ordem}}{Fator numérico inteiro que representa o número do ordenamento aleatório.} - -\item{\code{e}}{Fator numérico inteiro que representa o erro - experimental.} - -\item{\code{rc}}{Fator numérico inteiro que representa o valor da - reação em função da concentração.} - -\item{\code{rexp}}{Fator numérico inteiro que representa o valor da - reação em função temperatura.} - + +\item{\code{e}}{[ESSE NÃO É UM RESULTADO? POIS DIFICILMENTE SE + OBSERVA UM SINAL ALEATÓRIO] Fator numérico inteiro que representa + o erro experimental.} + +\item{\code{rc}}{[PROVAVELMENTE É RESULTADO DE ALGUMA FORMULA FÍSICA + TEORICA] Fator numérico inteiro que representa o valor da reação + em função da concentração.} + +\item{\code{rexp}}{[VERIFICAR SE HÁ INDICAÇÃO DE UNIDADE DE MEDIDA, + POIS ESTÁ DIFERENTE DE SCHWAABTB2.2.7] Fator numérico inteiro + que representa o valor da reação em função temperatura.} + \item{\code{rpm}}{Fator numérico inteiro que representa o valor da - velocidade de agitação.} + velocidade de agitação.} }} \source{ Schwaab, M., Pinto, J. C.(2007). Análise de Dados @@ -50,16 +57,16 @@ str(SchwaabTb2.3.5) library(lattice) - xyplot(rc ~ rexp, - groups = rpm , - auto.key = list(title = "rpm", - cex.title = 1, - columns = 3), - type = c("a", "p"), - data = SchwaabTb2.3.5, - ylab = "Função temperatura ", - xlab = "Função da concentração") +xyplot(rc ~ rexp, + groups = rpm , + auto.key = list(title = "rpm", + cex.title = 1, + columns = 3), + type = c("a", "p"), + data = SchwaabTb2.3.5, + ylab = "Função temperatura ", + xlab = "Função da concentração") } -\keyword{FAT} +\keyword{FAT3K} diff --git a/man/SchwaabTb5.5.3.Rd b/man/SchwaabTb5.5.3.Rd index 7ced453431cdc953c02c4330cf31f1e0442c99bd..4cae08a34a9ffcf88f9f94816029fdd6cd3c9c0d 100644 --- a/man/SchwaabTb5.5.3.Rd +++ b/man/SchwaabTb5.5.3.Rd @@ -2,44 +2,63 @@ % Please edit documentation in R/SchwaabTb5.5.3.R \name{SchwaabTb5.5.3} \alias{SchwaabTb5.5.3} -\title{Reação de Isomeração} +\title{Modelo Cinético de uma Reação de Isomerização} \format{Um \code{data.frame} com 41 observações e 3 variáveis, em que \describe{ -\item{\code{tempo(s)}}{Fator numérico que representa o tempo em que a - reação ocorreu.} +\item{\code{tempo}}{Tempo decorrido da reação, mensurado em segundos.} -\item{\code{temp(K)}}{Fator numérico que representa o valor da - temperatura.} - -\item{\code{reag}}{Fator numérico que representa o valor da - fração do reagente.} +\item{\code{temp}}{Temperatura sob a qual a reação ocorreu, mensurada + na escala kelvin.} + +\item{\code{reag}}{Fração que restou do reagente após a reação.} }} \source{ Schwaab, M., Pinto, J. C.(2007). Análise de Dados Experimentais I. Rio de Janeiro, RJ. COPPE UFRJ e-papers. } \description{ -Modelo cinético que representa uma reação de isomeração - por uma cinética de primeira ordem. +Dados resultantes de um estudo sobre uma reação de + isomerização. Foram registrados no estudo o tempo da reação + (\code{tempo}), a temperatura sob a qual ela ocorreu + (\code{temp})e a fração que restou do reagente após o tempo + decorrido (\code{reag}). O interesse no estudo é relacionar a + fração que restou do reagente com o tempo e temperatura da + reação. + + Para expressar tal relação existe um modelo cinético, descrito + por uma cinética de primeira ordem, que representa uma reação de + isomerização conforme descrito na equação abaixo + \deqn{\texttt{reag}_i = \exp \left ( -(\alpha 10^{17})\, + \texttt{tempo}_i \exp \left ( -\frac{\beta}{\texttt{temp}_i} + \right )\right )} + em que o índice \eqn{i} representa cada realização do experimento + e \eqn{\alpha} e \eqn{\beta} são parâmetros a serem estimados. } \examples{ data(SchwaabTb5.5.3) str(SchwaabTb5.5.3) + library(lattice) +splom(~SchwaabTb5.5.3) + +# Ajuste do modelo cinético +model <- nls(reag ~ exp(-(alpha * 1e+17) * tempo * exp(-beta/temp)), + start = list(alpha = 0.8, beta = 27500), + data = SchwaabTb5.5.3) +summary(model) + +# Observados vs preditos +xyplot(SchwaabTb5.5.3$reag ~ fitted(model), + type = c("g", "p", "r"), + panel = function(x, y, ...) { + panel.xyplot(x, y, ...) + panel.abline(0, 1) + }) -xyplot(reag ~ temp, - groups = tempo , - type = c("a", "p"), - auto.key = list(title = "Tempo(s)", - cex.title = 1, - columns = 6), - data = SchwaabTb5.5.3, - ylab = "Reação", - xlab = "Temperatura(k)") - } -\keyword{DBC} +\keyword{RM} +\keyword{RNL}