Distribuição Gamma

quarta-feira, 31 de outubro de 2007

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")

Distribuição Weibull

terça-feira, 30 de outubro de 2007

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")

Imagem de um Vulcão

segunda-feira, 29 de outubro de 2007

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)

Distribuição Empírica

domingo, 28 de outubro de 2007

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)

Distribuição Normal Acumulada

sábado, 27 de outubro de 2007

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%")

Box Plot de Percentuais

sexta-feira, 26 de outubro de 2007

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)

Matriz de Dispersão de Dados

quinta-feira, 25 de outubro de 2007

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)

Gráfico Balloon

quarta-feira, 24 de outubro de 2007

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)

Histograma (04)

terça-feira, 23 de outubro de 2007

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)

Histograma (03)

segunda-feira, 22 de outubro de 2007

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)

Histograma (02)

domingo, 21 de outubro de 2007

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
set.seed(123)
x <- sample(1:100, 50, rep = TRUE)
hist.meu <- hist(x)
text(hist.meu$mid, hist.meu$counts + 0.2, hist.meu$counts)

Histograma (01)

sábado, 20 de outubro de 2007

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
set.seed(123)
x <- sample(1:100, 50, rep = TRUE)
hist(x)

Gráfico de Barras (05)

sexta-feira, 19 de outubro de 2007

Esse é um exemplo encontrado no help do software R.
Comando: help (barplot)



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
require(grDevices)
tN <- table(Ni <- stats::rpois(100, lambda=5))
r <- barplot(tN, col=rainbow(20))
#- type = "h" plotting *is* 'bar'plot
lines(r, tN, type='h', col='red', lwd=2)

Gráfico de Barras (04)

quinta-feira, 18 de outubro de 2007

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(MASS)
data(cats)
attach(cats)
aptidao <- ifelse(Hwt >= 9.5, "apto", "inapto")
sex.t <- table(Sex, aptidao)
barplot(sex.t, beside = TRUE, legend.text = rownames(sex.t), col = c("white", "gray"))

Gráfico de Barras (03)

quarta-feira, 17 de outubro de 2007

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(MASS)
data(cats)
attach(cats)
aptidao <- ifelse(Hwt >= 9.5, "apto", "inapto")
sex.t <- table(Sex, aptidao)
barplot(sex.t)

Gráfico de Barras (02)

terça-feira, 16 de outubro de 2007

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
cores <- c(0.3, 0.25, 0.2, 0.15, 0.1)
names(cores) <- c("este", "ou", "aquele", "estat", "2007")
barplot(cores, col = c("green", "red", "blue", "orange", "black"))

Gráfico de Barras (01)

segunda-feira, 15 de outubro de 2007

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
cores <- c(0.3, 0.25, 0.2, 0.15, 0.1)
names(cores) <- c("este", "ou", "aquele", "estat", "2007")
barplot(cores)

Gráfico de Setores ou Pizza (02)

domingo, 14 de outubro de 2007

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
cores <- c(0.3, 0.25, 0.2, 0.15, 0.1)
names(cores) <- c("aquele", "ou", "este", "2007", " estatística")
pie(cores)

Gráfico de Setores ou Pizza (01)

sábado, 13 de outubro de 2007

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(MASS)
data(cats)
attach(cats)
aptidao <- ifelse(Hwt >= 9.5, "apto", "inapto")
sex.t <- table(Sex, aptidao) pie(margin.table(sex.t, 2))

Cara de Boneco

sexta-feira, 12 de outubro de 2007

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
plot(1:10, 1:10, type = "n", ann = FALSE, axes = FALSE)
symbols(5, 5, 2, add = TRUE, inc = FALSE)
symbols(c(4, 6), c(6, 6), c(0.3, 0.3), add = TRUE, inc = FALSE)
points(c(4.041275, 6.086176), c(6.006826, 6.006826), pch = c(44,
44), cex = 4, col = "blue")
points(c(7.022771, 2.997524), c(4.974744, 5.011604), pch = c(41,
40), cex = 7)
segments(4.467296, 3.426621, 5.617553, 3.297611, col = "red")
segments(4.22, 7.1, 5.04, 9, 1)
segments(5.04, 9.01, 5.96, 7.02)
symbols(5.16, 8.97, 0.2, add = TRUE, inc = FALSE, bg = "yellow")
points(4.989172, 4.827304, pch = 94, cex = 3)
points(c(3.49, 3.54, 3.84, 4.07), c(4.71, 4.31, 4.21, 4.6),
pch = 42, col = 5:8)
points(c(5.81, 6.05, 6.41, 6.13), c(4.53, 4.01, 4.18, 4.66),
pch = 42, col = 5:8)

Bandeira do Brasil

