-
Notifications
You must be signed in to change notification settings - Fork 211
/
Copy pathconfig.yaml
80 lines (72 loc) · 2.95 KB
/
config.yaml
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
# path or URL to sample sheet (TSV format, columns: sample, condition, ...)
samples: config_complex/samples.tsv
# path or URL to sequencing unit sheet (TSV format, columns: sample, unit, fq1, fq2)
# Units are technical replicates (e.g. lanes, or resequencing of the same biological
# sample).
units: config_complex/units.tsv
ref:
# Ensembl species name
species: saccharomyces_cerevisiae
# Ensembl release
release: 100
# Genome build
build: R64-1-1
trimming:
# If you activate trimming by setting this to `True`, you will have to
# specify the respective cutadapt adapter trimming flag for each unit
# in the `units.tsv` file's `adapters` column
activate: False
mergeReads:
activate: False
pca:
activate: True
# Per default, a separate PCA plot is generated for each of the
# `variables_of_interest` and the `batch_effects`, coloring according to
# that variables groups.
# If you want PCA plots for further columns in the samples.tsv sheet, you
# can request them under labels as a list, for example:
# - relatively_uninteresting_variable_X
# - possible_batch_effect_Y
labels:
# columns of sample sheet to use for PCA
- jointly_handled
diffexp:
# variables where you are interested in whether they have
# an effect on expression levels
variables_of_interest:
treatment_1:
# any fold change will be relative to this factor level
base_level: untreated
treatment_2:
# any fold change will be relative to this factor level
base_level: untreated
batch_effects:
- jointly_handled
# contrasts for the deseq2 results method to determine fold changes
contrasts:
treatment_1_alone:
# must be one of the variables_of_interest
variable_of_interest: treatment_1
# the variable's level to test against the base_level
level_of_interest: treated
treatment_2_alone:
# must be one of the variables_of_interest
variable_of_interest: treatment_2
# the variable's level to test against the base_level
level_of_interest: treated
# Must be a valid expression for option two in the contrasts description
# of ?results in the DESeq2 package. For a more detailed intro, also see:
# https://github.com/tavareshugo/tutorial_DESeq2_contrasts/blob/main/DESeq2_contrasts.md
both_treatments: 'list(c("treatment_1_treated_vs_untreated", "treatment_2_treated_vs_untreated", "treatment_1treated.treatment_2treated"))'
# The default model includes all interactions among variables_of_interest,
# and batch_effects added on. For the example above this implicitly is:
# model: ~jointly_handled + treatment_1 * treatment_2
# For the default model to be used, simply specify an empty `model: ""` below.
# If you want to introduce different assumptions into your model, you can
# specify a different model to use, for example skipping the interaction:
# model: ~jointly_handled + treatment_1 + treatment_2
model: ""
params:
cutadapt-pe: ""
cutadapt-se: ""
star: ""