Skip to content

Commit

Permalink
now actually fix
Browse files Browse the repository at this point in the history
  • Loading branch information
3mmaRand committed Oct 24, 2023
1 parent d1f8448 commit 9d5bacd
Showing 1 changed file with 11 additions and 21 deletions.
32 changes: 11 additions & 21 deletions omics/week-5/workshop.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -110,9 +110,8 @@ s30_results <- read_csv("results/S30_results.csv")
the dataframes (perhaps using `glimpse()`)

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
glimpse(s30_count_norm)
```

Expand All @@ -135,9 +134,8 @@ glimpse(s30_count_norm)
<!-- xenbase_gene_id -->

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
glimpse(s30_results)
```

Expand Down Expand Up @@ -293,9 +291,8 @@ s30_results_sig0.01 <- s30_results |>
🎬 Write the dataframe to file

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
# write to csv file
write_csv(s30_results_sig0.01,
file = "results/s30_results_sig0.01.csv")
Expand All @@ -305,9 +302,8 @@ write_csv(s30_results_sig0.01,
write to file:

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
s30_results_sig0.05 <- s30_results |>
filter(padj <= 0.05)
Expand Down Expand Up @@ -560,7 +556,7 @@ interpretation.
the names `prog` and `hspc` for the dataframes.

```{r}
#| echo: false
#| include: false
#---CODING ANSWER---
# import the normalised counts
prog <- read_csv("data-raw/surfaceome_prog.csv")
Expand All @@ -580,7 +576,7 @@ prog_hspc <- bind_cols(prog, hspc[-1])
used the name `prog_hspc_results` for the dataframe.

```{r}
#| echo: false
#| include: false
#---CODING ANSWER---
# import the DE results
prog_hspc_results <- read_csv("results/prog_hspc_results.csv")
Expand All @@ -590,9 +586,8 @@ prog_hspc_results <- read_csv("results/prog_hspc_results.csv")
the dataframes (perhaps using `glimpse()`)

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
glimpse(prog_hspc)
```

Expand All @@ -601,9 +596,8 @@ glimpse(prog_hspc)
<!-- The normalised counts dataframe has 280 rows and 1,500 columns: which are the cell id names and the gene ids. -->

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
glimpse(prog_hspc_results)
```

Expand Down Expand Up @@ -688,8 +682,7 @@ Oh, **`conflicted`** has flagged a conflict for us.


```{r}
#| echo: false
#| output: false
#| include: false
#---CODING ANSWER---
prog_hspc_results |> dplyr::select(ensembl_gene_id) |>
filter(!ensembl_gene_id %in% gene_info$ensembl_gene_id)
Expand Down Expand Up @@ -729,19 +722,16 @@ prog_hspc_results_sig0.01 <- prog_hspc_results |>
🎬 Write the dataframe to file

```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
write_csv(prog_hspc_results_sig0.01,
file = "results/prog_hspc_results_sig0.01.csv")
```

🎬 Create a dataframe of the genes significant at the 0.05 level and write to file:
```{r}
#| include: false
#---CODING ANSWER---
#| echo: false
#| output: false
prog_hspc_results_sig0.05 <- prog_hspc_results |>
filter(FDR <= 0.05)
Expand Down

0 comments on commit 9d5bacd

Please sign in to comment.