From af83a4f17ffa19f8ac1935cbcdc8e8e2d46b858b Mon Sep 17 00:00:00 2001
From: Sajjad Ghaffarinasabsharabiani <ghaffas@hhu.de>
Date: Tue, 9 Aug 2022 12:06:26 +0000
Subject: [PATCH] Upload New File
---
.../organismInformation/taxInfo.R | 35 +++++++++++++++++++
1 file changed, 35 insertions(+)
create mode 100644 networkComplexityBigg/organismInformation/taxInfo.R
diff --git a/networkComplexityBigg/organismInformation/taxInfo.R b/networkComplexityBigg/organismInformation/taxInfo.R
new file mode 100644
index 0000000..891c42c
--- /dev/null
+++ b/networkComplexityBigg/organismInformation/taxInfo.R
@@ -0,0 +1,35 @@
+#!/usr/bin/Rscript
+library(methods)
+library(dplyr)
+library(XML)
+
+modelTable <- read.csv("../envirDist/modelTable.csv")
+
+
+xmlURL <- paste0("https://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=taxonomy&retmode=xml&id=",
+ paste(modelTable$taxid, collapse=","))
+xmlText <- paste(scan(xmlURL, what="",sep="\n"),"\n",collapse="\n")
+doc <- xmlParse(xmlText)
+
+
+getNodeSet(doc, "/TaxaSet/Taxon/Lineage")
+
+
+modelTableTaxonInformation <- data.frame(bigg_id=modelTable$bigg_id,
+taxid=xpathSApply(doc, "/TaxaSet/Taxon/TaxId", xmlValue),
+name=xpathSApply(doc, "/TaxaSet/Taxon/ScientificName", xmlValue),
+division=xpathSApply(doc, "/TaxaSet/Taxon/Division", xmlValue),
+lineage=xpathSApply(doc, "/TaxaSet/Taxon/Lineage", xmlValue)
+)
+
+
+
+
+
+
+
+
+
+
+
+
--
GitLab