Untitled

#install.packages("palmerpenguins")
library(palmerpenguins)
data("penguins")

penguin_df <- na.omit(penguins[, c("species", "island")])

head(penguin_df)

## # A tibble: 6 × 2
##   species island   
##   <fct>   <fct>    
## 1 Adelie  Torgersen
## 2 Adelie  Torgersen
## 3 Adelie  Torgersen
## 4 Adelie  Torgersen
## 5 Adelie  Torgersen
## 6 Adelie  Torgersen

table_species_island <- table(penguin_df$species, penguin_df$island)

print(table_species_island)

##            
##             Biscoe Dream Torgersen
##   Adelie        44    56        52
##   Chinstrap      0    68         0
##   Gentoo       124     0         0

chisq_test <- chisq.test(table_species_island)

print(chisq_test)

## 
##  Pearson's Chi-squared test
## 
## data:  table_species_island
## X-squared = 299.55, df = 4, p-value < 2.2e-16

p_value <- format(chisq_test$p.value, scientific = FALSE)

print(p_value)

## [1] "0.000000000000000000000000000000000000000000000000000000000000001354574"