annotate sm_Tax4Fun.R @ 10:444f1e7fe177 draft

Uploaded
author sarahinraauzeville
date Thu, 21 Dec 2017 09:09:09 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
1 # The function -------------------------------------
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
2 #---------------------------------------------------
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
3 GalaxyFrogsTax <- function(input.file,
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
4 reference,
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
5 output.file = "Tax4FunProfile.txt"
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
6 ) {
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
7 library(Tax4Fun)
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
8 folderReferenceData<- reference
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
9
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
10 Tax4FunOutput<-Tax4Fun(importSilvaNgsData(input.file), folderReferenceData) # optionnal: fctProfiling = TRUE, refProfile = "UProc", shortReadMode = TRUE, normCopyNo = TRUE)
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
11 #folderReferenceData (required): a character vector with one character string indicating the folder location of the unzipped reference data. The reference data can be obtained from the Tax4Fun website http://tax4fun.gobics.de/ ("SILVA Reference data").
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
12 #reference = "/galaxydata/galaxy-preprod/my_tools/sm_Tax4Fun/SILVA123/"
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
13
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
14 Tax4FunProfile <- Tax4FunOutput$Tax4FunProfile
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
15 Tax4FunProfile <- data.frame(t(Tax4FunOutput$Tax4FunProfile))
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
16 #View(Tax4FunProfile)
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
17
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
18 write.table(Tax4FunProfile, "Tax4FunProfile.txt", sep="\t")
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
19
444f1e7fe177 Uploaded
sarahinraauzeville
parents:
diff changeset
20 } #fin de la fonction