covid_comp_script0_table.Rnw

@@ -56,6 +56,43 @@ names(dginnT)<-c("File", "Name", "Gene.name", "GeneSize",
   "dginn-primate_codemlM7M8.PSS")
 @
 
+Add SELENOS
+
+<<selenos>>=
+selenos<-read.delim(paste0(workdir,
+      "data/resSELENOS.tab"))
+
+# liste of colonne
+
+colonnes<-c("File", "Name", "Gene", "GeneSize",
+  "NbSpecies", "omegaM0Bpp", "omegaM0codeml", "BUSTED",
+  "BUSTED_p.value",  "MEME_NbSites", "MEME_PSS", "BppM1M2",
+  "BppM1M2_p.value", "BppM1M2_NbSites", "BppM1M2_PSS", "BppM7M8", 
+  "BppM7M8_p.value", "BppM7M8_NbSites", "BppM7M8_PSS","codemlM1M2",
+  "codemlM1M2_p.value", "codemlM1M2_NbSites",  "codemlM1M2_PSS", "codemlM7M8",   
+  "codemlM7M8_p.value", "codemlM7M8_NbSites", "codemlM7M8_PSS")         
+
+selenos<-selenos[,colonnes]
+
+
+
+@
+
+<<>>=
+names(selenos)<-names(dginnT)
+selenos[,6]<-as.factor(selenos[,6])
+selenos[,9]<-as.factor(selenos[,9])
+selenos[,11]<-as.factor(selenos[,11])
+
+selenos[,13]<-as.factor(selenos[,13])
+selenos[,17]<-as.factor(selenos[,17])
+selenos[,21]<-as.factor(selenos[,21])
+selenos[,25]<-as.factor(selenos[,25])
+
+## convertir les pvalues
+dginnT<-rbind(dginnT, selenos)
+@
+    
 \subsection{Bats}
 
 <<>>=