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.
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
# 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
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