This project will demonstrate your ability to do exploratory data analysis on single variables of data in R and RStudio. The entire project will use the NSCC Student Dataset, which you will need to load into R in question 1.
Download the “nscc_student_data.csv” file from MyOpenMath and use the read.csv() function to store it into an object called “nscc_student_data”. Print the first few lines of the dataset using the head() function. Also print the structure of the dataset using the str() function.
# Load dataset in and store as object "nscc_student_data"
nscc_student_data<-read.csv("nscc_student_data.csv")
# Preview first 6 lines of dataset
head(nscc_student_data)
## Gender PulseRate CoinFlip1 CoinFlip2 Height ShoeLength Age Siblings
## 1 Female 64 5 5 62 11.00 19 4
## 2 Female 75 4 6 62 11.00 21 3
## 3 Female 74 6 1 60 10.00 25 2
## 4 Female 65 4 4 62 10.75 19 1
## 5 Female NA NA NA 66 NA 26 6
## 6 Female 72 6 5 67 9.75 21 1
## RandomNum HoursWorking Credits Birthday ProfsAge Coffee VoterReg
## 1 797 35 13 July 5 31 No Yes
## 2 749 25 12 December 27 30 Yes Yes
## 3 13 30 6 January 31 29 Yes No
## 4 613 18 9 6-13 31 Yes Yes
## 5 53 24 15 02-15 32 No Yes
## 6 836 15 9 april 14 32 No Yes
# Structure of dataset
str(nscc_student_data)
## 'data.frame': 40 obs. of 15 variables:
## $ Gender : Factor w/ 2 levels "Female","Male": 1 1 1 1 1 1 2 2 1 2 ...
## $ PulseRate : int 64 75 74 65 NA 72 72 60 66 60 ...
## $ CoinFlip1 : int 5 4 6 4 NA 6 6 3 7 6 ...
## $ CoinFlip2 : int 5 6 1 4 NA 5 6 5 8 5 ...
## $ Height : num 62 62 60 62 66 ...
## $ ShoeLength : num 11 11 10 10.8 NA ...
## $ Age : int 19 21 25 19 26 21 19 24 24 20 ...
## $ Siblings : int 4 3 2 1 6 1 2 2 3 1 ...
## $ RandomNum : int 797 749 13 613 53 836 423 16 12 543 ...
## $ HoursWorking: int 35 25 30 18 24 15 20 0 40 30 ...
## $ Credits : int 13 12 6 9 15 9 15 15 13 16 ...
## $ Birthday : Factor w/ 40 levels "02-15","03.14.1984",..: 32 25 30 18 1 21 19 27 35 31 ...
## $ ProfsAge : int 31 30 29 31 32 32 28 28 31 28 ...
## $ Coffee : Factor w/ 2 levels "No","Yes": 1 2 2 2 1 1 2 2 2 1 ...
## $ VoterReg : Factor w/ 2 levels "No","Yes": 2 2 1 2 2 2 2 2 2 1 ...
a.) What are the dimensions of the nscc_student_data dataframe?
# Find the dimensions of the nscc_student_data dataframe
dim(nscc_student_data)
## [1] 40 15
The dimensions are 40 obs. by 15 variables
b.) The chunk of code below will tell you how many values in the PulseRate variable exist (FALSE) and how many are NA (TRUE). How many values in the variable are missing?
# How many values in PulseRate variable are missing
table(is.na(nscc_student_data$PulseRate))
##
## FALSE TRUE
## 38 2
There are 2 values missing in the PulseRate variable
Use an r chunk to calculate the mean and median of the pulse rate variable. Do they differ by much? If yes, explain why and which would be a better choice as the “center” or “average” of this variable.
#Calculate the mean and median of the pulse rate variable
mean(nscc_student_data$PulseRate, na.rm = TRUE)
## [1] 73.47368
median(nscc_student_data$PulseRate, na.rm = TRUE)
## [1] 70.5
** Yes, both the mean and median differ by 3. The better choice for the “average” would be the median because this variable has a few outliers. When there are outliers, it could possibly affect the overall “average”**
Use an r chunk to calculate the sample standard deviation of the pulse rate variable.
#Standard Deviation of the pulse rate variable
sd(nscc_student_data$PulseRate,na.rm = TRUE)
## [1] 12.51105
The standard deviation of the ppulse rate variable is 12.51105
Use an r chunk to calculate the Five Number Summary and IQR of the pulse rate variable. Clearly state your answer for each below the r chunk. Based on the definition of outliers being more than 1.5 IQRs below Q1 or above Q3, are there any values in the pulse rate variable that are considered outliers? Clearly state below the thresholds for a data to be considered an outlier.
#Calculate the Five number summary and IQR
summary(nscc_student_data$PulseRate, na.rm = TRUE)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 50.00 64.25 70.50 73.47 83.75 98.00 2
IQR(nscc_student_data$PulseRate,na.rm = TRUE)
## [1] 19.5
**The five number summary of the PulseRate variable is Min=50, Q1=64.25, Median=70.50, Q3=83.75, Max=98.00. The IQR of the PulseRate variable is 19.5.
Any data values that are below 64.25 or above 83.75 would be considered to be outliers. Based on those values, there are outliers and those outliers are 50, 56, 60, 61, 62, 85, 87, 88, 89, 92, 96, 98 .**
The Gender variable gives whether students identified as male or female. Create a table and a barplot of that variable.
#Create a table and barplot for the gender variable
table(nscc_student_data$Gender)
##
## Female Male
## 27 13
barplot(table(nscc_student_data$Gender))
Split the dataframe into two subsets – one that has all the males and another that has all the females. Store them into objects called “NSCC_males” and “NSCC_females”. The first one has been done for you as a template.
# Create males subset
NSCC_males <- subset(nscc_student_data, nscc_student_data$Gender == "Male")
# Create females subset
NSCC_females <- subset(nscc_student_data, nscc_student_data$Gender == "Female")
What is the Five Number Summary for the pulse rate variable for each of the male and female subsets.
# Five Number Summary of males subset
summary(NSCC_males$PulseRate)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.00 60.00 71.00 70.85 80.00 96.00
# Five Number Summary of females subset
summary(NSCC_females$PulseRate)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 56.00 65.00 70.00 74.84 88.00 98.00 2
The Five Number Summary of each dataset are:
Males: Females: Min=50 Min=56 Q1=60 Q1=65 Median=71 Median=70 Q3=80 Q3=88 Max=96 Max=98
Create side-by-side boxplots for the pulse rate variable each of the male and female subsets. Is there any noticeable difference between the spread of the variables?
# Create side-by-side boxplots for each subset
boxplot(NSCC_males$PulseRate, NSCC_females$PulseRate)
Yes there is a noticeable difference, and it is that feales have a higher pulse rate.
Create a frequency distribution for how many males and females answered “Yes” or “No” to the variable “Coffee” by using the table() function. What percent of this sample of NSCC students drink coffee? Is there any noticeable difference in coffee drinking based on gender?
# Male Coffee Drinkers
table(NSCC_males$Coffee)
##
## No Yes
## 3 10
# Females Coffee Drinkers
table(NSCC_females$Coffee)
##
## No Yes
## 7 20
# Percent of Males that Drink Coffee
10/13
## [1] 0.7692308
# Percent of Females that Drink Coffee
20/27
## [1] 0.7407407
Overall of the sample of students in this data, 75% of these students drink coffee, seperated by gender, 74.04% of females drink coffe and 76.92% of males drink coffee. There is a noticeable difference of 2% between bothsubsets