%Crime rate Case
>> f=[334 372 327 328]
f =
   334   372   327   328

>> e=1361*[1/4 1/4 1/4 1/4]
e =
  340.2500  340.2500  340.2500  340.2500

>> Q=sum((f-e).^2./e)
Q =
    4.0345

>> 1-chi2cdf(Q,4-1)
ans =
    0.2578 %can't reject H_0

%Education vs. Marital Statue
>>f=[550 60
     681 144]

>>e= [611*1231 611*205
    825*1231 825*205]/1436
e =
  523.7751   87.2249
  707.2249  117.7751

>>(f-e).^2./e
ans =
    1.3131    8.4975
    0.9725    5.8395

>>Q=sum(sum((f-e).^2./e))
Q =
    16.6225

>>1-chi2cdf(Q,(3-1)*(3-1))
ans =
    0.0023 %reject H_0

%Austen
f=[83 29 15 22 43 4]
>>sum(f)
ans = 
     196

austen = [434 62 86 236 161 38]
>> sum(austen)
ans =
        1017

>>e =austen/1017*196
e =
   83.6421   11.9489   16.5742   45.4828   31.0285    7.3235

>>Q=sum((f-e).^2./e) 
Q =
   42.7378

>>1-chi2cdf(Q,6-1)
ans =

   4.1759e-08 %reject H_0