- R-br - DINF UFPR

Re: [R-br] Converter Shape de Pontos para Shape de Polígonos (Yury Duarte)
by Maurício Sangiogo 23 Ago '18

23 Ago '18

Não seria mais fácil tu fazer essa conversão pelo SIG? Eu faço essa conversão no QGis, mas outros programas da área também fazem essa conversão facilmente. Att, Maurício Sangiogo Obter o Outlook para Android<https://aka.ms/ghei36> ________________________________ From: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> on behalf of r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Sent: Wednesday, August 22, 2018 2:23:48 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 92, assunto 14 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. Converter Shape de Pontos para Shape de Polígonos (Yury Duarte) ---------------------------------------------------------------------- Message: 1 Date: Wed, 22 Aug 2018 14:23:33 -0300 From: Yury Duarte <yurynepomuceno(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Converter Shape de Pontos para Shape de Polígonos Message-ID: <CADFS6E5bHfa1JhUawhruRXnRYCFB_e0UEdJwdey-WPLLv_DVyA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Boa tarde colegas listeiros, como vão? Estou começando a trabalhar com dados georreferenciados e geoestatistica no R e me deparei com uma dificuldade que não estou conseguindo contornar. Tenho um shapefile de pontos que eu gostaria de transformar em um shape de polígonos e não estou sabendo como fazer essa conversão. Depois de fazer essa transformação de pontos para polígonos, queria saber se também é possível eu unir os polígonos adjacentes quando seus valores forem iguais, afim de reduzir o numero total de polígonos do meu shapefile. Desde já, agradeço pela ajuda e atenção de todos! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180822/648d2750/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: X.rar Tipo: application/octet-stream Tamanho: 142643 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180822/648d2750/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 92, assunto 14 *****************************************

1 0

1. Variáveis intervalares nos modelos mistos do R
by Thiago Wendling gonçalves de oliveira 20 Ago '18

20 Ago '18

Assunto: Variáveis intervalares nos modelos mistos do R Bom dia, meu nome é Thiago e estou precisando de ajuda na avaliação de alguns dados. Quero modelar a densidade básica da madeira de Eucalyptus de acordo com variáveis como: Espécies (3 espécies) Idade (contínua) Tipo de Solo (variável categórica - 7 classes) Textura do Solo (variável intervalar - 3 intervalos)Classes de Altitude (variável intervalar - 5 intervalos)Precipitação (variável inteira)Município (variável categórica - várias classes) e Estado (variável categórica - 3 classes). Eu estudei e penso em utilizar modelos mistos não lineares (do pacote nlme) visto que o comportamento da densidade básica perante as variáveis independentes não segue uma função linear, contudo me deparo com alguns problemas: 1) O primeiro é as variáveis intervalares pois não sei como trabalhar com elas no modelo.2) O segundo é que variáveis escolher para os efeitos aleatórios, eu primeiramente pensei no Estado e dentro dele o Município. Alguém poderia me auxiliar em algum caminho que poderia tomar para resolver esse problema? Quem sabe modelos mistos não sejam o caso e eu poderia utilizar alguma técnica de inteligência artificial, visto que, atualmente, a presença de variáveis intervalares e categóricas é comum em muitas bases de dados. Obrigado. Thiago Wendling Gonçalves de Oliveira Engenheiro FlorestalMestre em Recursos Florestais (ESALQ/USP)Doutorando em Engenharia Florestal (UFPR)Skype: Thiago W.G. de Oliveira+ 55 41 98858-2031

2 1

Ler arquivo JSON com diversos outros JSONs dentro dele
by Yury Duarte 17 Ago '18

17 Ago '18