quinta-feira, 11 de outubro de 2007

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
plot(1:10, 1:10, type = "n", ann = FALSE, axes = FALSE, main = "Bandeira do BRASIL")
rect(1, 2, 10, 9, col = "darkgreen")
polygon(c(5.5, 1.2, 5.5, 9.8), c(2.2, 5.5, 8.8, 5.5), col = "yellow")
symbols(5.5, 5.5, 2, bg = rgb(0, 0, 1), add = TRUE)
rect(3.54, 5.3, 7.46, 5.8, col = "white")
text(5.5, 5.5, "ORDEM E PROGRESSO", col = "darkgreen")
points(5.92, 5.98, , pch = 8, col = "white")
points(c(6.06, 5.91, 6.07, 6.13, 6.24, 5.68, 5.29, 5.19,
5.03, 4.87, 5.13, 5.18, 5.38), c(5.14, 4.82, 4.67, 4.99,
5.1, 4.9, 4.55, 4.86, 5.06, 5.06, 4.69, 4.47, 4.31),
pch = 8, col = "white")
points(c(5.73, 5.88, 5.51, 5.79, 5.96, 6.08, 5.34, 5.5, 5.23,
5.12, 5.55, 5.87, 6.04), c(5.01, 4.58, 4.66, 4.84, 5.1,
4.71, 4.88, 4.56, 4.55, 4.95, 4.42, 4.49, 4.64), pch = 8,
col = "white")

Mapa do Brasil

quarta-feira, 10 de outubro de 2007

Essa foi a 4ª questão da 3ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
require(sp)
require(maptools)
require(maps)

map("world","brazil",fill=TRUE,col="gray")
#map.axes()
box()
title("Mapa do Brasil")
text(-49,-24.6,"Curitiba",cex=.7)
points(-49.16,-25.25,pch=17,cex=1)

Gráfico de Bolhas

terça-feira, 9 de outubro de 2007

Essa foi a 3ª questão da 3ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
set.seed(123)
xc = round(runif(10), 2)
yc = round(runif(10), 2)
xy = cbind(xc, yc)

xy.sp = SpatialPoints(xy)
df <- data.frame(Bolhas = seq(10,50,length=10))
xy.spdf <- SpatialPointsDataFrame(xy.sp, df)

print(bubble(xy.spdf, "Bolhas", col="black",fill=FALSE, maxsize=2))

Curvas de Nível

segunda-feira, 8 de outubro de 2007

Essa foi a 2ª questão da 3ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.






CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
z.f<-function(x,y){y<--15 + 1 * x + 0.5 * y}
z <- outer(x, y, z.f)
contour(x, y, z, xlab = "X", ylab = "Y")

Gráfico 3D

domingo, 7 de outubro de 2007

Essa foi a 1ª questão da 3ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
x <- seq(0, 1.5, length = 30)
y <- seq(0, 1.5, length = 30)

#z.f<-function(x,y){y<--15 - 1 * x - 0.5 * y + .5*x^2 + .5*y^2 -.3*x*y}
z.f<-function(x,y){y<--15 - 1 * x - 0.5 * y}
z <- outer(x, y, z.f)

persp(x,y,z, theta = -35, phi = 10, expand = 0.5, col = "gray",
xlab = "X", ylab = "Y", zlab = "Resposta",scale = T, ticktype = "detailed",xlim=c(0,1.5))

Box Plot

sábado, 6 de outubro de 2007

Essa foi a 2ª questão da 2ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
set.seed(123)
x3<-rnorm(50) boxplot(x3, boxwex = 0.15, ylab = "Resposta") rug(x3, side = 2)

Distribuição Exponencial

sexta-feira, 5 de outubro de 2007

Essa foi a 1ª questão da 2ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
curve(exp(x),from=10, to=15)

Curva da Distribuição Normal

quinta-feira, 4 de outubro de 2007

Essa foi a 2ª questão da 1ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.




CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
a <- -4
b <- -1
c <- 2
d <- 4
a.n <- "-4"
b.n <- "-1"
c.n <- "2"
d.n <- "4"
plot(dnorm, -4, 4, xlab = "", ylab = "", axes = F, bty = "n")
polygon(x = c(a, seq(a, b, l = 20), b), y = c(0, dnorm(seq(a,b, l = 20)), 0), col = "gray")
text(c(a, b), c(-0.011, -0.011), c(a.n, b.n), cex = 0.5)
polygon(x = c(c, seq(c, d, l = 20), d), y = c(0, dnorm(seq(c,d, l = 20)), 0), col = "gray")
text(c(c, d), c(-0.011, -0.011), c(c.n, d.n), cex = 0.5)
segments(-4, -0.004, 4, -0.004)

Triângulo, Círculo, Fórmula e Texto

quarta-feira, 3 de outubro de 2007

Essa foi a 1ª questão da 1ª prova, aplicada pelo Prof Adilson, no segundo semestre de 2006.
Gabarito.

CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
plot(1:10,1:10,ann=FALSE,type='n')
polygon(c(1.5,2.5,2),c(8,8,9))
text(8,3,"Este ou Aquele?",srt=-45)
text(8,8,expression(hat(y)==hat(beta)[0]+hat(beta)[1]*x),cex=1.5)
symbols(2,2,, circles = 0.8, inches = FALSE, add = TRUE)

