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

Calculo do Volume a partir de imagens MDT
by julio cesar oliveira 23 Fev '22

23 Fev '22

Prezados Pretendo realizar o cálculo do volume de um MDT (no formato raster) para os pixels que estão acima ou abaixo de um determinado valor. Por acaso existe algum pacote / função que faça esse cálculo? Grato Julio

1 0

Dúvida reestruturar variáveis.
by Diogo Jerônimo 16 Fev '22

16 Fev '22

Prezados, bom dia de novo, lá vou perturbar os amigos... Segue outro exemplo reproduzível (tem a ver com minha dúvida de ontem): library(dplyr) x01<-c(1001,"Potássio",3.49,3.62,3.63,3.49,3.44,2) x02<-c(2000,"Potássio",4.18,4.21,4.15,4.86,4.36,3)x03<-c(1002,"Potássio",3.49,3.62,3.63,3.49,3.44,2)x04<-c(1003,"Potássio",3.49,3.62,3.63,3.49,3.44,2)x05<-c(2001,"Potássio",4.18,4.21,4.15,4.86,4.36,3)x06<-c(2000,"Albumina",2.35,2.36,2.21,2.15,2.27,10)x07<-c(3000,"Albumina",2.69,2.78,2.88,2.91,2.84,12)x08<-c(2001,"Albumina",2.35,2.36,2.21,2.15,2.27,10)x09<-c(3001,"Albumina",2.69,2.78,2.88,2.91,2.84,12) exemplo<-data.frame(rbind(x01,x02,x03,x04,x05,x06,x07,x08,x09))names(exemplo)<-c("CodCli","Substância","v1","v2","v3","v4","v5","Cluster")row.names(exemplo)<-NULL exemplo<-exemplo %>%mutate(v1=as.numeric(v1),v2=as.numeric(v2),v3=as.numeric(v3),v4=as.numeric(v4), v5=as.numeric(v5)) exemplo$Média<-rowMeans(exemplo[,3:7],na.rm = TRUE) Como o de ontem, esse banco tem variáveis que identificam o código do cliente, a média dos dados (iguais entre laboratórios) e o Cluster, que foi obtido previamente. A substância é detalhe adicional. O que eu gostaria é: reestruturar esse banco, colocando como variáveis de identificação a substância, o cluster e a média, e o código do cliente indo para as colunas, como está nessa tabela abaixo: Subst Cluster Média CodCli.1 CodCli.2 CodCli.3Albumina 10 2,2680 2000 2001 Albumina 12 2,8200 3000 3001 Potássio 2 3,5340 1001 1002 1003Potássio 3 4,3520 2000 2001 Eu tentei usar o comando reshape, com esse código abaixo, mas ele não funcionou comigo (tentei também o spread, do tidyr): exemplo<-exemplo %>% select(CodCli,Substância,Cluster,Média) exemplo_quebra<-reshape(exemplo, direction="wide", idvar=c("Substância","Cluster","Média"), timevar="Cod.Cli") Algum amigo poderia indicar onde estou errando? Ou se teria uma rotina melhor para realizar esse procedimento? Novamente obrigado!!! Diogo Jerônimo Bacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

3 4

Contagem de fatores
by Sergio Pinto 15 Fev '22

15 Fev '22

Prezados Tenho um data-frame (dfp21) com duas variáveis, a primeira variável, chamada Deficiências, é categórica com sete níveis ('TEA', 'SCe', 'NE', 'DV', 'DF', 'DA/S', 'AH/S') a segunda, chamada p21, também é categórica com dois níveis ("8" e "9"). Coloquei os dados em anexo. Não consegui usar o DATAFILEHOST. > dfp21 # A tibble: 127 × 2 p21 DEFICIÊNCIAS <fct> <fct> 1 9 TEA 2 9 TEA 3 8 TEA 4 9 TEA 5 9 TEA 6 9 TEA 7 9 TEA 8 9 TEA 9 9 TEA 10 9 TEA # … with 117 more rows *A dúvida:* como criar uma nova coluna com a contagem da segunda variável em cada nível da primeira? ex: TEA 8 20 TEA 9 13 NE 8 17 NE 9 15 ... Estou tentando usar o dplyr para isso. Desde já agradeço a atenção. -- *_______________________________________* *Prof. Dr. SERGIO MEDEIROS PINTO* LABICOM/IEFD-UERJ FAMERC-RJ UNESA-RJ http://lattes.cnpq.br/5897044784217523

3 3

Dúvida lógica do código R.
by Diogo Jerônimo 15 Fev '22

15 Fev '22

