Dezembro 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

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

Salvar como Geotiff objeto krige.conv do geoR
by ASANTOS 31 Dez '13

31 Dez '13

Boa tarde Pessoal, Gostaria de pegar o resultado de uma krigagem realizada com a função krige.conv() do pacote geoR e salvar como uma imagem Geotiff, abaixo fiz um CRM com a função writeGDAL(), sem sucesso, alguém teria alguma ideia para resolver isto, obrigado, require(geoR) loci <- expand.grid(seq(0,1,l=21), seq(0,1,l=21)) # Predição por krigagem ordinária kc <- krige.conv(s100, loc=loci, krige=krige.control(cov.pars=c(1, .25))) require(rgdal) geoTiff <- tempfile() writeGDAL(kc["predict"], geoTiff , drivername="GTiff") GDALinfo(tf) -- ====================================================================== 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 ======================================================================

2 2

Qual a melhor opção de Editor, GUI, IDE para o R?
by Jose Claudio Faria 29 Dez '13

29 Dez '13

Gostaria de convidar a Comunidade Brasileira de Usuários do R para uma lida crítica no link: http://nbcgib.uesc.br/lec/llec/avale-es/editor-gui-ide Minha intenção não é criar mais polêmica do que a existente sobre o assunto. Mas sim, prover esclarecimentos para quem está em busca de informações nesse sentido, baseado na experiência adquirida ao longo da minha vida de usuário e desenvolvedor. Inicialmente havia escrito o texto para postar diretamente na lista, mas acabei optando por colocar em uma página do LEC (Lab. Est. Computacional/UESC) devido a: - Melhor possibilidade de formatação; - Links; - Para ficar didático. As possíveis manifestações podem se dar aqui mesmo na lista, pois o LEC ainda não possui um forum de discussão. Boa leitura, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

10 13

Como produzir gráficos com diferente símbolos e cores in PCA?
by Jackson Rodrigues 29 Dez '13

29 Dez '13

Olá pessoal meu nome é Jackson Rodrigues e estou utilizando o pacote vegan para fazer análise multivariada, PCA 'Principal Component Analysis"para ser mais preciso. Direto ao ponto. Tenho um conjunto de dados polínicos com 54 linhas (representando os sites) e 145 colunas (representando os táxons). Eu quero plotar os 2 primeiros PC´s sendo de maneira que seja um grupo com 27 sites com um único símbolo e cores diferentes e os outro grupo com 27 sites utilizando outro símbolo repetindo as mesmas cores anteriores. Alguém pode me indicar como fazê-lo. Estou usando o script abaixo, no entanto não sei como utilizar símbolos para cada grupo. Pre_euro_veg.Pollen<-read.csv("E:/Modern data/Vegetation changes/Veg chang 1/sem Pinus/pca test.csv", header=TRUE)#Apenas trees/shrubs e herbs rownames(Pre_euro_veg.Pollen) <- c("Conf", "Pires", "Silv", "Capa", "SDO", "SDB", "Ita", "Botu", "SCG", "Arac", "VoVe_1", "PoGd", "Misi", "SBV", "Ciama", "Tabu", "MDI", "CaPu", "Mate", "SJdA", "CDS", "Pint", "SMar", "RdCa", "SFA", "Sant", "StMon", "Conf1", " Pires1", "Silv1", " Capa1", " SDO1", " SDB1", "Ita1", "Botu1", "SCG1", "Arac1", "VoVe_11", "PoGd1", "Misi1", "SBV1", "Ciama1", "Tabu1", "MDI1", "CaPu1", "Mate1", "SJdA1", "CDS1", "Pint1", "SMar1", "RdCa1", "SFA1", "Sant1", "StMon1") colvec <- c("antiquewhite", "antiquewhite4", "aquamarine", "aquamarine4", "bisque1", "bisque4", "black", "blue4", "blueviolet", "brown", "burlywood", "chartreuse", "chocolate", "cornflowerblue", "cornsilk", "darkgoldenrod1", "darksalmon", "darkorchid4", "darkslategrey", "gold", "gold4", "lavenderblush2", "lemonchiffon", "lemonchiffon3", "lightblue", "lightpink", "navyblue", "antiquewhite", "antiquewhite4", "aquamarine", "aquamarine4", "bisque1", "bisque4", "black", "blue4", "blueviolet", "brown", "burlywood", "chartreuse", "chocolate", "cornflowerblue", "cornsilk", "darkgoldenrod1", "darksalmon", "darkorchid4", "darkslategrey", "gold", "gold4", "lavenderblush2", "lemonchiffon", "lemonchiffon3", "lightblue", "lightpink", "navyblue") Pre_euro.colour.pca<-plot(Pre_euro_veg.1.all.pca, display = "sites", type = "t") with(Veg_chng_clim, points(Pre_euro_veg.1.all.pca, display = "sites", col = colvec[Sites], scaling = 2, pch = 21, bg = colvec[Sites])) head(with(Veg_chng_clim, colvec[Sites])) with(Veg_chng_clim, legend("bottomleft", legend = levels(Sites), bty = "o", col = colvec, pch = 21, pt.bg = colvec)) title("PCA of All sites") Se não fui claro o suficiente, por favor pergunte. Obrigado. Jackson -- Jackson M. Rodrigues Department of Palynology and Climate Dynamics Albrecht-von-Haller-Institute for Plant Sciences Georg-August-University Göttingen Untere Karspuele 2 37073 Göttingen/Germany Tel.: 0049 (0) 176 6133 1991 Web: http://www.uni-goettingen.de/en/306700.html

