Setembro 2014 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

2 2

RCurl: passando parâmetros e executando função javascript
by Éder Comunello 04 Dez '14

04 Dez '14

Colegas, boa tarde! Estou precisando de ajuda com relação ao pacote *RCurl*. Tomando por exemplo o endereço < http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001>, ao acessá-lo via *browser *pode-se notar dois botões que executam trechos de código em *javascript *e geram links para o *download *de arquivos de dados: * Botão [Arquivo Access] remete ao comando "javascript:criarArq(2254001,2)" * Botão [Arquivo Texto] remete ao comando "javascript:criarArq(2254001,1)" Estou tentando usar o *RCurl *para passar diretamente os parâmetros para criação do *link *para o "Arquivo Texto". Isto feito, posteriormente pretendo usar o *link *gerado para *download* dos dados. Contudo, não estou tendo sucesso e ficarei muito grato por qualquer ajuda. Abaixo segue o código que tentei aplicar: ### <BEGIN> ### library(RCurl) myURL <- 'http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001' ### URL resultante da operação via browser: ### myURL2 <- ' http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001&CriaArq=true&TipoArq=1 ' myCurl <- getCurlHandle(cookiefile = "myCookies.txt", useragent="Mozilla/5.0", verbose = T, autoreferer = T, followlocation = T) myParams <- list(Codigo=2254001, CriaArq='true', TipoArq=1) ### Codigo: ID da estação ### CriaArq: true (executa função?), ### TipoArq: 1 (tipo texto) e 2 (Access) postForm(myURL, .params=myParams, curl=myCurl) getURL(myURL, curl=myCurl) ### Operação falha, retornando página original! ### Em caso de sucesso embutiria no código um link no formato: ### <a href="ARQ/A20131008-164720-14/CHUVAS.ZIP"> ### <END> ### -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 3

PNAD
by Vanessa Felix 23 Out '14

23 Out '14

Pessoal, Eu tô querendo usar os dados da PNAD 2013 mas não consigo usar o pacote IBGEPesq. Fica aparecendo esse erro aqui: library(IBGEPesq) Error: package ‘IBGEPesq’ was built before R 3.0.0: please re-install it Tem alguma opção que não seja reinstalar o R? Se alguém puder ajudar... att, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ (32) 8854-4881 / 9142-8428 (21) 97656-3448 -------------------------------------------------------------------------------------------

4 4

mensagem de erro usando função glmer
by Ana Maria Nievas 16 Out '14

16 Out '14

Oi pessoal! Tentando reproduzir um script do livro do Zuur, me deparei com a seguinte mensagem de erro: O1.glmer <-glmer(NCalls~offset(LBroodSize)+SexParent*FoodTreatment+SexParent*ArrivalTime+(1 | fNest),data = an1, family = poisson) Mensagens de aviso perdidas: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, : Model is nearly unidentifiable: very large eigenvalue - Rescale variables? Alguém sabe como solucionar? -- Ana Maria Nievas Bióloga e Mestre em Ciências, com ênfase em Ecologia Aplicada (16) 98843-7744

3 3

Variáveis contínuas explicando Variáveis categóricas
by Jefferson Ferreira-Ferreira 08 Out '14

08 Out '14

Prezados; Uma dúvida básica de um iniciante em análises estatísticas. Tenho uma série de variáveis categóricas binárias e duas variáveis contínuas. Eu gostaria de saber o quando minhas variáveis categóricas explicam minhas duas variáveis contínuas. Por exemplo: c1 c2 c3 c4 x y 0 1 1 1 931 623 0 1 0 0 163 259 1 0 1 0 690 306 1 0 1 0 690 306 A pergunta é: o quanto a combinação das variáveis c explicam os valores de x e y? Ou, qual a correlação entre as variáveis c e as variáveis x e y? Podem me dar alguma ideia de análises possíveis? Obrigado. -- *Jefferson Ferreira-Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> *Contatos particulares:* *(55) 9615-0100*

9 9

Criar função no R com if e for
by Rodolfo Souza 02 Out '14

02 Out '14

Olá pessoal, gostaria mais uma vez da ajuda de vocês. Eu estou precisando criar uma função no R que tem três condições e após testar todas as condições o programa me retorne um resultado na forma de vetor para cada valor (input) na função. Segue o código que consegui fazer até agora: x=c(-5,-4,-3,-2,-1,0,1,2,3,4,5);x fr=function(x){ r=numeric(length = length(x)); for(i in seq(along=x)){ if(x<-3){r[i]=0} else if(x>=-3 && x<=2){r[i]=3*x} else{r[i]=1} } return(r) } rs=fr(x);rs O resultado saí: 0 0 0 0 0 0 0 0 0 0 0 Se eu tenho só um valor no vetor e não utilizo o "for" consigo resolver, mas se o vetor tem mais de um valor dá erro. Desde já grato pela atenção/ajuda. -- Rodolfo Marcondes Silva Souza Engenheiro Agrônomo Mestre em Produção Vegetal - UFRPE/UAST Doutorando em Tecnologias Energéticas e Nucleares - UFPE/DEN Currículo Lattes <http://lattes.cnpq.br/4486320338778858>

