library(palmerpenguins)

data <- na.omit(palmerpenguins::penguins)

num_simulations <- 3200
sample_size <- 5

results <- numeric(num_simulations)

# Flipper length simulation:
for (i in 1:num_simulations) {
  sampled_data <- data[sample(nrow(data), sample_size, replace = TRUE), ]
  results[i] <- mean(sampled_data$flipper_length_mm)
}

# Calculate the overall average from all simulations:
estimated_average_flipper_length <- mean(results)
cat("Estimated Average Flipper Length (mm):", estimated_average_flipper_length, "\n")