A explicação do script a seguir está no artigo Matriz de Transição Litológica – Teste do chi-quadrado em Linguagem
#dados segundo Davis 2002 p. 169
MT=as.table(rbind(c(18,0,5,0),c(0,5,2,0),
c(5,2,18,3),c(0,0,3,2)))
dimnames(MT)=list(De=c("A","B","C","D"),
Para=c("A","B","C","D"))
MT
x2=chisq.test(MT)
round(x2$expected,3)
x2