Boa tarde colegas listeiros, como vão? Ultimamente venho encontrando certa dificuldade para trabalhar com arquivos JSON que contém outros JSON dentro deles. Sempre que precisei trabalhar com esse tipo de arquivo utilizei a expressão: library(jsonlite) raiz= ('C:\\Users\\Desktop\\') a = fromJSON(paste0(raiz,'document.json')) Dessa forma eu obtenho uma lista e consigo selecionar os valores dos campos que quero. Mas quando tento aplicar a mesma lógica para um JSON que contém diversos outros JSONs, o código quebra e retorna: Error in parse_con(txt, bigint_as_char) : parse error: trailing garbage 99809265137 } }, { "_index": "ch (right here) ------^ Meu interesse é acessar apenas as informações de alguns campos desses JSONs e construir uma tabela, onde nas colunas eu terei os campos que desejo e nas linhas eu terei os valores de cada um dos JSONs. A seguir, estou colando no corpo da mensagem um pedaço de um desses arquivos que estou tentando lidar no R. Desde já, agradeço pela ajuda e atenção. Att { "_index": "busca", "_type": "data_point", "_id": "interpolation:699646|mgper-0-20|2018-06-22|5x21sxq1km3b", "_score": 1, "_source": { "name": "Mg% 0-20", "data_type": "interpolation", "attribute": "Mg%", "farm_id": 22706, "depth": "0-20", "location": [ -48.91500982855756, -22.293821593345108 ], "inceres_id": 699646, "date_created": "2018-06-22T14:06:13.587631-03:00", "slug": "mgper-0-20", "quantity": 10.565964698791504 } }, { "_index": "busca", "_type": "data_point", "_id": "interpolation:701809|ph-20-40|2018-06-22|5x21uc94cc55", "_score": 1, "_source": { "name": "ph 20-40", "data_type": "interpolation", "attribute": "ph", "farm_id": 22706, "depth": "20-40", "location": [ -48.90241501215151, -22.28980716917012 ], "inceres_id": 701809, "date_created": "2018-06-22T14:06:35.183784-03:00", "slug": "ph-20-40", "quantity": 4.381259918212891 } }, { "_index": "busca", "_type": "data_point", "_id": "interpolation:699646|mgper-0-20|2018-06-22|5x21sz1rcjp5", "_score": 1, "_source": { "name": "Mg% 0-20", "data_type": "interpolation", "attribute": "Mg%", "farm_id": 22706, "depth": "0-20", "location": [ -48.915268161911925, -22.289490533842073 ], "inceres_id": 699646, "date_created": "2018-06-22T14:06:13.587631-03:00", "slug": "mgper-0-20", "quantity": 8.669662475585938 } }, { "_index": "busca", "_type": "data_point", "_id": "interpolation:701809|ph-20-40|2018-06-22|5x21uc3097d7", "_score": 1, "_source": { "name": "ph 20-40", "data_type": "interpolation", "attribute": "ph", "farm_id": 22706, "depth": "20-40", "location": [ -48.904161243141765, -22.289829735373146 ], "inceres_id": 701809, "date_created": "2018-06-22T14:06:35.183784-03:00", "slug": "ph-20-40", "quantity": 4.352087497711182 } } Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 8

Como aplicar um modelo em novos dados
by Edmar Caldas 13 Ago '18

13 Ago '18

Boa Tarde! Preciso aplicar um modelo em novos dados. Esse modelo é supervisionado. Novos dados não tem variável target, como proceder nesse caso, fazendo merge dos arquivos? Se fizer merge vai somar cerca de 10 milhões de registros. Edmar

2 1

Fatorial duplo e teste de tukey
by Geraldo Cabral Gouveia 13 Ago '18

13 Ago '18

Bom dia pessoal, estou precisando de uma ajuda Tenho que rodar a estatística de um fatorial duplo em tukey e para as interações quando significativas preciso do desdobramento Meus dados estão tabulados no excel Não tenho idéia de como proceder pois estou começando a trabalhar agora com o R. O fatorial é presença e ausência de um produto aplicado em 4 genótipos de arroz com 3 repetições para cada tratamento e 8 variáveis Obrigado Obter o Outlook para Android<https://aka.ms/ghei36>

2 1

Curvas não lineares - qual o melhor modelo?
by Thallita Grande 11 Ago '18

11 Ago '18

Thallita Grande compartilhou arquivos do OneDrive com você. Para exibi-los, clique nos links abaixo. <https://1drv.ms/u/s!AppBtJtj9_6eoi_LWHEfO01CSJ2u> [https://r1.res.office365.com/owa/prem/images/dc-png_20.png]<https://1drv.ms/u/s!AppBtJtj9_6eoi_LWHEfO01CSJ2u> IICxnativa_p1_2017 1.png<https://1drv.ms/u/s!AppBtJtj9_6eoi_LWHEfO01CSJ2u> <https://1drv.ms/u/s!AppBtJtj9_6eojB0ldbGtpxKiffC> [https://r1.res.office365.com/owa/prem/images/dc-png_20.png]<https://1drv.ms/u/s!AppBtJtj9_6eojB0ldbGtpxKiffC> IIC_fracoes_p1_2017_area 2.png<https://1drv.ms/u/s!AppBtJtj9_6eojB0ldbGtpxKiffC> Tenho essas 4 curvas não lineares entre os índices de conectividade e a porcentagem de área nativa remanescente. Quais são os modelos mais adequados para o ajuste das curvas nas imagens? Meu palpite é que para os 3 gráficos da primeira figura, seria um modelo logístico ou quadrático-platô e na outra figura, sera exponencial ou Michaelis-Manten. Thallita Grande

3 2

Dados de plano altimetria
by Ze Henrique 10 Ago '18

10 Ago '18

Bom dia Senhores, estou tendo que construir uma uma mapa com as altitudes dos municípios brasileiros e não estou conseguindo onde encontrar está informação para que consiga utilizar no R. Alguém conhece onde posso conseguir está informação. Desde já agradeço a todos. -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

3 2

ANOVA com medidas repetidas no tempo
by Amanda Martini 09 Ago '18

09 Ago '18

Olá pessoal, Sou relativamento nova no e R e estou trabalhando em análises com medidas repetidas, porém estou insegura com o código que estou utilizando. Nós trabalhamos com avaliação de desempenho animal em diferentes períodos, sendo que nossas repetições são os animais que são aleatorizados nos tratamento e coletadas as informações em diferentes períodos. Para montar o código utilizei o pacote lme4 considerei como fixos: trat e período , aleatório: animais, ficando desse forma: a1 = lme(desempenho ~ trat * per, random = ~1|animal, method="REML",data=desempenho). Estou no caminho certo? Alguém tem alguma sugestão ? Outra duvida é: como eu indico que esses meus animais é dentro do tratamento? Montei o código para indicar essa questão da seguinte forma: a2= lme(desempenho ~ trat * per, random = ~1|animal/trat, method="REML",data=desempenho). Está certo ou estou seguindo um raciocínio errado? Além disso, como descubro qual a matriz de covariância melhor para os meus dados e onde indico no codigo? No SAS sempre realizamos essa analise para posteriormente fazer a analise de variância. Agradeço desde já, Att, Amanda Martini Mestranda em Produção Animal – Forragicultura Laboratório Pastos e Suplementos

1 0

Formatação de summary table no RStudio apresentando &nbsp
by Elias Carvalho 05 Ago '18

05 Ago '18

*Estou testando o pacote "qwraps2", então estou rodando o exemplo em: https://github.com/r-lib/xml2/issues/151 <https://github.com/r-lib/xml2/issues/151>, conforme abaixo.* set.seed(42) library(dplyr) ## ## Attaching package: 'dplyr' ## The following objects are masked from 'package:stats': ## ## filter, lag ## The following objects are masked from 'package:base': ## ## intersect, setdiff, setequal, union library(qwraps2) # define the markup language we are working in. # options(qwraps2_markup = "latex") is also supported. options(qwraps2_markup = "markdown") data(mtcars) mtcars2 <-dplyr::mutate(mtcars, cyl_factor = factor(cyl,levels = c(6, 4, 8), labels = paste(c(6, 4, 8), "cylinders")), cyl_character = paste(cyl, "cylinders")) our_summary1 <- list("Miles Per Gallon" = list("min" = ~ min(mpg), "max" = ~ max(mpg), "mean (sd)" = ~ qwraps2::mean_sd(mpg)), "Displacement" = list("min" = ~ min(disp), "max" = ~ max(disp), "mean (sd)" = ~ qwraps2::mean_sd(disp)), "Weight (1000 lbs)" = list("min" = ~ min(wt), "max" = ~ max(wt), "mean (sd)" = ~ qwraps2::mean_sd(wt)), "Forward Gears" = list("Three" = ~ qwraps2::n_perc0(gear == 3), "Four" = ~ qwraps2::n_perc0(gear == 4), "Five" = ~ qwraps2::n_perc0(gear == 5)) ) our_summary1 ### Overall summary_table(mtcars2, our_summary1) *O código foi executado com sucesso no RStudio, porém a tabela aparece como:* [image: Screenshot from 2018-08-04 18-29-38.png] *Como faço para aparecer a tabela formatada como no exemplo:* [image: Screenshot from 2018-08-04 18-31-38.png] -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

3 2

Dados ausaentes (NA)
by Edson Lira 03 Ago '18

03 Ago '18

Boa tarde, a todos, estou com um problema com essa tabela, dentre outras.Vejam que há 38 observações em branco (seriam NAs). $DispSecNao Qtde Perc 38 79.17DispAcionamento, DispParadaEmergencia, DispInterTrav 1 2.08DispParadaEmergencia 1 2.08DispParadaEmergencia, DispInterTrav 7 14.58DispPartida, DispParadaEmergencia, DispInterTrav 1 2.08Total 48 100.00 Usei a rotina abaixo para eliminar os NAs. b<-read.csv2("C:\\Users\\Edson Lira\\Google Drive\\CONSULTA\\CLAUDIO_SILVIO\\banco1.csv")b[b==""]<-NA Para outras tabelas que existem registro em branco, o NA foi eliminado das frequências, exceto na tabela acima. Alguém tem alguma sugestão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 3