Boa tarde a todos!!! Dúvida boba, mas tenho de explicar antes: estou fazendo uma análise de cluster, para poder identificar resultados iguais de laboratórios analíticos. O resultado foi um banco parecido com esse reproduzível: library(dplyr) x01<-c(8.49,8.62,8.63,8.49,8.44,2) x02<-c(9.18,9.21,9.15,8.86,9.36,3)x03<-c(8.49,8.62,8.63,8.49,8.44,2)x04<-c(9.74,9.69,9.81,9.25,9.98,3)x05<-c(9.74,9.69,9.81,9.25,9.98,3)x06<-c(8.15,9.52,9.89,9.95,8.33,2)x07<-c(8.32,8.73,9.60,9.04,9.97,3)x08<-c(7.42,7.01,7.31,7.20,7.61,1)x09<-c(9.79,8.76,9.68,8.31,9.37,3)x10<-c(7.42,7.01,7.31,7.20,7.61,1) dados<-data.frame(rbind(x01,x02,x03,x04,x05,x06,x07,x08,x09,x10)) names(dados)<-c("v1","v2","v3","v4","v5","Cluster") Para facilitar essa identificação, ordenei os resultados por cluster, pelos dados das variáveis, e então calculei a média em cada linha: dados<-dados %>% arrange(Cluster, v1,v2,v3,v4,v5) for (i in (1:nrow(dados))){dados$Média[i]<-mean(dados$v1[i],dados$v2[i],dados$v3[i],dados$v4[i],dados$v5[i],na.rm = TRUE)} MEU PROBLEMA: para finalmente identificar os resultados iguais, eu apliquei esse código abaixo, onde em resumo, se o cluster for igual e a média do resultado do laboratório da linha i for igual a média do resultado do laboratório da linha i+1, então ambos os laboratórios responderam resultados iguais. Apliquei o código abaixo, e ele está resultando em erro: i=1for (i in (1:nrow(dados))){if ((dados$Cluster[i]==dados$Cluster[i+1]) & (dados$Média[i]==dados$Média[i+1])){dados$sit[i] == "Resultado Igual"dados$sit[i+1] == "Resultado Igual"}else {}i=i+1} Alguém poderia dizer qual o meu erro? Ou se faltou algo no código? Muito agradecido pela atenção, e desculpa a extensão!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

3 4

Problema de reconhecimento de projeções no RGDAL
by Helder Gramacho 09 Fev '22

09 Fev '22

Boa tarde pessoal, Estou tendo problemas em reconhecer o sistema de coordenadas de imagens de satélite no R, provavelmente devido a algum problema com o rgdal, quando executo o comando crs aparece o seguinte erro: crs(imagem_de_satelite) Coordinate Reference System: Deprecated Proj.4 representation: NA Warning message: In wkt(x) : CRS object has no comment Alguém já teve este problema e conseguiu resolver? Desde já agradeço, *Hélder G. dos Santos*

1 0

Texto na margem do ggplot
by Pedro Emmanuel Alvarenga Americano do Brasil 03 Fev '22

03 Fev '22

Amigos de R, Saudações a todos. Gostaria de colocar um texto nas margens superior e direita do gráfico de tal forma que eu consiga diferenciar duas variáveis que têm as mesmas categorias. No graphics seria equivalente a função mtext(). Por exemplo... plot(bla bla bla) mtext("Esquema inicial", side = 3) mtext("Reforço", side = 4) Mas no ggplot não estou sabendo fazer. Alguma luz? # Simulando os dados x <- rbind( data.frame(id = 1:20, Visita = "Visita 1", Inicial = "Vacina 1", Reforço = "Vacina 1", Titulos = rnorm(20, 20, 5)), data.frame(id = 1:20, Visita = "Visita 2", Inicial = "Vacina 1", Reforço = "Vacina 1", Titulos = rnorm(20, 30, 5)), # data.frame(id = 21:40, Visita = "Visita 1", Inicial = "Vacina 2", Reforço = "Vacina 1", Titulos = rnorm(20, 40, 5)), data.frame(id = 21:40, Visita = "Visita 2", Inicial = "Vacina 2", Reforço = "Vacina 1", Titulos = rnorm(20, 50, 5)), # data.frame(id = 41:60, Visita = "Visita 1", Inicial = "Vacina 1", Reforço = "Vacina 2", Titulos = rnorm(20, 15, 5)), data.frame(id = 41:60, Visita = "Visita 2", Inicial = "Vacina 1", Reforço = "Vacina 2", Titulos = rnorm(20, 30, 5)), # data.frame(id = 61:80, Visita = "Visita 1", Inicial = "Vacina 2", Reforço = "Vacina 2", Titulos = rnorm(20, 15, 5)), data.frame(id = 61:80, Visita = "Visita 2", Inicial = "Vacina 2", Reforço = "Vacina 2", Titulos = rnorm(60, 30, 5))) head(x) # Fazendo o grafico... ggplot(x, aes(x = Visita, y = Titulos, group = id)) + geom_line(colour = "lightgreen", show.legend = FALSE) + geom_smooth(aes(group = 1), method = lm, colour = "red", na.rm = TRUE) + facet_grid(rows = vars(Inicial), cols = vars(Reforço))+ theme_light() + ylab("Titulos de anticorpos") + xlab("Visitas de seguimento") Abraços fraternos, Pedro Brasil

3 4