- R-br - DINF UFPR

Assunto: Digest R-br, volume 98, assunto 1
by Edmar Caldas 01 Fev '19

01 Fev '19

instala o pacote ("plyr" uniao = rbind.fill(arquivo_a, arquivo_b) Enviado do Yahoo Mail no Android <div>Em sex, 1 1e fev 1e 2019 às 12:00, r-br-request(a)listas.c3sl.ufpr.br</div><div><r-br-request(a)listas.c3sl.ufpr.br> escreveu:</div> Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Assunto: Re: JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES (Edson Lira) ---------------------------------------------------------------------- Message: 1 Date: Thu, 31 Jan 2019 22:19:36 +0000 (UTC) From: Edson Lira <edinhoestat(a)yahoo.com.br> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Assunto: Re: JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES Message-ID: <1661977038.1186878.1548973176805(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Muito obrigado. Enviado do Yahoo Mail no Android <div>Em qua, 30 30e jan 30e 2019 às 22:56, Marcus Nunes por (R-br)</div><div><r-br(a)listas.c3sl.ufpr.br> escreveu:</div> _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190131/bb012f5a/atta…> -------------- Próxima Parte ---------- Um texto embutido e sem conjunto de caracteres especificado foi limpo... Nome: Untitled Url: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190131/bb012f5a/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 98, assunto 1 ****************************************

1 0

JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES
by Edson Lira 31 Jan '19

31 Jan '19

Boa tarde a todos!!! Tenho dois arquivos: um com os CEP's de pacientes de um hospital (57000 registros) sem as coordenadas, o outro arquivo que consegui no link abixo http://results.openaddresses.io/?runs=all#runs No link acima tem POSTCODES (CEP"S) de todo o mundo. Baixei o do Amazonas que tem o CEP DE MANAUS (SEPAREI ESSE ARQUIVO com 554000 registros) Quero fazer uma rotina aonde eu possa juntar o arquivo de consultas sem as coordenadas com o arquivo de CEP'S com as coordendas. (LAT LON), ou seja, quero fazer um arquivo mantendo a dimensão do arquivo das consultas. Espero ter sido claro. Abaixo o link dos arquivos (um com CEP e o outro com CEP e Coordenadas) https://www.datafilehost.com/d/b93da0f3 https://www.datafilehost.com/d/b8f65984 [ ]'s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 5

Problema para plotar IC para modelo glm Gamma usando stat_smooth() no ggplot2
by ASANTOS 30 Jan '19

30 Jan '19

Prezados Membros do R-br, Estou tentando plotar um modelo glm Gamma com stat_smooth() no ggplot2 sem sucesso. Minha ideia é representar o intervalo de confiança para um modelo onde realizei a junção de dois níveis para o fator Feature. Segundo exemplo fictício do meu problema: #Simulação de um banco de dados, onde a variável resposta Production tem distribuição Gamma set.seed(123) d<-NULL N<-50 d$Production <- rgamma(N,10) d$Feature <- ifelse(d$Production >7 & d$Production<10, c("green"),ifelse(d$Production>11, c("red"), c("blue"))) d$Temp<-rnorm(N,20,5) d<-as.data.frame(d) # # Ajusto um modelo de Gamma completo mG<- glm(Production~ Feature + Temp, family= Gamma, data = d) summary(mG) anova(mG,test="Chi") # Comparo os níveis do fator Feature library(multcomp) PW<-summary(glht(mG, linfct = mcp(Feature= "Tukey"))) PW cld(PW) # Faço de conta que os níveis green = blue e realizo um novo ajuste Feature2<-d$Feature levels(Feature2) levels(Feature2)[1]<-"blue&green" levels(Feature2)[2]<-"blue&green" levels(Feature2) d$Feature2<-Feature2 mG2<- glm(Production~ Feature2 + Temp, family= Gamma, data = d) # Crio os valores de predição pred.data = data.frame( Feature2<-d$Feature2, Temp<-d$Temp ) pred.data$Production = predict(mG2, newdata=pred.data, type="response") #Realizo o plot library("ggplot2") ggplot(d, aes(Temp, Production, colour = Feature2)) + geom_point() + geom_line(data=pred.data) + stat_smooth(method = "glm", formula = y ~ x, family = Gamma) # E tenho uma representação errada do intervalo de confiança do meu modelo mG2 porque não consigo fazer a função stat_smooth() considerar o ajuste mG2. Como posso resolver isso? 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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

2 2

análise da estrutura de populações
by Ze Henrique 29 Jan '19

29 Jan '19

Boa noite Srs. gostaria de fazer uma figura de análise da estrutura de populações. Sei que o programa STRUCTURE, faz estas análises, porém gostaria de saber se existe algum pacote implementado no R que faz. Alguém poderia me recomendar algum pacote. Desde já agradeço pela atenção. -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

3 2

Como Salvar Shapefile
by Yury Duarte 22 Jan '19

22 Jan '19

Bom dia colegas listeiros! Ultimamente venho tentando manusear arquivos espaciais utilizando o R. Em alguns casos, objetos do tipo 'large SpatialPolygons' são gerados dentro do código. Quando isso ocorre (acredito que seja algo que acontece por padrão, devido a grande quantidade de informações armazenadas no objeto), a função writeOGR, que utilizo para salvar os shapefiles, quebra, por só tratar de objetos do tipo 'Spatial' e não os 'large Spatial'. Sendo assim, gostaria de saber se existe a possibilidade de salvar um shapefile a partir de um objeto do tipo 'large Spatial' ou ainda se existe a possibilidade de transformar esse num objeto do tipo SpatialPolygonsDataFrame, para que a função writeOGR possa funcionar normalmente. Segue código desenvolvido seguido do erro gerado. Desde já, agradeço pela ajuda de todos! rodar_bibliotecas = function(necessarias_para_o_projeto){ if(necessarias_para_o_projeto){ library(rgrass7) library(spatstat) library(maptools) library(shapefiles) library(foreign) library(magrittr) library(formattable) library(e1071) library(rlang) library(rgdal) library(rgeos) library(raster) library(sp) library(sf) library(RcppCNPy) library(deldir) library(dismo) library(dplyr) library(ggplot2) library(gstat) library(tidyverse) library(smoothr) library(viridisLite) } return('Ok') } rodar_bibliotecas(TRUE) raiz = '/home/yury/pesquisa/arquivos/areas/' clientes = dir(raiz) for (cliente in 1:length(clientes)) { analises = dir(paste0(raiz, clientes[cliente])) q = list() for (analise in 1:length(analises)) { shapefile_contorno = readOGR(paste0(raiz, clientes[cliente], '/', analises[analise], '/', 'vetores.shp')) shapefile_contorno$id = 1 dissolvido = unionSpatialPolygons(shapefile_contorno, IDs = shapefile_contorno$id) #plot(dissolvido) q[analise] = dissolvido } fim = q[[1]] for (x in 2:length(q)){ fim = bind(fim, q[[x]]) } plot(fim) path_to_save = paste0(raiz, clientes[cliente]) writeOGR(fim, dsn = path_to_save, layer = paste0('contorno_', clientes[cliente]), driver = 'ESRI Shapefile', overwrite_layer = TRUE) } Error in writeOGR(fim, dsn = path_to_save, layer = paste0("contorno_", : obj must be a SpatialPointsDataFrame, SpatialLinesDataFrame or SpatialPolygonsDataFrame Yury Duarte Engenheiro Agrônomo - ESALQ/USP

4 7

Tinn-R: novo website sob protocolo https
by Jose Claudio Faria 19 Jan '19

19 Jan '19

Prezados, O projeto Tinn-R está com novo website sob o protocolo https: Dears, The Tinn-R projetc has a new website under the https protocol: https://nbcgib.uesc.br/tinnr/en/ Best, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)99966.9100 - VIVO ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ If you have software to deal with statistics, you have arms; if you have good software, you have arms and legs; if you have software like R, you have arms, legs and wings... the height of your flight depends only on you!

2 1

Ajuda com a passagem incorreta de variável string para uma função do R
by Bruno César Pasquini 18 Jan '19

18 Jan '19

PessoALL, bom dia! Estou com dificuldade em passar os valores de uma variável string para uma função do R de forma que esta função a interprete corretamente. Por exemplo, para indicar qual coluna eu quero que a função selecione para um cruzamento de 2 data.tables (acho melhor que usar data.frames pela melhor performance do data.table)... Eu poderia levar a tabela inteira para o cruzamento, mas depois teria que deletar as demais colunas manualmente, gostaria de resolver com a função conforme o exemplo abaixo. Construí este exemplo baseado nos resultados do primeiro turno das últimas eleições presidenciais. Normalmente eu importaria as tabelas, mas eu as declaro aqui para a devida reprodução do código. Fonte: https://pt.wikipedia.org/wiki/Resultados_da_eleição_presidencial_no_Brasil_… Meu propósito é entender o que estou fazendo de errado (ou deixando de fazer...), não estou propondo um debate político-ideológico que foge ao escopo desta lista... Especificamente, eu não consigo fazer a função "*consolida*" (na verdade as versões 1, 2 e 3 desta função) interprete(m) o parâmetro "*Candidato*" corretamente. Eu tento usar o parâmetro sozinho e com as funções "*substitute*" e "*eval*", mas não consigo levar a coluna com os números de votação do candidato por região para o cruzamento. Sim, eu sei que o cruzamento como o desenhei vai gerar duplicidade, mas o que preciso é que o cruzamento "funcione"... Vocês podem me ajudar a entender o que estou fazendo de errado? Muito grato pela atenção de cada um(a), Bruno César Pasquini Estatístico ####### Codigo Reprodutivel... library(data.table) ### Percentual de Votos Validos nos estados e no exterior... perc_validos <- data.table(Regiao = c('Centro-Oeste', 'Centro-Oeste', 'Centro-Oeste', 'Centro-Oeste', 'Exterior', 'Nordeste', 'Nordeste', 'Nordeste', 'Nordeste', 'Nordeste', 'Nordeste', 'Nordeste', 'Nordeste', 'Nordeste', 'Norte', 'Norte', 'Norte', 'Norte', 'Norte', 'Norte', 'Norte', 'Sudeste', 'Sudeste', 'Sudeste', 'Sudeste', 'Sul', 'Sul', 'Sul'), Estado = c('Distrito Federal', 'Goias', 'Mato Grosso', 'Mato Grosso do Sul', '', 'Alagoas', 'Bahia', 'Ceara', 'Maranhao', 'Paraiba', 'Pernambuco', 'Piaui', 'Rio Grande do Norte', 'Sergipe', 'Acre', 'Amapa', 'Amazonas', 'Para', 'Rondonia', 'Roraima', 'Tocantins', 'Espirito Santo', 'Minas Gerais', 'Rio de Janeiro', 'Sao Paulo', 'Parana', 'Rio Grande do Sul', 'Santa Catarina'), Bolsonaro = c(0.5837, 0.5724, 0.6004, 0.5506, 0.5879, 0.3440, 0.2341, 0.2174, 0.2428, 0.3130, 0.3057, 0.1876, 0.3021, 0.2721, 0.6224, 0.4074, 0.4348, 0.3619, 0.6224, 0.6297, 0.4464, 0.5476, 0.4831, 0.5979, 0.5300, 0.5689, 0.5263, 0.6582), Haddad = c(0.1187, 0.2186, 0.2476, 0.2387, 0.1010, 0.4475, 0.6028, 0.3312, 0.6126, 0.4546, 0.4887, 0.6340, 0.4119, 0.5009, 0.1853, 0.3277, 0.4030, 0.4139, 0.2036, 0.1785, 0.4112, 0.2420, 0.2765, 0.1469, 0.1642, 0.1970, 0.2281, 0.1513), Ciro = c(0.1660, 0.0860, 0.0559, 0.0804, 0.1452, 0.1012, 0.0941, 0.4095, 0.0839, 0.1675, 0.1356, 0.1142, 0.2231, 0.1302, 0.0517, 0.1234, 0.0750, 0.1003, 0.0603, 0.0536, 0.0717, 0.0954, 0.1164, 0.1522, 0.1135, 0.0831, 0.1137, 0.0668), Outros = c(0.1316, 0.1230, 0.0961, 0.1304, 0.1659, 0.1072, 0.0690, 0.0419, 0.0606, 0.0650, 0.0700, 0.0642, 0.0631, 0.0967, 0.1405, 0.1417, 0.0873, 0.1237, 0.1739, 0.1383, 0.0707, 0.1149, 0.1241, 0.1030, 0.1922, 0.1510, 0.1318, 0.1236) ) ### Analisando o data.table perc_validos... perc_validos; class(perc_validos); str(perc_validos); dim(perc_validos) ### Total de Votos Validos por regiao e no exterior... votos_regiao <- data.table(Regiao = c('Centro-Oeste', 'Exterior', 'Nordeste', 'Norte', 'Sudeste', 'Sul'), 'Jair Bolsonaro' = c(4555415, 113690, 7453186, 3785038, 23915925, 9453736), 'Fernando Haddad' = c(1642054, 19540, 14583334, 3211002, 8623232, 3262843), 'Ciro Gomes' = c(750776, 28073, 4892961, 748170, 5425104, 1499282), 'Demais Candidatos' = c(952142, 32078, 1894092, 982103, 6968773, 2258047) ) ### Analisando o data.table votos_regiao... votos_regiao; class(votos_regiao); str(votos_regiao); dim(votos_regiao) ### Fazendo o cruzamento... Ou pelo menos tentando... consolida1 <- function(Candidato) { temporario <- merge(x = perc_validos, y = votos_regiao[, .(Regiao, Candidato)], by = "Regiao", all.x = TRUE, sort = TRUE) return(temporario) } consolida1(Candidato = "Jair Bolsonaro") consolida2 <- function(Candidato) { temporario <- merge(x = perc_validos, y = votos_regiao[, .(Regiao, substitute(Candidato))], by = "Regiao", all.x = TRUE, sort = TRUE) return(temporario) } consolida2(Candidato = "Ciro Gomes") consolida3 <- function(Candidato) { temporario <- merge(x = perc_validos, y = votos_regiao[, .(Regiao, eval(Candidato))], by = "Regiao", all.x = TRUE, sort = TRUE) return(temporario) } consolida3(Candidato = "Fernando Haddad")

3 6

Construir um laço que modifique o conteúdo do meu contador
by Yury Duarte 18 Jan '19

18 Jan '19

Bom dia a todos! Há algum tempo venho tentando encontrar uma lógica de laço aonde eu consiga eliminar determinadas posições de uma lista (que uso de contador), de acordo com uma condição, e seguir com laço. Em suma, dentro de uma lista de polígonos, preciso identificar os polígonos que se sobrepõe (mesmo que apenas se encostem), uni-los num novo polígono (sobrescrevendo o primeiro?), eliminando o que restou da união e seguindo o fluxo do laço com o resultado da união. Assim, um polígono que foi unido não mais seria considerado no laço. Segue um código como exemplo, aonde testei algumas possibilidades para contornar a questão. Acredito que a leitura do código elucide melhor meu problema. Caso alguém tenha uma ideia de logica/abordagem diferente, ficarei feliz em testar! Desde já, agradeço o auxilio de todos! ############################################################################# rodar_bibliotecas = function(necessarias){ if(necessarias){ library(rgrass7) library(spatstat) library(maptools) library(shapefiles) library(foreign) library(magrittr) library(formattable) library(e1071) library(rlang) library(rgdal) library(rgeos) library(raster) library(sp) library(sf) library(RcppCNPy) library(deldir) library(dismo) library(dplyr) library(ggplot2) library(gstat) library(tidyverse) library(smoothr) library(viridisLite) library(gpclib) } return('Ok') } ############################################################################# rodar_bibliotecas(TRUE) xpol = c(0,1,1,0,0) ypol = c(0,0,1,1,0) pol = SpatialPolygons(list( Polygons(list(Polygon(cbind(xpol-1.05,ypol))), ID="x1"), Polygons(list(Polygon(cbind(xpol,ypol))), ID="x2"), Polygons(list(Polygon(cbind(xpol,ypol - 1.0))), ID="x3"), Polygons(list(Polygon(cbind(xpol + 1.0, ypol))), ID="x4"), Polygons(list(Polygon(cbind(xpol+.4, ypol+.1))), ID="x5"))) plot(pol) ############################################################################# #nao precisava ter feito isso, mas fiz para deixar o exemplo no mesmo formato com que meus dados estao vindo areas = length(pol) lista_vazia = list() for (area in 1:length(areas)) { lista_vazia[area] = pol[area] } lista = lista_vazia lista.1 = lista_vazia #tentativa de loop para unir poligonos que se encostam # for (i in 1:length(lista)) { while (i < length(lista)) { if(i >= length(lista.1)){ break() }else{ j = 2 while (j < length(lista.1)) { # for (j in 2:length(lista.1)) { if(!is.na((over(lista[[i]], lista.1[[j]], returnList = F) == 1))) { lista[[i]] = gUnion(lista.1[[i]], lista.1[[j]]) lista.1[[i]] = gUnion(lista.1[[i]], lista.1[[j]]) #lista[[j]] = NULL lista.1 = lista.1[-j] lista = lista[-j] }else{ j = j+1 } } i = i+1 } } ############################################################################# fim = lista.1[[1]] plot(fim) for (x in 2:length(lista)){ if(x==2){ if(is.na((over(fim, lista[[x]], returnList = T) == 1))){ fim = bind(fim, lista[[x]]) } }else{ a = is.na((over(fim, lista[[x]], returnList = T) == 1)) if(a[length(a)] == TRUE){ fim = bind(fim, lista[[x]]) } } } #o objetivo final seria ter todas as areas da lista unidas num unico objeto #dessa forma, o contorno do polt(pol) devera ser igual ao contorno do plot(fim) plot(pol) plot(fim) Yury Duarte Engenheiro Agrônomo - ESALQ/USP

1 0

Tinn-R: nova página
by Jose Claudio Faria 12 Jan '19

12 Jan '19

Prezados, O projeto Tinn-R está com nova web page: http://nbcgib.uesc.br/tinnr/en/index.php Abraços, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)99966.9100 - VIVO ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ If you have software to deal with statistics, you have arms; if you have good software, you have arms and legs; if you have software like R, you have arms, legs and wings... the height of your flight depends only on you!

1 0

Linha de tendência em xyplot (lattice): melhor ajuste
by Marcelo Laia 12 Jan '19

12 Jan '19

Colegas, Fiz o seguinte: xyplot(Clorofila~Tempo|Isolado, groups=Genotipo, data=a, type = c("p", "spline", "spline"), distribute.type = FALSE, auto.key=TRUE) Na ajuda do panel.xyplot consta que: "smooth"’ adds a loess fit "spline"’ adds a cubic smoothing spline fit "a"’ has the effect of calling ‘panel.average’, which can be useful for creating interaction plots Minha pergunta é: Nos dois primeiros casos, tem como eu saber, diretamente no gráfico, qual seria o melhor ajuste? Plotar algo como um R2, por exemplo? (Se falei besteira, aqui, perdão!) Obrigado -- Marcelo

2 1