- Loading ggfortify
- Plotting matrix
- Plotting diagnostics for LM and GLM
- Plotting time series
- Plotting PCA (Principal Component Analysis)
- Plotting K-means
- Plotting cluster package
- Plotting Local Fisher Discriminant Analysis
- Plotting survival curves
- Learn more
ggfortify extends ggplot2 for plotting some popular R packages using a standardized approach, included in the function autoplot().
This article describes how to draw:
- a matrix,
- a scatter plot,
- diagnostic plots for linear model,
- time series,
- the results of principal component analysis,
- the results of clustering analysis,
- and survival curves
The following R packages and functions are covered in the package ggfortify:
|base||matrix and table|
|cluster||clara, fanny and pam|
|dlm||dlmFilter and dlmSmooth|
|forecast||bats, forecast, ets and nnetar|
|KFAS||KFS and signal|
|lfda||klfda and self|
|MASS||isoMDS and sammon|
|stats||acf, ar, Arima, smdscale, decomposed.ts, density, fractanal, glm, HoltWinters, kmeans, lm, prcomp, princomp, spec, stepfun, stl and ts|
|survival||survfit and survfit.cox|
|strucchange||breakpoints and breakpointsfull|
ggfortify can be installed from GitHub or CRAN:
# Github if(!require(devtools)) install.packages("devtools") devtools::install_github("sinhrks/ggfortify")
# CRAN install.packages("ggfortify")
The function autoplot.matrix() is used:
autoplot(object, geom = "tile")
- object: an object of class matrix
- geom: allowed values are “tile” (for heatmap) or “point” (for scatter plot)
The mtcars data set is used in the example below.
df <- mtcars[, c("mpg", "disp", "hp", "drat", "wt")] df <- as.matrix(df)
Plot a heatmap:
# Heatmap autoplot(scale(df))
Plot a scatter plot: The data should be a matrix with 2 columns named V1 and V2. The R code below plots mpg by wt. We start by renaming column names.
# Extract the data df2 <- df[, c("wt", "mpg")] colnames(df2) <- c("V1", "V2") # Scatter plot autoplot(df2, geom = 'point') + labs(x = "mpg", y = "wt")
Plotting diagnostics for LM and GLM
The function autoplot.lm() is used to draw diagnostic plots for LM and GLM [in stats package].
autoplot(object, which = c(1:3, 5))
- object: stats::lm instance
- which: If a subset of the plots is required, specify a subset of the numbers 1:6.
- ncol and nrow allows you to specify the number of subplot columns and rows.
Diagnostic plots for Linear Models (LM)
iris data set is used for computing the linear model
# Compute a linear model m <- lm(Petal.Width ~ Petal.Length, data = iris) # Create the plot autoplot(m, which = 1:6, ncol = 2, label.size = 3)
# Change the color by groups (species) autoplot(m, which = 1:6, label.size = 3, data = iris, colour = 'Species')
Diagnostic plots with Generalized Linear Models (GLM)
USArrests data set is used.
# Compute a generalized linear model m <- glm(Murder ~ Assault + UrbanPop + Rape, family = gaussian, data = USArrests) # Create the plot # Change the theme and colour autoplot(m, which = 1:6, ncol = 2, label.size = 3, colour = "steelblue") + theme_bw()
Plotting time series
Plotting ts objects
- Data set: AirPassengers
- R Function: autoplot.ts()
The function autoplot() can handle also other time-series-likes packages, including:
Plotting with changepoint package
The changepoint package provides a simple approach for identifying shifts in mean and/or variance in a time series.
ggfortify supports cpt object in changepoint package.
Plotting with strucchange package
strucchange is an R package for detecting jumps in data.
Data set: Nile
library(strucchange) autoplot(breakpoints(Nile ~ 1))
Plotting PCA (Principal Component Analysis)
- Data set: iris
- Function: autoplot.prcomp()
# Prepare the data df <- iris[, -5] # Principal component analysis pca <- prcomp(df, scale. = TRUE) # Plot autoplot(pca, loadings = TRUE, loadings.label = TRUE, data = iris, colour = 'Species')
- Data set: USArrests
- Function: autoplot.kmeans()
The original data is required as kmeans object doesn’t store original data. Samples will be colored by groups (clusters).
autoplot(kmeans(USArrests, 3), data = USArrests, label = TRUE, label.size = 3, frame = TRUE)
Plotting cluster package
ggfortify supports cluster::clara, cluster::fanny and cluster::pam classes. These functions return object containing original data, so there is no need to pass original data explicitly.
The R code below shows an example for pam() function:
library(cluster) autoplot(pam(iris[-5], 3), frame = TRUE, frame.type = 'norm')
Plotting Local Fisher Discriminant Analysis
library(lfda) # Local Fisher Discriminant Analysis (LFDA) model <- lfda(iris[,-5], iris[, 5], 4, metric="plain") autoplot(model, data = iris, frame = TRUE, frame.colour = 'Species')
Plotting survival curves
library(survival) fit <- survfit(Surv(time, status) ~ sex, data = lung) autoplot(fit)
Read more on ggfortify.
This analysis has been performed using R software (ver. 3.2.1) and ggplot2 (ver. 1.0.1)
Show me some love with the like buttons below... Thank you and please don't forget to share and comment below!!
Montrez-moi un peu d'amour avec les like ci-dessous ... Merci et n'oubliez pas, s'il vous plaît, de partager et de commenter ci-dessous!
Recommended for You!
Recommended for you
This section contains best data science and self-development resources to help you on your path.
Coursera - Online Courses and Specialization
- Course: Machine Learning: Master the Fundamentals by Standford
- Specialization: Data Science by Johns Hopkins University
- Specialization: Python for Everybody by University of Michigan
- Courses: Build Skills for a Top Job in any Industry by Coursera
- Specialization: Master Machine Learning Fundamentals by University of Washington
- Specialization: Statistics with R by Duke University
- Specialization: Software Development in R by Johns Hopkins University
- Specialization: Genomic Data Science by Johns Hopkins University
Popular Courses Launched in 2020
- Google IT Automation with Python by Google
- AI for Medicine by deeplearning.ai
- Epidemiology in Public Health Practice by Johns Hopkins University
- AWS Fundamentals by Amazon Web Services
- The Science of Well-Being by Yale University
- Google IT Support Professional by Google
- Python for Everybody by University of Michigan
- IBM Data Science Professional Certificate by IBM
- Business Foundations by University of Pennsylvania
- Introduction to Psychology by Yale University
- Excel Skills for Business by Macquarie University
- Psychological First Aid by Johns Hopkins University
- Graphic Design by Cal Arts
Books - Data Science
- Practical Guide to Cluster Analysis in R by A. Kassambara (Datanovia)
- Practical Guide To Principal Component Methods in R by A. Kassambara (Datanovia)
- Machine Learning Essentials: Practical Guide in R by A. Kassambara (Datanovia)
- R Graphics Essentials for Great Data Visualization by A. Kassambara (Datanovia)
- GGPlot2 Essentials for Great Data Visualization in R by A. Kassambara (Datanovia)
- Network Analysis and Visualization in R by A. Kassambara (Datanovia)
- Practical Statistics in R for Comparing Groups: Numerical Variables by A. Kassambara (Datanovia)
- Inter-Rater Reliability Essentials: Practical Guide in R by A. Kassambara (Datanovia)
- R for Data Science: Import, Tidy, Transform, Visualize, and Model Data by Hadley Wickham & Garrett Grolemund
- Hands-On Machine Learning with Scikit-Learn, Keras, and TensorFlow: Concepts, Tools, and Techniques to Build Intelligent Systems by Aurelien Géron
- Practical Statistics for Data Scientists: 50 Essential Concepts by Peter Bruce & Andrew Bruce
- Hands-On Programming with R: Write Your Own Functions And Simulations by Garrett Grolemund & Hadley Wickham
- An Introduction to Statistical Learning: with Applications in R by Gareth James et al.
- Deep Learning with R by François Chollet & J.J. Allaire
- Deep Learning with Python by François Chollet
Want to Learn More on R Programming and Data Science?
Follow us by Email On Social Networks: