- R-br - DINF UFPR

Superfície resposta
by Adriano Carvalho Costa 18 Mai '21

18 Mai '21

Boa noite pessoal, estou com um artigo de aquaponia em que estamos fazendo a viabilidade econômica de acordo com o preço de venda dos produtos, sendo as variáveis respostas TIR, VPL entre outros. Estou precisando de um parceiro que interesse em participar como co-autor para fazer o gráfico. O trabalho faz parte de uma dissertação de mestrado. Quem tiver interesse meu whats 64 981186749.

1 0

Fatorial 2^3 + ponto central
by ana paula coelho madeira 23 Abr '21

23 Abr '21

Prezados, Boa tarde! Estou analisando um experimento fatorial , 3 fatores e 2 níveis, adicionando 3 repetições no ponto central. Utilizando o pacote FrF2, a função library(FrF2) plan.pc<-FrF2(nruns=8, nfactors=3, ncenter=3, factor.names=c("AM", "pH", "T"), randomize=F) summary(plan.pc) lm1<-lm(y~A*B*C, iscube(plan.pc), data=plan.pc) # iscube: termo para avaliar a curvatura está apresentando na saída: Residuals: ALL 8 residuals are 0: no residual degrees of freedom! Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.155375 NA NA NA AM 0.009875 NA NA NA pH 0.050875 NA NA NA T -0.103625 NA NA NA AM:pH 0.048875 NA NA NA AM:T -0.006625 NA NA NA pH:T -0.033625 NA NA NA AM:pH:T -0.045125 NA NA NA Residual standard error: NaN on 0 degrees of freedom Multiple R-squared: 1, Adjusted R-squared: NaN F-statistic: NaN on 7 and 0 DF, p-value: NA Ao que parece, por mais que eu tenha acrescentado os 3 pontos centrais, eles não estão sendo considerados. Alguma sugestão? Usava anteriormente o pacote qualityTools, o pacote foi desabilitado. Desde já agradeço Ana Paula

2 1

Análise ordenamento (rank) de dados respostas
by Marcelo Laia 19 Abr '21

19 Abr '21

Colega, Tenho uma tabela com dados na seguinte forma: Tesp1 Tesp2 Tesp3 Ind1 3 1 2 Ind2 1 3 2 Ind3 3 2 1 Ind4 2 1 3 Ind5 3 1 2 Pediu-se aos Ind que ordenasse de acordo com suas preferências. Preciso saber qual é o ranqueamento de Tesp1, Tesp2 e Tesp3. Por exemplo, pelas respostas dos indivíduos (1 a 5), o ranqueamento seria: Tesp2 > Tesp3 > Tesp1 Eu busquei na internet e achei: https://rpubs.com/kyle_macdonald/433105 https://www.programmingr.com/tutorial/rank-in-r/ https://towardsdatascience.com/r-rank-vs-order-753cc7665951 https://stackoverflow.com/questions/31617299/ranking-data-based-on-subcateg… https://www.r-bloggers.com/2019/02/introducing-rank-data-analysis-with-arkh… Você conhece alguma outra maneira de analisar esse tipo de dado? Eu preciso saber a ordem preferida e alguma estatística associada que me permita afirmar algo. Além disso, eu gostaria de saber quem se relaciona com quem. Por exemplo, Tesp1 se relaciona mais com Tesp2 ou com Tesp3? Ou seja, quem escolhe Tesp1, teria mais chances de escolher Tesp2 ou Tesp3? Isso é só um exemplo. O banco de dados é maior. Obrigado -- Marcelo

2 2

Regressão não linear
by Chiara Lubich 18 Abr '21

18 Abr '21

Bom dia, Pessoal, fiz a plotagem dos meu dados: ano n 1994 1 1996 1 2001 1 2002 1 2004 1 2005 3 2006 3 2007 2 2008 4 2009 1 2010 3 2011 2 2012 5 2013 4 2014 5 2015 3 2016 8 2017 2 2018 4 2019 3 2020 12 2021 2 [image: image.png] E vi na internet alguns vídeos em que gráficos que tinham esse tipo de distribuição, era feito o uso de Regressão Não Linear, por meio da função "nls". No entanto, tentei começar a digitar os comandos, mas não entendi como são feitas as escolhas para valores de a=?, b=? e c=?. Segue o script abaixo: plot(n ~ ano, data = chiara, xlab = "Ano da publicação", ylab = "Quantidade de artigos (N)") a_maximovalor=max(n) modelo<-nls(y~x*(1-exp(-b*x))^c, data=chiara, start = list(a= a_maximovalor, b=?, c=?)) Agradeço desde já pela colaboração

3 7

Re: [R-br] Digest R-br, volume 124, assunto 1
by Andrea Vita Reis Mendonca 16 Abr '21

16 Abr '21

