Skip to content

Commit

Permalink
Merge pull request #14 from Christina-Straham/main
Browse files Browse the repository at this point in the history
I only changed the font size for the heat map and added some comments!
  • Loading branch information
YanqingLiu824 authored Nov 29, 2021
2 parents 956dafd + da94736 commit 02c5fc9
Showing 1 changed file with 7 additions and 1 deletion.
8 changes: 7 additions & 1 deletion vignettes/TET2.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ knitr::opts_chunk$set(echo = TRUE)
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")
library(devtools)
load_all("./")
# not sure what this is for but it is important. When this chunk is ran, a pretty Rmd is pulled up of the code.
```

# Installation
Expand All @@ -34,6 +35,8 @@ To extract just the R code, you can use knitr::knit(input, tangle=TRUE):
```{r, tangle, eval = FALSE, message = FALSE, echo = FALSE}
# knitr::knit("TET2.Rmd", tangle = TRUE)
# [1] "TET2.R"
#What is the point of this chunk?
```

# Introduction
Expand Down Expand Up @@ -93,7 +96,7 @@ Is it the case that TET2, IDH1, and IDH2 mutations are exclusive?
library(ComplexHeatmap)
mutations <- t(as.matrix(pData(DNAme)[, c("TET2", "IDH")]))
Heatmap(mutations, col=c("lightgray","darkred"), name="mutant", column_km=4,
column_names_gp = gpar(fontsize = 7))
column_names_gp = gpar(fontsize = 0))
```

Expand All @@ -112,6 +115,8 @@ fit1 <- eBayes(lmFit(exprs(DNAme), design1))
coef=grep("IDH", colnames(design1)),
p.value=0.05, # change if you like
number=Inf)))
#if you're having trouble running this entire chunk or receiving error "'pData': object 'DNAme' not found", run design1 chunk of code first and the rest of the chunk should follow suit. I don't know why that is or if it is just me.
# 6513 probes for IDH
(TET_diffmeth_probes_fit1 <- nrow(topTable(fit1,
Expand All @@ -127,6 +132,7 @@ fit2 <- eBayes(lmFit(exprs(DNAme), design2))
coef=grep("IDH", colnames(design2)),
p.value=0.05, # change if you like
number=Inf)))
# 6651 probes for IDH
(TET2_diffmeth_probes_fit2 <- nrow(topTable(fit2,
Expand Down

0 comments on commit 02c5fc9

Please sign in to comment.