Outubro 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

Excluir valores da matriz de correlação
by Elias Carvalho 31 Out '17

31 Out '17

Boa noite pessoal Como posso excluir de uma matriz de correlação valores muito baixo como *esse em negrito* como abaixo G001 G002 G008 C051 C071 C101 C151 R001 R02A R02C R2B1 G001 1.00000000 0.14842780 -0.114680730 0.302947649 0.32931462 0.37740045 0.30480517 -0.34884430 -0.39049631 -0.34041558 0.056306588 G002 0.14842780 1.00000000 0.070276175 0.170691108 0.09668893 *0.02550203* 0.13058268 -0.15025183 -0.15148344 -0.18614711 0.118041910 G008 -0.11468073 0.07027618 1.000000000 *0.007789258* 0.06977162 0.01623666 -0.14927575 0.01368586 *0.03586060 0.00202448* 0.486717281 C051 0.30294765 0.17069111 0.007789258 1.000000000 0.37174418 0.26435392 0.22513174 -0.47614134 -0.48203009 -0.42861271 0.372495377 C071 0.32931462 0.09668893 0.069771618 0.371744175 1.00000000 0.45945064 0.40883404 -0.24972945 -0.22757192 -0.17468398 0.312496407 -- *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 1

Re: [R-br] Digest R-br, volume 82, assunto 19
by Magnun Rochel 29 Out '17

29 Out '17

Oi Glenio, Vc ja tentou com a função bwplot() do lattice? Ficaria algo assim: bwplot(weight ~ Litter | Treatment + sex, data= RatPupWeight) []s, Magnun 2017-10-28 16:12 GMT-02:00 <r-br-request(a)listas.c3sl.ufpr.br>: > 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. Ajuda para elaboração de boxplot em painel > (Gilenio Borges Fernandes) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 28 Oct 2017 15:11:59 -0300 > From: Gilenio Borges Fernandes <gilenio(a)ufba.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Ajuda para elaboração de boxplot em painel > Message-ID: > <CA+7i11aPG1Oi69uA8-4uz0aL9ZjbYQpOMkWou_uuUp_+ > GR3Ung(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Prezad@s > > Estou precisando de auxilio para elaborar boxplots em painel pelo menos > aproximadamente iguais aos dos arquivos pdf anexos. > > O banco de dados é do pacote nlme do R; pode ser acessado com o CMR > > library(nlme) > > data(RatPupWeight) > > head(RatPupWeight) > > names(RatPupWeight) > > Grato > > > Gilênio Borges Fernandes > Professor Associado IV (Aposentado) > Professor Adjunto A (Substituto, 03/2015 a 03/2017) > Universidade Federal da Bahia > Instituto de Matemática e Estatística > Departamento de Estatística > Av. Adhemar de Barros, s/n – Ondina. > 40.170-110 - Salvador - BA, Brasil > Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 > Skype: gilenio.fernandes > Lattes: http://lattes.cnpq.br/6764860618464860 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171028/e953defd/attachment.html> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: rat-pup-grafico1.pdf > Tipo: application/pdf > Tamanho: 80064 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171028/e953defd/attachment.pdf> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: rat-pup-Grafico2.pdf > Tipo: application/pdf > Tamanho: 87655 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171028/e953defd/attachment-0001.pdf> > > ------------------------------ > > 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 82, assunto 19 > ***************************************** > -- Magnun R. Madruga Linguista +55 11 953956550

1 0

Ajuda para elaboração de boxplot em painel
by Gilenio Borges Fernandes 28 Out '17

28 Out '17

Prezad@s Estou precisando de auxilio para elaborar boxplots em painel pelo menos aproximadamente iguais aos dos arquivos pdf anexos. O banco de dados é do pacote nlme do R; pode ser acessado com o CMR library(nlme) data(RatPupWeight) head(RatPupWeight) names(RatPupWeight) Grato Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto, 03/2015 a 03/2017) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860

1 0

Identificação de variáveis
by Elias Carvalho 25 Out '17

25 Out '17

Como posso identificar se uma variável é 'ordered factor nas seguintes condições: - o read.csv leu ela como int: - o read.csv leu ela como factor e não ordered factor: Obrigado

2 1

ajuda
by Edmar Caldas 25 Out '17

25 Out '17

Turma tem hora que funciona e hora não. alguém pode me ajudar? a variável datanasc no Rstudio está como <dttm> o que significa isso? Edmar

4 3

ANOVA in R using the function lapply ans terms
by Euriann Yamamoto 24 Out '17

24 Out '17

Hi everyone. In my dataset I´ve got 36 genotypes in 5 locations in partially balanced incomplete design 6x6, simple. We are analyzing productivity. I need to use adjusted means to do the orders analysis. TO obtain ANOVA individual to the 5 locations I run command: mi <- lapply(split(da, f=da$local), aov, formula=producao ~ repeticao/bloco + cultivar) lapply(mi, summary)#not adjusted for blocks But to transform cultivar adjusted for blocks (to obtain adjusted means) I need use command terms(formula,keep.order=TRUE) I DO NOT KNOW HOW TO USE THIS IN THE CORRECT ORDER. Something like this? m0 <- lapply(split(da, f=da$local), aov(terms(producao ~ repeticao/bloco + cultivar,keep.order = TRUE),da)) lapply(m0, summary)#adjusted to blocks Thanks Euriann

2 1