- R-br - DINF UFPR

Leitura de valores no csv
by Elias Carvalho 28 Ago '23

28 Ago '23

Bom dia amigos Tenho um arquivo .csv com uma variável chamada patrimônio que possui valores acima de 100.000 reais, assim: Patrimonio 100.000,00 110.000,00 120.000,00 120.000,00 120.000,00 Leio no RStudio com o comando: read.csv("dados.csv", sep = ',', dec = ",", header = TRUE), porém a variável vem sem casas decimais [image: Captura de tela de 2023-08-25 09-58-01.png] Tentei converter com: data.to.work$Patrimonio <- as.double(data.to.work$Patrimonio) Porém só mudou para num e sem casas decimais. Onde estou errando? -- *In Jesu et Maria* *Obrigado / Thank you* *Prof. Dr. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

3 4

Problemas atalhos
by Alfredo Diaz 28 Fev '23

28 Fev '23

Cordiais cumprimentos, Preciso de sua ajuda porque, quando estou trabalhando no Rstudio, os atalhos *CTRL + 1, CTRL + 2* não estão funcionando, também quando uso o atalho *CTRL + SHIFT + M *o operador que recebo é: |>. Atalhos como *CTRL + ENTER, ALT + -, ALT + SHIFT + K* estão funcionando perfeitamente, então eu não entendo o que está acontecendo, por que está acontecendo, o que devo fazer para resolver este problema? Muito obrigado de antemão por suas respostas. Fraternalmente, --- José Diaz E.

3 8

RMarkdown - Aplicação knitr e render
by Diogo Jerônimo 28 Fev '23

28 Fev '23

Bom dia amigos, minha dúvida é específica do RMarkdown hoje. Seguinte, irei automar relatórios que serão formatados no Mark, sendo salvos em html. Como reproduzível, coloco abaixo esse exemplo, que é executado em arquivo *.Rmd. ---title: "Relatório"author: "Diogo"output: html_document--- ```{r informações do participante, echo = FALSE, include = FALSE} library(tidyr) library(dplyr)library(rmarkdown)library(kableExtra)cli<-rep(1,10);x<-round(rnorm(10,10,1),2);e<-round(rnorm(10,0,1),2);y<-x+e;dados_ini<-data.frame(cbind(cli,x,y)) dados<-dados_ini %>% filter(cli==1) ``` ```{r tabela, echo = FALSE}tabela<-data.frame(dados[,c(2,3)]) tabela %>% kbl()``` ```{r gráfico, echo = FALSE}plot(x,y) ``` Aplicando o famoso botão "knitr", ele me retorna o relatório em html COM o gráfico. Acontece que para a automação, terei de usar o comando render dentro de uma function, para poder formatar o relatório. Ao aplicar o mesmo comando nessa estrutura, eu gero tabela e gráfico dentro de um arquivo R, e depois chamo o comando render. Da seguinte forma: library(tidyr)library(dplyr)library(rmarkdown)library(kableExtra)cli<-rep(1,10);x<-round(rnorm(10,10,1),2);e<-round(rnorm(10,0,1),2);y<-x+e;dados_ini<-data.frame(cbind(cli,x,y)) dados<-dados_ini %>% filter(cli==1) tabela<-data.frame(dados[,c(2,3)]) graf<-plot(x,y) render(input="c:/users/diogo/desktop/rmarkdown_via render_listaR.Rmd", output_format = "html_document", output_file = "relatório.html", output_dir = "c:/users/diogo/desktop/") Sendo que o arquivo Rmarkdown "rmarkdown_via render_listaR.Rmd" tem essa estrutura dentro: ---title: "Relatório"author: "Diogo"output: html_document--- ```{r tabela, echo = FALSE}tabela %>% kbl() ``` ```{r gráfico, echo = FALSE}print(graf)``` O meu problema é que em ambos os casos, eu até obtenho as mesmas tabelas, mas nesse segundo modo (que é o que viabiliza a function, para formatar vários relatórios), ele me retorna o relatório em html dessa vez SEM o gráfico. Enfim a dúvida: alguém sabe se tem algum subcomando dentro do render, ou opção, para forçar o gráfico a aparecer na segunda forma (usando o comando render)? Ou se estou aplicando alguma funcionalidade erradamente? Obrigado e bom dia aos amigos!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

1 0

[ OFF TOPIC ] Qual abordagem utilizar - associação de compostos com resistência a praga/doença
by Marcelo Laia 23 Fev '23

23 Fev '23

Caro(a) listeiro(a), Possuo um grupo de dados com valores representando a proporção de dado composto na planta. Esses valores variam de 0 a 100%. Assim, o composto Y pode apresentar valor 41 e o Z valor 0,02, por exemplo. Há em torno de 70 compostos. Nem todas as plantas contém todos os compostos. Hipótese: plantas resistentes possuem compostos exclusivos e/ou em proporção significativamente diferente das plantas suscetíveis. Para os compostos exclusivos, é tranquilo, pois, se o Z está presente somente no resistente, é muito possível que ele tenha influência sobre a resistência. O problema são aqueles que aparecem tanto no resistente quanto no suscetível em proporção diferente. Outra questão que pode acontecer é a interação entre compostos. Ou seja, dois ou mais compostos explicarem o efeito de resistência ou suscetibilidade. Neste caso, um composto pode agir como facilitador para a ação de outro, que seria o composto tóxico. Logo, sem o primeiro, o composto tóxico, sozinho, não teria efeito. Pensei em usar abordagem genômica, tratando cada composto como um gene. Não deu certo. Resultado estatístico não explica o biológico. Outra opção seria pensar na análise de QTLs. Não tentei essa. Mas, lendo trabalhos na área médica, vi que alguns utilizam vários sintomas para explicar a doença. Por exemplo, associam faixas de valores de hemoglobina, ácido úrico e proteína X para explicar a ocorrência de determinado câncer. Pergunto: você já viu ou já fez análises estatísticas no R para esse último caso? Explicar a doença com base em um roll de sintomas? Qual pacote usou? Poderia sugerir alguma literatura? Tutorial? Site? Muito obrigado Marcelo

2 1

Pacote grDevices - função plotmath
by Luiz Paulo Serrano 17 Fev '23

17 Fev '23

Prezadas (os), bom dia! Estava usando normalmente em legendas nos gráficos as notações %+-%, %~~% da função plotmath. Atualizei a versão do R para 4.2.2 e não consigo mais usar essas notações. Algum dos colegas está enfrentando esse problema? Como posso resolver? Solicito a habitual atenção. Forte abraço, Luiz Paulo Serrano 31 99899-0866

2 1

Pesquisa online
by Glaucia Ferreira 14 Fev '23

14 Fev '23

Gente desculpa incomodar com um assunto diferente do escopo do grupo mas eu estou com essa pesquisa em andamento se puderem responder e me ajudar a divulgar agradeço imensamente 📢 Galera gostaria de 2 minutos do seu tempo pra responder essa pesquisa sobre mercado de trabalho 👩🏽‍💻👨🏽‍💻 Para acessar, clique aqui: https://pt.surveymonkey.com/r/trabalhobr👆🏾🗂 As informações coletadas terão o sigilo preservado. Em caso de dúvidas, nos escreva no endereço de email br.gender.pep(a)gmail.com .Muito Obrigada!

1 0

GEOquery package
by Michele Claire Breton 03 Fev '23

03 Fev '23

