Setembro 2019 - 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

Ajuda com análise conjunta de dois experimentos com parcelas subdividas
by José Lucas Safanelli 02 Out '19

02 Out '19

Caros colegas, Peço ajuda a vocês para conseguir entender e analisar o seguinte experimento: - são dois experimentos feitos lado a lado: um com plantio direto, outro com convencional. - em cada local, existe um experimento de parcelas subdividades, sendo o fator manejo na parcela e dose de nitrogenio na subparcela. Diante disso, gostaria de saber como expressar o modelo da ANOVA no R, a fim de fazer comparação de médias, e mais pra frente, análise multivariada via discriminante canonica. Sei que são necessárias duas coisas: - Dividir o erro dentro da abordagem da análise conjunta. - Dividir o erro na parcela subdividida. Tenho utilizado a seguinte expressão, mas não estou confiante quanto a expressão. Infelizmente nunca trabalhei com experimento mais complexos como este. aov_N <- aov(N_kg_ha ~ Manejo * Rotacao * Nitrogenio + Bloco + Error(Nitrogenio/Bloco/Manejo), data = dados) summary(aov_N) Tanto o croqui quanto os dados estão neste link: https://drive.google.com/drive/folders/1rjVv9KwtE5A7y-VPXTV9Y2GlilKLXK2F?us… Agradeço se alguém puder me ajudar com alguma orientação. José Lucas.

3 4

Ajuda: Escolha de distribuição de probabilidade (GLM)
by Chiara Lubich 30 Set '19

30 Set '19

Oi, Estou tendo dificuldade para achar uma distribuição certa para minha análise (GLM)... Mesmo dados não são normais e também não apresentam homoscedasticidade, mesmo após realizar transformação (log). Qual família devo escolher? Estava trabalhando com binomial negativa, é a correta? Obrigada e aguardo retorno -- *Chiara Lubich Cardoso Furtado* Mestranda em Biologia de Água Doce e Pesca Interior - INPA Engenheira de Pesca - UFAM Telefone: (92) 99110-6657 e-mail: lubichchiara(a)gmail.com Lattes: http://lattes.cnpq.br/0673892253719514

2 1

dic parcela subdividida no tempo
by ANDREA VITA REIS MENDONCA 20 Set '19

20 Set '19

Trata-se de um experimento deduzas temperaturas e seis períodos de avaliação, em DIC. Resolvi conforme mostrado a seguir. Gostaria de saber se está correto. As analises ficaram boas. pkg= c("lattice", "latticeExtra", "doBy", "multcomp", "reshape", "plyr", "nlme") sapply(pkg, require, character.only=TRUE) source("http://dl.dropboxusercontent.com/u/48140237/bandas.R") sessionInfo() trellis.device(color=FALSE) ang=read.csv("ang.csv") str(ang) dh=melt(subset(ang, select=1:8), id.vars=1:2) d=diff(as.Date(c("2017-01-01","2017-01-16", "2017-01-31","2017-02-15", "2017-03-02","2017-03-17"))); d d =cumsum(c(15, d)); d dh$dia=d[as.integer(dh$variable)] levels(dh$variable) xyplot(value~dia|Temperatura, groups=Temperatura, data=dh, type=c("p","a"), auto.key=TRUE) dh1=transform(dh, D=factor(dia), parc=interaction(Rep, Temperatura), subp=interaction(Rep, Temperatura)) str(dh1) m2=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="ML") m3=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="REML") -- Andrea Vita Reis Mendonça Profa. Associado I UFRB/CCAAB

2 1

Não carrega o biblioteca - lmer4
by Andre Oliveira 13 Set '19

13 Set '19

boa noite,de uma hora para outra não consigo carregar a biblioteca lme4.Alguém tem alguma sugestão de ajuda? (estou usando win10) require(lmer4)Carregando pacotes exigidos: lmer4Warning message:In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package called ‘lmer4’ att., André

5 5

gráficos com bandas de confiança no R - com sombra
by Nilson Guiotoku 12 Set '19

12 Set '19

