# Goodness-of-fit test
# Example data is poll from 2003 California gubernatorial recall election, with
# proportions in p representing actual election results.

x=c(315,197,141,39,16,79)
p=c(.486,.315,.125,.028,.006,.04)
chisq.test(x,p=p)

# Calculate everything instead.
n=sum(x)
e=n*p
Xsq=sum((x-e)^2/e)
df=length(x)-1
Pval=1-pchisq(Xsq,df)
Xsq
df
Pval