Boa tarde caros colegas, Espero que todos estejam bem. Eu preciso de ajuda de quem tenha experiência com o GEOquery package. Eu gostaria de fazer a anotação das sondas e entender porque alguns genes possuem 2 ou mais valores de expressão (representados cada um por uma sonda) e com valores muito diferentes entre si, dentro de uma mesma amostra. E saber se este pacote que escolhi para trabalhar é o ideal para este fim, pois não tenho familiaridade com a técnica. Como resultado final, quero valores únicos de expressão de cada gene, em cada amostra, para plotar heatmaps. Por favor, há alguém disponível para me ajudar? Muito obrigada, Michele -- ------------------------------------------------------------------------ *Dra. Michele Claire Breton* Researcher in Bioinformatics PhD in Cellular and Molecular Biology CICS - Health Sciences Research Center Faculty of Health Sciences University of Beira Interior Covilhã - Castelo Branco - Portugal

1 0

(sem assunto)
by Ana Rita Saraiva 27 Jan '23

27 Jan '23

1 0

Simular sobrevivência por KM
by Pedro Emmanuel Alvarenga Americano do Brasil 27 Jan '23

27 Jan '23

Salve Amigos de R, Estou analisando dados intermediários de uma pesquisa que ainda está recrutando participantes sobrevivência por Kaplan-Meier. Esses dados são dados incompletos de 3 anos de seguimento de um estudo com um dispositivo de saúde. Um dos questionamentos é se é razoável interromper o recrutamento agora e somente seguir os já em seguimento. Pra isso, o grupo acho interessante que houvesse uma projeção do riscos e IC até 10 anos de seguimento. Um exemplo de tabela de resultado de interesse saída de um survival::survfit Time N at risk Events 1-Surv 12 186 1 0.004 24 85 1 0.010 36 15 1 0.023 Outro exemplo da saída de outro desfecho Time N at risk Events Cum Hazard 12 213 3 0.014 24 139 5 0.043 36 11 1 0.051 A pergunta é se há uma função que a partir desses resultados ou dos dados de onde vieram esses resultado que consigam projetar a sobrevivência e o Hazard até 10 anos de seguimento. Olhei por ai mas não achei nada que me ajudasse. A maioria dos que achei simula comparações de curvas KM. Mas não uma única curva e os ICs. Alguma luz? Abraço forte, Pedro Brasil

1 0

Cluster - agrupar grupos de pessoas
by Chiara Lubich 20 Jan '23

20 Jan '23

Boa tarde, pessoal Estou com dificuldade para entender como fazer para juntar os dados e formar grupos na análise e mostrar graficamente isso. Estou usando o seguinte banco de dados: encr.pw/Gs0ZK Estou usando o seguinte comando: dac<-read.table("clipboard",sep="\t", header=T, dec=".", row.names=1) dac attach(dac) # Ward Hierarchical Clustering d<-vegdist(dac[,-1],distance="jaccard")# distance matrix fit <- hclust(d, method="ward.D2") plot(fit) # display dendogram groups <- cutree(fit, k=5) # cut tree into 5 clusters # draw dendogram with red borders around the 5 clusters rect.hclust(fit, k=5, border="red") que gera esse gráfico: [image: image.png] No entanto, queria que formasse de acordo como nome dos grupos, semelhante a esse gráfico, mas em função dos meus dois atores avaliados (pescador e comunidade). [image: image.png] Só que no meu caso, ao invés de ser as cidades, são os grupos de pessoas. Aos que puderem ajudar eu agradeço Muito obrigada Abraços,

3 5

Salário médio por raça
by Elias Carvalho 27 Out '22

27 Out '22

Alguem saberia me informar por obséquio onde consigo salário médio por raça e estado dos estados unidos? -- *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

Somar linhas de colunas específicas
by Rogério Santos 25 Out '22

25 Out '22

Boa tarde. Estou tentando somar valores nas linhas mas de colunas específicas tentei usar o comando starts_with( ) para selecionar colunas que começam com letras específicas (ex: A1, A2, A3) sendo que essas colunas não estão nessa ordem, estão misturadas com outras colunas(ex: A1, B1, B2, A2, C1, C2, C3, A3...) e quando tento somar as colunas específicas que selecionei usando o comando starts_with dá erro dizendo que as colunas não são numéricas. Alguém sabe como resolver esse problema? Desde já mto obrigado! -- ------------------------------------------------------------ *Rogério Santos* *Estatístico* ------------------------------------------------------------

2 1

dúvida gráfico ggplot2
by Diogo Jerônimo 17 Out '22

17 Out '22

Boa tarde grupo, acho que é uma dúvida boba, mas não estou conseguindo resolver. Segue o exemplo desse gráfico que sai via ggplot2: grupo<-rep(c(1,2,3),8)nível<-rep(c(rep("I",3),rep("II",3),rep("III",3),rep("IV",3)),2)lote<-c(rep("L1",12),rep("L2",12))res<-rnorm(24,100,30)datos<-data.frame(cbind(grupo,nível,lote,res))datos$res<-as.numeric(datos$res) ggplot(datos,aes(y=res,x=nível))+ geom_point(aes(colour=grupo))+ facet_wrap(~lote) O resultado gráfico é algo assim: Eu gostaria de passar uma linha que ligasse as médias dos dados de cada nível, segundo cada lote, ou seja ligando a média do nível I a média do nível II... Teriam como me dar uma ajuda? Muito obrigado!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

2 1

Adicionar média no gráfico box plot
by ana paula coelho madeira 04 Out '22

04 Out '22

Prezados, Bom dia! Estou analisando um experimento com dois fatores e fiz um box plot para verificar o comportamento de um fator dentro do outro usando a função bwplot(peso.seco~F2|F1...) do pacote lattice. Gostaria de indicar a média nesse gráfico, porém não estou conseguindo. Encontrei muitos materiais com o ggplot. Alguma sugestão? Desde já agradeço. Atenciosamente, Ana Paula

1 0

Adicionar a média a um gráfico boxplot - lattice
by ana paula coelho madeira 03 Out '22

03 Out '22

Prezados, Bom dia! Estou analisando um experimento com dois fatores e fiz um box plot para verificar o comportamento de um fator dentro do outro usando a função bwplot(peso.seco~F2|F1...) do pacote lattice. Gostaria de indicar a média nesse gráfico, porém não estou conseguindo. Encontrei muitos materiais com o ggplot. Alguma sugestão? Desde já agradeço. Atenciosamente, Ana Paula

1 0

Importar formato json para o R
by Diogo Jerônimo 28 Set '22

28 Set '22

Prezados, bom dia!!! Novamente o chato rs. De forma direta, eu tenho essa consulta que eu fiz com essas packages: library("httr")library("jsonlite")library("data.table") url<-str_c("https://servicodados.ibge.gov.br/api/v1/localidades/estados/RJ/regioes-metr…")r<-GET(url) r_text<-content(r, as="text") data_json<-fromJSON(r_text, flatten=TRUE)data_json<-as.data.frame(data_json) Em resumo, dessa API do IBGE, busquei as informações das regiões metropolitanas, e as importei para dentro do R em dataframe, para realizar uma análise. Problema: quando ele vêm, ele está com a seguinte cara: Sendo que em municipios, está a informação dos municipios. Eu até consigo (via comando rbindlist) "quebrar" esse info de municipios assim: data_json_municipios<-rbindlist(data_json$municipios) Mas eu queria, no mesmo banco de origem (chamado aqui de data_json), expandir para o banco essa informações dos municípios. Alguém poderia me ajudar nessa? Obrigado!!! Diogo Jerônimo Bacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

2 1

Re: [R-br] Ajudar com organização de dados horários
by sznelwar＠uol.com.br 18 Set '22

