Dezembro 2021 - 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

ggplot + geom_bar
by Pedro Emmanuel Alvarenga Americano do Brasil 21 Dez '21

21 Dez '21

Amigos de R, Estou num esforço de tentar aprender a usar o tidyverse. Gostaria de manter a consistência de um relatório que estou trabalhando com gráficos do ggplot. Mas há um particular gráfico de barras que não estou conseguindo fazer. Alguém pode me ajudar a fazer o gráfico no ggplot + geom_bar? Segue abaixo um exemplo como eu faria com o graphics. # Preparando os dados simulados tab <- array(c(92,74,11,0,0,94,98,10,1,0, 83,56,13,2,0,96,68,9,7,0, 46,28,7,1,1,54,59,16,7,0, 26,20,3,0,0,36,32,16,4,0, 8,2,0,0,0,3,3,1,0,0), c(5, 2, 5)) dimnames(tab) <- list(c("Definitively satisfied", "Satisfied", "Little satisfied", "Little dissatisfied", "Definitively dissatisfied"), c("Surface 1","Surface 2"), c(1:5)) names(dimnames(tab)) <- c("Satisfaction","Intervention","Visit") tab a <- epitools::expand.table(tab) a$Visit <- paste0("Visit ", a$Visit) head(a) # Fazendo as tabelas ja com os percentuais x <- sapply(unique(a$Visit), function(i) prop.table(table(a$Satisfaction[a$Visit==i],a$Intervention[a$Visit==i]),2), simplify = FALSE) # Criando o objeto height do barplot y <- x[[1]] ; for(i in 2:length(x)) {y <- cbind(y,x[[i]])} ramp <- colorRampPalette(c("darkgreen","red")) # Criando o gráfico z <- barplot(y, beside = T, col = ramp(5), legend = T, args.legend = list(x="top",xpd=NA,inset=-.15,bty="n",cex=.8, ncol = 2), axisnames = FALSE) axis(2, at = seq(0,1,.025),labels = FALSE, tck = -.01) text(z[3,], -.005, colnames(y), cex = .5, srt = 35, xpd = NA, adj = 1, offset = 5) axis(1, at = sapply(seq(0,ncol(z)-2,2),function(i) mean(z[,c(1:2) + i])),labels = unique(a$Visit), line = 1.5, tick = FALSE) sapply(1:ncol(y), function(i) text(x = z[1:5,i], y = y[1:5,i], labels = round(y[1:5,i],2), pos = 3, cex = .5, xpd = NA)) Abraço forte a todos, Pedro Brasil

3 3

Help Teste de Dunn
by ana paula coelho madeira 19 Dez '21

19 Dez '21

Boa noite a todos Apliquei um teste de Kruskal Wallis visando comparar grupos independentes e no posto hoc apliquei o teste de Dunn. Utilizei os pacotes FSA e dunn.test e nenhum dos dois apresenta na saída as tradicionais letras usadas nos testes paramétricos. Encontrei um material que faz referência a um pacote "rcompanion" porém não está mais disponível. Desta forma, gostaria de saber se há alguma indicação de pacotes que fazer o teste de Dunn e apresentem na saída as letras para indicar diferenças entre os grupos. Desde já agradeço, Ana Paula OBS: Peço desculpas caso o mesmo email tenha sido enviado mais de uma vez, achei que não estava mais inscrita

2 1

k means em variável 2d
by Stefane Fonseca Freitas 18 Dez '21

18 Dez '21

Boa tarde usuários, Sou iniciante no uso do R e tenho dúvidas em relação ao uso da função de clusterização kmeans. Meu banco de dados é um perfil de refletividade, logo tenho observações para cada nível na vertical, neste caso, 144 níveis. Meu objetivo é aplicar o kmeans nesse banco de dados, preservando a estrutura inicial dele, que é 2d. Já realizei testes no R para outra variável 1d e li que era necessário preparar o dado em formato data frame para então usar o kmeans. No entanto, fiquei em dúvida se para esse caso específico, é possível realizar a clusterização sem modificar essa estrutura que é importante que se mantenha. Desde já agradeço,

1 0

Teste de Dunn
by ana paula coelho madeira 18 Dez '21

18 Dez '21

Boa tarde a todos, Apliquei um teste de Kruskal Wallis visando comparar grupos independentes e no posto hoc apliquei o teste de Dunn. Utilizei os pacotes FSA e dunn.test e nenhum dos dois apresenta na saída as tradicionais letras usadas nos testes parametricos. Encontrei um material que faz referência a um pacote "rcompanion" porém não está mais disponível. Desta forma, gostaria de saber se há alguma indicação de pacotes que fazer o teste de Dunn e apresentem na saída as letras para indicar diferenças entre os grupos. Desde já agradeço Ana Paula Obter o Outlook para Android<https://aka.ms/AAb9ysg>

2 1

Extraindo o "cumulative hazard" de objeto survfit
by Pedro Emmanuel Alvarenga Americano do Brasil 16 Dez '21

