This paper aims at studying the elementary chemical composition of 70 ceramic fragments from the São Paulo II archaeological site by instrumental neutron activation analysis. The concentration data was studied next using multivariate statistical methods, such as cluster analysis, principal component analysis and discriminant analysis. The results showed three different chemical groups of samples regarding the similarity/dissimilarity between the samples. Ceramics from each group have been selected and dated using thermoluminescence. The firing temperature of the ceramics was determined by electron paramagnetic resonance.