Skip to content

Commit

Permalink
Merge pull request #15 from 3mmaRand/fix-omics-2
Browse files Browse the repository at this point in the history
Fix omics 2
  • Loading branch information
3mmaRand authored Oct 19, 2023
2 parents 4a7dd22 + 0bc3400 commit 74ff924
Show file tree
Hide file tree
Showing 29 changed files with 123,949 additions and 6 deletions.
2 changes: 1 addition & 1 deletion omics/week-3/workshop.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -63,7 +63,7 @@ These are the three datasets. Each set compromises several files.
- [xlaevis_counts_S14.csv](data-raw/xlaevis_counts_S14.csv)
- [xlaevis_counts_S20.csv](data-raw/xlaevis_counts_S20.csv)
- [xlaevis_counts_S30.csv](data-raw/xlaevis_counts_S30.csv)
- [meta_data.txt](data-raw/meta_data.txt)


🐭 Stem cell data:

Expand Down
15 changes: 10 additions & 5 deletions omics/week-4/workshop.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -523,6 +523,8 @@ prog_hspc |>
```

Notice that we have summed across all the columns.

❓ What do you conclude?

<!-- #---THINKING ANSWER--- -->
Expand Down Expand Up @@ -552,7 +554,7 @@ rowSums(prog_hspc) |> sort() |> head(10)
To find the genes that are expressed in only one cell type, we can use
the same approach as above but only sum the columns for one cell type.

🎬 Find the genes that are 0 in every column for the Prog cells:
🎬 Find the genes that are 0 in every column for the HSPC cells:

```{r}
prog_hspc |>
Expand All @@ -561,18 +563,21 @@ prog_hspc |>
```

We have summed across the HSPC cells only.
Note that if we knew there were some rows that were all zero across both
cerll types, we would need to add
cell types, we would need to add
`|> filter(sum(c_across(Prog_001:Prog_852)) != 0)`

🎬 Now you find the genes that are 0 in every column for the HSPC cells:
meaning zero in all the HSPC but not zero in all the Prog

🎬 Now you find the genes that are 0 in every column for the Prog cells:

```{r}
#| echo: false
#---CODING ANSWER---
prog_hspc |>
rowwise() |>
filter(sum(c_across(HSPC_001:HSPC_852)) == 0)
filter(sum(c_across(Prog_001:Prog_852)) == 0)
```

Expand Down Expand Up @@ -694,7 +699,7 @@ practical. Qmd stands for Quarto markdown. It allows code and ordinary
text to be interleaved to produce well-formatted reports including
webpages. Right-click on the link and choose Save-As to download. You
will be able to open the Qmd file in RStudio. Alternatively, [View in
Browser](https://github.com/3mmaRand/). Coding and thinking answers are
Browser](https://github.com/3mmaRand/BIO00088H-data/blob/main/omics/week-4/workshop.qmd). Coding and thinking answers are
marked with `#---CODING ANSWER---` and `#---THINKING ANSWER---`

Pages made with R [@R-core], Quarto [@allaire2022], `knitr` [@knitr],
Expand Down
281 changes: 281 additions & 0 deletions omics/week-5/data-processed/hspc_summary_gene.csv

Large diffs are not rendered by default.

702 changes: 702 additions & 0 deletions omics/week-5/data-processed/hspc_summary_samp.csv

Large diffs are not rendered by default.

281 changes: 281 additions & 0 deletions omics/week-5/data-processed/prog_summary_gene.csv

Large diffs are not rendered by default.

799 changes: 799 additions & 0 deletions omics/week-5/data-processed/prog_summary_samp.csv

Large diffs are not rendered by default.

Loading

0 comments on commit 74ff924

Please sign in to comment.