TAO<-c(63,14,44) PAOD<-c(375,235,440) freq_TAO<-121/1050 expected_TAO <- freq_TAO*PAOD chisq.test(TAO, p=expected_TAO, rescale.p=TRUE) chisq.test(TAO, p=PAOD, rescale.p=TRUE) chisq.test(TAO, p=PAOD/sum(PAOD))