# Create a data frame 
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) 

# Calculate the mean of ice cream sales 
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")


# Calculate the standard deviation of ice cream sales
sd_sales <- sd(data$ice_cream_sales) 
cat("Standard Deviation of Sales:", sd_sales, "\n")


# Calculate the correlation coefficient between temperature and sales 
correlation_coefficient <- cor(data$temperature, data$ice_cream_sales) 
cat("Correlation between Temperature and Sales:", correlation_coefficient, "\n")