16 Dez '21

Saudações amigos do R, Estou as voltas de estimar taxas de eventos e estou batendo cabeça. Antes de escrever uma função eu mesmo para fazer uma tabela com alguns valores gostaria de uma luz dos amigos de R. O banco aml está no pacote survival então bastaria carregar o pacote pra reproduzir o exemplo. Eu gostaria de ver em formato de tabela alguns momentos específicos da tabela de sobrevivência. Só que o summary.survfit só retorna a sobrevivência e não a taxa cumulativa. library("survival") fit <- survfit(Surv(time, status) ~ x, data = aml) > fit Call: survfit(formula = Surv(time, status) ~ x, data = aml) n events median 0.95LCL 0.95UCL x=Maintained 11 7 31 18 NA x=Nonmaintained 12 11 23 8 NA # Summary com todos os momentos de eventos > summary(fit) Call: survfit(formula = Surv(time, status) ~ x, data = aml) x=Maintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 9 11 1 0.909 0.0867 0.7541 1.000 13 10 1 0.818 0.1163 0.6192 1.000 18 8 1 0.716 0.1397 0.4884 1.000 23 7 1 0.614 0.1526 0.3769 0.999 31 5 1 0.491 0.1642 0.2549 0.946 34 4 1 0.368 0.1627 0.1549 0.875 48 2 1 0.184 0.1535 0.0359 0.944 x=Nonmaintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 5 12 2 0.8333 0.1076 0.6470 1.000 8 10 2 0.6667 0.1361 0.4468 0.995 12 8 1 0.5833 0.1423 0.3616 0.941 23 6 1 0.4861 0.1481 0.2675 0.883 27 5 1 0.3889 0.1470 0.1854 0.816 30 4 1 0.2917 0.1387 0.1148 0.741 33 3 1 0.1944 0.1219 0.0569 0.664 43 2 1 0.0972 0.0919 0.0153 0.620 45 1 1 0.0000 NaN NA NA # Summary com os momentos desejados > summary(fit, times = c(14,28,35)) Call: survfit(formula = Surv(time, status) ~ x, data = aml) x=Maintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 14 8 2 0.818 0.116 0.619 1.000 28 6 2 0.614 0.153 0.377 0.999 35 3 2 0.368 0.163 0.155 0.875 x=Nonmaintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 14 7 5 0.583 0.142 0.3616 0.941 28 4 2 0.389 0.147 0.1854 0.816 35 2 2 0.194 0.122 0.0569 0.664 > plot(fit) > plot(fit, cumhaz = T) > Reparem que há uma opção para o gráfico cumhaz = T, o que significa que o cumhaz está depositado no objeto, inclusive dentro do summary também. Tipo summary(fit)$cumhaz. Só que não há uma opção summary(fit, cumhaz = T) que retorne o cumhaz ao invés da sobrevivência. Alguém tem algum bizu pra fazer isso organizado, extrair a mesma tabela só que com o cumhaz, como summary(fit, times = c(14,28,35)) sem muito trabalho? Abraço forte, Pedro Brasil

4 11

Download de dados públicos do DATASUS
by Pedro Emmanuel Alvarenga Americano do Brasil 07 Dez '21

07 Dez '21

Saudações turma, Talvez seja meio off topic mas não estou sabendo a quem recorrer. Eu estou procurando uma maneira inteligente de acessar dados públicos que estariam disponíveis pelo DATASUS. 1. Tem alguém na lista que trabalha com dados secundário do PNI (programa nacional de imunizações)? Alguém que poderia me ajudar a encontrar dados públicos dos CRIE? 2. Há alguns lugares meio escondidos do ministério da saúde e do datasus que dão acesso a algumas coisas a partir de um ftp como aqui ftp://ftp.datasus.gov.br/ Mas aqui os dados que me interessariam só vão até 2019, ou seja, não há continuidade na disponibilização. Neste exemplo eu acho que usar um utils::download.file() daria certo pra baixar por dentro do R. Alguma sugestão mais inteligente para baixar arquivos por dentro do R? 3.Como baixar dados que me interessam de um sítio que não tem um ftp visível e eu deveria dar um monte de cliques pra ver a o botão de baixar em csv? http://tabnet.datasus.gov.br/cgi/dhdat.exe?bd_pni/cpnibr.def por exemplo nesse endereço eu deveria marcar na caixa "periodo" todos os anos, na caixa "linha" marcar ano, na caixa "coluna" marcar imuno, depois na seleções Imuno+ marcar "Sarampo" e clicar em pesquisa e só depois aparece o botão salvar como csv ou excel. Não há um jeito inteligente de fazer isso por dentro do R? Aceito sugestões de leitura!!! Eu tentei umas paradas no RCurl, downloader e remotes mas sem sucesso. Abraço forte a todos, Pedro Brasil

3 3

Análise de experimento fatorial com parcela subdividida no tempo
by Maurício Sangiogo 04 Dez '21

04 Dez '21

1 0