Fevereiro 2018 - 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

Diferença resultados agricolae::HSD.test e lsmeas
by Fernando Antonio de souza 05 Nov '18

05 Nov '18

Pessoal, poque o resultado dos testes HSD.test e lsmeans diferem neste exemplo? O resultado fornecido pelo lsmeas está em acordo com o material de onde tirei esse exemplo e que havia sido rodado no SAS. Já o HSD.test difere tanto no valor das médias estimadas como na significancia. install.packages(c("agricolae",lsmeans")) library(agricolae) library(lsmeans) COVARIANCIA<-read.csv(" https://www.dropbox.com/s/0uvwqnnvsdsvphc/Covariancia.csv?raw=1",head=TRUE) modeloCOVAR<-lm(GANHO~PESOINIC+DIETA,data=COVARIANCIA) anova(modeloCOVAR) HSD.test(modeloCOVAR,"DIETA",console=TRUE) lsmeans(modeloCOVAR,pairwise~DIETA,adjust="tukey") -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

3 2

[OFF - TOPIC] - YouAgro, a rede social do mundo agro
by Guilherme Moraes Ferraudo 28 Fev '18

28 Fev '18

Pessoal bom dia, ultimamente venho trabalhando no desenvolvimento de um app especializado em agro e gostaria de compartilhar com vcs. Segue abaixo um press release. Obrigado pela atencao. Projeto com participação de egressos de universidades do Brasil lança novo aplicativo social para mundo Agro, disponivel para Android e iOS. A tecnologia tem conectado pessoas de diferentes lugares e formações, integrando experiências. É nesse contexto colaborativo, participativo e gratuito que surgiu a comunidade YouAgro. Com apenas 3 meses de existência, lançada em dezembro de 2017, a proposta já tem indicadores digitais interessantes. Até o momento são mais de 40 fóruns de discussão e mais de 1000 participantes, integrando toda a cadeia do agronegócio. Participam dos fóruns desde gestores corporativos de grandes empresas, consultores, técnicos agrícolas, estudantes e produtores. A proposta do YouAgro é compartilhar experiências de forma simples, objetiva e dinâmica com foco em AgroTech. Os integrantes da equipe do YouAgro, convidam você para participar dessa comunidade e conhecer novas alternativas para aumentar o valor de seu trabalho com menos recursos, compartilhando experiências similares e informações imparciais. A conversão do conhecimento em tecnologia para o campo é uma das alternativas para gerar vantagem competitiva no agronegócio brasileiro. Saiba mais sobre YouAgro: youagro.com facebook.com/youagrosocial/ instagram.com/youagro Segue algumas matérias publicadas do YouAgro em diferentes midias: http://sna.agr.br/agronegocio-ganha-aplicativo-social-inedito/ https://www.youtube.com/watch?v=KHLviOBxeMA&feature=youtu.be&app=desktop http://podcast.unesp.br/radiorelease-19022018-egressos-da-unesp-desenvolvem- aplicativo-gratuito-dedicado-ao-agronegocio http://www.mtitecnologia.com.br/novo-aplicativo-social-para-mundo-agro/ http://www.ambientelegal.com.br/novo-aplicativo-social-para-mundo-agro/ http://www.unesp.br/portal#!/noticia/31757/novo-aplicativo- social-para-mundo-agro/ -- Guilherme Moraes Ferraudo São Paulo - SP http://lattes.cnpq.br/2096118558794430

1 0

Empilhar arquivos de um diretório em arquivo único
by Yury Duarte 26 Fev '18

26 Fev '18

Boa tarde colegas listeiros, como vão? Estou buscando unir diversos arquivos .txt de determinados diretórios num arquivo único. Todos os diretórios tem o mesmo número de arquivos (10) e cada um deles tem o mesmo numero de linhas e de colunas (36 x 23), de forma que o arquivo final gerado deveria ser um objeto de dimensão 360 x 23. Entretanto, em algum ponto, o comando não funciona de modo correto e o arquivo final gerado nunca bate com o desejado. Segue o código: (*em anexo, seguem os arquivos de um dos diretórios para teste*) rm(list = ls()) raiz = "C:\\Users\\Yury\\Desktop\\" x = "NASA" setwd(paste0(raiz, x)) #--------------------Juntar Arquivos num Unico .txt--------------------# for(i in 1:length(dir(paste0(raiz, x)))){ A1 = read.table(dir(paste0(raiz, x))[i], header = T, sep = "\t"); head(A1) if(i == 1){ write.table(A1, paste0(raiz, x, '\\', x,'.txt'), sep = "\t", append = T, col.names = T, row.names = F) } else{ write.table(A1, paste0(raiz, x, '\\', x,'.txt'), sep = "\t", append = T, col.names = F, row.names = F) } print (paste0(i, ' ok')) } Desde já, agradeço pela ajuda! Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 4

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