Esse foi um exercício dado em sala de pelo Prof Adilson no segundo semestre de 2006.
Disciplina: Análise Gráfica.
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
curve(dgamma(x, scale = 1.5, shape = 2), from = 0, to = 15, main = “Distribuição Gamma”)
Esse foi um exercício dado em sala de pelo Prof Adilson no segundo semestre de 2006.
Disciplina: Análise Gráfica.
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
curve(dweibull(x, scale = 2, shape = 1), from = 0, to = 15, main = “Distribuição Weibull”)
Esse é um exemplo encontrado no help do software R.
Comando: help (image)
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
require(grDevices)
x <- 10*(1:nrow(volcano))
y <- 10*(1:ncol(volcano))
image(x, y, volcano, col = terrain.colors(100), axes = FALSE)
contour(x, y, volcano, levels = seq(90, 200, by = 5),
add = TRUE, col = “peru”)
axis(1, at = seq(100, 800, by = 100))
axis(2, at = seq(100, 600, by = 100))
box()
title(main = “Maunga Whau Volcano”, font.main = 4)
Esse é um exemplo encontrado no help do software R.
Comando: help (ecdf)
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
Fn <- ecdf(rnorm(12))
tt <- seq(-2,2, by = 0.1)
y <- round(rnorm(12),1); y[3] <- y[1]
Fn12 <- ecdf(y)
require(graphics)
op <- par(mfrow=c(3,1), mgp=c(1.5, 0.8,0), mar= .1+c(3,3,2,1))
F10 <- ecdf(rnorm(10))
plot(F10)
plot(F10, verticals= TRUE, do.points = FALSE)
plot(Fn12 , lwd = 2) ; mtext(”lwd = 2″, adj=1)
Esse foi um exercício dado em sala de pelo Prof Adilson no segundo semestre de 2006.
Disciplina: Análise Gráfica.
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
curve(pnorm(x), from = -4, to = 4, ylab = “P%”)
Esse foi um exercício dado em sala de pelo Prof Adilson no segundo semestre de 2006.
Disciplina: Análise Gráfica.
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
require(Hmisc)
set.seed(12)
x1 <- rnorm(500)
x2 <- runif(500, -2, 2)
x3 <- abs(rnorm(500)) - 2
bpplot(x1, x2, x3)
Esse foi um exercício dado em sala de pelo Prof Adilson no segundo semestre de 2006.
Disciplina: Análise Gráfica.
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
require(”gplots”)
set.seed(12)
vendas.mar <- sample(100:1000, 50)
vendas.abr <- sample(100:800, 50)
marca <- rep(LETTERS[1:5], each = 10)
cidade <- c("Curitiba", "Maringá", "Londrina", "Cascavel", "Foz do Iguaçu")
cidades <- rep(cidade, 10)
dados <- data.frame(cidades, marca, vendas.mar, vendas.abr)
pairs(dados)
Esse foi um exercício dado em sala de pelo Prof Adilson no segundo semestre de 2006.
Disciplina: Análise Gráfica.
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
require(”gplots”)
set.seed(12)
vendas.mar <- sample(100:1000, 50)
vendas.abr <- sample(100:800, 50)
marca <- rep(LETTERS[1:5], each = 10)
cidade <- c("Curitiba", "Maringá", "Londrina", "Cascavel", "Foz do Iguaçu")
cidades <- rep(cidade, 10)
dados <- data.frame(cidades, marca, vendas.mar, vendas.abr)
balloonplot(dados$marca, dados$cidades, dados$vendas.mar)
Esse é um exemplo encontrado no help do software R.
Comando: help (hist)
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
set.seed(14)
x <- rchisq(100, df = 4) hist(x, freq = FALSE, ylim = c(0, 0.2)) curve(dchisq(x, df = 4), col = 2, lty = 2, lwd = 2, add = TRUE)
Esse é um exemplo encontrado no help do software R.
Comando: help (hist)
CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
op <- par(mfrow=c(2, 2))
hist(islands)
utils::str(hist(islands, col=”gray”, labels = TRUE))
hist(sqrt(islands), breaks = 12, col=”lightblue”, border=”pink”)
##– For non-equidistant breaks, counts should NOT be graphed unscaled:
r <- hist(sqrt(islands), breaks = c(4*0:5, 10*3:5, 70, 100, 140), col='blue1')
text(r$mids, r$density, r$counts, adj=c(.5, -.5), col=’blue3′)
sapply(r[2:3], sum)
sum(r$density * diff(r$breaks)) # == 1
lines(r, lty = 3, border = “purple”) # -> lines.histogram(*)
par(op)
