R Markdown

Create an RMarkdown that computes the measures of central tendency and measures of variability and the relationships for each of the seven variables in the attitude dataset. Use the functions below:

mean, median, mode, max, min, range, quantile, IQR, var( ), sd( ), and cor( )

Rating

paste("The Mean is:", mean(attitude$rating, na.rm=TRUE))
## [1] "The Mean is: 64.6333333333333"
paste("The Median is:", median(attitude$rating))
## [1] "The Median is: 65.5"
#Need a table to calculate mode or else it outputs Numeric
table(attitude$rating) 
## 
## 40 43 48 50 53 58 61 63 64 65 66 67 68 69 71 72 74 77 78 81 82 85 
##  1  2  1  2  1  1  1  2  2  2  1  2  1  1  2  1  1  1  1  2  1  1
modeforrating <- table(attitude$rating)
as.numeric(names(modeforrating)[modeforrating==max(modeforrating)])
## [1] 43 50 63 64 65 67 71 81
paste("The Max is:", max(attitude$rating))
## [1] "The Max is: 85"
paste("The Min is:", min(attitude$rating))
## [1] "The Min is: 40"
paste("The Range is:", range(attitude$rating))
## [1] "The Range is: 40" "The Range is: 85"
quantile(attitude$rating)
##    0%   25%   50%   75%  100% 
## 40.00 58.75 65.50 71.75 85.00
paste("The IQR is:", IQR(attitude$rating))
## [1] "The IQR is: 13"
paste("The Variance is:", var(attitude$rating))
## [1] "The Variance is: 148.171264367816"
paste("The Standard deviation is:", sd(attitude$rating, na.rm = TRUE))
## [1] "The Standard deviation is: 12.1725619475859"
summary(attitude$rating)                                                
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   40.00   58.75   65.50   64.63   71.75   85.00

Complaints

paste("The Mean is:", mean(attitude$complaints, na.rm=TRUE))             
## [1] "The Mean is: 66.6"
paste("The Median is:", median(attitude$complaints))                     
## [1] "The Median is: 65"
#Need a table to calculate mode or else it outputs Numeric               
table(attitude$complaints)                                
## 
## 37 40 51 53 54 55 57 58 60 61 62 63 64 66 67 70 75 77 78 82 83 85 90 
##  1  1  1  1  1  1  1  1  2  2  1  1  1  1  1  2  2  2  1  2  1  2  1
modeforcomplaints <- table(attitude$complaints)                         
as.numeric(names(modeforcomplaints)[modeforcomplaints==max(modeforcomplaints)])                                                                   
## [1] 60 61 70 75 77 82 85
paste("The Max is:", max(attitude$complaints))                           
## [1] "The Max is: 90"
paste("The Min is:", min(attitude$complaints))            
## [1] "The Min is: 37"
paste("The Range is:", range(attitude$complaints))
## [1] "The Range is: 37" "The Range is: 90"
quantile(attitude$complaints)
##   0%  25%  50%  75% 100% 
## 37.0 58.5 65.0 77.0 90.0
paste("The IQR is:", IQR(attitude$complaints))                           
## [1] "The IQR is: 18.5"
paste("The Variance is:", var(attitude$complaints))                       
## [1] "The Variance is: 177.28275862069"
paste("The Standard deviation is:", sd(attitude$complaints, na.rm =TRUE))
## [1] "The Standard deviation is: 13.3147571746799"
summary(attitude$complaints)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    37.0    58.5    65.0    66.6    77.0    90.0

Privileges