Pictograma

terça-feira, 2 de outubro de 2007

Esse foi um trabalho realizado pelos alunos Admilson Buzinaro e Marcos Tadeu Andrade Cordeiro, para a disciplina de Análise Gráfica, ministrada pelo Profº Adilson dos Anjos, no segundo semestre de 2006. Foi tão criativo que o professor disponibilizou o código no site.


EXEMPLO DE PICTOGRAMA

Em meados de 1969, uma região situada ao sul do estado do Texas, no Estados Unidos, tinha sua economia baseada quase que exclusivamente no comércio varejista e no turismo. No entanto, em 1970, descobriu-se que nesta região havia petróleo. Estudos geológicos indicavam que ali, encontrava-se a maior reserva petrolífera já encontrada em solo norte-americano. Desde então, várias empresas começaram a explorar este petróleo. O gráfico abaixo, mostra a evolução da produção, entre os anos de 1970 e 1976 (em milhares de barris).

OBS: As linha verde, indica a produção máxima, e a linha vermelha, a produção mínima.



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
set.seed(123) #semente utilizada
(ano=seq(1970,1975,1))
(prod=rbinom(6,50,.4))

pictograma.petroleo = function(ano,prod){
tit="Produção de Petróleo entre 1970 e 1975"
plot(ano,prod,las=1,cex.main=1,main=tit,bty="u",xlab="Ano",cex.axis=.65,ylim=c(0,max(prod)+8),xlim=c(min(ano)-.5,max(ano)+.5),type="h",ylab="Produtividade \n (em milhares de barris)",cex=4,col="black",cex.lab=.7)
for(i in 1:length(ano)){
segments(ano[i],0,ano[i],prod[i],lwd=2)


v=prod[i]/150

polygon(c(ano[i]-v,ano[i],ano[i]+
v,ano[i]+.25*v,ano[i]-.25*v),c(0,40*v,0,prod[i],prod[i]),density=10,angle=45)
polygon(c(ano[i]-v,ano[i],ano[i]+
v,ano[i]+.25*v,ano[i]-.25*v),c(0,40*v,0,prod[i],prod[i]),density=10,angle=-45)
polygon(c(ano[i]-v,ano[i],ano[i]+
v,ano[i]+.25*v,ano[i]-.25*v),c(0,40*v,0,prod[i],prod[i]),density=10,angle=0)

s=.7 # tamanho
d=.95 # localização
lines(c(ano[i]-v*s,ano[i]+s*v),rep(prod[i]*d,2))
polygon(c(ano[i]-v*s,ano[i]+s*v,ano[i]+v*s,ano[i]-s*v),c(rep(prod[i]*(d-.02),2),rep(prod[i]*(d+.02),2)))

s=1 # tamanho
d=.6 # localização
lines(c(ano[i]-v*s,ano[i]+s*v),rep(prod[i]*d,2))
polygon(c(ano[i]-v*s,ano[i]+s*v,ano[i]+v*s,ano[i]-s*v),c(rep(prod[i]*(d-.02),2),rep(prod[i]*(d+.02),2)))

points(ano[i]-v,0,pch=15,cex=1)
points(ano[i]+v,0,pch=15,cex=1)
}
segments(min(ano)-5,min(prod),max(ano)+5,min(prod),lty=4,col="red",lwd=1.5)
segments(min(ano)-5,max(prod),max(ano)+5,max(prod),lty=4,col="darkgreen",lwd=1.5)
segments(min(ano)-5,mean(prod),max(ano)+5,mean(prod),lty=2,col="darkblue",lwd=1)
}

pictograma.petroleo(ano,prod)

legend("topright",title=c("Produção"),c("Máxima", "Mínima","Média"), lty =
c(4,4,2),col=c("darkgreen","red","darkblue"),cex=.7)

Software R

segunda-feira, 1 de outubro de 2007

O programa computacional R é um software gratuito e código aberto que permite manipular dados, fazer cálculos e gráficos estatísticos. A criação desse novo marcador tem como objetivo divulgar sua utilização.




Para ter o R instalado em seu computador é só fazer o download: R-Project

==> para os que tem Windons é só escolher base e depois a versão do programa.





Caso deseje aprofundar seus conhecimentos nesse software, aconselho entrar para esse grupo de discussão: R-STAT

==> Esta lista de discussão tem como objetivo a troca de informações, perguntas e respostas sobre o uso da linguagem R.





Apostilas sobre a utilização do R em português:



Introdução ao sistema estatístico R - Mini-curso EMBRAPA
Profº. Dr. Paulo Justiniano Ribeiro Jr (e-mail) PDF HTML

Análise Gráfica
Profº. Me. Adilson dos Anjos (e-mail) PDF HTML

==> a página do professor Adilson tem senha, dessa forma havendo interesse no material, sugiro solicitar para ele por e-mail.