Janeiro 2012 - R-br

Python para R
by Daniel Marcelino 04 Jun '12

04 Jun '12

Caros, alguem da lista conhece programação python? Estou há dias tentando desenvolver uma rotina em R para explorar dados eleitorais, mas não tenho conseguido ir muito longe. Eu tenho um modelo do que mais ou menos eu preciso, mas a rotina está em python. Alguem encara me ajudar nisso? Daniel ----------------------------------------------------- My contacts: Phone 1+514 343 6111 Ext 3799 Skype: d.marcelino 3200 Jean Brillant, Office C5075 Montreal, QC; H3T 1N8 Canada -----------------------------------------------------

3 2

WebScraping - Como passar informações ao site antes do readLines
by Paulo Nogueira 23 Fev '12

23 Fev '12

Boa tarde R-Br. Mais uma vez peço a vocês uma orientação de caminho para descobrir como solucionar problemas relacionados a web scraping. Hoje consigo ler códigos HTML de ṕáginas atráves do readLines e depois extrair os dados que me interessam utilizando funções como grep, gsub e regexpr. Evolui nisso e tenho conseguido resultados bastante interessantes. O problema é que nem sempre o que é exibido no navegador, está disponível no código fonte da página. Por exemplo, a página inicial do site Webmotors ( http://www.webmotors.com.br/index.html) apresenta um box de pesquisa. Escolhendo, na página, os campos marca = Chevrolet, modelo = Astra, estado = AM e cidade = Manaus e clicando em BUSCAR, encontro esse resultado: http://www.webmotors.com.br/Webmotors/Compra/carrosResultado/carros-resulta… E assim consigo extrair com o R as informações dos veículos Astra a venda em Manaus. Porém, como posso fazer para o R informar ao site que desejo fazer essa busca? Gostaria de automatizar essa parte. Sei que na URL resultado basta trocar os campos como &modelo=1042 ou marca=2, etc... mas também não sei quais são os possíveis valores para esses campos. Ressalto que esse problema é um exemplo. Gostaria de aprender a fazer chamada às funções (acredito que são as .js) da página para construir as páginas com os resultados e então importar com o readLines para guardar os dados de meu interesse. Espero ter sido claro no questionamento. Agradeço qualquer ajuda. Abraços, Paulo Nogueira Starzynski

3 8

[NEWS] Mudanças nas versões do R e melhora de performance
by Fernando Mayer 16 Fev '12

16 Fev '12

Não sei se todos já estão sabendo, mas li isso hoje e achei importante repassar para a lista: Sobre as mudanças de releases do R: R 2.14 to be released on October 31; R 2.13 patch on September 13 http://blog.revolutionanalytics.com/2011/09/r-214-to-be-released-on-october… Brian Ripley on The R Development Process http://blog.revolutionanalytics.com/2011/08/brian-ripley-on-the-r-developme… Sobre a melhora de performance: With byte compiler, R 2.14 will be even faster http://blog.revolutionanalytics.com/2011/08/with-byte-compiler-r-214-will-b… []s, --- Fernando Mayer Dept. Ecologia e Zoologia - ECZ/CCB Universidade Federal de Santa Catarina - UFSC URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com

4 4

Re: [R-br] arredondar soma para 100%
by Ivan Bezerra Allaman 06 Fev '12

06 Fev '12

Se você arredondar o seu objeto "porc" para 3 digitos, terás 20,455 para a letra a. Não sei porque cargas d'água o R arredonda para 20,45, pois aprendi na escolinha e segundo o wikipedia também, o correto seria arredondar este número para 20,46. Fazendo isso, o seu problema fica resolvido. (S,f,P) Allaman \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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

7 7

ajuste com nls
by Gilson Sanchez 31 Jan '12

31 Jan '12

Pessoal mais uma duvida, eu fiz um ajuste com o nls, mas eu quero saber como utilizar a e b ou seja se tem algum comando para que eu possa utilizar os valores que foram estimadas pela equação Segue o CMR que estou trabalhando # Planta 1 planta1=foto[1:90,-3:-6] names(planta1) [1] "especie" "esp" "PAR" "Tleaf" "W" "Pn" "E" [8] "C" "int_CO2" colnames(planta1)<-c("especiea", "espa", "PARa","Tleafa", "Wa", "Pna", "Ea", + "Ca", "int_CO2a") # Obtendo PAR planta1.luz =tapply(planta1$PARa, INDEX=planta1$PARa, FUN=mean) planta1.luz 0 25 50 75 100 500 1000 1500 2000 0 25 50 75 100 500 1000 1500 2000 planta1.me.Pna=tapply(planta1$Pna, INDEX=planta1$PARa, FUN=mean) planta1.me.Pna 0 25 50 75 100 500 1000 1500 2000 -0.428 0.144 0.411 1.275 2.052 4.349 5.082 4.078 3.794 planta1.ajuste=nls((planta1.me.Pna~(a+0.428)*(1-exp(-b*planta1.luz/(a+0.428)))-0.428), + start=list(a=200, b=0.5)) planta1.ajuste Nonlinear regression model model: planta1.me.Pna ~ (a + 0.428) * (1 - exp(-b * planta1.luz/(a + 0.428))) - 0.428 data: parent.frame() a b 4.37564 0.02943 residual sum-of-squares: 1.258 Number of iterations to convergence: 5 Achieved convergence tolerance: 3.175e-06 max=c(planta1.ajuste$a) = ????? alfa=c(plnat1.ajuste$b) = ??????? Desde já muito obrigado. -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 http://lattes.cnpq.br/3684553005529634 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *

3 4

(sem assunto)
by Paulo Ricardo Silva Moreira 31 Jan '12

31 Jan '12

Boa tarde Galera, Gostaria de saber como faço para postar várias séries no mesmo gráfico, e como faço para colocar uma cor diferente para cada série, de forma a diferenciá-las. Att, Paulo Ricardo.

3 2

Inserir legenda em Boxplot
by Marcelo Claro de Souza 31 Jan '12

31 Jan '12

Tentei inserir uma legenda dentro de um boxplot utilizando o comando abaixo e sempre recebo a informção de erro inesperado. legend("topleft","(x,y)",pch=1,tittle=" A",inset=0.02) Teria alguma outra opção? Abraço e muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

4 3

Dúvida sobre teste
by Katia Alves Campos 31 Jan '12

31 Jan '12

2 1

Importar Microdados - arquivos .txt sem separação entre as colunas
by Roney Fraga 31 Jan '12

31 Jan '12

Caros, Diante da necessidade de importar microdados, principalmente do IBGE, por meio do pacote IBGEPesq()* *que vem apresentando alguns problemas, pergunto se algum usuário tem alguma sugestão para importar microdados alternativa ao pacote IBGEPesq(). Microdados consiste geralmente em grandes pesquisas que são disponibilizadas em arquivos .txt apenas com números sem nenhuma separação entre as colunas, sendo o arquivo dicionário disponibilizado a parte em outro arquivo .txt, nesse arquivo dicionário consta em que coluna começa e termina cada variável. No link http://www.datafilehost.com/download-2565fe79.html segue exemplo do arquivo .txt com os dados da PNAD 2009 banco domicílios. E no link http://www.datafilehost.com/download-fe9c31a8.html segue exemplos de dicionários para o SAS, disponibilizado pelo IBGE, e para o Stata que fiz com base no dicionário anterior. Se necessário o mesmo dicionário pode ser encontrado no R através dos comandos install.packages('dicionariosIBGE')library(dicionariosIBGE) data(dicPNAD2009) dicdom2009 No Stata uma vez que tenha o dicionário, na extensão .dct, e o arquivo .txt com os dados a serem lidos basta executar o comando infile using “/Caminho/Arquivo/dicionario.dct” no R não sei se é possível seguir o mesmo raciocínio para importar dados em arquivos .txt, caso algum usuário possa contribuir será uma grande ajuda. Os arquivos em .txt geralmente são grandes, por exemplo, a PNAD 2009 arquivo pessoas tem 320MB. Uma vez ultrapassada a dificuldade de importar os dados, será mais fácil utilizar o R para trabalhar com os dados públicos das seguintes pesquisas: PNAD, POF, PME, Enade, Censo Escolar, Censo da Educação Superior, Censo dos Profissionais do Magistério, Saeb, Enem, Provão, PNERA, Prova Brasil, Pesquisa de Ações Discriminatórias no Âmbito Escolar, etc... Atenciosamente Roney

6 6

Comparando series temporais com uma variavel externa, como analisar?
by Augusto Ribas 30 Jan '12

30 Jan '12

Ola pessoal. Bem pra economizar tempo, antes de ficar falando minha historinha eu vou fazer a minha pergunta pra não tomar muito o tempo de vocês. Como correlacionar series temporais com um atributo externo de cada serie temporal e estimar o quanto esse atributo deixa a serie temporal estável. Acho que essa é a pergunta que eu tenho que responder e ainda não sei como. Estava tentando aprender sobre series temporais mas estou apanhando aqui com o R. Basicamente eu gostaria de comparar varias series temporais com uma variável externa. A ideia é mais ou menos o seguinte, tenho dados da fenologia de varias especies de planta. Por exemplo, todo mes eu vou la no campo e para varias arvores de cada especie eu olho se ha brotos e estimo qual a porcentagem da copa que tem esses brotos. Ai fico com uma matriz onde colunas são especies e linhas são o estado fenológico da planta. Além disso, cada especie tem uma medida caracteristica dela, da especie, que é a densidade de madeira. A principio, me parece que especies com densidade de madeira alta so dão brotos no período de seca, enquanto as especies que tem densidade de madeira baixa são bagunçadas, brotam na chuva, fora da chuva, estão nem ai pra paçoca. Eu estava lendo o livro Analyzing Ecological Data (2007 by Zuur, Ieno & Smith), e pra esse tipo de coisa eles sugerem o uso de Min/max auto-correlation factor analysis (MAFA) e Dynamic factor analysis (DFA). Mas nos scripts do livro ele fala que ambos os métodos não eram implementados no R na época e ainda não achei no R eles. Eu estou garimpando atrás de algo nesse sentido aqui http://cran.r-project.org/web/views/TimeSeries.html Mas Time Series parece ser um campo muito amplo e eu to ficando perdido. Alguem tem alguma sugestão de literatura para Time Series mais voltado pra ecologia, pra fazer relação de series temporais com variáveis externas, mas não tão matematico ja que eu não tenho muita facilidade com matematica? Eu fiz um exemplo de dados de mais ou menos como estão as coisas. O exemplo é uma matriz de series temporais chamada brota, um vetor chamado chuva que é a pluviosidade no mês e um vetor que são valores de densidade de madeira das especies: Eu estava a principio explorando graficamente tudo , se alguem tiver alguma sugestão. Segue o CMR. #Exemplo: #Criando dados #Densidade de madeira set.seed(69) dmad<-round(c(rnorm(5,mean=50),rnorm(5,mean=90))) names(dmad)<-c(paste("sp",1:10,sep="")) dmad #Chuva chuva<-c(150,100,50,40,20,20,30,50,70,100,150,200,150,100,50,30,20,20,40,50,70,100,150,200) #Padrão Anual ydeb<-c(10,10,10,10,20,40,50,40,20,10,10,10) ydea<-rnorm(12,mean=mean(ydeb),sd=sd(ydeb)) #Amostras por 2 anos para especies com baixa e alta densidade de madeira brota<-matrix(NA,nrow=24, ncol=10,dimnames=list(paste("mes",1:24,sep=""),paste("sp",1:10,sep=""))) for (i in 1:5) { brota[,i]<-round(c(ydeb+rnorm(12,mean=5,sd=5),ydeb+rnorm(12,mean=5,sd=5)),digits=2) } for (i in 6:10) { brota[,i]<-round(c(ydea+rnorm(12,mean=5,sd=5),ydea+rnorm(12,mean=5,sd=5)),digits=2) } #TRansformando em serie temporal brota<-ts(brota) #Observando as Series plot(brota) #Observando as series juntas, separando as densidade de madeira por cores #Mais vermelho = mais densa, mais claro = menos densa cores<-heat.colors(10,alpha=1) matplot(c(1:24),brota,type="l",lwd=2,lty=1,xlab="meses",ylab="brota",col=cores[order(dmad)]) #Observando a chuva junto no grafico anterior graficobarras<-barplot(chuva/max(chuva)*100) matlines(graficobarras,brota,type="l",lwd=2,lty=1,xlab="meses",ylab="brota",col=cores[order(dmad)]) #Obsevando a partir de uma pca como são, mas acredito que isso esta errado fazer. brota.pca<-princomp(brota) biplot(brota.pca) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 0