- R-br - DINF UFPR

FAOSTAT
by Alisson Lucrécio 26 Fev '18

26 Fev '18

Caro Colegas, Bom dia. Tentei usar o pacote FAOSTAT para baixar dados da FAO, mas não funcionar. Eu acho que tem um erro no código arableLand ou domainCode. Vocês saberiam o que eu poderia fazer para resolver esse erro? Obrigado. -- Alisson Lucrecio da Costa -------------------------- library(FAOSTAT) FAOquery.df = data.frame(varName = c("arableLand", "cerealExp", "cerealProd"), domainCode = c("RL", "TP", "QC"), itemCode = c(6621, 1944, 1717), elementCode = c(5110, 5922, 5510), stringsAsFactors = FALSE) ## Download the data from FAOSTAT FAO.lst = with(FAOquery.df, getFAOtoSYB(name = varName, domainCode = domainCode, itemCode = itemCode, elementCode = elementCode, useCHMT = TRUE, outputFormat = "wide")) FAO.lst$entity[, "arableLand"] = as.numeric(FAO.lst$entity[, "arableLand"]) -------------------------- library(FAOSTAT) > > FAOquery.df = data.frame(varName = c("arableLand", "cerealExp", "cerealProd"), + domainCode = c("RL", "TP", "QC"), + itemCode = c(6621, 1944, 1717), + elementCode = c(5110, 5922, 5510), + stringsAsFactors = FALSE) > > ## Download the data from FAOSTAT > FAO.lst = with(FAOquery.df, + getFAOtoSYB(name = varName, domainCode = domainCode, + itemCode = itemCode, elementCode = elementCode, + useCHMT = TRUE, outputFormat = "wide")) ---------------------------------------- ** FAOSTAT Data Download (3 in Total) ** ---------------------------------------- (1): Downloading variable arableLand ... The specified query has no data, consult FAOSTAT (2): Downloading variable cerealExp ... NOTE: Multiple China detected in 'Value' sanitization is performed OK (3): Downloading variable cerealProd ... The specified query has no data, consult FAOSTAT Number of variables successfully downloaded: 1 out of 3 > FAO.lst$entity[, "arableLand"] = as.numeric(FAO.lst$entity[, "arableLand"]) Error in `[.data.frame`(FAO.lst$entity, , "arableLand") : undefined columns selected

2 1

Dúvida Básica
by Maurício Soares 20 Fev '18

20 Fev '18

Eu estou revisando uma programação em R, e fiquei com essa dúvida básica. dados$Nome é uma variável com nomes de pessoas. n=length(dados$Nome) num=NA for(i in 1:n){ num[i]=sum(dados$Nome==dados$Nome[i]) } Não entendi essa linha -> num[i]=sum(dados$Nome==dados$Nome[i]). Alguém poderia me dar auxílio? Obrigado.

2 1

Experimento em parcelas subdivididas no tempo
by ana paula coelho madeira 16 Fev '18

16 Fev '18

Prezados colegas, boa tarde! Estou com um experimento em DIC no esquema de parcelas subdividas no tempo. No pacote ExpDes, temos implementado o esquema de parcelas subdivididas, porém, no espaço. Peço, por favor, sugestões de materiais e/ou scripts que possam me auxiliar. Desde já, agradeço

3 2

Experimento dados de contagem
by Maurício Lordêlo 15 Fev '18

15 Fev '18

