====== Aula semana 4 ======

==== Comandos do R digitados na aula em 21/03/2007 ====

<code R>
a1 <- -10:10
a1
a2 <- seq(-10,10,by=1)
a2
a3 <- seq(-10,10,len=21)
a3
a4 <- rep(c(5,7,8), 3)
a4
a5 <- rep(c(7,9,11), c(2,3,4))
a5
a6 <- rep(c(5,7,8), each=3)
a6
a7 <- rep(seq(1:10, by=2), each=5)
a7
a8 <- rep(seq(1:10, l=6), 4)
a8 
a9 <- rep(paste("grupo", 1:5, sep="-"), each=4)
a9

a10 <- rep(paste(c("Tarde", "Noite"), 1:6, sep=""), each=4)
a10

## "Semana 3 Quarta","Semana 3 Sexta","Semana 4 Quarta","Semana 4 Sexta", 
rep("Semana", 4)
rep(3:4, each=2)
rep(c("quarta", "sexta"), 2)

paste(rep("Semana", 4), rep(3:4, each=2),rep(c("quarta", "sexta"), 2))
paste("Semana", rep(3:4, each=2),rep(c("quarta", "sexta"), 2))
paste("Semana", rep(3:4, each=2),c("quarta", "sexta"))

a12 <- rep(c("Azul", "Vermelho", "Branco"), 3)
a12
class(a12)
attributes(a12)

a13 <- factor(a12)
a13
class(a13)
attributes(a13)
a13[-2]
a13[-c(2,5,8)]

matrix(1:20, nc=5)
matrix(1:23, nc=5)
matrix(1:20, nc=6)
matrix(1:10, nc=5, nr=6)
matrix(1:10, nr=5, nc=6)
t(matrix(1:10, nr=5, nc=6))
matrix(1:10, nc=5, nr=6, byrow=TRUE)

m <- matrix(sample(1:20), nc=5)
m
m[3,2]
m[7]

m[,2]
m[,2:3]
m[c(1,3,4),c(2,4)]
m[-2,c(2,4)]
m[1:2,]
m[2,]
m[2,,drop=FALSE]

dim(m)
length(m)
m
apply(m,1,sum)
apply(m,2,sum)
apply(m,1,mean)
apply(m,2,mean)
</code>

==== Comandos do R digitados na aula em 21/03/2007 ====

<code R>
# dat <- scan()
#1: 58
#2: 37
#3: 43
#4: 25
#5: 39
#6: 26
#7: 57
#8: 38
#9: 
#Read 8 items
dat
is.vector(dat)
is.matrix(dat)
is.numeric(dat)
is.character(dat)
class(dat)

datm <- matrix(dat, nc=2)
datm

dim(datm)
dimnames(datm)
dimnames(datm) <-  list(c("Péssimo", "Regular", "Bom", "Ótimo"), c("Masc", "Fem"))
datm
dimnames(datm)
apply(datm,1,sum)
apply(datm,2,sum)
sum(datm[,1])
sum(datm[,2])

sc <- apply(datm,2,sum)
sc
sc[1]

#tab3 <- scan()
#1: 25
#2: 12
#3: 5
#4: 13
#5: 39
#6: 14
#7: 7
#8: 20
#9: 43
#10: 27
#11: 8
#12: 11
#13: 15
#14: 13
#15: 13
#16: 5
#17: 18
#18: 27
#19: 23
#20: 8
#21: 14
#22: 32
#23: 39
#24: 12

tab3a <- array(tab3, dim=c(4,3,2))
tab3a

## corrigindo um dado
tab3a[2,2,1]
tab3a[2,2,1] <- 15
tab3a
dimnames(tab3a) <-  list(c("Péssimo", "Regular", "Bom", "Ótimo"), c("Baixa","Média","Alta"),c("Masc", "Fem"))
tab3a
dim(tab3a)

## quantos homens?
sum(tab3a[,,1])

## quantas mulheres?
sum(tab3a[,,1])

## ou ....
apply(tab3a,3,sum)
apply(tab3a,3,sum)[1]
apply(tab3a,3,sum)[2]

## proporção de homens e de mulheres
apply(tab3a,3,sum)/sum(tab3a)

## quantos de baixa renda?
tab3a[,1,]
sum(tab3a[,1,])
apply(tab3a,2,sum)
apply(tab3a,2,sum)[1]
## quantos de opinião Bom?
sum(tab3a[3,,])
## número de homens com bom
tab3a[3,,1]
sum(tab3a[3,,1])
apply(tab3a,c(1,3), sum)
apply(tab3a,c(1,3), sum)[3,1]
## proporção de homens com bom
sum(tab3a[3,,1])/sum(tab3a)

## proporção de homens com bom ou ótimo
## proporção de homens com bom ou ótimo
## proporção de mulheres com bom ou ótimo
## proporção de ótimo para cada faixa de renda
## proporção de ótimo para cada combinaçõa de faixa de renda e sexo
## proporção de cada opinião para cada combinação de faixa de renda e sexo
</code>