Skip to content

use peppy #22

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
stolarczyk wants to merge 5 commits into from
Closed

use peppy #22

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
750e481
use peppy
stolarczyk Mar 2, 2020
4ebd07d
add PEP config
stolarczyk Mar 2, 2020
e893958
reformat cfg
stolarczyk Mar 2, 2020
25459cb
correct config file path
stolarczyk Mar 3, 2020
adb285c
add pep_version to PEP config
stolarczyk May 21, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions pep_config.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,4 @@
pep_version: 2.0.0
sample_table: samples.tsv
subsample_table: units.tsv
snake_config: "config.yaml"
11 changes: 7 additions & 4 deletions rules/common.smk
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,20 +1,23 @@
import pandas as pd
from peppy import Project
from snakemake.utils import validate
from snakemake.utils import min_version

min_version("5.7.1")

report: "../report/workflow.rst"

p = Project(cfg="pep_config.yaml", sample_table_index="sample",
subsample_table_index=["sample", "unit"])

###### Config file and sample sheets #####
configfile: "config.yaml"
configfile: p.config.snake_config
validate(config, schema="../schemas/config.schema.yaml")

samples = pd.read_table(config["samples"]).set_index("sample", drop=False)
samples = p.sample_table
validate(samples, schema="../schemas/samples.schema.yaml")

units = pd.read_table(config["units"], dtype=str).set_index(["sample", "unit"], drop=False)
units.index = units.index.set_levels([i.astype(str) for i in units.index.levels]) # enforce str in index
units = p.subsample_table
validate(units, schema="../schemas/units.schema.yaml")


Expand Down