Dataset used: Airline Stats

data<-read.csv("airline_stats.csv")
View(data)
boxplot(data$pct_carrier_delay ~ data$airline ,
        xlab="Airlines" , 
        ylab="Percentage of Carrier Delay (%)" , 
        main="(%) Carrier Delay by Airlines" , 
        col="#c4dec1")
mtext("Prepared by Nabeeha", side=3, line=0.5, font=1)
mtext("Source: Airline Stats Dataset", side=1, line=4, font=1)

summary(data$pct_carrier_delay, data$airline)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
##   0.000   4.145   6.357   7.041   9.140 100.000      28