-
Notifications
You must be signed in to change notification settings - Fork 1
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
1 changed file
with
11 additions
and
1 deletion.
There are no files selected for viewing
12 changes: 11 additions & 1 deletion
12
Data Processing/Psychencode_June2019/Single Label-based exp/ReadMe.md
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -1,4 +1,14 @@ | ||
# Single-Sample Data | ||
## Peak filtering | ||
We selected rows that have at least 50% samples have a signal "1" (defined as peak height > 0) using [```peak_filtering_single_Sample.R```](https://github.com/Akmazad/deepPsych/blob/master/Data%20Processing/Psychencode_June2019/Single%20Label-based%20exp/peak_filtering_single_Sample.R)) script. It outputs filtered data in BED format along with binarized peak value at threshold 0. | ||
|
||
|Name|nSample (only Controls)|nPeaks (total)| nPeaks (>50% samples have non-zero peak-height) | | ||
|---|---|---|---| | ||
|HumanFC|137 (out of 288)|197,263|15,867| | ||
|
||
## Bin overlapping | ||
Next, we need to intersect each peaks with chromosomal bins of fixed-width. Following commands were used to get it done. | ||
|
||
```sh | ||
intersectBed -wao -f 0.05 -a hg19_bins_200bp.bed -b mergedPeakHeightMatrix_HumanFC_filtered_signle_label.BED > mergedPeakHeightMatrix_HumanFC_filtered_signle_label.overlaps.bed | ||
``` | ||
|