Module 4 Chi-Squared - Goodness of Fit

#install.packages("readxl")
library(readxl)

## Warning: package 'readxl' was built under R version 4.2.3

DatasetName <- read_excel("C:/Users/whitn/OneDrive - Saint Louis University/Applied Analytics Fall 2026/AA 5221/Week 04/Whitney_Moss_Assignment_4/A4Q1.xlsx")
table(DatasetName$flavor)

## 
##  Chocolate      Mango Strawberry    Vanilla 
##         87         32         57         74

observed <- table(DatasetName$flavor)
observed

## 
##  Chocolate      Mango Strawberry    Vanilla 
##         87         32         57         74

barplot(observed,
  main = "Summer Icecream Flavor",
  xlab = "flavor",
  ylab = "Count",
  col = rainbow(length(observed)))

expected <- c(.348, .128, .228, .296)



chi_result <- chisq.test(x = observed, p = expected)
chi_resultchi_result <- chisq.test(x = observed, p = expected)
chi_result

## 
##  Chi-squared test for given probabilities
## 
## data:  observed
## X-squared = 0, df = 3, p-value = 1

w <- sqrt(as.numeric(chi_result$statistic) / sum(observed))
w

## [1] 0