Fig.6A.EPICON.DOC.2019.04.06.Rmd

---
title: "Untitled"
author: "Cheng Gao"
date: "4/6/2019"
output: html_document
---

```{r, message=FALSE, warning=FALSE, fig.height = 3, fig.width = 4}
wd<- "/Users/chengg/Google Drive/EPICON/Mycobiome/Fungal ITS/statistic/Total.fungi"
setwd(wd)
rm(list = ls())
library(DOC)
load("EPICON.DOC.2018.09.22.Rdata")
da<-doc.endophyte
plot(da);
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)


da<-doc.pathogen
plot(da);
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)


da<-doc.yeast
plot(da);
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.saprotroph
plot(da);
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

load("EPICON.DOC.2018.06.19.Rdata")
da<-doc.fung
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

load("EPICON.DOC.2018.05.24.Rdata")

da<-doc.amf
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.fung.L1.7
plot(da)

median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.fung.L8.17
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.fung.L
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.fung.R
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.fung.Z
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)

da<-doc.fung.S
plot(da)
median(da$FNS$Fns)
(sum(da$LME$Slope >= 0) + 1) / (length(da$LME$Slope) + 1)
median(da$NEG$Neg.Slope)
```