18 Set '22

1 0

Re: [R-br] Ajudar com organização de dados horários
by sznelwar＠uol.com.br 17 Set '22

17 Set '22

2 1

Re: [R-br] Ajudar com organização de dados horários
by sznelwar＠uol.com.br 16 Set '22

16 Set '22

2 1

Ajudar com organização de dados horários
by Bruce Kelly 14 Set '22

14 Set '22

Boa noite pessoal, Estou com dúvida de como acrescentar a coluna horaria de algumas. Estou usando as seguintes scritp: rm(list=ls()) # library(lubridate) library(tidyverse) library(dplyr) library(hrbrthemes) library(GGally) library(ggplot2) # Dados ------------------------------------------------------------------- # Chuva input<- read.table('Rain_1.txt', header = T, sep = ',', dec = '.', na.strings = c('NAN', 'NaN', 'Inf', 'NA')) # dec = '.', na.strings = c('NAN', 'NaN', 'Inf')) # Organizar os dados ------------------------------------------------------ # converter o TIMESTEP em data input$TIMESTAMP<- as.Date(ymd_hms(input$TIMESTAMP)) #tentei input$TIMESTAMP<- as.Date(input$TIMESTAMP)) # # #Adicionado os meses e anos nas observações #### a saída é em um caracter converter para numerico input$Year<-format(input$TIMESTAMP, '%Y') input$Month<-format(input$TIMESTAMP, '%m') input$Day<-format(input$TIMESTAMP, '%d') input$Doy<-format(input$TIMESTAMP, '%J') input$Hour.s<-format(input$TIMESTAMP, '%H') Quando executo a coluna Hour.s fica zerada como mostra a figura. E preciso que essa coluna horaria com mediação a cada 30min, fique na variação de 0-23h (0.5, 1, 1.5,...23,5) [cid:66d78bc1-f1a2-49bc-a0c7-85bb36bcd032] estrutura do TIMESTAMP antes de ser convertido em data [cid:b17558fa-d648-4fbd-afb0-24424c5586e8] Obrigado pela ajuda e atenção de todo (@s)

3 4

Ajuda com manipulação de dados horários
by Bruce Kelly 13 Set '22

13 Set '22

Boa tarde listeiros, Estou uma base de dados que tem a coluna TIMESTAMP em ymd_hms, estou criando as colunas meses, dia e hora. Da seguinte forma: input$Year<-format(input2$TIMESTAMP, '%Y') input$Month<-format(input2$TIMESTAMP, '%m') input$Day<-format(input2$TIMESTAMP, '%d') input$Hour.S<-formart(input$TIMESTAMP, ‘%H’) No entanto, a coluna da hora aparece zerada. Lembrando que fiz a conversão de string para data: input$TIMESTAMP<- as.Date(ymd_hms(input$TIMESTAMP)) input2$TIMESTAMP<- as.Date(input$TIMESTAMP) Ambas meu resultado é 00. Eu preciso que minha saída seja de 00-24, por exemplo, 0.5;1;1.5;...23.5 Muito Obrigados a todos(@s)

2 1

apache arrow erro na instacao ubuntu 22
by marco aurelio kasmin 26 Ago '22

26 Ago '22

Por favor, alguém sabe resolver o erro de instalação do arrow no ubuntu 22? Eu olhei o site do apache arrow que fala sobre o erro mas não sei como resolver para o ubuntu. installing *source* package ‘arrow’ ... ** package ‘arrow’ successfully unpacked and MD5 sums checked ** using staged installation *** Arrow C++ libraries found via pkg-config at /usr/lib/R/bin/config: 1: eval: make: not found ------------------------- NOTE --------------------------- There was an issue preparing the Arrow C++ libraries. See https://arrow.apache.org/docs/r/articles/install.html --------------------------------------------------------- ERROR: configuration failed for package ‘arrow’ * removing ‘/home/mak/R/x86_64-pc-linux-gnu-library/4.2/arrow’ * restoring previous ‘/home/mak/R/x86_64-pc-linux-gnu-library/4.2/arrow’ Warning in install.packages : installation of package ‘arrow’ had non-zero exit status

2 1

Exportação de imagem no RStudio
by Marília Teresinha 12 Ago '22

12 Ago '22

Pessoal, boa tarde. Não estou conseguindo exportar um gráfico como imagem no RStudio. Aparece a seguinte mensagem: "TIFFOpen: C://Figura4.tiff: Cannot open." O mesmo acontece quando tento pelo comando tiff: tiff("Figura 4.tiff", units="in", width=10, height=5, res=150)dev.off() Aparece a mensagem: "Warning message:In dev.off() : unable to open TIFF file 'C://Figura4.tiff' O mesmo vale para qualquer formato de imagem (tiff, jpeg, png...).Esse problema começou depois da instalação da versão atual do RStudio. Alguém sabe o está acontecendo e alguma dica de como resolver?Obrigada. Marilia T. de S. Machado, Dra. Pesquisadora Colaboradora, Centro de Desenvolvimento Sustentável - CDS/UnB.Universidade de Brasília, Campus Universitário Darcy Ribeiro. CV Lattes: http://lattes.cnpq.br/4190451945124694

3 4

Como instalar o edgeR no Rstudio versão 3.6.3
by Michele Claire Breton 05 Ago '22

05 Ago '22

Caríssimos(as) utilizadores do R, bom dia Espero que todos estejam bem. Preciso converter dados de expressão gênica de FPKM para CPM, a fim de comparar estes dados com outros bancos que estão em CPM. Para fazer a conversão, eu deveria usar o pacote edgeR, entretanto, não consigo instalá-lo no RStudio V.3.6.3. a seguinte mensagem aparece: "Warning in install.packages : package ‘edgeR’ is not available (for R version 3.6.3)". Por favor, como posso instalar o pacote? Ou então, há outro pacote que eu possa usar para fazer a conversão que seja compatível com esta versão do R? Muitíssimo obrigada! Michele -- ------------------------------------------------------------------------ *Dra. Michele Claire Breton* Researcher in Bioinformatics PhD in Cellular and Molecular Biology CICS - Health Sciences Research Center Faculty of Health Sciences University of Beira Interior Covilhã - Castelo Branco - Portugal

2 2

Não reconhecimento do data.frame como numérico
by Michele Claire Breton 05 Ago '22

05 Ago '22

Caríssimos(as) utilizadores(as) do R, Estou em outro nível de problema. Quero converter FPKM em log2CPM. Para tanto, fiz o seguinte: dados<-read.csv("/home/Genes interesse NePCa.csv") head(dados) class(dados) *[1] "data.frame"* geneLength <- rowMeans(dados) *Error in rowMeans(dados) : 'x' must be numeric* dadoslog2cpm <- convertCounts(dados, unit = "cpm", log = TRUE, normalize = "tmm") *Error: countsMatrix must be a numeric matrix or dataframe of N genes x M Samples. All columns must be numeric.* Não entendo o que possa estar errado. A tabela de dados está em anexo. Vocês podem me ajudar, por favor? Obrigada, Michele -- ------------------------------------------------------------------------ *Dra. Michele Claire Breton* Researcher in Bioinformatics PhD in Cellular and Molecular Biology CICS - Health Sciences Research Center Faculty of Health Sciences University of Beira Interior Covilhã - Castelo Branco - Portugal

3 6

Cancelar assinatura
by Paulo Amorim 04 Ago '22

04 Ago '22

Olá, tudo bem? Como faço para cancelar a assinatura? já tentei pelo site, mas não tive sucesso. Abraços

1 0

