Julho 2020 - R-br - DINF UFPR

[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

Duvidas Análise de Componentes Principais
by Fernando Souza 09 Dez '20

09 Dez '20

Olá pessoal tudo bem Estou realizando uma análise de componentes principais, com o objetivo de selecionar as melhores variáveis para uma análise multivariada. Gostaria de uma orientação sobre a interpretação. Eu interpretei de forma correta? Variáveis que possuem alta correlação com um CP também estão correlacionadas entre si. Isso significa que posso utilizar somente uma das principais variáveis do CP1 (por exemplo) para incluir no modelo multivariado? No caso a variável que apresentar significado para o trabalho. Os círculos de correlação, screem plot e corrplot geradas pelo Pca estão disponíveis no link https://drive.google.com/open?id=1jS5tCX-2-qHo2Gbej00DkFQBpLztwXks A interpretação que fiz foi a seguinte: Considerei as variáveis PC,IC e Prod leite (Azul no gráfico) como variável suplementar quantitativa. Baseado no ScreenPlot selecionei somente 5 CP para avaliação que explicaram 67,80% da variancia total. Pela minha interpretação cheguei as seguintes conclusões Principais variáveis para cada componente: Consumos (CP1): CMS,CMO,CPB,CEE,CCNF,CFDNcp,CFDA,CNDT,CEB,CED,CEM,CEL Energia (CP2) : EB,ED,EM,EL Peso (CP3): PV,PCJ,PCVZ,PCVZ075 Composição do Leite(CP4) : Solido e DCNF Variaveis sanguíneas (CP5): PT e ALBUMINA Segue uma descrição do trabalho para compreensão. Aceito sugestões sobre a melhor forma trabalhar esses dados O banco de dados é oriundo dois estudos realizados pelo mesmo grupo de pesquisa nas mesmas condições experimentais (animais,instalação,dieta,manejo) e avaliaram 48 variáveis. Os fatores experimentais foram: *Raça*: Estudo 1 avaliou 12 animais (6 raça1 e 6 da raça 2). O estudo 2 avaliou 13 animais (6 raça 1 e 7 raça 2) *Plano nutricional:* O estudo 1 avaliou os níveis (0,15,30) o estudo 2 avaliou os níveis (0,5,10,20) *Dias em lactação-DEL* (medidas repetidas no tempo): Estudo 1, dias dispersos entre os períodos de 30 a 120 dias. No estudo 2 os DEL foram avaliados em dias entre 150 a 250 dias. Eu realizei uma análise gráfica e não há indicações de haver diferenças entre estudos, o que era de certo modo esperado. Raça e DEL parecem ser os fatores que afetam as variáveis respostas. -- ========================================= 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/ ==========================================

6 14

Explicação gráfico e resultados da PCA
by Chiara Lubich 05 Ago '20

05 Ago '20

Gerei os resultados da PCA, porém estou com dificuldades para explicação dos resultados e do gráfico. Alguém poderia me dar um norte ou sugerir algum livro/artigo que me ajude na interpretação? Desde já agradeço pela atenção Abraços

3 3

Manipulação de dados
by Cid Póvoas 31 Jul '20

31 Jul '20

Boa tarde, Eu tenho um banco de dados com 19 colunas sendo uma de data e 18 colunas com 67308 linhas e gostaria de saber se alguém sabe como editar esse dataframe quando valor de uma linha pra outra variar mais q X, repetir o valor da última linha, e continuar assim para todas as colunas. Na minha situação variar mais que 5 de uma linha pra outra seria uma anomalia do sensor. Banco de dados para exemplo df <- structure(list(a = c("A", "A","B","B","C","C","A", "A","B","B","C","C"), b = c(1,2,1,4,3,10,2,3,1,7,9,-10), c = c(5,1,2,3,3,1,2,3,9,7,-5,10), d = c(3,10,1,9,2,1,2,3,3,5,-5,9)), .Names = c("nome","b","c","d"), class = "data.frame",row.names = c(NA, -12L))

3 2

histograma no R, como fazer?
by Ana Maria 31 Jul '20

31 Jul '20

Bom dia! Alguém sabe como eu faço um histograma desse tipo abaixo no R? As classes que quero inserir são: [10,30), [30,50) e ≥50. Já tentei diversas vezes e não consigo criar as classes que eu desejo no eixo x. [image: image.png] Segue anexo os dados que eu tenho para criar o histograma, caso alguém possa me ajudar eu agradeço desde já. Att., Ana Maria

4 9

Re: [R-br] Digest R-br, volume 115, assunto 5
by Andrea Vita Reis Mendonca 29 Jul '20

29 Jul '20

7.082E-08 significa: 7,082X10^-8 Em qua., 29 de jul. de 2020 às 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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. Re: Fwd: =?utf-8?Q?D=C3=BAvidas_?=sobre leitura de "p". > (Fernando Souza) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 29 Jul 2020 09:39:54 -0300 > From: Fernando Souza <nandodesouza(a)gmail.com> > To: > > "=?utf-8?Q?a_lista_Brasileira_oficial_de_discuss=C3=83=C2=A3o_do_programa_R.?=" > > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Fwd: =?utf-8?Q?D=C3=BAvidas_?=sobre leitura de > "p". > Message-ID: <489AE422-ED50-4989-BFFE-B9C075CAF47B(a)getmailspring.com> > Content-Type: text/plain; charset="utf-8" > > p= 7.082E-08 é o mesmo que 7.082x(10^-08) ou seja 0.00000007082 > na prática p<0.001 > > On Jul 27 2020, at 5:49 pm, Ana Maria por (R-br) <r-br(a)listas.c3sl.ufpr.br> > wrote: > > > > ---------- Forwarded message --------- > > De: Ana Maria <moreirafernandesanamaria(a)gmail.com (mailto: > moreirafernandesanamaria(a)gmail.com)> > > Date: seg, 27 de jul de 2020 17:45 > > Subject: Dúvidas sobre leitura de "p". > > To: <r-br-request(a)listas.c3sl.ufpr.br (mailto: > r-br-request(a)listas.c3sl.ufpr.br)> > > > > > > > > Boa tarde. > > > > Alguém sabe me dizer como eu leio esse resultado de p= 7.082E-08? O > mesmo foi gerado após um teste de normalidade que eu fiz nos meus dados. O > que significa E-08? Esse valor de p é maior que 0.05? > > > > Agradeço desde já. > > Ana Maria > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200729/faeec36e/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 115, assunto 5 > ***************************************** >

2 1

Fwd: Dúvidas sobre leitura de "p".
by Ana Maria 29 Jul '20

29 Jul '20

---------- Forwarded message --------- De: Ana Maria <moreirafernandesanamaria(a)gmail.com> Date: seg, 27 de jul de 2020 17:45 Subject: Dúvidas sobre leitura de "p". To: <r-br-request(a)listas.c3sl.ufpr.br> Boa tarde. Alguém sabe me dizer como eu leio esse resultado de p= 7.082E-08? O mesmo foi gerado após um teste de normalidade que eu fiz nos meus dados. O que significa E-08? Esse valor de p é maior que 0.05? Agradeço desde já. Ana Maria

2 1

Re: [R-br] Fwd: Dúvidas sobre leitura de "p".
by sznelwar＠uol.com.br 27 Jul '20

27 Jul '20

2 1

GLM com alta proporção de zeros
by Emerson Cotta Bodevan 23 Jul '20

23 Jul '20

Prezados, boa tarde. Estou com dificuldade para cálculo de Razão de Prevalência, com respectivos IC95% para meus dados. Meus dados estão com a seguinte estrutura: Variável resposta: - ocorrência de AVC (0: Não, 1: Sim) Variáveis independentes: - Faixa Etária (18 a 28 - referência, 29 a 39, 40 a 59, 60 a 79 e 80 ou mais) - Escolaridade (Nunca estudou - referência, Ens Fun, Ens Med, Ens Sup) - Saúde (Excelente - referência, Muito Boa, Boa, Regular, Ruim) - Número Doenças Crônicas (Nenhuma - referência, Uma, Duas ou mais) - Automedicação (Não - referência, Sim) Importante dizer que tenho alta proporção de zeros na variável resposta: 371 zeros e 11 uns. *Tentei inicialmente a regressão de poisson da seguinte forma:* library(sandwich) library(lmtest) fit.poisson=glm(AVC~FE+Esc+Saude+NDC+Auto,data=dados,family=poisson) *Obtive a seguinte mensagem:* glm.fit: taxas ajustadas numericamente 0 ocorreu *Vi que o problema pode ser o excesso de zeros. Então, lendo algumas postagens aqui da lista tentei o seguinte:* library(pscl) fit.hurdle<-hurdle(AVC~FE+Esc+Saude+NDC+Auto,dist="poisson",zero.dist="poisson",data=dados) *Obtive a seguinte mensagem:* Warning messages: 1: glm.fit: fitted rates numerically 0 occurred 2: glm.fit: fitted rates numerically 0 occurred 3: In value[[3L]](cond) : Lapack routine dgesv: system is exactly singular: U[7,7] = 0FALSE *Tentei então* fit.zeroinflPoisson<-zeroinfl(AVC~FE+Esc+Saude+NDC+Auto,dist="poisson",link="logit",data=dados) *Obtive a seguinte mensagem:* Warning messages: 1: glm.fit: fitted probabilities numerically 0 or 1 occurred 2: In value[[3L]](cond) : sistema é computacionalmente singular: condição recíproca número = 8.9423e-20FALSE Não consigo encontrar uma saída. Alguém pode me dar uma dica? Agradeço imensamente. Atenciosamente, *Emerson*

4 6