paste("The Mean is:", mean(attitude$privileges, na.rm=TRUE))             
## [1] "The Mean is: 53.1333333333333"
paste("The Median is:", median(attitude$privileges))                     
## [1] "The Median is: 51.5"
#Need a table to calculate mode or else it outputs Numeric              
table(attitude$privileges)
## 
## 30 33 39 42 44 45 46 47 49 50 51 52 53 54 56 57 58 64 65 66 68 71 72 83 
##  1  1  1  3  1  2  1  1  1  2  1  2  1  1  1  1  1  1  1  1  2  1  1  1
modeforprivileges <- table(attitude$privileges)
as.numeric(names(modeforprivileges)[modeforprivileges==max(modeforprivileges)])
## [1] 42
paste("The Max is:", max(attitude$privileges))            
## [1] "The Max is: 83"
paste("The Min is:", min(attitude$privileges))                          
## [1] "The Min is: 30"
paste("The Range is:", range(attitude$privileges))
## [1] "The Range is: 30" "The Range is: 83"
quantile(attitude$privileges)                             
##   0%  25%  50%  75% 100% 
## 30.0 45.0 51.5 62.5 83.0
paste("The IQR is:", IQR(attitude$privileges))
## [1] "The IQR is: 17.5"
paste("The Variance is:", var(attitude$privileges))
## [1] "The Variance is: 149.705747126437"
paste("The Standard deviation is:", sd(attitude$privileges, na.rm =      TRUE))
## [1] "The Standard deviation is: 12.2354299935244"
summary(attitude$privileges)                              
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   30.00   45.00   51.50   53.13   62.50   83.00

Learning

paste("The Mean is:", mean(attitude$learning, na.rm=TRUE))
## [1] "The Mean is: 56.3666666666667"
paste("The Median is:", median(attitude$learning))
## [1] "The Median is: 56.5"
#Need a table to calculate mode or else it outputs Numeric               
table(attitude$learning)
## 
## 34 39 42 44 45 47 48 50 54 55 56 57 58 59 62 63 66 67 69 71 72 74 75 
##  1  2  1  1  2  2  1  1  2  1  1  1  2  1  1  1  1  1  2  1  2  1  1
modeforlearning <- table(attitude$learning)                              
as.numeric(names(modeforlearning)[modeforlearning==max(modeforlearning)])
## [1] 39 45 47 54 58 69 72
paste("The Max is:", max(attitude$learning))                             
## [1] "The Max is: 75"
paste("The Min is:", min(attitude$learning))
## [1] "The Min is: 34"
paste("The Range is:", range(attitude$learning))          
## [1] "The Range is: 34" "The Range is: 75"
quantile(attitude$learning)
##    0%   25%   50%   75%  100% 
## 34.00 47.00 56.50 66.75 75.00
paste("The IQR is:", IQR(attitude$learning))
## [1] "The IQR is: 19.75"
paste("The Variance is:", var(attitude$learning))                         
## [1] "The Variance is: 137.757471264368"
paste("The Standard deviation is:", sd(attitude$learning, na.rm = TRUE)) 
## [1] "The Standard deviation is: 11.7370128765529"
summary(attitude$learning)                                              
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   34.00   47.00   56.50   56.37   66.75   75.00

Raises

paste("The Mean is:", mean(attitude$raises, na.rm=TRUE))
## [1] "The Mean is: 64.6333333333333"
paste("The Median is:", median(attitude$raises))
## [1] "The Median is: 63.5"
#Need a table to calculate mode or else it outputs Numeric
table(attitude$raises)                                                  
## 
## 43 50 51 54 55 58 59 60 61 62 63 64 66 70 71 75 76 77 79 80 88 
##  1  1  1  2  2  1  2  1  1  1  2  2  3  1  2  1  1  1  2  1  1
modeforraises <- table(attitude$raises)
as.numeric(names(modeforraises)[modeforraises==max(modeforraises)])         
## [1] 66
paste("The Max is:", max(attitude$raises))                               
## [1] "The Max is: 88"
paste("The Min is:", min(attitude$raises))                
## [1] "The Min is: 43"
paste("The Range is:", range(attitude$raises))                           
## [1] "The Range is: 43" "The Range is: 88"
quantile(attitude$raises)                                                
##    0%   25%   50%   75%  100% 
## 43.00 58.25 63.50 71.00 88.00
paste("The IQR is:", IQR(attitude$raises))                              
## [1] "The IQR is: 12.75"
paste("The Variance is:", var(attitude$raises))
## [1] "The Variance is: 108.102298850575"
paste("The Standard deviation is:", sd(attitude$raises, na.rm = TRUE))   
## [1] "The Standard deviation is: 10.3972255361983"
summary(attitude$raises)                                                 
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   43.00   58.25   63.50   64.63   71.00   88.00