Caros, Foi realizado um experimento para avaliar o efeito de oito diferentes tratamentos na mortalidade de insetos (foram duas variáveis respostas). O pesquisador colocou em cada unidade experimental 10 larvas e registrou o número de larvas mortas e de pupas inviáveis em quatro diferentes tempos sendo cinco repetições. Este registro no tempo não foi feito de forma independente e sim de forma acumulada. Ou seja, se o número de mortos no tempo t1 é igual 2 e no tempo t2 é igual a 5, indica entre t1 e t2 "morreram" 3 insetos (na mesma unidade experimental) . O objetivo é avaliar o efeito do tratamento e do tempo. Como não há independência do tempo, alguém tem alguma sugestão de como analisar estes dados? rm(list = ls(all=TRUE)) #Tratamentos trat = as.factor(c(rep("Trat1",20),rep("Trat2",20), rep("Trat3",20),rep("Trat4",20),rep("Trat5",20), rep("Trat6",20),rep("Trat7",20),rep("Trat8",20))) #Tempo tempo = as.factor(rep(c("t1","t2","t3","t4"), each = 5, times =8)) #Número total de larvas em cada unidade experimental n_total_larvas= rep(10,160) #Número de larvas mortas n_larvas_mortas = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 3, 0, 0, 1, 1, 3, 0, 3, 0, 0, 0, 0, 3, 1, 1, 1, 1, 5, 3, 1, 2, 2, 5, 3, 3, 4, 3, 0, 1, 0, 0, 0, 0, 2, 1, 2, 1, 0, 4, 1, 2, 1, 4, 4, 2, 2, 2, 0, 0, 1, 1, 0, 0, 0, 1, 2, 0, 2, 2, 1, 2, 0, 2, 3, 2, 2, 4, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 1, 1, 0, 0, 2, 1, 1, 0, 0, 1, 1, 1, 1, 0, 3, 1, 2, 1, 0, 3, 1, 2, 2, 1, 3, 1, 3, 2, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 2, 0, 1, 0, 1, 3, 0, 1, 8, 8, 7, 9, 8, 8, 8, 7, 9, 8, 8, 9, 7, 9, 8, 8, 9, 7, 9, 8) #Número de pupas inviáveis n_pupas_inv = c(3, 4, 3, 6, 1, 3, 8, 6, 6, 1, 10, 9, 9, 4, 6, 10, 9, 9, 7, 10, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 5, 7, 4, 6, 5, 5, 7, 7, 6, 7, 2, 1, 0, 2, 0, 2, 4, 2, 5, 1, 4, 4, 2, 6, 7, 6, 6, 8, 6, 8, 0, 3, 0, 2, 0, 4, 5, 0, 4, 3, 7, 6, 5, 5, 5, 8, 7, 7, 8, 5, 1, 0, 2, 6, 0, 3, 3, 2, 9, 9, 6, 7, 7, 9, 9, 10, 8, 9, 9, 10, 0, 1, 0, 0, 0, 5, 4, 3, 2, 3, 9, 7, 8, 4, 7, 9, 7, 9, 6, 7, 1, 2, 0, 0, 1, 6, 8, 2, 2, 1, 10, 8, 6, 9, 6, 10, 9, 7, 10, 8, 2, 1, 1, 0, 0, 2, 1, 3, 1, 1, 2, 1, 3, 1, 2, 2, 1, 3, 1, 2) p_larvas_mortas = n_larvas_mortas/n_total_larvas table(p_larvas_mortas) p_pupas_inv = n_pupas_inv/n_total_larvas table(p_pupas_inv) dados_larvas = data.frame(trat, tempo, n_total_larvas, n_larvas_mortas,n_pupas_inv,p_larvas_mortas,p_pupas_inv) str(dados_larvas) tapply(p_larvas_mortas, list(trat, tempo), mean, na.rm=T)

4 10

Salvar saídas de um loop
by Yury Duarte 10 Fev '18

10 Fev '18

