Novembro 2017 - 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

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

Peso espacial em X usando spreg- regressão espacial
by Edimeire Alexandra Pinto 21 Nov '17

21 Nov '17

Boa tarde, pessoal. Alguém por acaso sabe com se trabalha com modelos "parecidos" com CLIFF-ORD usando spreg do pacote sphet? Sei que o Stata faz, mas gostaria de fazer no R. Estou trabalhando com vários tipos de modelos espaciais (na realidade, quase todos do pacote sphet) e esse me gerou muita dúvida, pois preciso de que uma das variaveis explicativas com e sem peso espacial da matriz W. Por exemplo, acrescentando peso espacial a uma das variaveis explicativas do modelo. library(spdep)data(columbus)listw<-nb2listw(col.gal.nb)res1<- spreg(CRIME~HOVAL + INC+ INC, data=columbus , listw= listw, het = TRUE, verbose = FALSE,model="sarar")res2<- spreg(CRIME~HOVAL + INC+ WINC, data=columbus , listw= listw, het = TRUE, verbose = FALSE,model="sarar") ### WINC é o peso espacial multiplicado #pelo variavel INC, especie de W*INC ou mesmo res3<- spreg(CRIME~HOVAL + INC+WINC, data=columbus , listw= listw, het = TRUE, verbose = FALSE,model="error")

2 3

Reshape
by Edson Lira 17 Nov '17

17 Nov '17

Bom dia caros amigos, estou trabalhando com a base de dados abaixo(somente 6 linhas): id UF x2002 x2003 x2004 x2005 x2006 x2007 x2008 x2009 x2011 x2012 x2013 x2014 1 1 Acre 192 189 311 311 333 327 338 354 357 381 337 373 2 2 Alagoas 1251 1269 1292 1357 1413 1410 1355 1401 1344 1391 1429 1417 3 3 Amapá 204 198 239 262 252 294 317 299 302 311 332 355 4 4 Amazonas 958 1050 1433 1528 1502 1488 1531 1673 1689 1720 1771 1908 5 5 Bahia 6885 6901 7167 7336 7292 7338 7665 7792 7602 7416 7530 7879 6 6 Ceará 3705 3871 3944 4120 4143 4138 4314 4380 4121 4184 4215 4319 Estou usando o reshape para tentar transformar as linhas em colunas, ou seja, cada uma UF seria uma coluna, e assim com as demaisComo gostaria que ficasse: Acre Alagoas Amapá Amazonas Bahia Ceará ano 192 1251 3705 2002 189 1269 3871 2003 311 1293 3944 2004 ... .... Estou usando a rotina: pea<-reshape(pea1, varying=c("x2002","x2003","x2004","x2005","x2006","x2007", "x2008","x2009","x2011","x2012","x2013","x2014"), v.names="Medida", timevar="UF", times=c("Acre","Alagoas","Amapá","Amazonas","Bahia", "Ceará","Distrito Federal","Espírito Santo", "Goiás","Maranhão","Mato Grosso","Mato Grosso do Sul", "Minas Gerais","Pará","Paraíba","Paraná"," Pernambuco", "Piauí","Rio de Janeiro","Rio Grande do Norte", "Rio Grande do Sul","Rondônia","Roraima", "Santa Catarina","São Paulo","Sergipe","Tocantins"), new.row.names=1:27, direction="wide") que está me dando o erro: Erro em varying[, i] : número incorreto de dimensões Alguém tem alguma sugesstão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

5 5

usar separador de milhar
by Adriele Giaretta Biase 16 Nov '17

16 Nov '17

Bom dia a todos, tenho uma tabela, em que os números não estão separados por milhar, uma tabela com muitas linhas, mas para exemplificar, tomei como exemplo essa tabela pequena abaixo: Ex: teste metodo_1 metodo_2 A 30000000 530000000 B 40000034 550003600 C 30000020 570000002 D 50000000 830000000 E 20000000 730000000 F 30000400 530000500 Eu faço graficos de barras usando cada linha de forma que os valores apareçam em cima de cada barra. A dúvida é como fazer com que esses valores (da tabela) fiquem separados por milhar. A minha intenção é fazer tanto em padrão ingles como em padrão português. Exemplo, tomando o padrão português, o teste A da tabela acima, no metodo1 queria que o R lesse 3.000.0000 e no metodo2 53.000.0000 Se tomasse o padrão ingles, o teste A no metodo1 seria 3,000,0000 e no metodo2 53,000,0000 Porém, a tabela não tem nem ponto e nem virgula para poder trocar de uma nomenclatura para outra.

2 1

Corrigir ausência de dados
by Edson Lira 15 Nov '17

15 Nov '17

Caros amigos, tenho uma série de dados de pea e pnea que têm uma ausência de dados no ano de 2010. Como faço para corrigir essa ausência? TrêsDuas dúvidas: 1. Como organizar a planilha visto que devo corrigir a ausência por estado para corrigir a ausência de pea e pnea para o ano de 2010 por estado?2. Qual a técnica que me aconselham a usar?3. E como fazê-lo com uma única rotina? Abaixo o link da planilha Download pea pnea 02 14.csv | | | Download pea pnea 02 14.csv | | | [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

4 9

Gráfico com médias e erro agrupando var. preditora através de duas categorias
by Carlos Peixoto Neto 13 Nov '17

13 Nov '17

Boa tarde a todos. Peço um auxílio para o caso : uma espécie de peixe dividida em duas categorias de tamanho "grande" e "pequeno" e dentro dessas categorias obtenho variáveis de interação social métrica em segundos "Peixe isolado ", "peixe interação pequeno ", "peixe interação grande " estas variáveis tenho tanto para o grupo de peixes "grande" como para o "pequeno". Desejo ver o efeito na variável "Altura coluna d'água (metros)". Como representar estes dados graficamente no R, de forma que no gráfico o eixo x tenha as var. de interação social e para cada uma apareça a representação de média e erro padrão para categoria de tamanho "grande" e "pequeno" ? Segue link com médias simuladas <https://drive.google.com/open?id=1c7EkAl0s7abiHdUdbTWrmKTowsz9z0GV>. Grato.

2 1

Correlação entre variáveis categoricas
by Elias Carvalho 10 Nov '17

10 Nov '17

Boa Tarde Pessoal Tenho um banco de dados com 9 variáveis categóricas binárias e 12 variáveis categóricas ordinais. Preciso saber se há correlação entre elas, preciso de um valor que diga se há uma correlação ou associação baixa, média ou alta. O objetivo é pré-processar um banco de dados antes de criar um grafo de rede (path analysis e rede bayesiana) e eliminar os relacionamentos de variáveis com baixa corrrelação para evitar conexões causais espurias. Lendo alguns materiais foi recomendado o seguinte: - Para variáveis categoricas x categoricas usar o qui-quadrado para determinar se elas são independentes ou não.Sendo o valor P considerado a medida de correlação - Calcular também o v de Crammer (quanto menor maior é a força da correlação) Alguuém pode me confirmar se esse é o melhor método ? Ou sugere outra coisa ? -- *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"*

2 2