Critical

paste("The Mean is:", mean(attitude$critical, na.rm=TRUE))
## [1] "The Mean is: 74.7666666666667"
paste("The Median is:", median(attitude$critical))
## [1] "The Median is: 77.5"
#Need a table to calculate mode or else it outputs Numeric               
table(attitude$critical)
## 
## 49 54 57 63 64 66 67 68 73 74 75 76 77 78 79 80 83 84 85 86 92 
##  1  1  1  1  1  1  1  1  1  2  1  1  2  3  1  4  3  1  1  1  1
modeforcritical <- table(attitude$critical)                              
as.numeric(names(modeforcritical)[modeforcritical==max(modeforcritical)])
## [1] 80
paste("The Max is:", max(attitude$critical))                             
## [1] "The Max is: 92"
paste("The Min is:", min(attitude$critical))
## [1] "The Min is: 49"
paste("The Range is:", range(attitude$critical))          
## [1] "The Range is: 49" "The Range is: 92"
quantile(attitude$critical)
##    0%   25%   50%   75%  100% 
## 49.00 69.25 77.50 80.00 92.00
paste("The IQR is:", IQR(attitude$critical))
## [1] "The IQR is: 10.75"
paste("The Variance is:", var(attitude$critical))                         
## [1] "The Variance is: 97.9091954022988"
paste("The Standard deviation is:", sd(attitude$critical, na.rm = TRUE)) 
## [1] "The Standard deviation is: 9.89490754895157"
summary(attitude$critical)                                              
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   49.00   69.25   77.50   74.77   80.00   92.00

Advance

paste("The Mean is:", mean(attitude$advance, na.rm=TRUE))
## [1] "The Mean is: 42.9333333333333"
paste("The Median is:", median(attitude$advance))
## [1] "The Median is: 41"
#Need a table to calculate mode or else it outputs Numeric
table(attitude$advance)                                                 
## 
## 25 31 33 34 35 36 37 38 39 41 45 46 47 48 49 52 55 60 63 72 
##  1  1  2  2  3  1  1  1  1  4  1  2  2  1  2  1  1  1  1  1
modeforadvance <- table(attitude$advance)
as.numeric(names(modeforadvance)[modeforadvance==max(modeforadvance)])
## [1] 41
paste("The Max is:", max(attitude$advance))
## [1] "The Max is: 72"
paste("The Min is:", min(attitude$advance))                              
## [1] "The Min is: 25"
paste("The Range is:", range(attitude$advance))                         
## [1] "The Range is: 25" "The Range is: 72"
quantile(attitude$advance)
##    0%   25%   50%   75%  100% 
## 25.00 35.00 41.00 47.75 72.00
paste("The IQR is:", IQR(attitude$advance))                              
## [1] "The IQR is: 12.75"
paste("The Variance is:", var(attitude$advance))                          
## [1] "The Variance is: 105.857471264368"
paste("The Standard deviation is:", sd(attitude$advance, na.rm = TRUE)) 
## [1] "The Standard deviation is: 10.2887060053424"
summary(attitude$advance)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   25.00   35.00   41.00   42.93   47.75   72.00

Relationships

cor(attitude)                                                            
##               rating complaints privileges  learning    raises  critical
## rating     1.0000000  0.8254176  0.4261169 0.6236782 0.5901390 0.1564392
## complaints 0.8254176  1.0000000  0.5582882 0.5967358 0.6691975 0.1877143
## privileges 0.4261169  0.5582882  1.0000000 0.4933310 0.4454779 0.1472331
## learning   0.6236782  0.5967358  0.4933310 1.0000000 0.6403144 0.1159652
## raises     0.5901390  0.6691975  0.4454779 0.6403144 1.0000000 0.3768830
## critical   0.1564392  0.1877143  0.1472331 0.1159652 0.3768830 1.0000000
## advance    0.1550863  0.2245796  0.3432934 0.5316198 0.5741862 0.2833432
##              advance
## rating     0.1550863
## complaints 0.2245796
## privileges 0.3432934
## learning   0.5316198
## raises     0.5741862
## critical   0.2833432
## advance    1.0000000