Heat: Heatmap with Row Names Colored by Group

In statVisual: Statistical Visualization Tools

Description

Heatmap with row names colored by group.

Usage

Heat(data, 
     group = NULL, 
     fontsize_row=10,
     fontsize_col=10, 
     scale = "none",
     cluster_rows = TRUE,
     cluster_cols = TRUE,
     color = colorRampPalette(rev(brewer.pal(n = 7, name ="RdYlBu")))(100),
     angle_col = c("270", "0", "45", "90", "315"), 
     ...) 

Arguments

data	A data frame. Rows are subjects; Columns are variables describing the subjects. Except the column indicating subject group, all columns of data should be numeric.
group	character. The column name of data that indicates the subject groups. The row names of the heatmap will be colored based on group.
fontsize_col	x axis label font size
fontsize_row	y axis label font size
scale	character. Indicate how data will be scaled: “none” (i.e., no scaling), “row” (i.e., row scaled), “column” (i.e., column scaled).
cluster_rows	logic. Indicates if rows should be clustered.
cluster_cols	logic. Indicates if columns should be clustered.
color	vector indicating colors used in heatmap
angle_col	angle of the column labels. Please refer to the manual in pheatmap
...	other input parameters for facet & theme.

Value

A list with 10 elements: “rowInd”, “colInd”, “call”, “carpet”, “rowDendrogram”, “colDendrogram”, “breaks”, “col”, “colorTable”, “layout”.

Note

This function is based on the function pheatmap in pheatmap R package.

Author(s)

Wenfei Zhang <[e-mail]>, Weiliang Qiu <[e-mail]>, Xuan Lin <[e-mail]>, Donghui Zhang <[e-mail]>

Examples

data(esSim)
print(esSim)

# expression data
dat = exprs(esSim)
print(dim(dat))
print(dat[1:2,])

# phenotype data
pDat = pData(esSim)
print(dim(pDat))
print(pDat[1:2,])

# feature data
fDat = fData(esSim)
print(dim(fDat))
print(fDat[1:2,])

# choose the first 6 probes (3 OE probes, 2 UE probes, and 1 NE probe)
pDat$probe1 = dat[1,]
pDat$probe2 = dat[2,]
pDat$probe3 = dat[3,]
pDat$probe4 = dat[4,]
pDat$probe5 = dat[5,]
pDat$probe6 = dat[6,]

print(pDat[1:2, ])

# check histograms of probe 1 expression in cases and controls
print(table(pDat$grp, useNA = "ifany"))

statVisual(type = 'Heat', 
           data = pDat[, c(2:8)], 
           group = 'grp')

Heat(
     data = pDat[, c(2:8)], 
     group = 'grp')

statVisual documentation built on Feb. 21, 2020, 1:08 a.m.