Estatstica Bsica e Experimentao no RWalmes Marques Zeviani
Sumrio1 Introduo manipulao de objetos e funes 1.1 Instalao do R
. . . . . . . . . . . . . . . . . . . . . 1.2 Pedindo ajuda . . . .
. . . . . . . . . . . . . . . . . . 1.3 Criao, atribuio, acesso e
modicao de objetos 1.4 Informaes sobre objetos (atributos) . . . .
. . . . . 1.5 Operaes matemticas . . . . . . . . . . . . . . . .
1.6 Operaes estatsticas . . . . . . . . . . . . . . . . . . 1.7
Construindo funes . . . . . . . . . . . . . . . . . . Importao de
dados e anlise exploratria 2.1 Impotando dados . . . . . . . . . .
. . . 2.2 Exploraes grcas (1) . . . . . . . . . . 2.3 Exploraes
grcas (2) . . . . . . . . . . 2.4 Recursos grcos avanados . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 3 3 4 5
6 6 7 8 8 9 9 10 11 11 12 12 13 13 14 14 15 15 15 16 16 17 18 19 21
22
2
3
Regresso linear 3.1 Importando e manipulando dados . . . . . . .
. . . . . . . . 3.2 Regresso linear simples . . . . . . . . . . . .
. . . . . . . . . 3.3 Regresso linear mltipla . . . . . . . . . . .
. . . . . . . . . 3.4 Seleo de modelos/variveis . . . . . . . . . .
. . . . . . . . 3.5 Remoo de pontos discrepantes . . . . . . . . .
. . . . . . . 3.6 Predio de valores a partir do modelo escolhido .
. . . . . . 3.7 Representao grca do ajuste . . . . . . . . . . . .
. . . . . 3.8 Mais sobre anlise de resduos e R2 (quarteto de
Anscombe)
4
Regresso no linear 4.1 Motivao . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . 4.2 Denio . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . 4.3 Exemplo de
modelos no lineares . . . . . . . . . . . . . . . . . . . . . 4.4
Uso de recursos grcos para entender o signicado dos parmetros 4.5
Estimao de parmetros em modelos no lineares . . . . . . . . . . .
4.6 Ajuste de modelo no linear aos dados de DAP . . . . . . . . . .
. . . 4.7 Comparao de curvas ajustadas . . . . . . . . . . . . . .
. . . . . . . 4.8 Ajuste de modelos no lineares com a library{nlme}
. . . . . . . . . .
Documento concludo em 1 de dezembro de 2010 s 18:08:02 Centro
Politcnico Universidade Federal do Paran.
1
5
Anlise de experimento com um fator em DIC 5.1 Importando dados .
. . . . . . . . . . . . . . . . . . 5.2 Anlise de varincia . . . .
. . . . . . . . . . . . . . 5.3 Aplicando teste de Tukey para
comparar mdias . . 5.4 Aplicando teste de Scott-Knott para agrupar
mdias 5.5 Aplicando contrastes . . . . . . . . . . . . . . . . . .
5.6 Estudo das taxas de erro tipo I dos testes . . . . . . . Anlise
de experimentos de um fator em DBC 6.1 Entrada de dados . . . . . .
. . . . . . . . 6.2 Anlise de varincia . . . . . . . . . . . . 6.3
Teste de mdias . . . . . . . . . . . . . . . 6.4 Observaes perdidas
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
.
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
23 23 24 24 24 25 25 26 26 26 27 27 28 28 29 30 30 30 31 32 33
33 34 35 35 36 37 37 37 38 38 38 39 39 40 42 42 44
6
7
Anlise de experimento fatorial duplo em DIC 7.1 Anlise de
varincia . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . 7.2 Testes de mdias . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . Anlise de fatorial duplo em
DBC 8.1 Entrando com os dados . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . 8.2 Anlise de varincia e
desdobramento das somas de quadrados . . . . . . . . . . 8.3
Desdobramento da interao com testes de mdias . . . . . . . . . . .
. . . . . . . Anlise de experimento fatorial com um tratamento
adicional
8
9
10 Anlise de covarincia 10.1 Anlise de varincia . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . 10.2
Constraste entre nveis dos fatores . . . . . . . . . . . . . . . .
. . . . . . . . . . . 11 Experimento fatorial com fatores
qualitativos e quantitativos 11.1 Desdobramento da interao . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . 11.2 Obteno das
equaes de regresso e R2 . . . . . . . . . . . . . . . . . . . . . .
. 12 Fatorial com fatores quantitativos - superfcie de resposta
12.1 Anlise de varincia e obteno do modelo emprico . . . . . . . .
. . . . . . . . . 12.2 Grco do modelo nal . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . 13 Anlise de experimentos
em parcela subdividida 13.1 Anlise de varincia . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . 13.2 Teste de
mdias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . 14 Experimentos em parcelas subsubdivididas 14.1
Anlise de varincia . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . 14.2 Testes de mdias . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . 15 Recursos
grcos 15.1 Grcos do pacote graphics . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . 15.2 Grcos do pacote lattice . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2
11.1
Introduo manipulao de objetos e funesInstalao do R
.
#-----------------------------------------------------------------------------------------#
pgina do R browseURL(URLencode("http://www.r-project.org/")) #
#-----------------------------------------------------------------------------------------#
pgina de download
browseURL(URLencode("http://cran.stat.ucla.edu/bin/windows/base/"))
#
#-----------------------------------------------------------------------------------------#
documento com instrues de intalao e primeiros passos
browseURL(URLencode("http://cran.r-project.org/doc/contrib/Itano-installation.pdf"))
#
#-----------------------------------------------------------------------------------------#
curiosidades
browseURL(URLencode("http://www.nytimes.com/2009/01/07/technology/business-computing/07program.html"))
browseURL(URLencode("http://jeromyanglim.blogspot.com/2010/05/abbreviations-of-r-commands-explained.html"))
#
#-----------------------------------------------------------------------------------------#
quando no souber os links consulte o google
browseURL(URLencode("http://www.lmgtfy.com/?q=R+download+windows"))
browseURL(URLencode("http://www.lmgtfy.com/?q=R+reference+card")) #
#-----------------------------------------------------------------------------------------.
1.2
Pedindo ajuda
.
#-----------------------------------------------------------------------------------------#
quando voc s sabe algo sobre apropos("tukey") apropos("help") #
#-----------------------------------------------------------------------------------------#
fazendo a busca do termo help(TukeyHSD) help(TukeyHSD,
help_type="html") ?TukeyHSD #
#-----------------------------------------------------------------------------------------#
buscando em pacotes o termo help.search("Tukey") ??Tukey #
#-----------------------------------------------------------------------------------------#
fazendo a busca na web RSiteSearch("curve fitting") #
#-----------------------------------------------------------------------------------------.
3
1.3
Criao, atribuio, acesso e modicao de objetos
.
#-----------------------------------------------------------------------------------------#
vetores, sequncias e nmeros aleatrios c(2,4,7,3,8,9) 1:7 seq(0, 20,
by=3.1) seq(0, 20, length=4) rep(1:3, times=3) rep(1:3, each=3)
rnorm(5, 3, 2) rnorm(5, sd=2, mean=3) rnorm(5, mean=3, sd=2)
runif(5) #
#-----------------------------------------------------------------------------------------#
matrizes matrix(c(1,5,38,400), 2, 2) matrix(1:6, 2, 3)
matrix(rnorm(9), 3, 3) matrix(c("a","c","b","j"), 2, 2) #
#-----------------------------------------------------------------------------------------#
estrutura de dados (planilha) data.frame(A=1:4, B=runif(4),
C=letters[1:4]) data.frame(trat=c(1:2,1:2), bloc=rep(1:2, e=2))
expand.grid(cult=c("A","B"), bloc=c("I","II","III"), dose=1:3) #
#-----------------------------------------------------------------------------------------#
listas list(A=rnorm(4), B=matrix(1:4,2,2), C=data.frame(a=1:4,
b=runif(4), c=letters[1:4]), D="O R livre") #
#-----------------------------------------------------------------------------------------#
atribuio, acesso e modificao de vetores x
