Skip to content

Commit

Permalink
Update command line args for dlcpar
Browse files Browse the repository at this point in the history
  • Loading branch information
davidemms committed Oct 3, 2016
1 parent d861544 commit bb048b4
Show file tree
Hide file tree
Showing 4 changed files with 8 additions and 10 deletions.
2 changes: 1 addition & 1 deletion Tests/test_orthofinder.py
Original file line number Diff line number Diff line change
Expand Up @@ -34,7 +34,7 @@
goldResultsDir_smallExample = baseDir + "ExpectedOutput/SmallExampleDataset/"
goldPrepareBlastDir = baseDir + "ExpectedOutput/SmallExampleDataset_PreparedForBlast/"

version = "1.0.6"
version = "1.0.7"
requiredBlastVersion = "2.2.28+"

citation = """When publishing work that uses OrthoFinder please cite:
Expand Down
Binary file removed orthofinder/Orthogroups_Orthologues.png
Binary file not shown.
14 changes: 6 additions & 8 deletions orthofinder/scripts/get_orthologues.py
Original file line number Diff line number Diff line change
Expand Up @@ -448,16 +448,14 @@ def RunDlcpar(treesPat, ogSet, nOGs, speciesTreeFN, workingDir):
- run
"""
rootedTreeDir = workingDir + "Trees_ids_arbitraryRoot/"
if not os.path.exists(rootedTreeDir): os.mkdir(rootedTreeDir)
RootGeneTreesArbitrarily(treesPat, nOGs, rootedTreeDir)
geneMapFN = WriteGeneSpeciesMap(rootedTreeDir, ogSet)

dlcparResultsDir = workingDir + 'dlcpar/'
if not os.path.exists(dlcparResultsDir): os.mkdir(dlcparResultsDir)
filenames = [rootedTreeDir + os.path.split(treesPat % i)[1] for i in xrange(nOGs)]
RootGeneTreesArbitrarily(treesPat, nOGs, dlcparResultsDir)
geneMapFN = WriteGeneSpeciesMap(dlcparResultsDir, ogSet)

filenames = [dlcparResultsDir + os.path.split(treesPat % i)[1] for i in xrange(nOGs)]

dlcCommands = ['dlcpar_search -s %s -S %s -D 1 -C 0.125 %s -O %s' % (speciesTreeFN, geneMapFN, fn, dlcparResultsDir + os.path.splitext(os.path.split(fn)[1])[0]) for fn in filenames]
dlcCommands = ['dlcpar_search -s %s -S %s -D 1 -C 0.125 %s -I .txt' % (speciesTreeFN, geneMapFN, fn) for fn in filenames]
# print(dlcCommands[0])
# use this to run in parallel
util.RunParallelOrderedCommandLists(nThreads, [[c] for c in dlcCommands], qHideStdout = True)
Expand Down Expand Up @@ -531,7 +529,7 @@ def GetResultsFilesString(rootedSpeciesTreeFN):

def CleanWorkingDir(dendroBlast):
dendroBlast.DeleteBlastMatrices()
dirs = ['Distances/', "matrices_orthologues/", "Trees_ids_arbitraryRoot/"]
dirs = ['Distances/', "matrices_orthologues/"]
for d in dirs:
dFull = dendroBlast.workingDir + d
if os.path.exists(dFull):
Expand Down
2 changes: 1 addition & 1 deletion orthofinder/scripts/util.py
Original file line number Diff line number Diff line change
Expand Up @@ -47,7 +47,7 @@
FileInfo = namedtuple("FileInfo", "inputDir outputDir graphFilename")

picProtocol = 1
version = "1.0.6"
version = "1.0.7"

def PrintNoNewLine(text):
sys.stdout.write(text)
Expand Down

0 comments on commit bb048b4

Please sign in to comment.