Boa tarde colegas listeiros! Estou tendo um pequeno problema para salvar as saídas de um loop que fiz para preencher uma matriz com dados de diferentes arquivos. Meus arquivos base possuem informações anuais de produtividade de 15 locais para um determinado mês. Construí o loop com intenção de gerar um arquivo por local, contemplando as produtividades de cada mês (nas linhas) e de cada ano (nas colunas). Conferi o output do único arquivo salvo pelo script e os valores correspondem ao desejado, então acredito que a lógica para o preenchimento da matriz em branco esteja correta. Os arquivos para executar a simulação estão em anexo e o script está no corpo do email. Desde já, agradeço pela colaboração de todos! rm(list = ls()) raiz = 'C:\\Users\\Yury\\Desktop\\Mestrado\\1_TESE\\Model_Data\\Maize\\PREVISAO\\PA\\' estrategia = c('E1', 'E2', 'E3', 'E4', 'E5') meses = c('setembro', 'outubro', 'novembro', 'dezembro', 'janeiro', 'fevereiro', 'marco', 'abril') cidade = c("abelardoluz", "altamira", "brasilia", "caarapo", "catalao", "lagarto", "machado", "mateiros", "muitoscapoes", "piracicaba", "pontagrossa", "primaveradoleste", "srmangab", "sdesiderio", "urucui") anos = c(1980:2012) previsao = matrix(NA, nrow = length(meses), ncol = length(1980:2012)) colnames(previsao) = 1980:2012 #for(e in 1:length(estrategias)){ for(m in 1:length(meses)){ a = read.table(paste0(raiz, 'E1\\', meses[m],'\\', 'final.txt'), header = T, sep = '\t') for(city in 1:length(cidade)){ for(i in 1:length(anos)){ previsao[m,i] = a[city,i] } } write.table(previsao, paste0(raiz, 'E1\\', cidade[city],'.txt'), row.names = F, col.names = T, sep = '\t') } #} Yury Duarte Engenheiro Agrônomo - ESALQ/USP

2 1

Re: [R-br] [WordCloud] Problema com encoding
by Rafael Tieppo 09 Fev '18

09 Fev '18

