> layout(matrix(c(1, 1, 2, 3), 2, 2, byrow = TRUE))

> boxplot(free.sulfur.dioxide  ~ quality,

+         data= RED.Wine.14MARCH ,

+         main=" Variable #1 with Quality",

+         sub=" Boxplot #1",

+         col.main="orange",

+         ylab="Values of Variable #1 ",

+         xlab="Quality Scores Range - MIN=3 and MAX =8",

+         col.axis="black",

+         col.lab="blue",

+         cex.lab=1,

+         cex.axis=1,

+         cex.main=2,

+         cex.sub=1,

+         col="pink",

+         notch=FALSE,

+         pch=23,       

+         bg="green",

+         las=1)

> boxplot(residual.sugar ~ quality,

+         data= RED.Wine.14MARCH ,

+         main=" Variable #2 with Quality",

+         sub=" Boxplot #2",

+         col.main="orange",

+         ylab="Values of Variable #2 ",

+         xlab="Quality Scores Range - MIN=3 and MAX =8",

+         col.axis="black",

+         col.lab="blue",

+         cex.lab=1,

+         cex.axis=1,

+         cex.main=2,

+         cex.sub=1,

+         col="grey",

+         notch=FALSE,

+         pch=23,       

+         bg="green",

+         las=1)

> boxplot(citric.acid  ~ quality,

+         data= RED.Wine.14MARCH ,

+         main=" Variable #3 with Quality",

+         sub="Boxplot #3",

+         col.main="orange",

+         ylab="Values of Variable #3",

+         xlab="Quality Scores Range - MIN=3 and MAX =8",

+         col.axis="black",

+         col.lab="blue",

+         cex.lab=1,

+         cex.axis=1,

+         cex.main=2,

+         cex.sub=1,

+         col="yellow",

+         notch=FALSE,

+         pch=23,       

+         bg="green",

+         las=1)

>