Tfe3 and hyperLOPIT plot

hyperlopit-tfe.R

## load packages and data
library("pRoloc")
library("pRolocdata")
data(hyperLOPIT2015)

## produce PCA plot
plot2D(hyperLOPIT2015)
addLegend(hyperLOPIT2015, where = "bottomleft", cex = .7)

## Find Tfe3's identifier and define it as a
## feature of interest that can be overlaid on the PCA plot
id <- featureNames(hyperLOPIT2015)[grep("TFE3", fData(hyperLOPIT2015)[, 1])]
tfe3 <- FeaturesOfInterest(id, description = "Tfe3")
highlightOnPlot(hyperLOPIT2015, tfe3, pch = 19, cex = 2)
highlightOnPlot(hyperLOPIT2015, tfe3, pch = 19, cex = 1, labels = TRUE, pos = 3)