- R-br - DINF UFPR

Problemas com a paleta de cores
by Michele Claire Breton 31 Mai '22

31 Mai '22

Caríssimos(as), boa tarde! Espero que todos se encontrem bem. Estou tentando plotar um gráfico de bolhas e estou com problemas com as cores. A mensagem de erro é a seguinte: [image: image.png] O comando foi: bubble_plot <- ggplot(m3,aes(replicates_list, gene_per)) + geom_point(aes(size=gene_per, fill=gene_col),shape=21,color="black") + theme(panel.grid.major=element_line(linetype=1,color="grey"), axis.text.x=element_text(angle=90,hjust=1,vjust=0), panel.background = element_blank()) + ylab("Gene ID") + xlab("Prostate Adenocarcinome cell types") + scale_fill_discrete(name="Expression values")+ scale_size(name="% of cells with expression") bubble_plot Minha intenção é que haja um gradiente de cores (tipo heatmap) para valores que vão de 0.01 a 5. e o resultado que espero ver é algo semelhante a isto: [image: image.png] Por favor, me ajudem! Falta só isto para plotar o gráfico! Aproveito a oportunidade para perdir-lhes indicações de cursos online brasileiros de R. Muitíssimo obrigada pela atenção, Michele -------------------------------------------------------------------------- *Dra. Michele Claire Breton* Técnica Superior em Bioinformática Doutora em Biologia Celular e Molecular CICS - Centro de Investigação em Ciências da Saúde Faculdade de Ciências da Saúde Universidade da Beira Interior Covilhã - Castelo Branco - Portugal

2 1

Trifatorial com trat adicionais
by Alessandro Dal Col Lúcio 31 Mai '22

31 Mai '22

Bom dia. Alguma dica de pacote para analisar experimentos trifatoriais com mais de um tratamento adicional? Abraços. Alessandro -- Enviado via Windows Phone.

2 2

(sem assunto)
by Elias Carvalho 30 Mai '22

30 Mai '22

Boa noite Digamos que extrai dados sobre mortalidade de uma população de certa cidade. Daí quero avaliar a taxa de mortalidade de brancos na faixa etária de 25-34 anos. Se a cidade tem 100.000 habitantes, porém brancos de 25-34 anos são apenas 50.000, eu devo normalizar a taxa de mortalidade por 100.000 ou 50.00?

2 1

função geom_oint
by Michele Claire Breton 25 Mai '22

25 Mai '22

Caríssimos(as) colegas(as), muito bom dia! Espero que todos se encontrem bem. Estou com uma dificuldade. Quero fazer um bubble plot usando as funções ggplot2 e dplyr, onde pretendo plotar no eixo x o gene_ID (são 84 linhas), no y a % de ocorrência da expressão na célula e no color, o valor da expressão. O meu primeiro problema é que no y são 9 tipos celulares diferentes e cada um tem, na tabela de entrada de dados, uma coluna com os dados de expressão e outra com os dados de % de ocorrência (amostra abaixo). Como eu escrevo a função geom_point para plotar o gráfico? Amostra da tabela de entrada gene BE-exp BE% CE_exp CE% E_exp E% F_exp F% HE_exp HE% L-exp L% LE-exp LE% N-exp N% SM-exp SM% ACTB 4.006 99.8 3.690 99.6 2.691 95.2 2.377 80.1 2.901 97.7 3.182 94.1 2.864 94.0 3.243 96.0 3.214 93.5 ACTG1 3.850 5.9 3.592 7.1 2.469 1.5 2.252 1.1 3.129 6.4 2.451 0.6 2.979 1.1 3.104 0 2.500 2.4 AJUBA 1.185 5.9 0.972 7.1 1.042 1.5 0.912 1.1 1.043 6.4 0.865 0.6 1.024 1.1 0 0 1.149 2.4 Muitíssimo obrigada pela ajuda! Michele -- ------------------------------------------------------------------------ *Dra. Michele Claire Breton* Técnica Superior em Bioinformática CICS - Centro de Investigação em Ciências da Saúde Faculdade de Ciências da Saúde Universidade da Beira Interior Covilhã - Castelo Branco - Portugal

4 4

Anaconda/Jupyter erro :Graphics API version mismatch -
by Fernando Souza 25 Mai '22

25 Mai '22

Boa noite pessoal, Após atualizar para a versão 4.1.3 do R no Anaconda passei a receber esta mensagem: "Error in Cairo::Cairo(width, height, tf, "png", pointsize, bg, "transparent", : Graphics API version mismatch Em pesquisa na rede vi que tem a ver com uma mudança na api Cairo do R. Vi sugestões de desinstalar e instalar o Cairo entre outras, porém nehuma funcionou. No Rstudio, tudo funciona bem, então acredito ser um problema no Anaconda/Jupyter. Alguém passou por isso? Sabe como resolver? Att

1 1

Heatmap com dois conjuntos de dados
by Michele Claire Breton 20 Mai '22

20 Mai '22

