Diagrama de dispersão

# o diagrama de dispersão representa duas (ou mais) variáveis quantitativas em um gráfico bidimensional

# é utilizado para observar possível relação entre duas (ou mais) variáveis

# vamos simular quatro variáveis x, y1, y2, y3

x=1:12

y1=c(10,15,22,26,33,38,42,41,43,44,42,43)

y2=c(50,51,46,42,39,36,29,24,20,15,10,6)

y3=rnorm(12,20,12)

# y1 em função de x (relação linear positiva até um platô)

plot(y1~x)

# y2 em função de x (relação linear negativa)

plot(y2~x)

# y3 em função de x (não tem relação)

plot(y3~x)

# adicionando algumas modificações e juntado em um mesmo gráfico as respostas y1 e y2

plot(y1~x, type="o", lty=2, col="dark red", xlab="Variável X", ylab="Variável Y", ylim=c(min(y2),max(y2)), lwd=3)

points(y2~x, type="o", lty=3, col="dark green", lwd=3, pch=2)

legend(6,50, c("y1","y2"), lty=c(2,3), col=c("dark red", "dark green"), pch=c(1,2), bty="n", cex=1.5, lwd=3)

# vamos simular a degradabilidade (%) de três alimentos (deg1, deg2 e deg3) em função do tempo em horas de incubação em rúmen bovino

tempo=c(6,9,24,48,56,72,96)

deg1=c(10,12,22,35,45,65,78)

deg2=c(23,45,66,82,88,92,93)

deg3=c(9,11,22,36,49,58,66)

# para observar em um único gráfico

dados=data.frame(tempo, deg1, deg2, deg3) # x deve ser a primeira coluna

# carregar o pacote ds

require(ds)

# usar a função dplot() do pacote ds

dplot(dados)

# com algumas alterações

dplot(dados, xlab = "Tempo (horas)", ylab = "Degradabilidade (%)", position = 3)

# agora um exemplo com dados carregados do pacote epr

require(epr)

data(data5)

dplot(data5, position = 9)

# para informações sobre parâmetros gráficos da função plot faça

?plot

?par

# para informações sobre parâmetros gráficos da função dplot faça

?dplot