> # Perform chi-squared test > chisq_result <- chisq.test(cont_table) > chisq_result Pearson's Chi-squared test data: cont_table X-squared = 105.87, df = 6, p-value < 2.2e-16 please conduct post-hoc test in this data

### Contingency table 
cont_table <- table(iris$Species, iris$Petal.Width)
cont_table

### Perform chi-squared test
chisq_result <- chisq.test(cont_table)
chisq_result

### Post-hoc test
# step 1: install the package
install.packages("vcd")

# step2: load the package
library(vcd)

# step3: test
assocstats(cont_table)