Caríssimos(as) mestres(as), bom dia! Estou com uma dúvida sobre como deve ser meu ficheiro de entrada para plotar um gráfico como o terceiro deste link: https://stackoverflow.com/questions/67746044/r-heatmap-with-circles Eu criei 2 ficheiros .csv, um contento os valores de expressão gênica nos diferentes tipos celulares (col->genes, lin->tipos de células); e o outro, contendo os valores de freqüência em que a expressão do gene ocorre em cada célula, com a mesma estrutura do primeiro. A representação destes dados no gráfico deve ser: a área do círculo construído a partir dos dados de frequencia e a cor do círculo deve corresponder aos valores de expressão, como em um heatmap. Eu tenho que trabalhar com os ficheiros em separado ou tenho que montar um único ficheiro para ter este resultado? Muitíssimo obrigada pela ajuda!! Michele -- ------------------------------------------------------------------------ *Dra. Michele Claire Breton* Técnica Superior em Bioinformática CICS - Centro de Investigação em Ciências da Saúde Faculdade de Ciências da Saúde Universidade da Beira Interior Covilhã - Castelo Branco - Portugal

2 1

Ajuda sobre licença
by Alexandre Loures 19 Mai '22

19 Mai '22

Bom dia, Listeiros! Com relação ao package covid19italy alterei apenas os endereços, no arquivo data_refresh.R, para busca da base para meu github bem como os nomes, nos arquivos data.R e data_refresh.R, de covid19italy para covid19brazil, italy_total para brazil_total, etc. Gostaria de subir este pacote para o Cran, mas minha dúvida é sobre direitos autorais. A licença do covid19italy está definida como: MIT + file LICENSE. Com essa licença do covid19italy posso subir meu pacote? Package: covid19italy Type: Package Title: The 2019 Novel Coronavirus COVID-19 (2019-nCoV) Italy Dataset Version: 0.3.1 Authors@R: person("Rami", "Krispin", email = "rami.krispin(a)gmail.com", role = c("aut", "cre")) Maintainer: Rami Krispin <rami.krispin(a)gmail.com> Description: Provides a daily summary of the Coronavirus (COVID-19) cases in Italy by country, region and province level. Data source: Presidenza del Consiglio dei Ministri - Dipartimento della Protezione Civile <https://www.protezionecivile.it/>. License: MIT + file LICENSE Encoding: UTF-8 LazyData: true Depends: R (>= 3.0.2) Imports: devtools Suggests: knitr (>= 1.28), readr (>= 1.3.1), remotes (>= 2.1.1), rmarkdown (>= 2.1), testthat (>= 2.1.0) URL: https://github.com/RamiKrispin/covid19italy BugReports: https://github.com/RamiKrispin/covid19italy/issues RoxygenNote: 7.1.1 VignetteBuilder: knitr NeedsCompilation: no Packaged: 2021-07-28 05:39:41 UTC; ramikrispin Author: Rami Krispin [aut, cre] Repository: CRANDate/Publication: 2021-07-28 06:00:02 UTC Meu pacote no github: https://github.com/AlexandreLoures/covid19Brazil -------- Alexandre Loures Programa de Pós-Graduação em Organizações e Mercados Universidade Federal de Pelotashttp://www.wp.ufpel.edu.br/economia/ www.rodriguesloures.com Sent with [ProtonMail](https://protonmail.com/) secure email.

2 4

Analise Geoespacial
by Elias Carvalho 18 Mai '22

18 Mai '22

Boa tarde Alguém poderia por favor indicar um bom livro com exemplos de análise geoespacial? -- *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"*

1 0

Cores do heatmap
by Michele Claire Breton 18 Mai '22

18 Mai '22

Caríssimos conterrâneos, bom dia! Espero que todos se encontrem bem. Eu estou trabalhando em um heatmap e resolvi mudar o esquema de cores, entretanto, após a mudança, aparece a mensagem de erro: Error in image.default(1:nc, 1:nr, x, xlim = 0.5 + c(0, nc), ylim = 0.5 + : must have one more break than colour. A princípio, eu usei: my_palette <- colorRampPalette(c("green","yellow","red"))(n=299) col_breaks = c(seq(-2,-1,length=100), seq(-0.9,0.9, length=100), seq(1,2,length=100)) E quis mudar para: my_palette <- colorRampPalette(c("blue","white","red"))(n=20) col_breaks = c(seq(-2,-1,length=100), seq(-0.9,0.9, length=100), seq(1,2,length=100)) ... e enttão surgiu a mensagem de erro. O que fiz de errado? Muito obrigada por me ajudarem Michele -- ------------------------------------------------------------------------ *Dra. Michele Claire Breton* Técnica Superior em Bioinformática CICS - Centro de Investigação em Ciências da Saúde Faculdade de Ciências da Saúde Universidade da Beira Interior Covilhã - Castelo Branco - Portugal

2 2

Mudar script encoding
by Diogo Jerônimo 28 Abr '22

28 Abr '22

Procurei nas ridículas e não achei... Tem como mudar, no RStudio, a opção de default text encoding de "UTF-8" para "ISO8859-1"? Nas opções em que eu entrei, ora desconfigura meu código, ora apaga o mesmo... Obrigado!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

2 1