Agrupamento - Listar elementos
by Toni 25 Jul '22

25 Jul '22

Olá Pessoal! Alguém sabe como faço para listar os elementos de um determinado agrupamento? Comando utiizado: fviz_cluster(objeto, geom="point", data = df_utilizado) Nesse caso o objeto tem 23 clusters; Como resultado ele me fornece o agrupamento dos pontos mas eu gostaria de listar quem são os pontos de cada um desses 23 agrupamentos. Tem como fazer isso? OBS: Se trocar o geom para "text" ele mostra os pontos mas é impossível ler isso. Grato, Toni Borges

4 7

Geração de vetor, segundo linha da tabela
by Diogo Jerônimo 05 Jul '22

05 Jul '22

Boa tarde prezados, uma vez mais, solicitando a ajuda dos amigos!! Estou necessitando gerar valores de uma normal, condicionada a parâmetros de interesse, como no exemplo reproduzível abaixo: x<-c("X1","X2","X3","X4")média<-c(100,120,150,180)dp<-c(20,15,10,25)exemplo<-data.frame(cbind(x,média,dp)) O meu objetivo é gerar de forma automática, para cada linha, uma distribuição normal com os parâmetros nessas linhas. No exemplo, para a linha que contêm X1, gerar 1000 valores de uma normal com média=100 e dp=20, para a linha que contêm X2, gerar 1000 valores de uma normal de média=120 e dp=15, e assim sucessivamente, até o último vetor, gerado pelos parâmetros da última linha (média=180 e dp=25). Algo similar a isso: for (i in 1:nrow(exemplo)){ simu[i]<-rnorm(100,exemplo[i,2],exemplo[i,3]) } Alguém teria uma solução para isso? É que o número de linhas dessa tabela "exemplo" poderá variar. Muito obrigado!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

3 3

Re: [R-br] Digest R-br, volume 137, assunto 6
by Alexandre Dos Santos 04 Jul '22

04 Jul '22