regressão não linear Bom dia Não conheço este modelo. Mas tenho a impressão que deveria ser um logístico, porque coloca o A entre os parâmetros estimados e ele não aparece no modelo, acho que o primeiro X é o A. Além deste equívoco, os parâmetros em parte depende do conhecimento do comportamento dos dados pelo pesquisador, A é a assíntota, equivale ao valor de Y referente ao ponto que ele tende a estabilizar. Procure na literatura trabalhos com o mesmo tema do seu e pode lhe ajudar no chute inicial dos parâmetros. Leia teses e dissertações sobre regressão não linear para entender melhor este tema. Em qui., 15 de abr. de 2021 à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. Regressão não linear (Chiara Lubich) > 2. Re: Regressão não linear (Cesar Rabak) > 3. Sugestão de material R (Emerson Cotta Bodevan) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 14 Apr 2021 11:00:14 -0400 > From: Chiara Lubich <lubichchiara(a)gmail.com> > To: a lista Brasileira oficial de discussÃ£o do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Regressão não linear > Message-ID: > <CAH-ZE+DswkC8wt= > 1vpDxM_95YuWHtXgwase8d0574c5HCfp8Rg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Bom dia, > Pessoal, fiz a plotagem dos meu dados: > ano n > 1994 1 > 1996 1 > 2001 1 > 2002 1 > 2004 1 > 2005 3 > 2006 3 > 2007 2 > 2008 4 > 2009 1 > 2010 3 > 2011 2 > 2012 5 > 2013 4 > 2014 5 > 2015 3 > 2016 8 > 2017 2 > 2018 4 > 2019 3 > 2020 12 > 2021 2 > [image: image.png] > E vi na internet alguns vídeos em que gráficos que tinham esse tipo de > distribuição, era feito o uso de Regressão Não Linear, por meio da função > "nls". No entanto, tentei começar a digitar os comandos, mas não entendi > como são feitas as escolhas para valores de a=?, b=? e c=?. > Segue o script abaixo: > plot(n ~ ano, data = chiara, xlab = "Ano da publicação", ylab = "Quantidade > de artigos (N)") > a_maximovalor=max(n) > modelo<-nls(y~x*(1-exp(-b*x))^c, data=chiara, start = list(a= > a_maximovalor, b=?, c=?)) > > Agradeço desde já pela colaboração > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20210414/13050f49/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9748 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20210414/13050f49/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 14 Apr 2021 14:48:24 -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> > Cc: Chiara Lubich <lubichchiara(a)gmail.com> > Subject: Re: [R-br] Regressão não linear > Message-ID: > <CAKrF98= > xscBmoCWH4kVnn-+NfhohtA4wQ8m5Hmy+XhRapnGj6A(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Chiara, > > Você menciona « ? esse tipo de distribuição ? ». > > Qual tipo você acha que é? > > Pelos nomes que deste aos eixos, parece que a variável dependente é > inteira¹ ("Qtde. de artigos [N]"), por outro lado, a variável independente > é um ano calendário, que por definição é uma variável ordinal arbitrária (a > faixa 1994 -- 2021 poderia ser noutros calendários diferente, por exemplo > no Bahá'i de 157 a 178, etc.). > > Para decidir por uma regressão (de qualquer tipo) antes de mais nada > deve-se refletir sobre o fenômeno e a modelagem que o *processo* que gera > esses dados. > > As regressões são abordagens que nos permitem fazer ajustes dos dados > experimentais à essas hipóteses modeladas e enfrenta os dois problemas dos > dados obtidos por meio de observação: a) amostragem; b) perturbações que > geram erros nas observações, modeladas como "ruído" ou "erros" com > distribuição gaussiana de média zero e variância em função da dispersão dos > dados vis-à-vis à abstração matemática (equação da regressão). > > Daí a sua confusão sobre « ? como são feitas as escolhas para valores de > a=?, b=? e c=? » ser esperada. > > A equação da regressão ("não linear") propõe uma complexa equação > exponencia multiplicada pela variável independente, etc. > > O quê você precisa responder é: essa curva descreve um processo que explica > a geração dos dados que está em análise? > > OBS.: Faz sentido incluir um ano que ainda não acabou numa regressão que > conta eventos por unidade de ano calendário? > > HTH > -- > Cesar Rabak > > [1] Mais importante, parece ser o tipo de dados "de contagem" que tem > outras restrições, como não poder ser menor que zero, etc. > > On Wed, Apr 14, 2021 at 12:01 PM Chiara Lubich por (R-br) < > r-br(a)listas.c3sl.ufpr.br> wrote: > > > Bom dia, > > Pessoal, fiz a plotagem dos meu dados: > > ano n > > 1994 1 > > 1996 1 > > 2001 1 > > 2002 1 > > 2004 1 > > 2005 3 > > 2006 3 > > 2007 2 > > 2008 4 > > 2009 1 > > 2010 3 > > 2011 2 > > 2012 5 > > 2013 4 > > 2014 5 > > 2015 3 > > 2016 8 > > 2017 2 > > 2018 4 > > 2019 3 > > 2020 12 > > 2021 2 > > [image: image.png] > > E vi na internet alguns vídeos em que gráficos que tinham esse tipo de > > distribuição, era feito o uso de Regressão Não Linear, por meio da função > > "nls". No entanto, tentei começar a digitar os comandos, mas não entendi > > como são feitas as escolhas para valores de a=?, b=? e c=?. > > Segue o script abaixo: > > plot(n ~ ano, data = chiara, xlab = "Ano da publicação", ylab = > > "Quantidade de artigos (N)") > > a_maximovalor=max(n) > > modelo<-nls(y~x*(1-exp(-b*x))^c, data=chiara, start = list(a= > > a_maximovalor, b=?, c=?)) > > > > Agradeço desde já pela colaboração > > > > _______________________________________________ > > 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/20210414/7e8accc9/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9748 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20210414/7e8accc9/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 14 Apr 2021 19:13:18 -0300 > From: Emerson Cotta Bodevan <bodevan.ec(a)gmail.com> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Sugestão de material R > Message-ID: > < > CAP5i-swAprQOdhWxrH2CoWDnvS+a-8iJdraS3G8W4C-2f6uWOQ(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Prezados, boa noite. > > Estou interessado em trabalhar com o R para gerar mapas temáticos com: > > - Indice de Moran Local e Global e > - Kernel Density > > Gostaria de alguma sugestão de material para leitura dos assuntos acima, > bem como para uso do R para gerar os índices e mapas temáticos. > > Agradeço qualquer sugestão. > > Atenciosamente, > > *Emerson* > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20210414/e0c99258/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 124, assunto 1 > ***************************************** >

1 0

Sugestão de material R
by Emerson Cotta Bodevan 15 Abr '21

15 Abr '21

Prezados, boa noite. Estou interessado em trabalhar com o R para gerar mapas temáticos com: - Indice de Moran Local e Global e - Kernel Density Gostaria de alguma sugestão de material para leitura dos assuntos acima, bem como para uso do R para gerar os índices e mapas temáticos. Agradeço qualquer sugestão. Atenciosamente, *Emerson*

2 2

Erro no gap.boxplot
by Michele Claire Breton 30 Mar '21

30 Mar '21

Caros(as) senhores(as), Estou com dificuldade para usar o gap.boxplot. Eu fiz o dowload do pacote Plotrix e não consigo rodar nem o exemplo da documentação da sintaxe, volta como erro: Error in gap.boxplot(vec, gap = list(topo = c(20, 60)), principal = "Mostrar valores discrepantes separadamente") : could not find function "gap.boxplot" Para o gráfico de barras e o gráfico de linha, rodei o exemplo perfeitamente. Também gerei o boxplot sem o gap no eixo, aí sim, com meus dados. O que estou fazendo de errado? Muito obrigada pela ajuda. -- Com os melhores cumprimentos, *Dra. Michèle Claire Breton* Investigadora auxiliar Faculdade de Ciências da Saúde Universidade da Beira Interior Covilhã - Portugal

5 4

Limpar console!
by Nory 18 Mar '21

18 Mar '21

Olá pessoal, Uso o R dentro do Emacs no meu macbook. Gostaria de saber se alguém tem alguma dica para limpar o console? O "Ctrl + L" no RStudio funciona, mas não para o R dentro do Emacs. -- Nory Daniel de Carvalho Erazo Técnico em Sistemática Vegetal/Herbário INPA Coordenação de Biodiversidade - CBio Instituto Nacional de Pesquisas da Amazônia - INPA Ministério da Ciência, Tecnologia e Inovação - MCTI Av. André Araújo, 2.936 Bairro: Petrópolis, CEP: 69067-375 Manaus - AM, Brasil Telefone: (92) 3643-3125 Celular: (92) 99378-7592 nory(a)inpa.gov.br norydaniel(a)gmail.com

2 2

Recuperar mensagens
by sznelwar＠uol.com.br 15 Mar '21

15 Mar '21

2 1

Identificando e criando nova variável por uma expressão regular
by Felipe Barletta 15 Mar '21

15 Mar '21

Prezados, Preciso de uma ajuda para resolver a seguinte questão: Possuo uma lista de nomes de medicamentos. Preciso identificar seus nomes em um vetor em que há informações de medicamentos preescritos. Para isso estou usando a função grep(). Após identificar, quero incluir no data frame uma nova coluna contendo o nome do medicamento da primeira lista. Segue um código reproduzível ##### Lista de medicamentos que guardei em um objeto chamado patterns patterns <- structure(c(1L, 2L, 3L, 4L, 6L, 7L, 5L, 8L, 9L, 10L, 11L), .Label = c("Oritavancina", "Oxacilina", "Pefloxacino", "Penicilina", "Pexiganan", "Piperacilina", "Piperacilina-tazobactam", "Pirazinamida", "Plazomicina", "Polimixina B", "Posilozid"), class = "factor") ### Amostra do Data frame em que preciso encontrar os nomes da lista acima df <- structure(list(nome = structure(c(9L, 11L, 3L, 3L, 7L, 2L, 4L, 5L, 8L, 6L, 10L, 1L), .Label = c("CLORETO DE POTASSIO DRAGEA 600MG", "CLORETO DE SODIO 0,9% SERINGA PREENCHIDA 5ML", "CLORETO DE SODIO SOLUCAO INJETAVEL 0,9% 10ML", "CODEINA FOSFATO SOLUCAO ORAL 3MGML 10ML ISCMPA @", "CODEINA FOSFATO SOLUCAO ORAL 3MGML 5ML ISCMPA @", "DipiRONA SOLUCAO INJETAVEL 500MGML 2ML", "FUROSEMIDA SOLUCAO INJETAVEL 10MGML 2ML", "HIDROCORTISONA SUCCINATO SODICO PO LIOFILO INJETAVEL 100MG", "ONDANSETRONA CLORIDRATO SOLUCAO INJETAVEL 2MGML 4ML", "Penicilina G BENZATINA PO LIOFILO INJETAVEL 1200000UI", "PIPERACILINA SODICA 4G + TAZOBACTAM SODICA 0,5G PO LIOFILO INJETAVEL" ), class = "factor"), ID = c(13750598L, 13797248L, 13797248L, 13797248L, 13727853L, 13727853L, 13922033L, 13922033L, 13923030L, 13923030L, 13923030L, 13877978L), date = structure(c(1L, 2L, 2L, 2L, 3L, 3L, 4L, 4L, 5L, 5L, 5L, 6L), .Label = c("2021-01-01 00:08:14", "2021-01-01 00:20:39", "2021-01-01 00:22:18", "2021-01-31 23:53:48", "2021-01-31 23:58:50", "2021-01-31 23:59:47"), class = "factor")), row.names = c(1L, 2L, 3L, 4L, 5L, 6L, 583867L, 583868L, 583869L, 583870L, 583871L, 583872L), class = "data.frame") ### Usei esse código para identificar os medicamentos matches <- unlist(sapply(patterns, function(p) grep(p, df$nome, value = FALSE, ignore.case = TRUE) ) ) encontrados <- df[matches,] Porém além de identificar, quero incluir no df uma nova variável identificando o nome do medicamento que consta na lista (patterns) - como segue exemplo Abaixo: dados nome 1 ONDANSETRONA CLORIDRATO SOLUCAO INJETAVEL 2MGML 4ML 2 PIPERACILINA SODICA 4G + TAZOBACTAM SODICA 0,5G PO LIOFILO INJETAVEL 3 CLORETO DE SODIO SOLUCAO INJETAVEL 0,9% 10ML 4 CLORETO DE SODIO SOLUCAO INJETAVEL 0,9% 10ML 5 FUROSEMIDA SOLUCAO INJETAVEL 10MGML 2ML 6 CLORETO DE SODIO 0,9% SERINGA PREENCHIDA 5ML 583867 CODEINA FOSFATO SOLUCAO ORAL 3MGML 10ML ISCMPA @ 583868 CODEINA FOSFATO SOLUCAO ORAL 3MGML 5ML ISCMPA @ 583869 HIDROCORTISONA SUCCINATO SODICO PO LIOFILO INJETAVEL 100MG 583870 DipiRONA SOLUCAO INJETAVEL 500MGML 2ML 583871 Penicilina G BENZATINA PO LIOFILO INJETAVEL 1200000UI 583872 CLORETO DE POTASSIO DRAGEA 600MG ID date Medicamento 1 13750598 2021-01-01 00:08:14 Não consta na lista 2 13797248 2021-01-01 00:20:39 *Piperacilina* 3 13797248 2021-01-01 00:20:39 Não consta na lista 4 13797248 2021-01-01 00:20:39 Não consta na lista 5 13727853 2021-01-01 00:22:18 Não consta na lista 6 13727853 2021-01-01 00:22:18 Não consta na lista 583867 13922033 2021-01-31 23:53:48 Não consta na lista 583868 13922033 2021-01-31 23:53:48 Não consta na lista 583869 13923030 2021-01-31 23:58:50 Não consta na lista 583870 13923030 2021-01-31 23:58:50 Não consta na lista 583871 13923030 2021-01-31 23:58:50 *Penicilina* 583872 13877978 2021-01-31 23:59:47 Outro medicamento

3 3