3 2

Recortar raster com limite em data.frame
by ASANTOS 29 Dez '13

29 Dez '13

Boa tarde Pessoal, Gostaria de saber qual a melhor forma de recortar um rasterlayer se eu tenho o limite de uma área em data frame, pois segundo o exemplo abaixo não consigo fazer a operação com o uso de função overlay() para retirar os NA's que correspondem a malha fora do limite sem ficar fazendo inúmeras transformações dos objetos raster e data.frame, segue CRM: ## Raster----------------------------------------------------------- set.seed(20131226) r <- raster(nc=100, nr=100) r <- setValues(r, round(runif(ncell(r))* 255)) image(r,col=gray(seq(0,1,0.025))) ## Borda ----------------------------------------------------------- x<-c(-125,-100,0,100,-125) y<-c(-49,25,0,-50,-49) limite<-cbind(x,y) lines(limite, col="yellow") ## 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 ======================================================================

3 6

encoding (problema)
by Edson Lira 27 Dez '13

27 Dez '13

Caros amigos, estou lendo um arquivo em csv originado no odk agregate e não consegui qual o formato que usam, latin, UTF-8, … Já tentei vários formatos e não consegui. A rotina que estou usando está abaixo wj<-read.csv2("pasta/WJ_12_2013.csv",header=T,encoding="macroman") wj[wj==""]<-NA O link para o arquivo em parcial.csv está abaixo: http://www.datafilehost.com/d/ffb7270f Minha máquina é um Mac Book Pro. A versão do R é a 3.0.1 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

Selecionar coordenadas e valor de pixels contidos no interior de parcelas
by ASANTOS 27 Dez '13

27 Dez '13

Boa tarde Pessoal, Estou quebrando a cabeça em um desafio natalino e gostaria de compartilhar o desafio com vocês, tenho uma imagem simulada stackIm1, sobre a qual lancei pontos aleatórios position.parc que são o centro de parcelas circulares CIR de 5 de raio, sendo meu CRM: require(raster) require(spatstat) require(sp) ##Criando um raster ----------------------------------------------------------- r <- raster(nc=100, nr=100) set.seed(20130622) stackIm1 <- stack(lapply(1, function(x) setValues(r, round(runif(ncell(r))* 255)))) plot(stackIm1) # ## Criando 15 pontos amostrais ------------------------------------------------- att.Im1<-as(stackIm1, "SpatialGridDataFrame") df.Im1<-as.data.frame(att.Im1) colnames(df.Im1)<-c("px","x","y") gridded(df.Im1) = ~x+y position.parc<-spsample(df.Im1,15, type="random") points(position.parc, col="red",pch=20) # #Criar as parcelas circulares em torno do ponto amostral com 5 de raio # CIR <- vector("list",length(position.parc$x)) for(i in 1:length(position.parc$x)) CIR[[i]] <- disc(radius=5, centre=c(position.parc$x[i],position.parc$y[i])) # Gostaria de alguma luz para selecionar as informações("px","x","y") para os pixels contidos no interior de cada parcela circular CIR na imagem stackIm1. 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 ======================================================================

3 6

Substituir NA por 999
by Luciane Maria Pilotto 26 Dez '13

26 Dez '13

Olá pessoal, tem como substituir os valores NA por 999 diretamente no banco e salvar? E fazer o inverso tb, substituir 999 por NA? Tentei a função recode e não deu certo Usei ifelse e due certo, porém, cria novo vetor. Gostaria de substituir nomesmo vetor do banco. renddomd<-ifelse(rendadomexcp==999999999999,NA,rendadomexcp)cor<-cbind(cor,renddomd) Atenciosamente, Luciane Maria Pilotto Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS

3 2