Skip to content

Commit

Permalink
trying setup py for installations
Browse files Browse the repository at this point in the history
  • Loading branch information
mshakya committed Oct 26, 2018
1 parent bf8bd19 commit 80dd8da
Show file tree
Hide file tree
Showing 2 changed files with 19 additions and 5 deletions.
21 changes: 17 additions & 4 deletions scripts/gage_analysis.R
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,8 @@ option_list <- list(
help = "directory containing results of EdgeR or DESeq2"),
make_option(c("-m", "--method"), action = "store",
help = "output of which results to conduct gage analysis on"),
make_option(c("-c", "--org_code"), action = "store",
help = "kegg orgainsm code to download pathway information"),
make_option(c("-o", "--out_dir"), action = "store",
help = "output directory")
)
Expand All @@ -17,6 +19,7 @@ opt <- parse_args(OptionParser(option_list = option_list))
#
dir <- opt$dir
method <- opt$method
code <- opt$org_code

#==============================================================================#
# create the output directory
Expand All @@ -25,12 +28,22 @@ method <- opt$method

# grab fold change files
if ( method == "edgeR"){
temp = list.files(path=dir, pattern="gene_*_et.csv")
fcs = list.files(path = dir, pattern = "gene_*_et.csv", full.names = TRUE)

} else if (method == "DESeq2")
{
temp = list.files(pattern="gene_*_et.csv")
print(temp)
fcs = list.files(path = dir, pattern = "gene_*_et.csv", full.names = TRUE)

}


# perform gage analysis, one fc at a time

for (fc in fcs){
fc_df <- read.csv(fc, row.names = 1)
print(head(fc_df))

}
# read in all the


gene_sets <- kegg.gsets(species = code)
3 changes: 2 additions & 1 deletion setup.py
Original file line number Diff line number Diff line change
Expand Up @@ -8,9 +8,10 @@
description="RNAseq pipeline",
author="Migun Shakya",
author_email="[email protected]",
url="https://github.com/mshakya/pypiret",
url="https://github.com/mshakya/piret",
install_requires=open("requirements.txt").read().splitlines(),
packages=find_packages(),
scripts=['bin/piret'],
license="Apache License 2.0",
platforms="Posix; MacOS X",
classifiers=[
Expand Down

0 comments on commit 80dd8da

Please sign in to comment.