Guilherme

• Session settings

pkg <- c("ggplot2", "knitr","reshape2", "dplyr")
sapply(pkg, library, character.only=TRUE, logical.return=TRUE)

##  ggplot2    knitr reshape2    dplyr 
##     TRUE     TRUE     TRUE     TRUE

setwd("/home/epi/Dropbox/MyR/Análise otros/Gui")# lab

• Dataset

# Dados planilha laboratorio
dat = read.csv("gui.csv", dec=",", header=T, sep="\t", check.names=FALSE, na.strings=".")  
head(dat)

##   var dis rep fruto tam        ic pi
## 1 Ham  fm   1     1 1.5 -16.74434 57
## 2 Ham  fm   1     2 1.5 -17.55321 57
## 3 Ham  fm   1     3 1.5 -19.03078 57
## 4 Ham  fm   1     4 1.5 -17.93105 57
## 5 Ham  fm   1     5 1.5 -17.49463 57
## 6 Ham  fm   1     6 1.5 -17.83131 57

for (i in 1:3) dat[,i] <- as.factor(dat[,i])
for (i in 4:7) dat[,i] <- as.numeric(dat[,i])
str(dat)

## 'data.frame':    203 obs. of  7 variables:
##  $ var  : Factor w/ 3 levels "Ham","Pera","Val": 1 1 1 1 1 1 1 1 1 1 ...
##  $ dis  : Factor w/ 2 levels "fm","md": 1 1 1 1 1 1 1 1 1 1 ...
##  $ rep  : Factor w/ 3 levels "1","2","3": 1 1 1 1 1 1 2 2 2 2 ...
##  $ fruto: num  1 2 3 4 5 6 1 2 3 4 ...
##  $ tam  : num  1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 ...
##  $ ic   : num  -16.7 -17.6 -19 -17.9 -17.5 ...
##  $ pi   : num  57 57 57 57 57 57 65 65 65 65 ...

levels(dat$dis) <- c("Falsa melanose", "Mancha dura")

- Plot 1 » escala fixa

(plot1 <- qplot(pi, factor(tam), data=dat, colour=ic)
    + facet_wrap(~dis,  scales = "fixed") 
    + geom_point(size = 4) + scale_shape(solid = TRUE) 
    + scale_colour_gradient2(low='darkgreen',mid='yellow3',high='yellow',
                           limits=c(-20,6), expression("Indice de cor"))
    + scale_y_discrete(name="Diâmetro de fruto (cm)")
    + scale_x_continuous(name="Avaliação (dias)")
    + coord_cartesian(xlim = c(0, 380)) 
)

plot1 + theme_bw()

- Plot 2 » escala livre

(plot2 <- qplot(pi, factor(tam), data=dat, colour=ic)
    + facet_wrap(~dis,  scales = "free") 
    + geom_point(size = 4) + scale_shape(solid = TRUE) 
    + scale_colour_gradient2(low='darkgreen',mid='yellow3',high='yellow',
                           limits=c(-20,6), expression("Indice de cor"))
    + scale_y_discrete(name="Diâmetro de fruto (cm)")
    + scale_x_continuous(name="Avaliação (dias)")
)

plot2 + theme_bw()