data <- data.frame(
  temperature = c(15, 18, 22, 25, 28, 30, 32, 35, 38),
  ice_cream_sales = c(50, 70, 100, 120, 150, 180, 200, 220, 250),
)  
print(data) 

mean_sales <- mean(data$ice_cream_sales)
cat("Mean Ice Cream Sales:", mean_sales, "\n")


get_mode <- function(v) { 
   unique_v <- unique(v) 
   unique_v[which.max(tabulate(match(v, unique_v)))] 
}

mode_sales <- get_mode(data$ice_cream_sales) 
cat("Mode Ice Cream Sales:", mode_sales, "\n")


sd_sales <- sd(data$ice_cream_sales)
cat("Standard Deviation of Sales:", sd_sales, "\n")


correlation_coefficient <- cor(data$temperature, data$ice_cream_sales)
cat("Correlation between Temperature and Sales:", correlation_coefficient, "\n")
