Outubro 2012 - R-br

correlação de Pearson
by marcelo claro de souza 03 Out '12

03 Out '12

Boa tarde pessoal Montei uma matriz de correlação utilizando cor(dados, method='pearson') Gostaria de saber se existe alguma função que teste automaticamente o coeficiente de correlação entre as variáveis sem precisar rodar manualmente cor.test(x,y,method='pearson'), afinal são quase 30 variáveis que estou correlacionando. Muito obrigado pela ajuda. Abraço

5 4

obtenção de lat long para uma série de CEPs
by Mauricio Cardeal 03 Out '12

03 Out '12

Recentemente foi postado na lista uma solução para a obtenção de lat long a partir de um cep: require(XML) coordenadas<- function(cep) { url_lat_lon <- paste(sprintf(" http://maps.google.com/maps/api/geocode/xml?address=%s,", cep),"%20Brasil&sensor=false", sep="") lat_lon=xmlApply(xmlRoot(xmlTreeParse( readLines(url_lat_lon)))[['result']][['geometry']][['location']], "[[", 1) return(lat_lon) } cep=40110010 coordenadas(cep) > coordenadas(cep) $lat -12.9896162 $lng -38.4922117 Ficou excelente. Mas fiquei imaginando se ao invés de apenas 1 cep eu tivesse uma série, por exemplo: ceps <- c(40110010,01317000,17120000) Como modificaria a função acima rotina para ter algo do tipo (em um data frame), sem precisar rodar 1 a 1 cada cep: cep lat long 40110010 -12.98 -38.49 01317000 -25.29 -57.62 17120000 -22.52 -49.13 Valeu -- Mauricio Cardeal UFBA

2 2

Organizacao de dados Netcdf
by Humberto Hazin 03 Out '12

03 Out '12

Bom dia pessoal, Estou tentando extrair variáveis ambientais de arquivos formato netcdf. Até um certo momento consigo separtar em arquivos as Latitudes, Longitudes, time (mes/ano) e a dml (exemplo que mando aqui). Link para baixar os dados do exemplo https://www.dropbox.com/s/66dp4rgei791r52/DML05_12.nc nc <- open.nc("DML05_12.nc", write=TRUE) print.nc(nc) lat<-var.get.nc(nc, "LAT46_103") lon<-var.get.nc(nc, "LONN55_12") time<-var.get.nc(nc, "TIME") dml<-var.get.nc(nc, "MLD") Meu problema é que nao estou conseguindo organizar eles em colunas com seus respectivos valores para cada latitude, longitude, time e dml. o output final deveria ficar assim: lat lon time dml 1 -44.5 -55.5 2 -43.5 -55.5 3 -42.5 -55.5 4 -41.5 -55.5 5 -40.5 -55.5 6 -39.5 -55.5 Agradeço a ajuda desde já Humberto __________________________________________________ Prof. Dr. Humberto Gomes Hazin Universidade Federal Rural de Pernambuco - UFRPE Unidade Acadêmica de Serra Talhada - UAST Coordenador do Curso de Engenharia de Pesca Fazenda Saco, S/N - Caixa Postal 063 CEP 56900-000 Serra Talhada - PE Tel.(87) 3831-2206/3831-1927 Cel.(81) 9158 0564

2 3

Sweave absoleto?
by Andre Oliveira 03 Out '12

03 Out '12

Pessoas, boa tarde. Aqui em meu notbook, com ubuntu 12.04, uso R version 2.15.1 (2012-06-22) e o Sweave ainda está compilando normalmente. Mas em outras listras já dizem que o Sweave ficou absoleto. Isto procede? Alguma função nova com mesmas funções do Sweave? obrigado André Oliveira Souza André Oliveira Souza

3 2

Eixo fracionado
by Aira Nava 03 Out '12

03 Out '12

Boa tarde, Acredito que meu problema seja simples, mas tentei solucioná-lo de várias formas e não consegui. Estou começando a utilizar o R para análise de séries temporais, porém ao plotar os gráficos ts o eixo das abscissas aparece fracionado em semestres, por exemplo: 2008 e 2008.5 ou 2009 e 2009.5. Gostaria de saber se existe uma forma de modificar este formato para que as datas mostradas sejam correspondentes ao dia de coleta (31/mar/2008) ou que não sejam fracionas, mas apenas 2008, 2009 e 2010. Pta=read.table("p_ta.cvs", h=F, sep=";", dec=",") Pta Pta=ts(Pta[,2],st=c(2008,3),freq=12) plot.ts(Pta,yax.flip=12) #ou apenas plot(Pta) ou plot.ts(Pta) Obrigada pelo apoio, Aira Nava - UNESP Ourinhos.

