Skip to content
Snippets Groups Projects
Commit 702cbfbd authored by Alexis Mergez's avatar Alexis Mergez
Browse files

Added total length to GFAstats

parent f19a6e04
No related branches found
Tags v1.6.1
No related merge requests found
Pipeline #182811 passed
...@@ -164,6 +164,7 @@ for pathName in pathStats.keys(): ...@@ -164,6 +164,7 @@ for pathName in pathStats.keys():
if args.progress: if args.progress:
print(f"[GFAstats::{panname}] Computing general stats") print(f"[GFAstats::{panname}] Computing general stats")
# Mean length of paths # Mean length of paths
genStats["Path.length.mean"] = round(np.mean( genStats["Path.length.mean"] = round(np.mean(
[pathStats[pathname]["Path.length"] for pathname in pathStats.keys()] [pathStats[pathname]["Path.length"] for pathname in pathStats.keys()]
...@@ -184,6 +185,11 @@ genStats["Steps.count"] = np.sum( ...@@ -184,6 +185,11 @@ genStats["Steps.count"] = np.sum(
[pathStats[pathname]["Path.steps.count"] for pathname in pathStats.keys()] [pathStats[pathname]["Path.steps.count"] for pathname in pathStats.keys()]
) )
# Total length
genStats["Total.length"] = np.sum(
list(nodesLengthDict.values())
)
# Mean length of nodes # Mean length of nodes
genStats["Nodes.length.mean"] = round(np.mean( genStats["Nodes.length.mean"] = round(np.mean(
list(nodesLengthDict.values()) list(nodesLengthDict.values())
...@@ -201,7 +207,7 @@ genStats["Degree.mean"] = round(genStats["Edges.count"] / genStats["Nodes.count" ...@@ -201,7 +207,7 @@ genStats["Degree.mean"] = round(genStats["Edges.count"] / genStats["Nodes.count"
print(f"[GFAstats::{panname}] Statistics computed in {round(time.time() - start_time, 2)} s") print(f"[GFAstats::{panname}] Statistics computed in {round(time.time() - start_time, 2)} s")
## Exporting general data to tsv ## Exporting general data to tsv
# Tranforming keys to array like keys # Transforming keys to array-like keys
genStatsArr = {key: {panname : value} for key, value in genStats.items()} genStatsArr = {key: {panname : value} for key, value in genStats.items()}
# Converting to pandas dataframe # Converting to pandas dataframe
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment