Setembro 2013 - 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

Manual
by Maicon J. Galiazzi 13 Nov '13

13 Nov '13

Olá, estou desenvolvendo um manual de fácil utilização ao usuário, e no momento estou descrevendo a interfacedo software R Stuio 0.97.551, assim como a barra de ferramentas do mesmo. Podem me indicar algum material no qual posso conferir? Desde já agradeço! Att Maicon J. Galiazzi

2 2

Notação potências de dez (Gráficos)
by Eric Leandro 18 Out '13

18 Out '13

Boa tarde pessoal Tenho uma dúvida quanto a intervalo com potências de dez no R. Quando dou este comando : plot(areadrenagem,vazao, xlab="Área drenagem (m²)",ylab="Q(m³/s)",pch=16,xlim=range(0:5e+12),ylim=range(0:150000)), No eixo x aparece 1e+12..... 5e+12, como eu faria para aparecer neste formato : 2.1012 ..... 5.1012. ATT Eric Leandro Silva Pereira Mestrando em Clima e Ambiente Instituto Nacional de Pesquisas da Amazônia Programa de Grande Escala da Biosfera-Atmosfera na Amazônia Av. André Araújo, 2936 - Aleixo, CEP 69060-001, Manaus - AM.

4 14

Integração entre R e Calc
by Andre Oliveira 03 Out '13

03 Out '13

Pessoas boa noite, gostaria de saber se existe alguma pacote ou forma de integrar R + calc? André Oliveira Souza

3 2

Produção de Gráficos: PCA
by Vinícius Lionel Mateus 01 Out '13

01 Out '13

Prezados, Estou começando a aplicar PCA usando R, mas anteriormente o fazia no Statistica. Tenho utilizado os pacotes 'psych' e 'GPArotation' pois prefiro utilizar rotações em Varimax, porém os gráficos não ficaram ao meu gosto. Depois de algum esforço, eu consegui trocar as labels do meu gráfico, mas além disto gostaria de retirar o histograma, e plotar gráficos da dispersão das amostras(labels) e das variáveis entres os componentes, em 2D e 3D. Segue meu CMR. Desde já, obrigado. # instalando os pacotes install.packages("psych") install.packages("GPArotation") #carregando require (psych) install.packages (GPArotation) #entrada de dados (https://www.dropbox.com/s/l8v9g6qub27oig0/enqalog.xls) read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> d #entrada de dados1 read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> enqa enqa -> d attach(d) names(d) summary(d) #PCA(Varimax)/psych dpca<-principal(d[,2:16], nfactors=3, rotate="varimax", scores=T) #Scree plot plot (dpca$values, type="b", ylab="Eigenvalue", xlab="Component", lab=c(5,5,5)) #biplot biplot(dpca, labels=d[,1], cex=c(0.75,1)) Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

2 2

Desdobramento de interações em dados com distribuição de Poisson inflacionada por zeros
by ASANTOS 01 Out '13

01 Out '13

Bom dia Pessoal, Estou quebrando a cabeça para tentar desdobrar interações em dados com distribuição de Poisson inflacionada por zeros, na verdade é uma situação muito recorrente em entomologia, onde em algumas situações temos zero insetos e em outras muitos, considerando-se dados de contagens. Fiz uma pequena simulação abaixo para ilustrar o problema (Variáveis observadas Trat e Epoca e variável resposta y), se acho um modelo significativo que tem uma parte com distribuição de Poisson e outra com Binomial, como faço o desdobramento? Uma para cada distribuição ou tem alguma abordagem diferente? Segue CRM abaixo: ## Desdobramento de dados com distribuição de Poisson inflacionado por zeros --- require(pscl) y1<- c(mapply(rpois, lambda=c(5,20,45), MoreArgs=list(n=400)))##Criação da variável resposta Poisson y2<- c(mapply(rbinom, size=c(1,0,0), prob=c(0.5,1,1), MoreArgs=list(n=200)))##Criação da variável resposta Binomial y<-c(y1,y2) trat <- as.factor(gl(3,600)) ##Criação dos tratamentos tempo<- as.factor(rep(gl(6,100),3)) ### Criação da variável tempo dados<-as.data.frame(cbind(trat,tempo, y)) #------------------------------------------------------------------------------- # Análise de variância do dados inflacionados summary(m1 <- zeroinfl(y ~ trat | epoca, data = dados)) ## Modelo completo mnull <- update(m1, . ~ 1) ### Modelo nulo pchisq(2 * (logLik(m1) - logLik(mnull)), df = 2, lower.tail = FALSE) ## Teste de Chi o modelo completo foi significativo ## Comparando com GLM Poisson sem inflação por zeros --------------------------- summary(p1 <- glm(y ~ trat + epoca, family = poisson, data = dados)) ##Teste de Vuong---------------------------------------------------------------- vuong(p1, m1) ## O GLM Poisson estava tão mal ajustado que o Poisson iflacionado é a única opção # ## E agora, como desdobrar isto? Faço um desdobramento para parte de poisson e outra para parte binomial? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