2 1

ggplot2
by Alisson Lucrecio 03 Out '12

03 Out '12

Olá pessoal da r-br, Eu usava as funções ggplotGrob, grid.ls e getGrob para extrair os elementos do grid e renomear-los da forma que queria, mas o pacote ggplot2 teve uma alteração. Algum sabe como fazer as alterações do strip.elementes com o ggplot2? obrigado. Alisson Lucrécio da Costa

3 4

substituir valores nulos
by Samuel Carvalho 02 Out '12

02 Out '12

Caros (as) membros da r-br Seria possível substituir todos valores nulos ("NA") por algum valor específico? Como por exemplo pelo valor 1 ou somente dizendo se NA entao vazio? Exemplo de #RCMR #RCMR dados <- data.frame(parcela=rep(1:3, each=3), y=rnorm(9)) dados$substituir <- c(1,2,3,4,5,"NA",7, "NA",9) dados[2,1] <- "NA" dados[5,1] <- "NA" dados dados[which(dados=="NA")] <- 0 #nao rodou, sei que se especifico a coluna este comando funciona, # porem se tenho um data frame com varias colunas é um problema # ou para substituir os NA por vazios dados[which(dados=="NA")] <- "" Desde ja obrigado ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

3 3

(sem assunto)
by Edson Lira 02 Out '12

02 Out '12

Boa tarde caros amigos, com o comando table, tenho a tabela abaixo: table(ago$Grupo,ago$AvaliaTempo) Bom Otimo Regular Ruim Ambulatorio 0 14 9 32 35 Doador Ana Braga 0 58 49 5 3 Doador Coleta Externa 0 39 36 6 2 Doador Hemoam 0 61 48 14 10 Enfermaria 0 7 2 13 14 Laboratorio 0 61 26 25 20 Qual a rotina para ordenar a tabela pela avaliação na seguinte ordem: Otimo bom regular ruim Ambulatorio 9 14 32 35 Doador Ana Braga Doador Coleta Externa Doador Hemoam Enfermaria Laboratorio { ]'s. Edson Lira Estatístico Manaus-Amazonas

3 2

Re: [R-br] Mostrar valores no boxplot
by Walmes Zeviani 02 Out '12

02 Out '12

Edson, Submeta sua dúvida para a lista pois mais pessoas podem aproveitar-se da mesma, mais pessoas podem colaborar com você, você pode ser respondido em menos tempo, você pode contemplar diferentes soluções para o seu problema e sua mensagem fica no histórico da lista para consulta futura. A mediana já está representada pelo traço que divide a caixa no boxplot. Para adicionar as médias, após o boxplot() aplique um points() por exemplo. Envie um CMR. À disposição. Walmes. ========================================================================== 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/10/2 Edson Lira <edinhoestat(a)yahoo.com.br> > Valeu Walmes! Esqueci o assunto. > > Quando monto o boxplot qual a rotina para mostrar no gráfico por exemplo, > o valor da mediana e/ou média. > > boxplot(x) > > Edson Lira > Estatístico > Manaus-Amazonas > ------------------------------ > *De:* Walmes Zeviani <walmeszeviani(a)gmail.com> > *Para:* r-br(a)listas.c3sl.ufpr.br; Edson Lira <edinhoestat(a)yahoo.com.br> > *Enviadas:* Terça-feira, 2 de Outubro de 2012 11:46 > *Assunto:* Re: [R-br] (sem assunto) > > Submeta mensagens à lista sempre contendo título apropriado e envie CMR. > Seu problema é resolvido mexendo o índice das colunas, algo como > > x <- table(...) > x[,c(4,3,1,2)] > > À disposição. > Walmes. > > ========================================================================== > 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 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > >

1 0

Re: [R-br] correlação de Pearson
by Ivan Bezerra Allaman 02 Out '12

02 Out '12

cor.prob <- function(X, dfr=nrow(X)-2){ R <- cor(X,use='na.or.complete') above <- row(R) < col(R) r2 <- R[above]^2 Fstat <- r2*dfr/(1-r2) R[above] <- 1-pf(Fstat, 1, dfr) R } \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0