Skip to content

Commit d8d636e

Browse files
authored
Slides and sample code posted
1 parent ff22a3f commit d8d636e

File tree

2 files changed

+49
-0
lines changed

2 files changed

+49
-0
lines changed

Lab12/HW6_lab12_sample_code.Rmd

+49
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,49 @@
1+
---
2+
title: "Lab12_sample_code"
3+
author: "Qian Xiao"
4+
date: '2020-04-15'
5+
output: html_document
6+
---
7+
8+
#PartIV
9+
10+
##Sample Slurm Script
11+
```{r, engine='bash', eval=FALSE}
12+
#!/bin/bash
13+
#SBATCH -n 1 # Number of cores
14+
#SBATCH -N 1 # Ensure that all cores are on one machine
15+
#SBATCH -t 0-06:00 # Runtime in D-HH:MM
16+
#SBATCH -p serial_requeue # Partition to submit to
17+
#SBATCH --mem=1000 # Memory pool for all cores (see also --mem-per-cpu)
18+
#SBATCH -o mageck.out # File to which STDOUT will be written
19+
#SBATCH -e mageck.err # File to which STDERR will be written
20+
#SBATCH --mail-type=ALL
21+
#SBATCH --mail-user="[email protected]"
22+
23+
module load Anaconda/5.0.1-fasrc01
24+
25+
# your code here
26+
27+
```
28+
29+
##Pooled replicates
30+
```{r, engine='bash', eval=FALSE}
31+
#Generate read counts
32+
mageck count -l library.csv -n OUT --sample-label Day0,Day23 \
33+
--fastq Day0_Rep1.fastq.gz,Day0_Rep2.fastq.gz Day23_Rep1.fastq.gz,Day23_Rep2.fastq.gz
34+
35+
#Generate summary statistics
36+
mageck test -k OUT.count.txt -t Day23 -c Day0 -n OUT
37+
```
38+
39+
##Consider replicates separately
40+
```{r, engine='bash', eval=FALSE}
41+
42+
#Get the read counts
43+
44+
mageck count -l library.csv -n OUT_SEPARATE --sample-label Day0_Rep1,Day0_Rep2,Day23_Rep1,Day23_Rep2 \
45+
--fastq Day0_Rep1.fastq.gz Day0_Rep2.fastq.gz Day23_Rep1.fastq.gz Day23_Rep2.fastq.gz
46+
```
47+
48+
49+

Lab12/Lab12.pdf

3.65 MB
Binary file not shown.

0 commit comments

Comments
 (0)