4 13

[OFF] Gerenciador de referencias
by Andre Oliveira 01 Out '13

01 Out '13

Alguém me indica um bom gerenciador de referencias. Sem ser o JabRef, pois, ele não abre aqui no meu ubuntu? André Oliveira Souza

9 15

[OFF-TOPIC] regex
by Roney Fraga Souza 30 Set '13

30 Set '13

Caros, Tenho um arquivo em formato txt com abstracts de artigos científicos, segue exemplo: We examine the productive efficiency of 70 Indian commercial banks during the early stages (1986-1991) of the ongoing period of liberalization. We use data envelopment analysis to calculate radial technical efficiency scores. We then use stochastic frontier analysis to attribute variation in the calculated efficiency scores to three sources: a temporal component, an ownership component, and a random noise component. We find publicly-owned Indian banks to have been the most efficient, followed by foreign-owned banks and privately-owned Indian banks. We also find a temporal improvement in the performance of foreign-owned banks, virtually no trend in the performance of privately-owned Indian banks, and a temporal decline in the performance of publicly-owned Indian banks. We attempt to explain these patterns in terms of the government's evolving regulatory policies. (C) 1997 Elsevier Science B.V Recently, considerable attention has been focused on the performance of various airlines and air carriers in terms of efficiency. Although it is obvious that air carriers use airports, few studies have focused on airport operational efficiency. This empirical study evaluates the operational efficiencies of 44 major U.S. airports using data envelopment analysis and some of its recent developments. Various airport characteristics are evaluated to determine their relationship to an airport's efficiency. Efficiency measures are based on four resource input measures including airport operational costs, number of airport employees, gates and runways, and five output measures including operational revenue, passenger flow, commercial and general aviation movement, and total cargo transportation. The results of this study have operational as well as public policy implications. (C) 2000 Elsevier Science B.V. All rights reserved considerando que cada parágrafo esta em uma linha, ou seja, a quebra de linha existe apenas no final de cada abstract, preciso apagar tudo que estiver a partir do '(C)' até a quebra de linha. No primeiro parágrafo preciso eliminar '(C) 1997 Elsevier Science B.V', no segundo parágrafo '(C) 2000 Elsevier Science B.V. All rights reserved'. Uso um macbook, qualquer sugestão via R ou terminal é bem vinda. Atenciosamente Roney

2 2

Re: [R-br] ajuda normatização de grid
by Eder Comunello 27 Set '13

27 Set '13

Prezado Elias, Fornecendo dados junto com seu código, permitindo torná-lo reproduzível (CMR), você obtém ajuda mais facilmente. Suponho que o 17 seja em função da variável que você quer representar no grid que está exportando (predição). Estou supondo que 17 seja o índice da última coluna adicionada (cbind). Um objeto SpatialPixelsDataFrame permite multivariáveis, mas pode ser que o formato escolhido exija a seleção de uma única variável pra compor um grid. ================================================ Éder Comunello Em 27 de setembro de 2013 07:43, ELIAS MENDES COSTA MENDES escreveu: > > Obrigado Éder, era isso mesmo, agora que rodei um comando para tirar todos > os NAs deu certo. > Já aproveitando o email talvez pudesse me ajudar mais uma vez. No script > que estou usando o comando que uso para salvar o mapa predito em um formato > que > abra no SagaGis é: writeGDAL(grid_V[,17], "V_rna37.sdat", > drivername="SAGA", mvFlag=-99999) Você sabe explicar o porque do 17? > > Os comandos anteriores a esse foram. > > # fazer o mapa para V%, ou aplicar o modelo para a area de estudo > grid_V <- cbind(grid1, predito) # usar "predito" ou "test1" > names(grid_V) > gridded(grid_V) <- ~X+Y # transformar "grid_V" de DataFrame para > SpatialPixelsDataFrame > proj4string(grid_V) <- CRS(utm) # definir projecao > names(grid_V) > spplot(grid_V, "predito") # plotar no R apenas a coluna "predito" > >

1 0