#### Overview of this book

This cookbook offers a range of data analysis samples in simple and straightforward R code, providing step-by-step resources and time-saving methods to help you solve data problems efficiently. The first section deals with how to create R functions to avoid the unnecessary duplication of code. You will learn how to prepare, process, and perform sophisticated ETL for heterogeneous data sources with R packages. An example of data manipulation is provided, illustrating how to use the “dplyr” and “data.table” packages to efficiently process larger data structures. We also focus on “ggplot2” and show you how to create advanced figures for data exploration. In addition, you will learn how to build an interactive report using the “ggvis” package. Later chapters offer insight into time series analysis on financial data, while there is detailed information on the hot topic of machine learning, including data classification, regression, clustering, association rule mining, and dimension reduction. By the end of this book, you will understand how to resolve issues and will be able to comfortably offer solutions to problems encountered while performing data analysis.
R for Data Science Cookbook
Credits
www.PacktPub.com
Preface
Free Chapter
Functions in R
Data Preprocessing and Preparation
Visualizing Data with ggplot2
Making Interactive Reports
Simulation from Probability Distributions
Statistical Inference in R
Time Series Mining with R
Index

## Performing two-way ANOVA

Two-way ANOVA can be viewed as an extension of one-way ANOVA because the analysis covers more than two categorical variables rather than just one. In this recipe, we will discuss how to conduct two-way ANOVA in R.

### How to do it…

Perform the following steps to perform two-way ANOVA:

1. First, load the engineer's salary data from engineer.csv:

```>engineer<-read.csv("engineer.csv", header = TRUE)
```
2. Plot the two boxplots of the salary factor in regard to profession and region:

```>par(mfrow=c(1,2))
>boxplot(Salary~Profession, data = engineer,xlab='Profession', ylab = "Salary",main='Salary v.s. Profession')
>boxplot(Salary~Region, data = engineer,xlab='Region', ylab = "Salary",main='Salary v.s. Region')
```

Figure 14: A boxplot of Salary versus Profession and Salary versus Region

3. Also...