-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathData_exercise_template.Rmd
More file actions
87 lines (49 loc) · 1.9 KB
/
Data_exercise_template.Rmd
File metadata and controls
87 lines (49 loc) · 1.9 KB
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
81
82
83
84
85
86
87
---
title: "conception_intro_data"
output: html_document
date: '2023-08-04'
---
# Yoshino 2021 Mouse OSC Data Exercise
**Name:** Add your name here
## Instructions
Please clone this repo, create a branch for yourself, and then explore the scRNAseq data from [this](https://www.science.org/doi/epdf/10.1126/science.abe0237) paper.
<br>
Data is provided in the data directory as an RDS for a monocle3 object holding unprocessed matrices + metadata. Note: to reduce the RDS to within github size limits, we've filtered some non/lowly expressed genes and downsampled the data to 65% of the original cellcount per sample.
<br>
Please perform the following tasks and present your work product as a knitted RMD (knit to html). <br>
- Preprocess and cluster the data - annotate decisions in the RMD, we'd like to see your process.
- Roughly label celltypes -- this should include at least granulosa, soma, and germ cells. Can you get any more specific? Are there other additional celltypes? Show plots supporting your celltype labels.
- How do celltypes compare across ages or across in vivo/in vitro samples? (open ended)
<br>
Submit by pushing to your code & html files to your branch.
<br>
We suggest to not spend more than ~3hrs total on this. Feel free to explore beyond the starter questions here. Find anything particularly interesting? Let's discuss!
<br>
# Load Data and set up the RMD
```{r setup, include=FALSE}
library(knitr)
opts_chunk$set(echo=FALSE, message=FALSE, warning=FALSE)
```
```{r load_packages, include=FALSE}
library(monocle3)
library(ggplot2)
```
```{r load_data, include=FALSE}
cds<-readRDS("data/Yoshino_2021_unprocessed_downsampled.RDS")
```
# QC checks & Filtering
```{r QC_and_filter}
```
# Preprocess & Cluster
```{r preprocess}
```
# Celltype Identification
```{r celltypeID}
```
# Anything else of interest?
```{r play}
```
## SessionInfo
```{r sessionInfo}
sessionInfo()
```