#========================================================================================== # Estatística Computacional I - ce083-2012-01 # Aula 6 - Exercícios de fixação sobre seleção e agrupamento em classes (21/03/2012) # Prof. Walmes Zeviani - Dpto de Estatística - UFPR # #========================================================================================== #------------------------------------------------------------------------------------------ # exercício 16 str(rock) # estrutura do conjunto de dados subset(rock, subset=area>7500) # mostra valores de area>7500 rock[rock$area>7500,] # o mesmo usando o operador [] subset(rock, subset=area>6000 & area<8000) rock[rock$area>6000 & rock$area<8000,] rock$ap <- rock$area/rock$peri # cria uma nova coluna com a divisão de area por peri str(rock) rock <- subset(rock, select=c(5,1:4)) # altera a ordem das colunas no data.frame str(rock) #------------------------------------------------------------------------------------------ # exercício 17 rock2 <- subset(rock, select=c("area","peri")) # seleção pelo nome da coluna rock2 <- subset(rock, select=2:3) # seleção pelo índice da coluna rock3 <- subset(rock, subset=perm<800) rock3 rock4 <- data.frame("rarea"=sqrt(rock$area), "lperm"=log(rock$perm)) str(rock4) summary(rock) # summary retorna o extremos (min, max), os quartis (1,2,3 quartil) e a média #------------------------------------------------------------------------------------------ # exercício 18 help(cut, help_type="html") rock$carea <- cut(rock$area, seq(1000, 13000, 2000), dig.lab=5) rock$cperi <- cut(rock$peri, seq(from=0, by=1000, length.out=6)) min(rock$shape) # mínimo max(rock$shape) # máximo range(rock$shape) # mínimo e máximo length(rock$shape) # tamanho do vetor # algumas possíveis escolhas para os limites das classes seq(0.05, 0.5, 0.05) seq(0.0, 0.5, 0.1) seq(0.08, length.out=6, to=0.48) # agrupa os valores em classes rock$cshape <- cut(rock$shape, seq(0.08, length.out=6, to=0.48)) str(rock4) range(rock4$lperm) seq(1.5, 7.5, 1) rock$clperm <- cut(rock4$lperm, seq(1.5, 7.5, 1)) #------------------------------------------------------------------------------------------ # exercício 19 str(rock) table(rock$carea) # tabela de distribuição de frequências absoluta table(rock$carea)/nrow(rock) # tabela de distribuição de frequências relativas plot(rock$carea) # gráfico de barras das frequências das classes table(rock$cperi) table(rock$cperi)/nrow(rock) plot(rock$cperi) table(rock$cshape) table(rock$cshape)/nrow(rock) plot(rock$cshape) table(rock$clperm) table(rock$clperm)/nrow(rock) plot(rock$clperm) #------------------------------------------------------------------------------------------ # exercício 20 hist(rock$area) # histograma hist(rock$peri) hist(rock$shape) hist(rock$perm) hist(log(rock$perm)) pie(table(rock$carea)) # gráfico de setores pie(table(rock$cperi)) pie(table(rock$cshape)) pie(table(rock$clperm)) #------------------------------------------------------------------------------------------