Load Data and Prepare Analysis
# Load the data
url <- "https://raw.githubusercontent.com/tmatis12/datafiles/main/normtemp.csv"
data <- read.csv(url)
# Separate data by gender
males <- data[data$Sex == 1, ]
females <- data[data$Sex == 2, ]
# Extract heart rate data
male_heartrate <- males$Beats
female_heartrate <- females$Beats
# Check sample sizes
cat("Number of males:", length(male_heartrate), "\n")
## Number of males: 65
cat("Number of females:", length(female_heartrate), "\n")
## Number of females: 65
cat("Descriptive Statistics for Male Heart Rate\n")
## Descriptive Statistics for Male Heart Rate
cat("==========================================\n")
## ==========================================
cat("Minimum:", min(male_heartrate), "\n")
## Minimum: 58
cat("Maximum:", max(male_heartrate), "\n")
## Maximum: 86
cat("Mean:", round(mean(male_heartrate), 2), "\n")
## Mean: 73.37
cat("Standard Deviation:", round(sd(male_heartrate), 2), "\n")
## Standard Deviation: 5.88
cat("Median:", median(male_heartrate), "\n")
## Median: 73
cat("Q1 (25th percentile):", quantile(male_heartrate, 0.25), "\n")
## Q1 (25th percentile): 70
cat("Q3 (75th percentile):", quantile(male_heartrate, 0.75), "\n")
## Q3 (75th percentile): 78
hist(male_heartrate,
col = "blue",
main = "Histogram of Male Resting Heart Rate",
xlab = "Resting Heart Rate (beats per minute)",
ylab = "Frequency")
qqnorm(male_heartrate,
col = "blue",
main = "Normal Q-Q Plot for Male Resting Heart Rate")
qqline(male_heartrate, col = "red")
cat("Descriptive Statistics for Female Heart Rate\n")
## Descriptive Statistics for Female Heart Rate
cat("============================================\n")
## ============================================
cat("Minimum:", min(female_heartrate), "\n")
## Minimum: 57
cat("Maximum:", max(female_heartrate), "\n")
## Maximum: 89
cat("Mean:", round(mean(female_heartrate), 2), "\n")
## Mean: 74.15
cat("Standard Deviation:", round(sd(female_heartrate), 2), "\n")
## Standard Deviation: 8.11
cat("Median:", median(female_heartrate), "\n")
## Median: 76
cat("Q1 (25th percentile):", quantile(female_heartrate, 0.25), "\n")
## Q1 (25th percentile): 68
cat("Q3 (75th percentile):", quantile(female_heartrate, 0.75), "\n")
## Q3 (75th percentile): 80
hist(female_heartrate,
col = "pink",
main = "Histogram of Female Resting Heart Rate",
xlab = "Resting Heart Rate (beats per minute)",
ylab = "Frequency")
qqnorm(female_heartrate,
col = "pink",
main = "Normal Q-Q Plot for Female Resting Heart Rate")
qqline(female_heartrate, col = "red")
boxplot(male_heartrate, female_heartrate,
names = c("Male", "Female"),
col = c("blue", "pink"),
main = "Comparison of Resting Heart Rate by Gender",
ylab = "Resting Heart Rate (beats per minute)")
# All code used in this analysis
url <- "https://raw.githubusercontent.com/tmatis12/datafiles/main/normtemp.csv"
data <- read.csv(url)
males <- data[data$Sex == 1, ]
females <- data[data$Sex == 2, ]
male_heartrate <- males$Beats
female_heartrate <- females$Beats
# Male statistics
min(male_heartrate)
## [1] 58
max(male_heartrate)
## [1] 86
mean(male_heartrate)
## [1] 73.36923
sd(male_heartrate)
## [1] 5.875184
median(male_heartrate)
## [1] 73
quantile(male_heartrate, 0.25)
## 25%
## 70
quantile(male_heartrate, 0.75)
## 75%
## 78
# Female statistics
min(female_heartrate)
## [1] 57
max(female_heartrate)
## [1] 89
mean(female_heartrate)
## [1] 74.15385
sd(female_heartrate)
## [1] 8.105227
median(female_heartrate)
## [1] 76
quantile(female_heartrate, 0.25)
## 25%
## 68
quantile(female_heartrate, 0.75)
## 75%
## 80
# Plots
hist(male_heartrate, col = "blue", main = "Histogram of Male Resting Heart Rate")
hist(female_heartrate, col = "pink", main = "Histogram of Female Resting Heart Rate")
qqnorm(male_heartrate, col = "blue", main = "Normal Q-Q Plot for Male Resting Heart Rate")
qqline(male_heartrate, col = "red")
qqnorm(female_heartrate, col = "pink", main = "Normal Q-Q Plot for Female Resting Heart Rate")
qqline(female_heartrate, col = "red")
boxplot(male_heartrate, female_heartrate, names = c("Male", "Female"), col = c("blue", "pink"))
