Achei a idéia muito criativa! Seria possivel incentivar os alunos com os melhores scripts postar as análises no RPubs e publicar os links na lista para podermos aprender um pouco mais? Conheço colegas (eu inclusive) que também gostariam de ver os possíveis resultados Walmes

Abs

André Rutz

“We need to follow the data just like the yellow brick road, to the truth.” Samuel Klein

Em 04/10/2013 15:10, "walmes ." <walmeszeviani@gmail.com> escreveu:
Jackson,

Excelente. Com suas dicas eu consegui tirar os dados que eu precisava. Abaixo segue um CMR, para quem se sentir curioso, que lê as tabelas por meio da url correspondente aos dias da maratona e no final calcula algumas estatísticas canditadas à relacionar com o desempenho dos atletas.

Grato.
Walmes.

#-----------------------------------------------------------------------------
# Lê tabelas de clima nos dias da maratona de Paris.

ini <- "http://www.wunderground.com/history/airport/LFPB/"
fim <- "/DailyHistory.html??format=1&format=1"
dia <- c("2012/10/14")

# padrão do endereço das tabelas
url <- paste(ini, dia, fim, sep="")

da <- read.table(url, sep=",", header=TRUE,
                 na.string=c("N/A","-","","-9999"))
str(da)

#----------------------------------------------------------------------------
## Datas do evento.

## http://www.20kmparis.com/web/affiches.asp

dias <- c("2012/10/14","2011/10/9","2010/10/10",
          "2009/10/11","2008/10/12","2007/10/17",
          "2006/10/15","2005/10/16","2004/10/17")

lercada <- function(dia){
    ini <- "http://www.wunderground.com/history/airport/LFPB/"
    fim <- "/DailyHistory.html??format=1&format=1"
    url <- paste(ini, dia, fim, sep="")
    da <- read.table(url, sep=",", header=TRUE,
                     na.string=c("N/A","-","","-9999"))
    cat(dia, "\t ok!\n")
    return(da)
}

L <- lapply(dias, lercada)
sapply(L, class)
sapply(L, dim)
sapply(L, names)

names(L) <- 2012:2004

#-----------------------------------------------------------------------------

require(plyr)

M <- ldply(L)
str(M)
names(M)[1] <- "ano"

aggregate(cbind(TemperatureC,Humidity)~ano, data=M, max)
aggregate(cbind(TemperatureC,Humidity)~ano, data=M, mean)
aggregate(cbind(TemperatureC,Humidity)~ano, data=M, median)

#-----------------------------------------------------------------------------

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
skype: walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.