5 5

Gerar vetor de dados com assimetria negativa no R
by Andre Oliveira 01 Out '14

01 Out '14

Pessoal, como posso gerar um vetor de tamanho n de dados assimétricos com assimetria negativa? Pode ser qualquer distribuição! André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 3

Gráfico
by Antonio Armando Ulian do Lago Albuquerque 01 Out '14

01 Out '14

Caros amigos, Sou muito novato no R, mas tenho nos últimos 6 meses lido muito e tentando realizar algumas operações muito básicas. Infelizmente ainda estou gatinhando. Utilizo o UBUNTU 14.04, 64bits, RKward. A partir do banco de dados do TSE, baixei os arquivos referentes a raça nas eleições 2014. Como não consegui abri-lo na extensão csv transformei-os em xlsx e utilizei o pacote XLConnect. Quero elaborar 2 gráficos com os dados, mas não estou conseguindo. As tabelas são as seguintes. Quais comandos devo elaborar? Summary Raça Qtde Percent 1 Amarela 120 0.46 2 Branca 14379 54.96 3 Indígena 86 0.33 4 Parda 9158 35.00 5 Preta 2422 9.26 6 Total 26165 100.00 Raça Qtde Percent Length:6 Min. : 86.0 Min. : 0.33 Class :character 1st Qu.: 695.5 1st Qu.: 2.66 Mode :character Median : 5790.0 Median : 22.13 Mean : 8721.7 Mean : 33.34 3rd Qu.:13073.8 3rd Qu.: 49.97 Max. :26165.0 Max. :100.00 Muito obrigado e desculpe pela pergunta e dificuldade tão raza a qual vos exponho. Antonio

4 4

[R-br] Parâmetros_=?UTF-8?Q?modelo_de_regress=C3?==?UTF-8?Q?=A3o_din=C3=A2mico?=
by Fatima do Nascimento Silva 30 Set '14

30 Set '14

olá caros colegas, boa noite! estou ajustando um modelo de regressão dinâmico e preciso ponderar o vetor de parâmetros beta 2 (b2) por 7, que representaria um aumento de 7° na variável temperatura. E depois faria as previsões inserindo estes novos valores de b2 no calculo dos valores ajustados. Alguém poderia me dar uma dica de como alterar isso no R? segue os comandos utilizados. desde já agradeço a atenção. ponderar o beta 2 (L(Temperatura, 1) = -4.8725) por 7. caso possa me ajudar com alguma sugestão agradeço muito. #Comandos (CMR) rm(list=ls(all=TRUE)) #pacotes utilizados require(dynlm) require(tseries) require(car) require(verification) require(MASS) require(forecast) require(TSA) require(nortest) ################################################################## (dados=read.table(file="BDAC.txt",header=T,dec=",")) attach(dados) (dados1<-ts(dados,frequency=12, start=c(2000,1),end=c(2010,12))) #modelo ajustado MOD1 <- dynlm(BEAI~trend(BEAI)+season(BEAI)+L(Temperatura,1)+L(Poluente,3),data = dados1 ) summary(MOD1) par(mfrow = c(2,2)) plot(MOD1, main="") names(MOD1) #lista de objetos do modelo #previsões p<-(forecast(MOD1$fitted.values, h=12, level=c(80,95), fan=FALSE)) plot(p) #saída do modelo > summary(MOD1) Time series regression with "ts" data: Start = 2000(4), End = 2010(12) Call: dynlm(formula = BEAI ~ trend(BEAI) + season(BEAI) + L(Temperatura, 1) + L(Poluente, 3), data = dados1) Residuals: Min 1Q Median 3Q Max -42.946 -13.816 -1.482 12.327 80.586 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 263.2442 45.6925 5.761 7.21e-08 *** trend(BEAI) 0.4470 0.6201 0.721 0.47241 season(BEAI)Feb -19.3407 9.5712 -2.021 0.04565 * season(BEAI)Mar 9.0604 9.6612 0.938 0.35032 season(BEAI)Apr 5.1280 9.4526 0.543 0.58853 season(BEAI)May 25.2990 9.6410 2.624 0.00988 ** season(BEAI)Jun 35.7033 11.7395 3.041 0.00292 ** season(BEAI)Jul 30.0824 11.3742 2.645 0.00933 ** season(BEAI)Aug 19.1968 11.9829 1.602 0.11192 season(BEAI)Sep 4.6108 10.3265 0.447 0.65608 season(BEAI)Oct 3.7939 10.4920 0.362 0.71832 season(BEAI)Nov -1.3031 9.4690 -0.138 0.89078 season(BEAI)Dec -7.6315 9.4334 -0.809 0.42020 L(Temperatura, 1) -4.8725 1.5284 -3.188 0.00185 ** L(Poluente, 3) 2.6678 0.8455 3.155 0.00205 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 21.34 on 114 degrees of freedom Multiple R-squared: 0.618, Adjusted R-squared: 0.5711 F-statistic: 13.17 on 14 and 114 DF, p-value: < 2.2e-16 Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

3 5