Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior Próxima revisão Ambos lados da revisão seguinte | ||
disciplinas:ce223:comandos2008 [2008/03/13 23:45] paulojus |
disciplinas:ce223:comandos2008 [2008/03/17 16:51] ehlers |
||
---|---|---|---|
Linha 270: | Linha 270: | ||
Digitação e conversão de uma tabela de tripla entrada (dada no quadro durante a aula) em um objeto do tipo ''array'' | Digitação e conversão de uma tabela de tripla entrada (dada no quadro durante a aula) em um objeto do tipo ''array'' | ||
- | <code R> | + | | | PR || SC || RS || |
+ | | | Masculino | Feminino | Masculino | Feminino | Masculino | Feminino | | ||
+ | |Não Fuma ^ 45 ^ 16 ^ 21 ^ 33 ^ 40 ^ 45 ^ | ||
+ | |Fuma pouco ^ 28 ^ 22 ^ 34 ^ 21 ^ 50 ^ 37 ^ | ||
+ | |Fuma muito ^ 37 ^ 15 ^ 56 ^ 30 ^ 85 ^ 29 ^ | ||
+ | |||
+ | Comentários sobre ordem de entrada dos dados, cliclagem das variáveis e definição das dimensões do array<code R> | ||
freqs <- scan() | freqs <- scan() | ||
+ | 1: 45 | ||
+ | 2: 28 | ||
+ | 3: 37 | ||
+ | 4: 16 | ||
+ | 5: 22 | ||
+ | 6: 15 | ||
+ | 7: 21 | ||
+ | 8: 34 | ||
+ | 9: 56 | ||
+ | 10: 33 | ||
+ | 11: 21 | ||
+ | 12: 30 | ||
+ | 13: 40 | ||
+ | 14: 50 | ||
+ | 15: 85 | ||
+ | 16: 45 | ||
+ | 17: 37 | ||
+ | 18: 29 | ||
+ | 19: | ||
freqs | freqs | ||
Af <- array(freqs, dim=c(3,2,3)) | Af <- array(freqs, dim=c(3,2,3)) | ||
Linha 277: | Linha 302: | ||
</code> | </code> | ||
+ | ==== Semana 4 ==== | ||
+ | === 17/03/2008 === | ||
+ | <code R> | ||
+ | freqs = scan(file='http://leg.ufpr.br/~ehlers/CE223/fumo.dat') | ||
+ | |||
+ | array(freqs, dim=c(3,2,3)) | ||
+ | |||
+ | nomes = list(c('PR','SC','RS'), c('M','F'), c('nao fuma','fuma pouco','fuma muito')) | ||
+ | |||
+ | hf = array(freqs, dim=c(3,2,3), dimnames=nomes) | ||
+ | |||
+ | hf | ||
+ | |||
+ | m1 <- matrix(1:12, ncol = 3) | ||
+ | m1 | ||
+ | |||
+ | dimnames(m1) | ||
+ | |||
+ | dimnames(m1) <- list(c("L1", "L2", "L3", "L4"), c("C1", "C2", "C3")) | ||
+ | |||
+ | m1 | ||
+ | |||
+ | m2 <- cbind(1:5, 6:10) | ||
+ | m2 | ||
+ | |||
+ | m3 <- cbind(1:5, 6) | ||
+ | m3 | ||
+ | |||
+ | margin.table(m1, margin = 1) | ||
+ | |||
+ | apply(m1, 1, sum) | ||
+ | |||
+ | rowSums(m1) | ||
+ | |||
+ | margin.table(m1, margin = 2) | ||
+ | |||
+ | apply(m1, 2, sum) | ||
+ | |||
+ | colSums(m1) | ||
+ | |||
+ | #operacoes com matrizes | ||
+ | |||
+ | m4 <- matrix(1:6, nc = 3) | ||
+ | m5 <- matrix(10 * (1:6), nc = 3) | ||
+ | m4 | ||
+ | |||
+ | m5 | ||
+ | |||
+ | m4 + m5 | ||
+ | |||
+ | m4 * m5 | ||
+ | |||
+ | m5 - m4 | ||
+ | |||
+ | m5/m4 | ||
+ | |||
+ | m4 %*% m5 | ||
+ | |||
+ | t(m4) | ||
+ | |||
+ | m6 = t(m4)%*% m5 | ||
+ | |||
+ | solve(m6) | ||
+ | |||
+ | m6[3,3]=20 | ||
+ | |||
+ | solve(m6) | ||
+ | |||
+ | </code> |