Ao ler um arquivo tive os seguintes erros: A <- read.csv("tel.csv", sep = ",") Error in read.table(file = file, header = header, sep = sep, quote = quote, : empty beginning of file In addition: Warning messages: 1: In read.table(file = file, header = header, sep = sep, quote = quote, : invalid input found on input connection 'tel.csv' 2: In read.table(file = file, header = header, sep = sep, quote = quote, : line 1 appears to contain embedded nulls 3: In read.table(file = file, header = header, sep = sep, quote = quote, : incomplete final line found by readTableHeader on 'tel.csv' Com utf-16 resolveu A <- read.csv("tel.csv", sep = ",", fileEncoding="utf-16") Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: https://rafatieppo.github.io blog: https://fuidebicicleta.wordpress.com"Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". On Friday, February 9, 2018, 9:26:34 AM GMT-3, <r-br-request(a)listas.c3sl.ufpr.br> wrote: 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. [WordCloud] Problema com encoding (Max) ---------------------------------------------------------------------- Message: 1 Date: Tue, 26 Sep 2017 22:31:21 -0300 From: Max <max.lara(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] [WordCloud] Problema com encoding Message-ID: <CAK9SNW2=9j=pc3fYzFnM3ooYhN1eORufRONu-DpqN6jfR-ORhw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Pessoal, boa noite! Com base em um exemplo pego na internet, tentei fazer uma nuvem de palavras de um arquivo notepad (peguei uma pequena matéria na internet e salvei, em .txt, com o nome SaoBento). O código tem funcionado quase corretamente. O problema é que eu não tenho tido sucesso em corrigir o *encoding* do texto. - Tentei usar encoding = "UTF-8" na linha do readlines, mas sem sucesso. - Também tentei usar enc2native() na última linha do passo 7, mas ocorre erro argument is not a character vector. - Salvei o arquivo SaoBento.txt no Notepad++, usando UTF-8 e, também, a nuvem final acusou problemas de *encoding*. - Deixei o arquivo SaoBento.txt salvo sozinho em uma pasta e, também, não tive sucesso em acertar o *encoding*. Sendo assim, gostaria de saber se alguém poderia fornecer uma dica de como posso driblar esse problema. Agradeço pela atenção. Saudações, ****************************************

1 0

Assunto OFF
by Edson Lira 05 Fev '18

05 Fev '18

Caros amigos, boa noite, Estou com um problema na elaboração de um gráfico de controle. Estou analisando perdas em um período de 16 meses. Abaixo os dados. jan/16 fev/16 mar/16 abr/16 mai/16 jun/16 jul/16 2.096,23 3.532,09 1.995,22 2.294,34 6.471,51 6.161,56 3.687,96 ago/16 set/16 out/16 nov/16 dez/16 jan/17 30.469,19 52.931,60 39.768,36 21.855,01 16.399,77 17.521,40 fev/17 mar/17 abr/17 7.673,57 1.125,67 1.652,13 Vejam que de janeiro a julho de 2016, tenho uma grandeza, de agosto de 2016 a janeiro de 2017 a grandeza é bem maior e de fev de 2017 a abril a grandeza está na proporção do primeiro período. Tentei ajustar usando médias ponderadas pelas grandezas, porém ao usar a média encontrada no período divergente, ao calcular as diferenças para a média da amplitude tenho valores iguais e as diferenças zeram no período que uso a média ponderada e então fica difícil estabelecer os limites. Os limites são calculados da seguinte forma LSC=xbar +2*(MR/d2 ) d2=1,128 LSC=xbar - 2*(MR/d2) MR= soma(xi-xi-1)/m m= número de diferenças Alguém teria uma idéia de como ajustar os dados para fazer um gráfico de controle adequando? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 5

Ordenar os fatores do eixo X no plot.group()
by Rodrigo Simetti 03 Fev '18

03 Fev '18

Prezados, Alguém saberia como posso reordenar a sequência do eixo X ao utilizar o plot.group() do pacote agricolae? O padrão do comando é plotar da maior média para a menor, mas preciso plotar em sequência alfabética. Utilizando o exemplo do próprio pacotes temos como resultado a sequência oo, ff, cc e fc. Gostaria que fosse cc, fc, ff e oo. library(agricolae) data(sweetpotato) model<-aov(yield~virus,data=sweetpotato) comparison<- LSD.test(model,"virus",alpha=0.01,group=TRUE) plot(comparison) -- *Rodrigo Simetti* Engenheiro Industrial Madeireiro - UFPR Mestre em Engenharia Florestal - UFPR Doutorando em Ciência e Tecnologia da Madeira - UFLA (41) 9 8870 0587 (35) 9 9859 3966 Skype: rsimetti

4 4

Recode rotina
by Edson Lira 02 Fev '18

02 Fev '18

Bom dia caros amigos, tenho um pequeno problema que ainda não consegui encontrar a solução. É o seguinte: Estou aplicando um questionário para avaliar a satisfação do cliente em um hospital público. Na avaliação dos clientes estou usando uma escala de likert com as respostas para algumas questões. ótimo, bom, regular, ruim, péssimo e não utilizou o serviço Gostaria de fazer uma mudança nessas respostas usando uma rotina do R. Gostaria de mudar da seguinte forma para todo o banco de uma única vez. otimo e bom receberia o label "positiva" regular e não utilizou o serviço o label " "e ruim e péssimo receberia o label "negativa". Alguém tem alguma sugestão? Bom diae muito obrigado!!! Prof. Edson Lira, Me Estatístico Manaus-Amazonas

4 6

PERMDISP - interpretação dos resultados
by Carolina Lopes 31 Jan '18

31 Jan '18

Prezados. Rodei uma analise de PERMDISP (função betadisper) como pressuposto da PERMANOVA. Estou com dificuldades de interpretar os resultados. Algum material que possa ajudar no entendimento dessa analise?? Desde já agradeço a atenção Msc. Carolina A. Lopes Currículo Lattes: http://lattes.cnpq.br/9210875052648641

1 0