-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy path10-heatmaps.Rmd
More file actions
107 lines (73 loc) · 3.33 KB
/
10-heatmaps.Rmd
File metadata and controls
107 lines (73 loc) · 3.33 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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
# Heatmaps
Los heatmaps son una técnica de visualización de datos que mide la magnitud de un fenómeno en colores en dos dimensiones. La variación del color puede ser por tono o intensidad, haciendo obvia la lectura del fenómeno sobre el espacio que se trata.
Son el resultado obtenido al representar una matriz de valores en la que, en lugar de números, se muestra un gradiente de color proporcional al valor de cada variable en cada posición.
Los heatmaps son ampliamente utilizados para visualizar e interpretar datos de expresión génica, principalmente para experimentos con microarrays y RNA-seq. También se puede combinar con métodos de agrupación, que agrupan genes y/o muestras en función de la similitud de su patrón de expresión génica. Esto puede ser útil para identificar genes que están comúnmente regulados o firmas biológicas asociadas con una condición particular (por ejemplo, una enfermedad o una condición ambiental).
Para crear los heatmps usaremos `pheatmap()`. Instalamos la paquetería y después la cargamos.
```{r, eval=FALSE}
install.packages("pheatmap")
```
```{r}
library(pheatmap)
```
El set de datos que vamos a usar los puedes descargar de [aquí](https://drive.google.com/u/0/uc?id=1aS5-Zf03OeMRtZsCCLYFXAUQpesGy4IK&export=download). Recuerda guardarlos en tu carpeta de trabajo actual, como lo vimos en la sección 7.
```{r, eval=F}
data <- read.csv(file = "Datos.csv")
```
```{r, eval=T, echo=F}
data <- read.csv("~/Downloads/Datos.csv")
```
Vamos a restructurar los datos
```{r, eval=T}
rownames(data) <- data[,1]
samp2 <- data[,-1]
```
Transformaremos nuestros datos a una matriz para poder crear nuestro heatmap.
```{r, eval=T}
mat_data <- data.matrix(samp2[,1:ncol(samp2)])
mat_data
```
Construyendo un heatmap básico
```{r, eval=T, fig.show='hold'}
pheatmap(mat_data)
```
Con los comandos `fontsize_col`(columnas) y `fontsize_row` (filas), puedes cambiar el tamaño de letra y hacerlo más grande o más pequeño.
```{r, eval=T, fig.show='hold'}
pheatmap(mat_data,
fontsize_col=6,
fontsize_row=6)
```
`main` te permite agregar el título al Heatmap. Recuerda que el título que eligas debe estar entre comillas
```{r, eval=T, fig.show='hold'}
pheatmap(mat_data,
fontsize_col=10,
fontsize_row=10,
main = "Mi primer heatmap")
```
La escala de colores fríos generalmente se usa en los heatmaps. `color` te permite cambiar los colores de la gráfica.
[Aquí](http://www.stat.columbia.edu/~tzheng/files/Rcolor.pdf) y [aquí](https://htmlcolorcodes.com/) puedes revisar los colores disponibles en R que puedes usar en tus gráficas.
```{r, eval=T, fig.show='hold'}
pheatmap(mat_data,
fontsize_col=10,
fontsize_row=10,
main = "Mi primer Heatmap",
color = c("blue", "yellow","red"))
```
Agregar o cambiar las líneas de los clusters con `cluster_cols` y `cluster_rows`
```{r, eval=T, fig.show='hold'}
pheatmap(mat_data,
fontsize_col=10,
fontsize_row=10,
main = "Mi primer heatmap",
color = c("blue", "yellow","red"),
cluster_cols=F,
cluster_rows=T)
```
```{r, eval=T, fig.show='hold'}
pheatmap(mat_data,
fontsize_col=10,
fontsize_row=10,
main = "Mi primer heatmap",
color = c("blue", "yellow","red"),
cluster_cols=F,
cluster_rows=F)
```