Taxa de crscimento micelial in vitro de isolados…

O objetivo é testar se sao significativamente diferentes entre eles, e para aqueles que sao semelhantes fazer uma unica reta de regressão.

Session set

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/ISA/tese")# lab
#setwd("/media/juanchi/DATA/Dropbox/MyR/Análise otros/ISA/tese") # DELL

Dataset

# Dados planilha laboratorio
datlab = read.csv("micel_vitro.csv", dec=",", header=T, sep="\t", check.names=FALSE)  
head(datlab)
##       iso exp rep dia diam_col
## 1 SP05160   1   1   1     0.90
## 2 SP05160   1   2   1     0.75
## 3 SP05160   1   3   1     0.75
## 4 SP05160   1   4   1     0.80
## 5 SP05160   1   5   1     0.55
## 6 SP05160   1   6   1     0.55

Grafico exploratorio

explot = ggplot(datlab, aes(dia, diam_col)) + 
  geom_point(shape=20, size=2) + 
  geom_smooth(method="lm", color="red")+  
  facet_wrap(iso~exp)  
  #scale_y_continuous("Germinação (%)") 
explot