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 Ambos lados da revisão seguinte | ||
disciplinas:ce223:comandos2008 [2008/04/23 15:41] ehlers |
disciplinas:ce223:comandos2008 [2008/05/05 16:02] paulojus |
||
---|---|---|---|
Linha 519: | Linha 519: | ||
by(d4[,4],d4$sexo,function(x)as.character(x)) | by(d4[,4],d4$sexo,function(x)as.character(x)) | ||
</code> | </code> | ||
+ | |||
Listas | Listas | ||
Linha 937: | Linha 938: | ||
5 6 1 4 | 5 6 1 4 | ||
6 4 2 4 | 6 4 2 4 | ||
+ | </code> | ||
+ | |||
+ | ==== Semana 1 ==== | ||
+ | === 30/04/2008 === | ||
+ | |||
+ | Gerando 1000 amostras de tamanho n=20 de uma ormal padrão | ||
+ | <code R> | ||
+ | rnorm(20, m=70, sd=10) | ||
+ | ams <- matrix(rnorm(20*1000, m=70, sd=10), ncol=20) | ||
+ | dim(ams) | ||
+ | ams[1,] | ||
+ | ams[2,] | ||
+ | </code> | ||
+ | |||
+ | Calculando o valor da estatística de interesse para a primeira e segunda amostra | ||
+ | <code R> | ||
+ | max(ams[1,])/quantile(ams[1,], prob=0.75) | ||
+ | unname(max(ams[1,])/quantile(ams[1,], prob=0.75)) | ||
+ | unname(max(ams[2,])/quantile(ams[2,], prob=0.75)) | ||
+ | </code> | ||
+ | |||
+ | Escrevendo uma função que calcula o valor da estatística de interesse e calculando ovamente o valor para a primeira e segunda amostras. | ||
+ | </code R> | ||
+ | T.est <- function(x) unname(max(x)/quantile(x, prob=0.75)) | ||
+ | T.est(ams[1,]) | ||
+ | T.est(ams[2,]) | ||
+ | </code> | ||
+ | |||
+ | Calculando valor da estatística de interesse agora para todas as amostras de uma só vez | ||
+ | <code R> | ||
+ | ts <- apply(ams, 1, T.est) | ||
+ | length(ts) | ||
+ | ts | ||
+ | </code> | ||
+ | |||
+ | Explorando os resultados: medidas resumo, grafico de densidade estimada e IC (95%) | ||
+ | <code R> | ||
+ | summary(ts) | ||
+ | plot(density(ts)) | ||
+ | quantile(ts, prob=c(0.025, 0.975)) | ||
+ | </code> | ||
+ | |||
+ | Aumentando o número de amostras para 5000. | ||
+ | <code R> | ||
+ | ams <- matrix(rnorm(20*5000, m=70, sd=10), ncol=20) | ||
+ | ts <- apply(ams, 1, T.est) | ||
+ | plot(density(ts)) | ||
+ | medias <- apply(ams, 1, mean) | ||
+ | plot(density(medias)) | ||
+ | curve(dnorm, 60, 80, mean=70, sd=10/sqrt(20), col=2) | ||
</code> | </code> | ||