Pessoal, gostaria de fazer um gráfico conforme o anexo, cujo o intervalo de confiança fique com formato de sombra dentro do gráfico. Conforme figura anexada: URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190904/bfde02f0/atta… > Como tentativa usei o *polygon*, mas não sei se existe um outro recurso melhor. Se alguém souber e puder contribuir, a lista do R agradece. Segue o código: x <- seq(1:30) Y <- x + rnorm(30, 0,3) model <- lm(Y ~ x) pred <- data.frame(x1 = x) pred$y1 <- predict(model, newdata = pred, interval = "confidence") plot(consumo, type = "n", main = "Consumo médio, kg", xlab = "Consumo estimado, kg ", ylab = "Consumo observado, kg" ) icfit = pred$y1 yy = c(as.numeric(as.vector(icfit[,2])), rev(as.numeric(as.vector(icfit[,3])))) xx = c(x, rev(x)) polygon(xx, yy, col = "lightsteelblue2", border = NA) points(consumo) abline(model, col="blue", lwd= 2) points(seq(2:15), seq(2:15), type = "l", col = "green", lwd = 2) legend("topright", legend = c("modelo linear base y = x", "modelo ajustado", "Bandas de confiança, 95%" ), lty = c(1,1,1), col = c("green", "blue", "lightsteelblue2"), lwd = c(2,2,NA), bty = "n", pch = c(NA,NA, 15))

4 3

Re: [R-br] gráficos com bandas de confiança no R - com sombra
by sznelwar＠uol.com.br 05 Set '19

05 Set '19

1 0

gráficos com bandas de confiança no R - com sombra
by Nilson Guiotoku 05 Set '19

05 Set '19

Pessoal, gostaria de fazer um gráfico conforme o anexo, cujo o intervalo de confiança fique com formato de sombra dentro do gráfico. Conforme figura anexada: [image: image.png] Como tentativa usei o *polygon*, mas não sei se existe um outro recurso melhor. Se alguém souber e puder contribuir, a lista do R agradece. Segue o código: x <- seq(1:30) Y <- x + rnorm(30, 0,3) model <- lm(Y ~ x) pred <- data.frame(x1 = x) pred$y1 <- predict(model, newdata = pred, interval = "confidence") plot(consumo, type = "n", main = "Consumo médio, kg", xlab = "Consumo estimado, kg ", ylab = "Consumo observado, kg" ) icfit = pred$y1 yy = c(as.numeric(as.vector(icfit[,2])), rev(as.numeric(as.vector(icfit[,3])))) xx = c(x, rev(x)) polygon(xx, yy, col = "lightsteelblue2", border = NA) points(consumo) abline(model, col="blue", lwd= 2) points(seq(2:15), seq(2:15), type = "l", col = "green", lwd = 2) legend("topright", legend = c("modelo linear base y = x", "modelo ajustado", "Bandas de confiança, 95%" ), lty = c(1,1,1), col = c("green", "blue", "lightsteelblue2"), lwd = c(2,2,NA), bty = "n", pch = c(NA,NA, 15))

1 0

Problema criar colunas: socorro!
by Elerson 02 Set '19

02 Set '19

Prezados: Trabalho um banco de dados com mais de 10000 linhas, representado pelo exemplo abaixo. *PROCESSO* *DATA_DISTRIBUICAO* *NOME_PARTE* *DATA_CRIME* *DATA_DENUNCIA* *ENQUADRAMENTO* 180022121 02/01/2018 A_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022121 02/01/2018 B_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.35 - 11343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.33 - 11343/06 180022097 02/01/2018 D_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022097 02/01/2018 E_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR.4º INC.II E IV - CPB 180022097 02/01/2018 F_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022501 03/01/2018 G_NONIMATO 02/01/2018 03/01/2018 ART.157 PAR.2 INC.I E II - CP 180022501 03/01/2018 G_NONIMATO 02/01/2018 03/01/2018 ART.14 - CP Hipóteses: a) mesmo PROCESSO, NOME_PARTE diferentes, ENQUADRAMENTO igual (linhas 1 e 2) b) mesmo PROCESSO, mesmo NOME_PARTE, ENQUADRAMENTO diferentes (linhas 3 e 4) c) mesmo PROCESSO, mesmo NOME_PARTE, ENQUADRAMENTO diferentes (linhas 8 e 9) Minha questão se refere à hipótese c: preciso criar um coluna (enquadramento2) para transferir o art. 14 transformando o caso em apenas uma linha… sem que se mexa na hipótese b. Em todo o banco sempre estará escrito ART.14 – CP. Para ficar assim: *PROCESSO* *DATA_DISTRIBUICAO* *NOME_PARTE* *DATA_CRIME* *DATA_DENUNCIA* *ENQUADRAMENTO* *ENQUADRAMENTO2* 180022121 02/01/2018 A_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022121 02/01/2018 B_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.35 - 11343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.33 - 11343/06 180022097 02/01/2018 D_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022097 02/01/2018 E_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR.4º INC.II E IV - CPB 180022097 02/01/2018 F_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022501 03/01/2018 G_NONIMATO 02/01/2018 03/01/2018 ART.157 PAR.2 INC.I E II - CP ART.14 - CP Obrigado.

3 3