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/04/09 10:21] ehlers |
disciplinas:ce223:comandos2008 [2008/04/23 09:45] ehlers |
||
---|---|---|---|
Linha 306: | Linha 306: | ||
<code R> | <code R> | ||
freqs = scan(file='http://leg.ufpr.br/~ehlers/CE223/fumo.dat') | freqs = scan(file='http://leg.ufpr.br/~ehlers/CE223/fumo.dat') | ||
+ | freqs | ||
+ | [1] 45 16 21 33 40 45 28 22 34 21 50 37 37 15 56 30 85 29 | ||
- | array(freqs, dim=c(3,2,3)) | + | array(freqs, dim=c(2,3,3)) |
+ | , , 1 | ||
- | nomes = list(c('PR','SC','RS'), c('M','F'), c('nao fuma','fuma pouco','fuma muito')) | + | [,1] [,2] [,3] |
+ | [1,] 45 21 40 | ||
+ | [2,] 16 33 45 | ||
- | hf = array(freqs, dim=c(3,2,3), dimnames=nomes) | + | , , 2 |
+ | |||
+ | [,1] [,2] [,3] | ||
+ | [1,] 28 34 50 | ||
+ | [2,] 22 21 37 | ||
+ | |||
+ | , , 3 | ||
+ | |||
+ | [,1] [,2] [,3] | ||
+ | [1,] 37 56 85 | ||
+ | [2,] 15 30 29 | ||
+ | |||
+ | # Cada matrix 2x3 contem as contagens por sexo (linhas) e estado (colunas). | ||
+ | # A ultima dimensao refere-se ao habito de fumar. | ||
+ | |||
+ | nomes = list(c('M','F'),c('PR','SC','RS'),c('nao fuma','fuma pouco','fuma muito')) | ||
+ | |||
+ | hf = array(freqs, dim=c(2,3,3), dimnames=nomes) | ||
hf | hf | ||
+ | , , nao fuma | ||
+ | |||
+ | PR SC RS | ||
+ | M 45 21 40 | ||
+ | F 16 33 45 | ||
+ | |||
+ | , , fuma pouco | ||
+ | |||
+ | PR SC RS | ||
+ | M 28 34 50 | ||
+ | F 22 21 37 | ||
+ | |||
+ | , , fuma muito | ||
+ | |||
+ | PR SC RS | ||
+ | M 37 56 85 | ||
+ | F 15 30 29 | ||
m1 <- matrix(1:12, ncol = 3) | m1 <- matrix(1:12, ncol = 3) | ||
Linha 730: | Linha 769: | ||
25% 75% | 25% 75% | ||
30.66667 40.52083 | 30.66667 40.52083 | ||
+ | </code> | ||
Medidas robustas | Medidas robustas | ||
Linha 857: | Linha 897: | ||
</code> | </code> | ||
+ | |||
+ | |||
+ | ==== Semana 7 ==== | ||
+ | |||
+ | === 07/04/2008 e 09/04/2008 === | ||
+ | |||
+ | Analisar os dados do Exercicio 26, Capitulo 1 do livro NOÇÕES DE PROBABILIDADE E ESTATÍSTICA disponiveis em http://www.ime.usp.br/~noproest | ||
+ | |||
+ | Note que ha brancos no arquivo de dados (dados omissos). Uma forma de tratar este problema é abrir o arquivo Excel e salvar como um arquivo texto do tipo CSV (comma separated values). Posteriormente este arquivo pode ser lido como | ||
+ | <code R> | ||
+ | |||
+ | read.table('nome do arquivo', header=T, sep=',') | ||
+ | |||
+ | # ou | ||
+ | |||
+ | read.csv('nome do arquivo', header=T) | ||
+ | </code> | ||
+ | |||
+ | Uma alternativa melhor é utilizar a função read.xls do pacote gdata pois assim não precisamos abrir o arquivo Excel. Após salvar o arquivo aeusp.xls na sua area de trabalho execute | ||
+ | <code R> | ||
+ | library(gdata) ou require(gdata) | ||
+ | |||
+ | x = read.xls ('aeusp.xls') | ||
+ | |||
+ | head(x) | ||
+ | |||
+ | Num Comun Sexo Idade Ecivil X.Reproce X.Temposp X.Resid Trab Ttrab X.Itrab | ||
+ | 1 1 JdRaposo 2 4 4 Nordeste 21 9 3 NA 20 | ||
+ | 2 2 JdRaposo 2 1 1 Sudeste 24 9 1 1 14 | ||
+ | 3 3 JdRaposo 2 2 1 Nordeste 31 3 1 1 14 | ||
+ | 4 4 JdRaposo 1 2 2 Nordeste 10 3 1 4 10 | ||
+ | 5 5 JdRaposo 2 4 2 Nordeste 31 6 1 1 11 | ||
+ | 6 6 JdRaposo 2 4 2 Sudeste 24 4 2 NA 15 | ||
+ | X.Renda X.Acompu X.Serief | ||
+ | 1 1 2 1 | ||
+ | 2 2 2 7 | ||
+ | 3 5 2 7 | ||
+ | 4 5 2 11 | ||
+ | 5 6 1 4 | ||
+ | 6 4 2 4 | ||
+ | </code> | ||
+ | |||