Obrigado pela resposta Cesar, Então o maior entrave é a dificultade com a sintaxe e vou tentar trazer o resultado para o mais próximo que eu consigo: library(dplyr) library(tidyr) my.ds <- read.csv(" https://raw.githubusercontent.com/Leprechault/trash/main/test_ants.csv") my.ds$ClassType <- cut(my.ds$AT,breaks=c(-Inf,1,2.9,8.9,24.9,49.9,Inf), right=FALSE,labels=c("Class_0_1","Class_1_3","Class_3_9", "Class_9_25","Class_25_50","Class_50")) if (my.ds$nest<1000 & my.ds$ClassType!="Class_0_1") { #1) Se nest for < 1000 - considerando "Class_1_3" para "Class_50" counts <- my.ds%>%group_by(nest)%>%summarize(N=length(unique(ClassType))) my.ds.3 <- my.ds%>%group_by(nest)%>%dplyr::summarize(ClassType=ClassType, max_AT=max(AT)) my.ds.3$class_max<- as.numeric(gsub('Class_[^_]*[^_]_', '', my.ds.3$ ClassType)) # Extrai o maior valor da classe if(counts$N>1 | my.ds.3$max_AT>my.ds.3$class_max){ #ClassType não mudar por nest durante o tempo ou AT entrar e não sair da classe nests <- unique(rbind(counts$nest,my.ds.3$nest)) my.ds.2 <- my.ds[my.ds$nest %in% nests,] my.ds.2$velocity <- c(0,diff(my.ds.2$AT))/c(0,diff(my.ds.2$days)) # Procurando o último ClassType antes de AT e calcula a `velocity` entre o último ClassType e o atual ClassType para nest my.ds.2 <- %>% group_by(nest,ClassType)%>% summarize(velocity=mean(velocity)) %>% complete(ClassType, tidyr:fill = list(velocity = NA)) %>% fill(velocity, .direction = "downup") # Se atrasar qualquer ClassType em "Class_1_3" para "Class_50" em uma classe maior que o valor de AT no tempo 0, a classe não presente é o valor de velocity do último ClassType } } else { #Se o ninho for >= 1000 - apenas considerando "Class_0_1" my.ds$velocity <- c(0,diff(my.ds$AT))/c(0,diff(my.ds$days)) } Em qua., 29 de jun. de 2022 às 11: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: Criação de uma função complexa com dplyr (Cesar Rabak) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 28 Jun 2022 17:52:46 -0300 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Criação de uma função complexa com dplyr > Message-ID: > < > CAKrF98n75bODbfyc0O5wBg9VOn2PFx+HHExqUDkvJSOoHs6YAA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > A sua descrição está difícil de acompanhar porque é colocado entremeado nas > regras termos que necessitam melhor definição para entendê-las: > > > - considerando "Class_1_3" para "Class_50" ; > - ClassType não mudar por nest durante o tempo ; > - AT entrar mas não sair do ClassType ; > - o último ClassType antes de AT ; > - último ClassType e o atual ClassType para nest ; > - Se atrasar qualquer ClassType em ; > - "Class_1_3" para "Class_50" ; > - valor de AT no tempo 0 ; > - classe não presente ; > - último ClassType para nest ; > > Ademais, parece que sua formulação de velocidade e mesmo a saída desejável > indica, que existam velocidades diferentes para o mesmo nest e classe. > > Por certo esses termos todos acima devem fazer todo o sentido para vc mas > para podermos ajudá-lo, essas definições devem ser quebradas nas operações > que elas implicam, aí naturalmente convergir-se-á para a formulação final > desejada. > > Por exemplo, se ClassType implica em ordem, então necessita-se de listas > ordenadas dessas classe por nest? > > Parece que a formulação para o cálculo de velocidade é mais simples quando > Se AT for >= 50 && Se o ninho for >= 1000 -> a velocidade é > c(0,diff(AT))/ > c(0,diff(days)) por days para cada nest, Que se entendida corretamente sua > explicação, significa que apenas as linhas que tenham classe "Class_0_1", > seriam usadas para esses cálculos. > > Não entendo porque é catenado um zero no num e dem dos cálculos de > velocidade . . . > > Supondo as condições para esse cálculo para classe e ninho, você consegue > escrever essa spec na linguagem R ou tem problemas com a sintaxe para > produzir o resultado? > > HTH > > > On Tue, Jun 28, 2022 at 8:49 AM Alexandre Dos Santos por (R-br) < > r-br(a)listas.c3sl.ufpr.br> wrote: > > > Prezados Membros do R-br, > > > > Meu dados de exemplo são: > > > > library(dplyr) > > > > my.ds <- read.csv(" > > https://raw.githubusercontent.com/Leprechault/trash/main/test_ants.csv") > > > > my.ds > > > > # days nest AT > > > > # 1 0 2 10.9200 > > > > # 2 47 2 22.8600 > > > > # 3 76 2 23.2400 > > > > # 4 0 3 0.1400 > > > > # 5 47 3 0.4800 > > > > # 6 76 3 0.8300 > > > > # 7 118 3 0.8300 > > > > # 8 160 3 0.9400 > > > > # 9 193 3 0.9400 > > > > # 10 227 3 0.9400 > > > > # 11 262 3 0.9400 > > > > # 12 306 3 0.9400 > > > > # 13 355 3 11.9300 > > > > # 14 396 3 12.8100 > > > > # 15 450 3 29.3600 > > > > # 16 0 4 1.0000 > > > > # 17 76 4 1.5600 > > > > # 18 160 4 2.8800 > > > > # 19 193 4 2.8800 > > > > # 20 227 4 2.8800 > > > > # 21 262 4 2.8800 > > > > # 22 306 4 2.8800 > > > > # 23 355 4 17.0500 > > > > # 24 0 7 11.7100 > > > > # 25 47 7 24.7900 > > > > # 26 76 7 33.1200 > > > > # 27 0 10 4.1400 > > > > # 28 47 10 6.3000 > > > > # 29 76 10 13.9400 > > > > # 30 118 10 13.3000 > > > > # 31 160 10 29.4200 > > > > # 32 193 10 29.6000 > > > > # 33 227 10 45.1500 > > > > # 34 262 10 45.1500 > > > > # 35 306 10 45.1500 > > > > # 36 355 10 70.5300 > > > > # 37 396 10 109.2600 > > > > # 38 450 10 149.2000 > > > > # 39 0 1066 0.0289 > > > > # 40 29 1066 0.0870 > > > > # 41 70 1066 0.1254 > > > > # 42 112 1066 0.1254 > > > > # 43 146 1066 0.1254 > > > > # 44 180 1066 0.1254 > > > > # 45 215 1066 0.1350 > > > > # 46 259 1066 0.1350 > > > > # 47 0 1067 0.0520 > > > > # 48 29 1067 0.1254 > > > > # 49 70 1067 0.1440 > > > > # 50 112 1067 0.2160 > > > > # 51 146 1067 0.2538 > > > > # 52 180 1067 0.2538 > > > > # 53 215 1067 0.9600 > > > > # 54 259 1067 0.9600 > > > > # 55 349 1067 0.9600 > > > > # Class size > > > > my.ds$ClassType <- cut(my.ds$AT,breaks=c(-Inf,1,2.9,8.9,24.9,49.9,Inf), > > right=FALSE, > > > > labels=c("Class_0_1","Class_1_3","Class_3_9", > > "Class_9_25","Class_25_50","Class_50")) > > > > > > Eu gostaria de calular a variável velocity com as seguintes regras: > > > > Se nest for < 1000 - considerando "Class_1_3" para "Class_50" > > > > 1.1) Se AT for < 50 > > > > 1.1.1) Se ClassType não mudar por nest durante o tempo ou AT entrar mas > > não sair do ClassType, então não calcule velocity (no meu caso nest 2 e 7 > > )., > > > > 1.1.2) Procurando o último ClassType antes de AT e calcula a `velocity` > > entre o último ClassType e o atual ClassType para nest > > > > 1.1.3) Se atrasar qualquer ClassType em "Class_1_3" para "Class_50" em > uma > > classe maior que o valor de AT no tempo 0, a classe não presente é o > valor > > de velocity do último ClassType para nest. > > 1.2) Se AT for >= 50 > > > > 1.2.1) velocity é c(0,diff(AT))/c(0,diff(days)) por days para cada nest > > Se o ninho for >= 1000 - apenas considerando "Class_0_1" > > > > 2.1) a velocidade é c(0,diff(AT))/c(0,diff(days)) por days para cada nest > > também. > > > > Meu output desejável é: > > > > # nest ClassType velocity > > > > # 3 Classe_ 1_3 0,224285714 > > > > # 3 Classe_ 3_9 0,224285714 > > > > # 3 Classe_ 9_25 0,224285714 > > > > # 3 Classe_ 9_25 0,021463415 > > > > # 3 Classe_ 9_25 0,306481481 > > > > # 4 Classe_ 1_3 0,0075 > > > > # 4 Classe_ 1_3 0,015714286 > > > > # 4 Classe_ 1_3 0 > > > > # 4 Classe_ 1_3 0 > > > > # 4 Classe_ 1_3 0 > > > > # 4 Classe_ 1_3 0 > > > > # 4 Classe_ 1_3 0,289183673 > > > > # 4 Classe_ 3_9 0,289183673 > > > > # 10 Classe_9_25 0,263448276 > > > > # 10 Classe_9_25 -0,015238095 > > > > # 10 Classe_9_25 0,383809524 > > > > # 10 Class_25_50 0,383809524 > > > > # 10 Classe_25_50 0,005454545 > > > > # 10 Class_25_50 0,457352941 > > > > # 10 Classe_25_50 0 > > > > # 10 Classe_25_50 0 > > > > # 10 Class_25_50 0,517959184 > > > > # 10 Class_50 0,517959184 > > > > # 10 Class_50 0,944634146 > > > > # 10 Class_50 0,73962963 > > > > # 1066 Classe_0_1 0,002003448 > > > > # 1066 Classe_0_1 0,000936585 > > > > # 1066 Classe_0_1 0 > > > > # 1066 Classe_0_1 0 > > > > # 1066 Classe_0_1 0 > > > > # 1066 Classe_0_1 0,000274286 > > > > # 1066 Classe_0_1 0 > > > > # 1067 Classe_0_1 0,002531034 > > > > # 1067 Classe_0_1 0,000453659 > > > > # 1067 Classe_0_1 0,001714286 > > > > # 1067 Classe_0_1 0,001111765 > > > > # 1067 Classe_0_1 0 > > > > # 1067 Classe_0_1 0.020177143 > > > > # 1067 Classe_0_1 0 > > > > # 1067 Classe_0_1 0 > > > > Por favor, alguém poderia me ajudar com isso? > > > > > > > > -- > > > > -- > > Alexandre dos Santos > > Geotechnologies and Spatial Statistics applied to Forest Entomology > > Instituto Federal de Mato Grosso (IFMT) - Campus Caceres > > Caixa Postal 244 (PO Box) > > Avenida dos Ramires, s/n - Vila Real > > Caceres - MT - CEP 78201-380 (ZIP code) > > Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 > > Lattes CV: http://lattes.cnpq.br/1360403201088680 > > OrcID: orcid.org/0000-0001-8232-6722 > > ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 > > Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ > > -- > > > > _______________________________________________ > > 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/20220628/7cd50a50/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 137, assunto 6 > ***************************************** > -- -- Alexandre dos Santos Geotechnologies and Spatial Statistics applied to Forest Entomology Instituto Federal de Mato Grosso (IFMT) - Campus Caceres Caixa Postal 244 (PO Box) Avenida dos Ramires, s/n - Vila Real Caceres - MT - CEP 78201-380 (ZIP code) Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 Lattes CV: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ --

1 0

Criação de uma função complexa com dplyr
by Alexandre Dos Santos 28 Jun '22

28 Jun '22

Prezados Membros do R-br, Meu dados de exemplo são: library(dplyr) my.ds <- read.csv(" https://raw.githubusercontent.com/Leprechault/trash/main/test_ants.csv") my.ds # days nest AT # 1 0 2 10.9200 # 2 47 2 22.8600 # 3 76 2 23.2400 # 4 0 3 0.1400 # 5 47 3 0.4800 # 6 76 3 0.8300 # 7 118 3 0.8300 # 8 160 3 0.9400 # 9 193 3 0.9400 # 10 227 3 0.9400 # 11 262 3 0.9400 # 12 306 3 0.9400 # 13 355 3 11.9300 # 14 396 3 12.8100 # 15 450 3 29.3600 # 16 0 4 1.0000 # 17 76 4 1.5600 # 18 160 4 2.8800 # 19 193 4 2.8800 # 20 227 4 2.8800 # 21 262 4 2.8800 # 22 306 4 2.8800 # 23 355 4 17.0500 # 24 0 7 11.7100 # 25 47 7 24.7900 # 26 76 7 33.1200 # 27 0 10 4.1400 # 28 47 10 6.3000 # 29 76 10 13.9400 # 30 118 10 13.3000 # 31 160 10 29.4200 # 32 193 10 29.6000 # 33 227 10 45.1500 # 34 262 10 45.1500 # 35 306 10 45.1500 # 36 355 10 70.5300 # 37 396 10 109.2600 # 38 450 10 149.2000 # 39 0 1066 0.0289 # 40 29 1066 0.0870 # 41 70 1066 0.1254 # 42 112 1066 0.1254 # 43 146 1066 0.1254 # 44 180 1066 0.1254 # 45 215 1066 0.1350 # 46 259 1066 0.1350 # 47 0 1067 0.0520 # 48 29 1067 0.1254 # 49 70 1067 0.1440 # 50 112 1067 0.2160 # 51 146 1067 0.2538 # 52 180 1067 0.2538 # 53 215 1067 0.9600 # 54 259 1067 0.9600 # 55 349 1067 0.9600 # Class size my.ds$ClassType <- cut(my.ds$AT,breaks=c(-Inf,1,2.9,8.9,24.9,49.9,Inf),right =FALSE, labels=c("Class_0_1","Class_1_3","Class_3_9", "Class_9_25","Class_25_50","Class_50")) Eu gostaria de calular a variável velocity com as seguintes regras: Se nest for < 1000 - considerando "Class_1_3" para "Class_50" 1.1) Se AT for < 50 1.1.1) Se ClassType não mudar por nest durante o tempo ou AT entrar mas não sair do ClassType, então não calcule velocity (no meu caso nest 2 e 7)., 1.1.2) Procurando o último ClassType antes de AT e calcula a `velocity` entre o último ClassType e o atual ClassType para nest 1.1.3) Se atrasar qualquer ClassType em "Class_1_3" para "Class_50" em uma classe maior que o valor de AT no tempo 0, a classe não presente é o valor de velocity do último ClassType para nest. 1.2) Se AT for >= 50 1.2.1) velocity é c(0,diff(AT))/c(0,diff(days)) por days para cada nest Se o ninho for >= 1000 - apenas considerando "Class_0_1" 2.1) a velocidade é c(0,diff(AT))/c(0,diff(days)) por days para cada nest também. Meu output desejável é: # nest ClassType velocity # 3 Classe_ 1_3 0,224285714 # 3 Classe_ 3_9 0,224285714 # 3 Classe_ 9_25 0,224285714 # 3 Classe_ 9_25 0,021463415 # 3 Classe_ 9_25 0,306481481 # 4 Classe_ 1_3 0,0075 # 4 Classe_ 1_3 0,015714286 # 4 Classe_ 1_3 0 # 4 Classe_ 1_3 0 # 4 Classe_ 1_3 0 # 4 Classe_ 1_3 0 # 4 Classe_ 1_3 0,289183673 # 4 Classe_ 3_9 0,289183673 # 10 Classe_9_25 0,263448276 # 10 Classe_9_25 -0,015238095 # 10 Classe_9_25 0,383809524 # 10 Class_25_50 0,383809524 # 10 Classe_25_50 0,005454545 # 10 Class_25_50 0,457352941 # 10 Classe_25_50 0 # 10 Classe_25_50 0 # 10 Class_25_50 0,517959184 # 10 Class_50 0,517959184 # 10 Class_50 0,944634146 # 10 Class_50 0,73962963 # 1066 Classe_0_1 0,002003448 # 1066 Classe_0_1 0,000936585 # 1066 Classe_0_1 0 # 1066 Classe_0_1 0 # 1066 Classe_0_1 0 # 1066 Classe_0_1 0,000274286 # 1066 Classe_0_1 0 # 1067 Classe_0_1 0,002531034 # 1067 Classe_0_1 0,000453659 # 1067 Classe_0_1 0,001714286 # 1067 Classe_0_1 0,001111765 # 1067 Classe_0_1 0 # 1067 Classe_0_1 0.020177143 # 1067 Classe_0_1 0 # 1067 Classe_0_1 0 Por favor, alguém poderia me ajudar com isso? -- -- Alexandre dos Santos Geotechnologies and Spatial Statistics applied to Forest Entomology Instituto Federal de Mato Grosso (IFMT) - Campus Caceres Caixa Postal 244 (PO Box) Avenida dos Ramires, s/n - Vila Real Caceres - MT - CEP 78201-380 (ZIP code) Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 Lattes CV: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ --

2 1

Rcomander
by Emerson Cotta Bodevan 05 Jun '22

05 Jun '22

Prezados. Não estou conseguindo instalar o RComander no Linux ubuntu 20.04. Retorna erro de falta de repositórios e dependências. Alguém teria um passo a passo? A ideia é usar com alunos da educação física. Agradeço qualquer ajuda. *Emerson*

4 5

Loop para criar diferentes bancos de dados usando o mesmo script
by barbaramnascimento＠gmail.com 03 Jun '22

03 Jun '22

Boa tarde a todos! Estou com dificuldade em gerar um loop e espero que alguém possa me ajudar aqui. Estou trabalhando com dados de sobrevivência de animais até 5 idades determinadas e o efeito individual de 19 características dentro de um servidor. A partir de um banco inicial, eu preciso rodar o mesmo script para combinar todas as idades com todas as características, gerando 95 bancos de dados diferentes. Até agora, eu consegui fazer 5 bancos relativos a cada idade através do commandArg: #Loop ---------------------------------------------------------------------------- ------ ## Level 1: DATE OF BIRTH ==> z sdrumthreshold = c("2019-01-01", "2018-01-01", "2017-01-01", "2016-01-01", "2015-01-01") ### Level 2: SURVIVAL AGE ==> t sdactthreshold = c(36, 48, 60, 72, 84) ## Connection with the code in Linux. args = commandArgs(trailingOnly = TRUE) index.sdrum = as.numeric(args[1]) index.sdact = as.numeric(args[2]) Quando eu chamo as variáveis dentro do loop no Linux, o script roda certinho. Segue como ilustração a forma como eu chamei o loop no Linux. #!/bin/bash sdrumth_leader=(1 2 3 4 5) sdactth_leader=(1 2 3 4 5) for z in "${sdrumth_leader[@]}" do for t in "${sdactth_leader[@]}" do if [[ "$z" != "$t" ]] then continue else sbatch Survival.Rcode.bash $z $t #calling indexes. sleep 1 fi done done Porém, tentei fazer a mesma coisa, criando um vetor com o nome das colunas e não tive sucesso: ### Level 3: TYPE TRAITS ==> q sdtypetrait = c("Stature", "Strength", "FinalScore", "DairyForm", "RumpAngle", "RumpWidth", "RearLegs", "ForeUdder", "RearUdderHeight", "RearUdderWidth", "FootAngle", "UdderDepth", "UdderCleft", "TeatPlacement", "TeatLength", "RearTeatPlacementRv", "RearTeatPlacementSv", Pasterns", "RearTeatLength") ## Connection with the code in Linux. args = commandArgs(trailingOnly = TRUE) index.sdrum = as.numeric(args[1]) index.sdact = as.numeric(args[2]) index.typet = as.numeric(args[3]) Quando eu chamo a variável dentro do script, ela vem com as aspas e dá erro. As entradas da variável do vetor são: type13 <- subset(type12, sdtypetrait [index.typet] <= 50) count(type13, sdtypetrait [index.typet]) model <- lm(formula = sdtypetrait [index.typet] ~ CalvingAge + I(CalvingAge^2) + appraisalDIM + I(appraisalDIM^2), data = type15) dados_completo$sdtypetrait [index.typet]_adj <- dados_completo$residuals dados_completo$sdtypetrait [index.typet]_adj2 <- (dados_completo$sdtypetrait [index.typet]_adj)*2 dados_usar <- select(dados_completo, RegId, SireRegId, DamRegId, CG, Cross, AFC, AFC2, Inbreeding, Inbreed2, sdtypetrait [index.typet]_adj, sdtypetrait [index.typet]_adj2, Survival) write.table(dados_usar, file = paste0("/teste/", sdactthreshold[index.sdact],"/", sdtypetrait [index.typet],"/Survival.dat"), sep = " ", row.names = F, quote = F, col.names = F) Será que alguém tem alguma ideia de como eu posso chamar a variável sem entrar as aspas, ou se a forma como estou rodando tem algum erro? Obrigada! ____________________________________ Barbara Mazetti Nascimento Postdoctoral Researcher Department of Animal and Dairy Sciences University of Wisconsin-Madison

2 3

Re: [R-br] Digest R-br, volume 136, assunto 10
by Andrea Vita Reis Mendonca 02 Jun '22

02 Jun '22

Ferreira, E. B., Cavalcanti, P. P., Nogueira, D. A. (2021). ExpDes.pt: Experimental Designs pacakge (Portuguese). Rpackage version 1.2.2. Disponível em: https://cran.r-project.org/web/packages/ExpDes.pt/ExpDes.pt.pdf. Acesso em: 25/05/2022 Em ter., 31 de mai. de 2022 à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. Trifatorial com trat adicionais (Alessandro Dal Col Lúcio) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 31 May 2022 11:27:44 -0300 > From: Alessandro Dal Col Lúcio <adlucio(a)ufsm.br> > To: a lista Brasileira oficial de discussÃ£o do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Trifatorial com trat adicionais > Message-ID: > < > CAHVsyO4fQU4z_jBpT5cx2tTqmX2o0x4mB4Wxab09tMxyoSDojg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Bom dia. > Alguma dica de pacote para analisar experimentos trifatoriais com mais de > um tratamento adicional? > Abraços. > Alessandro > -- > Enviado via Windows Phone. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20220531/65ec63c9/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 136, assunto 10 > ****************************************** >

1 0

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

Otimizar código - executar em lote
by Marcelo Laia 24 Abr '22

24 Abr '22

Olá Saberia me ajudar com uma otimização (execução em lote) do código abaixo? A análise é a seguinte (exemplo): 1 - pego os dados da coluna 25 e obtenho o lambda (boxcox); 2 - crio uma nova coluna com os dados corrigidos 3 - aplico um lm 4 - checo as pressuposições 5 - calculo e jogo em um objeto o boxplot (função ggboxplot) 6 - ploto o boxplot 7 - armazeno no hd o boxplot com o nome da coluna lida no passo 1 acima 8 - vou para a coluna seguinte, 26 Estou fazendo isso com o código abaixo. Só que eu tenho 120 colunas. Por ser muito igual, já fiz confusão umas quatro ou cinco vezes com os valores, com os nomes dos arquivos ou com os nomes das colunas. Há como otimizar isso ou é melhor ir na "unha" mesmo? Obrigado! #### LINALOL #### colnames(dados01)[25] # retorna o nome da coluna boxcox(LINALOL ~ Genotipo * Tratamento, data=dados01, plotit = TRUE) boxcox(LINALOL ~ Genotipo * Tratamento, data=dados01, plotit = TRUE, lam = seq(-1, 1, 1/10)) boxcox(LINALOL ~ Genotipo * Tratamento, data=dados01, plotit = TRUE, lam = seq(-0.7, -0.65, 1/1000)) ## vou ajustando até achar o lambda mais próximo ## o valor aproximado de lambda é -0.655 dados01$LINALOL.boxcox <- (dados01$LINALOL^(-0.655) - 1)/(-0.655) fit25 <- lm(LINALOL.boxcox ~ Genotipo * Tratamento, data=dados01) summary(fit25) anova(fit25) ggqqplot(residuals(fit25)) # Os pontos devem estar ao longo da linha # de referência. Senão, não é normal shapiro_test(residuals(fit25)) # se p-value for significativo, # não é normal ### teste de homocedacidade das variâncias ### se p < 0.05 (significativo), as variâncias não são homogêneas dados01 %>% levene_test(LINALOL.boxcox ~ Genotipo * Tratamento) bxp <- ggboxplot( dados01, x = "Genotipo", y = "LINALOL.boxcox", color = "Tratamento", palette = "Set1", ylab = "LINALOL - boxcox", ) bxp pdf(file="LINALOL_boxcox.pdf") bxp dev.off() Obrigado -- Marcelo

2 1

Re: [R-br] Digest R-br, volume 135, assunto 1
by Samuel Carvalho 07 Abr '22

07 Abr '22

Obrigado Cid Edson,Estou quase lá. Desafio agora é exportar este resultado para MS Word ou Excel.Pela save_kable ele perde a formatação e não há opção de exportar para estes formatos.Outra tentativa sem sucesso foi forçar a transformação em flextable.Obrigado. ====================================Samuel P. C. Carvalho Engº Florestal [UFLA]Mestre em Ciências Florestais [UFLA]Doutor em Recursos Florestais [ESALQ/USP]============================================= Em quarta-feira, 6 de abril de 2022 12:00:52 BRT, <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. texto formatado em tabela (Samuel Carvalho) 2. Re: texto formatado em tabela (Cid Póvoas) ---------------------------------------------------------------------- Message: 1 Date: Tue, 5 Apr 2022 18:37:11 +0000 (UTC) From: Samuel Carvalho <samukajm(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] texto formatado em tabela Message-ID: <1064409059.950553.1649183831762(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Olá pessoal,Gostaria de uma ajuda de vocês para criar uma variável que tem algumas regras de formatação. No caso são nomes científicos de plantas.Envio o database bryophytes como exemplo. Nesta nova variável, as colunas Genus e Species são em itálico, a coluna Var normal e a coluna Intraspecific em itálico.Aqui vai um exemplo de uma tentativa sem sucesso: library(tidyverse); library(kableExtra); data("bryophytes")bryophytes %>% filter(!is.na(Intraspecific)) %>% kbl() %>% kable_styling() %>% column_spec(c(2,3,5), italic = T) %>% mutate(Scientific.name = paste(Genus, Species, Var, Intraspecific)) # Aqui esta o erro Obrigado.At.te. ====================================Samuel P. C. Carvalho Engº Florestal [UFLA]Mestre em Ciências Florestais [UFLA]Doutor em Recursos Florestais [ESALQ/USP]============================================= -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20220405/cb4290df/atta…> ------------------------------ Message: 2 Date: Tue, 5 Apr 2022 18:25:35 -0300 From: Cid Póvoas <cidedson(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] texto formatado em tabela Message-ID: <CAAmR27-X6ZODnUiL8GrA0944syvZnevyyp7kqjwJLJ8Vn39fmA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" library(tidyverse); library(kableExtra);library(Taxonstand) data("bryophytes") bryophytes %>% filter(Intraspecific!="") %>% mutate(Scientific.name = paste( paste0("$",Genus,"$$~$"), paste0("$",Species,"$$~$"), Var, paste0("$",Intraspecific,"$$~$","."))) %>% kbl(longtable = F, "html", escape = F) %>% kable_styling() *Cid Edson Mendonça Póvoas* *AnovAgro <http://www.anovagro.com/>* *Engenheiro Agrônomo - **Data Scientist* *CREA-BA: 051984991-4* *Técnico em Segurança do Trabalho * *Nº: **0012669/BA* *Tel: +55 73 99151-9565* *Lattes : *http://lattes.cnpq.br/2303498368142537 *LinkedIn :* http://br.linkedin.com/in/cidedson/ *Whatsapp :* https://wa.me/5573991519565 Em ter., 5 de abr. de 2022 às 15:37, Samuel Carvalho por (R-br) < r-br(a)listas.c3sl.ufpr.br> escreveu: > Olá pessoal, > Gostaria de uma ajuda de vocês para criar uma variável que tem algumas > regras de formatação. No caso são nomes científicos de plantas. > Envio o database bryophytes como exemplo. Nesta nova variável, as colunas > Genus e Species são em itálico, a coluna Var normal e a coluna > Intraspecific em itálico. > Aqui vai um exemplo de uma tentativa sem sucesso: > > library(tidyverse); library(kableExtra); > data("bryophytes") > bryophytes %>% > filter(!is.na(Intraspecific)) %>% > kbl() %>% > kable_styling() %>% > column_spec(c(2,3,5), italic = T) %>% > mutate(Scientific.name = paste(Genus, Species, Var, Intraspecific)) # > Aqui esta o erro > > Obrigado. > At.te. > > > > *====================================* > > *Samuel P. C. Carvalho*Engº Florestal [UFLA] > Mestre em Ciências Florestais [UFLA] > Doutor em Recursos Florestais [ESALQ/USP] > ============================================= > _______________________________________________ > 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/20220405/f0c768de/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 135, assunto 1 *****************************************

2 1

texto formatado em tabela
by Samuel Carvalho 05 Abr '22

05 Abr '22

Olá pessoal,Gostaria de uma ajuda de vocês para criar uma variável que tem algumas regras de formatação. No caso são nomes científicos de plantas.Envio o database bryophytes como exemplo. Nesta nova variável, as colunas Genus e Species são em itálico, a coluna Var normal e a coluna Intraspecific em itálico.Aqui vai um exemplo de uma tentativa sem sucesso: library(tidyverse); library(kableExtra); data("bryophytes")bryophytes %>% filter(!is.na(Intraspecific)) %>% kbl() %>% kable_styling() %>% column_spec(c(2,3,5), italic = T) %>% mutate(Scientific.name = paste(Genus, Species, Var, Intraspecific)) # Aqui esta o erro Obrigado.At.te. ====================================Samuel P. C. Carvalho Engº Florestal [UFLA]Mestre em Ciências Florestais [UFLA]Doutor em Recursos Florestais [ESALQ/USP]=============================================

2 1

Ajuda em análise de cluster p/ tese
by Suelen Takano 30 Mar '22

30 Mar '22

Bom dia pessoal, tudo bem? Gostaria de solicitar uma ajuda ao Luis Alberto, que está desenvolvendo uma tese e precisa de ajuda na replicação de algumas análises em R... O trabalho dele tem como base o objetivo o desenvolvimento de um índice composto para a Bioeconomia Circular, identificando de um conjunto aplicável de indicadores utilizando análise multicritério. É muito interessante, mas infelizmente não estou com disponibilidade para ajudá-lo. Alguém poderia auxiliá-lo? Ele segue em cópia aqui no email e tem total disponibilidade para falar esclarecer quaisquer dúvidas! Agradeço desde já! Att, Suelen Takano.

1 0

Buscar conjunto de códigos para executar rotina no R.
by Diogo Jerônimo 23 Mar '22

23 Mar '22

Boa tarde a todos!!! Novamente agradecido pela ajuda anterior, e já pedindo desculpas por nova perturbação, mas de novo apareceu uma situação que não sei resolver. Segue o código R parecido com o outro para baixar dados via URL, onde agora x1 representa o código de modulo de um serviço, x2 o ano e x3 o código da rodada trimestral de realização desse serviço:baixa_bancos<-function(x1,x2,x3){ dados<-read.csv2(paste('http://.../arquivo_m", x1 ,"_a",x2,"_s", x3,".csv'",sep=""),header=T, sep="|", dec =".", encoding = "UTF-8")} A chamada da macro fica algo assim: a<- vetor de código de módulos; b<-vetor de anos; c<-???? for(i in 1:length(a)); for(j in 1:length(b)); for(k in 1:length(c)){x1<-a[[i]]; x2<-b[[j]]; x3<-c[[k]]baixa_bancos(x1,x2,x3)} Coloquei interrogação no vetor c em vermelho por causa da dúvida. Os códigos de módulo e ano (óbvio) são fixos, basta vetorizar. Mas os códigos de realização de serviços trimestrais, identificados por código numérico, podem por conveniência de atraso, não ser iguais, embora sejam sempre quatro. Assim por exemplo: Mód:108;Ano:2020;codserv:1,2,3,5 Mód:108;Ano:2021;codserv:1,2,3,4 Mód:429;Ano:2020;codserv:1,2,3,4Mód:429;Ano:2021;codserv:1,3,4,5 Dessa forma o vetor c depende do código de serviço e do ano em que ele foi realizado. Para me ajudar nisso, o TI daqui me deu uma outra consulta que também baixo por URL, que após alguns ajustes, faço virar uma tabela dentro do RStudio, com essa cara: Resumindo agora: o que eu queria é adaptar meu código da seguinte forma: se viu que o módulo é 108 e o ano é 2020, então ele deve ler como parâmetros para k3 o valores 1,2,3,5; se for módulo 108 e ano 2021, ele deve ler para k3 o valores 1,2,3,4. E assim sucessivamente, até finalizar o baixar das bases de dados requeridas. Eu consegui ser claro? Se sim, alguém saberia como se faz isso? Obrigado novamente pela disponibilidade do time!!! Att, Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

2 1

Renomear arquivo dentro do R
by Diogo Jerônimo 22 Mar '22

22 Mar '22

Prezados boa tarde!!! Lá vou eu incomodar os amigos, rs... Estou fazendo uma macro para baixar dados "*.csv" de uma fonte URL, com um formato semelhante a esse (desculpa, não consegui fazer um reproduzível), onde k1 designa meses do ano e k2 é um conjunto de cinco (ou mais) substâncias: baixa_bancos<-function(k1,k2){ dados<-read.csv2(paste('http://.../arquivo_m", k1 ,"_s", k2 ,".csv'",sep=""), header=T, sep="|", dec =".", encoding = "UTF-8") } A chamada da macro fica algo assim: a<-seq(1:12);b<-seq(1:5) for(i in 1:length(a)) for(j in 1:length(b)) { k1<-a[[i]] k2<-b[[j]] baixa_bancos(k1,k2) } A dúvida: se eu executar assim a macro, dados vai sempre ser subscrito pelas informações de outro mês e outra substância. Assim, queria que o arquivo (tabela) fosse salvo no ambiente RStudio de forma similar abaixo (obs: sei que não está certo, é mais para entender o que estou precisando):paste("arquivo_m",k1,"_s",k2)<-dados E em seguida reunir todos eles via rbind ou outro comando. Consegui explicar certo? Se sim, alguém saberia como resolver? Obrigado e bom fim-de-semana!!! Diogo Jerônimo Bacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520

3 5

Correlação e Chi Squared
by Márcio Midon 09 Mar '22

09 Mar '22

Pessoal... Alguém pode me ajudar neste teste? Preciso ajustar dados para um artigo. Matrícula Aluno Especial Cordialmente, *Márcio Midon*Arquiteto e Urbanista UFPel Especialista Geoprocessamento UFPR Mestre Geomática UFSM Doutorando Geografia UFSM Tel 53 3305 7075 | Cel 53 98100 2020 SKYPE marcio_midon

2 3