> scatterplot(alcalinidade~hora | amostra, data=dados, smooth=F, legend.plot=F, reg.line=F,pch=c(16,17,15),col=c("blue","green","red"))

alcalinidade - y

hora = horário de coleta - x

amostra - fator

Gera um gráfico de "dispersão", com boxplots invés de pontos, isso pq minha variável hora está sendo considerada fator, já mudei mas não estou conseguindo. Alguém sabe como forçar a utilização de pontos?