- R-br - DINF UFPR

Re: [R-br] Digest R-br, volume 110, assunto 2
by juraci sampaio 04 Fev '20

04 Fev '20

Bom dia como posso aumentar a fonte dos números do correlograma? Usando a função: library(corrplot) library(Boruta) library(car) library(PerformanceAnalytics) library(caret) library(Rcmdr) library(corrgram) #Definindo os genótipos como nome dos casos: row.names(dados) <- as.character(dados$genotipos) dados$genotipos <- NULL #Gerando o correlograma dados<-as.matrix(dados) chart.Correlation(dados, histogram=TRUE, pch=19, main= "Moderately Resistante") # plot correlogram (com histograma) chart.Correlation(dados, histogram=FALSE, pch=19) # plot correlogram (sem histograma) At.te Juraci Souza Sampaio Filho Eng Agrônomo-UFRB Msc Ciências Agrárias -UFRB Doutorando em Ciências Agrárias- UFRB Tel: (75) 981722995 ________________________________ De: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> em nome de r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Enviado: terça-feira, 4 de fevereiro de 2020 12:00 Para: r-br(a)listas.c3sl.ufpr.br <r-br(a)listas.c3sl.ufpr.br> Assunto: Digest R-br, volume 110, assunto 2 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: Criar numeração única para dados iguais (Cesar Rabak) ---------------------------------------------------------------------- Message: 1 Date: Mon, 3 Feb 2020 17:29:49 -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] Criar numeração única para dados iguais Message-ID: <CAKrF98nwF0RPNK72uZXUCU96yyL5ZF4JZKWfALbg0qhFaiid_g(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Admitindo que você gerou a classificação na variável grupo por meios adequados, ou seja *você já disponha *dessa coluna, minha sugestão é que você crie a variável GC pela conversão em inteiro do grupo : Código não testado: seu_dataframe$GC <- as.integer(seu_dataframe$grupo) Dependendo como a foi feita a geração da variável *grupo* talvez para isso funcionar você tenha que primeiro converter para *factor*: seu_dataframe$GC <- as.integer(as.factor(seu_dataframe$grupo)) HTH -- Cesar Rabak On Mon, Feb 3, 2020 at 9:32 AM Bárbara Mazetti Nascimento por (R-br) < r-br(a)listas.c3sl.ufpr.br> wrote: > Boa tarde pessoal! > > Sou nova no R e estou me batendo um pouco com algumas edições de banco de > dados que era acostumada a fazer no SAS. > > Preciso criar uma coluna com uma nova identificação para uma variável que > eu criei agrupando outras variáveis. Eu trabalho com dados de bovinos e > criei o que chamamos de grupo de contemporâneos: animais que vivenciaram o > mesmo ambiente. Então geralmente são animais nascidos no mesmo ano, no > mesmo mês e na mesma propriedade. Então, por exemplo, temos os seguintes > animais: > > Animal Ano de nascimento Mês de nascimento Fazenda > A 2018 Oct > 1 > B 2018 Oct > 1 > C 2018 Oct > 1 > D 2018 Oct > 2 > E 2018 Oct > 2 > F 2018 Oct > 2 > G 2018 Nov > 1 > H 2018 Nov > 1 > I 2018 Nov > 1 > J 2018 Nov > 1 > > Eu agrupei essas variáveis em uma que eu chamei de Grupo: > > Animal Grupo > A 2018 Oct 1 > B 2018 Oct 1 > C 2018 Oct 1 > D 2018 Oct 2 > E 2018 Oct 2 > F 2018 Oct 2 > G 2018 Nov 1 > H 2018 Nov 1 > I 2018 Nov 1 > J 2018 Nov 1 > > Agora preciso criar essa nova numeração, para que os animais que tenham o > mesmo Grupo sejam chamados pelo mesmo número (GC). A saída que eu espero é > a seguinte: > > Animal Grupo GC > A 2018 Oct 1 1 > B 2018 Oct 1 1 > C 2018 Oct 1 1 > D 2018 Oct 2 2 > E 2018 Oct 2 2 > F 2018 Oct 2 2 > G 2018 Nov 1 3 > H 2018 Nov 1 3 > I 2018 Nov 1 3 > J 2018 Nov 1 3 > > Será que vocês poderiam me ajudar em como eu consigo criar essa nova > numeração? > > Muito obrigada!! > Bárbara > > -- > _____________________________________________________ > MSc. Bárbara Mazetti Nascimento > PhD student in Animal Science > Genetic Applied to Animal Breeding Group - GAMA > UFPR - Curitiba - PR - Brazil > +55 41 99142-9495 > > _______________________________________________ > 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/20200203/88cb6e34/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 110, assunto 2 *****************************************

1 0

Explicação Valores GLM
by Chiara Lubich 22 Jan '20

22 Jan '20

Boa tarde, Pessoal, que explicações eu poderia dar dos valores gerados por um modelo GLM, além de dizer que variáveis deram significativas? Por exemplo, abaixo: [image: image.png] Consigo ver que as variáveis "d", "n" e "intercepto" deram diferentes, mas o que isso quer dizer? Obrigada desde já Abraço

4 3

Re: [R-br] Digest R-br, volume 109, assunto 2
by Chiara Lubich 20 Jan '20

20 Jan '20

Oi, pessoal Muito obrigada pela ajuda Daniel e Cesar, os dois modos deram certo. Abraço Em seg., 20 de jan. de 2020 às 10: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. JUNÇÃO DE GRÁFICOS DE LINHA (Chiara Lubich) > 2. Re: JUNÇÃO DE GRÁFICOS DE LINHA (Daniel Guimarães Tiezzi) > 3. Re: JUNÇÃO DE GRÁFICOS DE LINHA (Cesar Rabak) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 19 Jan 2020 15:22:16 -0400 > From: Chiara Lubich <lubichchiara(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] JUNÇÃO DE GRÁFICOS DE LINHA > Message-ID: > < > CAH-ZE+C7RhyWZPGTu1jP7YRdhwZnfjxxtjV9RHF+323E+e7Xcg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Oi, gente > Estou com dificuldades para realizar a junção de gráficos de linha em > apenas um. > Minha tabela está da seguinte forma: > taxadecoleta quantidade especie > 0.000 81.873 m > 0.050 77.880 m > 0.100 74.082 m > 0.150 70.469 m > 0.200 67.032 m > 0.000 67.032 mo > 0.050 63.763 mo > 0.100 60.653 mo > 0.150 57.695 mo > 0.200 54.881 mo > 0.250 52.205 mo > 0.300 49.659 mo > 0.350 47.237 mo > 0.400 44.933 mo > 0.000 54.881 mor > 0.050 52.205 mor > 0.100 49.659 mor > 0.150 47.237 mor > Porém quando vou gerar o gráfico no r com esse comando abaixo: > > plot(taxadecoleta, quantidade, col=especie, ylim=c(0,100), t="l", > xlim=c(0,4), xlab = "Coleta", ylab = "Quantidade") > > O gráfico não fica da forma que eu gostaria, ele fica da seguinte forma: > *Ele junta todas as linhas, formando uma só! :(* > [image: image.png] > O máximo que consegui foi esse gráfico de pontos abaixo, porém eu precisava > que fosse em linha em diferente cores como o da figura abaixo. > [image: image.png] > > Agradeço desde já > Abraço > > Att., > Chiara Lubich > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/d8c5e046/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9073 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/d8c5e046/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9585 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/d8c5e046/atta… > > > > ------------------------------ > > Message: 2 > Date: Sun, 19 Jan 2020 16:46:42 -0300 > From: Daniel Guimarães Tiezzi <dtiezzi(a)usp.br> > To: a lista Brasileira oficial de discussÃ£o do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] JUNÇÃO DE GRÁFICOS DE LINHA > Message-ID: > <CA+RL3DPvO_+hd_ys3w= > B48JN7Du-xEeKiJgLJEwNm9+L1zA8Ug(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Pode usar a função lines() > > tx <- read.delim('tx.txt') > > plot(tx$taxadecoleta[tx$especie == 'm'], tx$quantidade[tx$especie == > 'm'], col='red', ylim=c(0,100), xlim=c(0,0.4), xlab = "Coleta", ylab > = "Quantidade", t='l')lines(tx$taxadecoleta[tx$especie == 'mo'], > tx$quantidade[tx$especie == 'mo'], > col='blue')lines(tx$taxadecoleta[tx$especie == 'mor'], > tx$quantidade[tx$especie == 'mor'], col='green') > > > Em dom., 19 de jan. de 2020 às 16:22, Chiara Lubich por (R-br) < > r-br(a)listas.c3sl.ufpr.br> escreveu: > > > Oi, gente > > Estou com dificuldades para realizar a junção de gráficos de linha em > > apenas um. > > Minha tabela está da seguinte forma: > > taxadecoleta quantidade especie > > 0.000 81.873 m > > 0.050 77.880 m > > 0.100 74.082 m > > 0.150 70.469 m > > 0.200 67.032 m > > 0.000 67.032 mo > > 0.050 63.763 mo > > 0.100 60.653 mo > > 0.150 57.695 mo > > 0.200 54.881 mo > > 0.250 52.205 mo > > 0.300 49.659 mo > > 0.350 47.237 mo > > 0.400 44.933 mo > > 0.000 54.881 mor > > 0.050 52.205 mor > > 0.100 49.659 mor > > 0.150 47.237 mor > > Porém quando vou gerar o gráfico no r com esse comando abaixo: > > > plot(taxadecoleta, quantidade, col=especie, ylim=c(0,100), t="l", > > xlim=c(0,4), xlab = "Coleta", ylab = "Quantidade") > > > > O gráfico não fica da forma que eu gostaria, ele fica da seguinte forma: > > *Ele junta todas as linhas, formando uma só! :(* > > [image: image.png] > > O máximo que consegui foi esse gráfico de pontos abaixo, porém eu > > precisava que fosse em linha em diferente cores como o da figura abaixo. > > [image: image.png] > > > > Agradeço desde já > > Abraço > > > > Att., > > Chiara Lubich > > _______________________________________________ > > 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. > > > > > -- > Daniel Tiezzi, MD, PhD > Breast Disease and Gynecologic Oncology Division > Department of Gynecology and Obstetrics > University of São Paulo > Brazil > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/031b1756/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9073 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/031b1756/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9585 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/031b1756/atta… > > > > ------------------------------ > > Message: 3 > Date: Sun, 19 Jan 2020 21:51:06 -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] JUNÇÃO DE GRÁFICOS DE LINHA > Message-ID: > < > CAKrF98knsthNV1wA-bPcuCSw5aYNW6wijmniepwg6rghu8toAw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Uma forma "mais automática" de fazer usando o pacote car, supondo seus > dados no dataframe dat. > > > library(car) > > sp(quantidade ~ taxadecoleta, dat, groups=dat$especie) > > HTH > -- > Cesar Rabak > > On Sun, Jan 19, 2020 at 4:22 PM Chiara Lubich por (R-br) < > r-br(a)listas.c3sl.ufpr.br> wrote: > > > Oi, gente > > Estou com dificuldades para realizar a junção de gráficos de linha em > > apenas um. > > Minha tabela está da seguinte forma: > > taxadecoleta quantidade especie > > 0.000 81.873 m > > 0.050 77.880 m > > 0.100 74.082 m > > 0.150 70.469 m > > 0.200 67.032 m > > 0.000 67.032 mo > > 0.050 63.763 mo > > 0.100 60.653 mo > > 0.150 57.695 mo > > 0.200 54.881 mo > > 0.250 52.205 mo > > 0.300 49.659 mo > > 0.350 47.237 mo > > 0.400 44.933 mo > > 0.000 54.881 mor > > 0.050 52.205 mor > > 0.100 49.659 mor > > 0.150 47.237 mor > > Porém quando vou gerar o gráfico no r com esse comando abaixo: > > > plot(taxadecoleta, quantidade, col=especie, ylim=c(0,100), t="l", > > xlim=c(0,4), xlab = "Coleta", ylab = "Quantidade") > > > > O gráfico não fica da forma que eu gostaria, ele fica da seguinte forma: > > *Ele junta todas as linhas, formando uma só! :(* > > [image: image.png] > > O máximo que consegui foi esse gráfico de pontos abaixo, porém eu > > precisava que fosse em linha em diferente cores como o da figura abaixo. > > [image: image.png] > > > > Agradeço desde já > > Abraço > > > > Att., > > Chiara Lubich > > _______________________________________________ > > 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/20200119/8a7fce1d/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9073 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/8a7fce1d/atta… > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9585 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200119/8a7fce1d/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 109, assunto 2 > ***************************************** > -- *Chiara Lubich Cardoso Furtado* Mestranda em Biologia de Água Doce e Pesca Interior - INPA Engenheira de Pesca - UFAM Telefone: (92) 99110-6657 e-mail: lubichchiara(a)gmail.com Lattes: http://lattes.cnpq.br/0673892253719514

1 0

JUNÇÃO DE GRÁFICOS DE LINHA
by Chiara Lubich 19 Jan '20

19 Jan '20

Oi, gente Estou com dificuldades para realizar a junção de gráficos de linha em apenas um. Minha tabela está da seguinte forma: taxadecoleta quantidade especie 0.000 81.873 m 0.050 77.880 m 0.100 74.082 m 0.150 70.469 m 0.200 67.032 m 0.000 67.032 mo 0.050 63.763 mo 0.100 60.653 mo 0.150 57.695 mo 0.200 54.881 mo 0.250 52.205 mo 0.300 49.659 mo 0.350 47.237 mo 0.400 44.933 mo 0.000 54.881 mor 0.050 52.205 mor 0.100 49.659 mor 0.150 47.237 mor Porém quando vou gerar o gráfico no r com esse comando abaixo: > plot(taxadecoleta, quantidade, col=especie, ylim=c(0,100), t="l", xlim=c(0,4), xlab = "Coleta", ylab = "Quantidade") O gráfico não fica da forma que eu gostaria, ele fica da seguinte forma: *Ele junta todas as linhas, formando uma só! :(* [image: image.png] O máximo que consegui foi esse gráfico de pontos abaixo, porém eu precisava que fosse em linha em diferente cores como o da figura abaixo. [image: image.png] Agradeço desde já Abraço Att., Chiara Lubich

3 2

componente da variância no modelo de efeitos aleatórios
by Luiz Leal 11 Jan '20

11 Jan '20

Prezados, boa tarde. No modelo de efeitos aleatórios estimamos a componente da variância conforme abaixo. y <-c(2370, 1687, 2592, 2283, 2910, 3020, 1282, 1527, 871, 1025, 825, 920, 562, 321, 636, 317, 485, 842, 173, 127, 132, 150, 129, 227, 193, 71, 82, 62, 96, 44)x <- as.factor(c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "D", "E", "E", "E", "E", "E", "E"))summary(aov(y~x)) Mas não sei como calcular quando utilizo mínimos quadrados generalizados. y <-c(2370, 1687, 2592, 2283, 2910, 3020, 1282, 1527, 871, 1025, 825, 920, 562, 321, 636, 317, 485, 842, 173, 127, 132, 150, 129, 227, 193, 71, 82, 62, 96, 44) x <- as.factor(c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "D", "E", "E", "E", "E", "E", "E"))library(nlme) modelo <- gls(y~x, weights = varExp() , method = "ML") Alguém sabe como calculo as componentes da variância neste caso? Desde já agradeço Luiz

2 4

Tinn-R project: new version (6.01.01.03) released
by Jose Claudio Faria 23 Dez '19

23 Dez '19

Hi, A new version of Tinn-R project (6.01.01.03) was released today. Download: - https://nbcgib.uesc.br/tinnr/en/download - https://sourceforge.net/projects/tinn-r/files/Tinn-R setup/ What is new: - https://nbcgib.uesc.br/tinnr/en/download#patch - https://sourceforge.net/projects/tinn-r/files/Tinn-R setup/6.1.1.3/ Best, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)99966.9100 - VIVO ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ If you have software to deal with statistics, you have arms; if you have good software, you have arms and legs; if you have software like R, you have arms, legs and wings... the height of your flight depends only on you!

1 0

Erro predição - Função predict
by Andre Oliveira 16 Dez '19

16 Dez '19

Boa noite a todos! Estou procurando obter a métrica RMSE e R2 entre os objetos pred e diameter no modelo abaixo, mas está ocorrendo um erro. Parece que reclama de missing, mas ainda não consegui resolver, mesmo depois de muitas pesquisas. O interesse é pegar o predito levando em consideração os dois efeitos do modelo. Alguém teria alguma dica e sugestão? # CMDR rm(list=ls(all=TRUE))require(lme4)data(Penicillin, package="lme4")attach(Penicillin) library(MCMCglmm)model <- MCMCglmm(diameter ~ 1, random = ~ plate + sample, verbose= T, nitt=200000,burn=10000,thin=10, rcov = ~ units, data=Penicillin) Penicillin$pred=predict(model, data=Penicillin, marginal = NULL) Error in object$Sol[, missing.random, FALSE] <- rnorm(nrow(object$Sol) * : incorrect number of subscripts Pedir os valores preditos como abaixo resolve, mas, demora muito tempo em caso de bases de dados grandes. Penicillin$pred=predict(model, data=Penicillin, marginal = NULL, interval="prediction") att,.André

1 0

Manter a estrutura do objeto da funcao randomForest
by Nilson Guiotoku 10 Dez '19

10 Dez '19

Pessoal estou com uma dúvida quanto o pacote randomForest do R. Usei a função randomForest para modelar o meu conjunto de dados. Como irei aplicar o modelo já treinado para meu conjunto teste, salvei como .rda. Na próxima etapa iria subir o treino (.rda) para rodar o banco de dados independente. Porém, a forma como salvei está dando erro quando aplico a função predict. É como se o load ou a extensão .rda não conservasse o formato do banco de dados. E o erro é que o objeto "model" deveria ser um randomForest e ele fica como "character". Eu sei que eu poderia usar o que foi salvo na memoria do R, o modelCT, esse se usar no predict não daria erro. Mas para cada novo conjunto de dados independente, não gostaria de ficar treinando a todo o momento, gostaria de passar somente os parâmetros já treinados. Como um código executável para exemplificar o problema usei o da Iris. library(randomForest) modelCT <- randomForest::randomForest(Sepal.Length ~ ., data = iris, importance = TRUE) # modelagem com randomForest str(modelCT) save(modelCT, file = "model.rda") # salvei como rda model = load("model.rda") # Faço a chamada do arquivo treinado .rda str(model) predValid <- randomForest::predict(model, iris, type = "class") # o banco de dados deveria ser independente, mas aqui é somente para mostrar o erro do formato como está o model predValid <- predict(model, iris, type = "class")

2 1

Comparar elementos em 2 dataframes
by sérgio marques 06 Dez '19

06 Dez '19

Boa tarde, Alguma dica para iniciante no R, para se obter o seguinte. Tenho 2 tabelas (dataframes) com várias linhas e várias colunas, em cada célula dessas tabelas está um valor (cada valor é único, mas pode estar repetido ou não na tabela 2), quero obter uma 3.ª tabela (dataframe) onde conste cada valor da tabela 1 que não esteja presente na Tabela 2. Obrigado! Sérgio ________________________________________________________________________________

5 5

Re: [R-br] Validação cruzada
by sznelwar＠uol.com.br 05 Dez '19

05 Dez '19

4 4

Tamanho do efeito para Wilcoxon
by Elias Carvalho 05 Dez '19

05 Dez '19

Boa tarde Tenho duas amostras independentes com distribuição nao-normal e usei o método abaixo para ver a diferença entre elas: coin::wilcoxsign_test(data.amostra1 ~ data.amostra2, distribution="exact") Alguem pode me ajudar a calcular o tamanho do efeito? -- *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"*

3 2

Ajuda com download de csv pela internet
by Olympio Neto 04 Dez '19

04 Dez '19

Pessoal, bom dia. Estou com o seguinte problema. Preciso fazer automaticamente um download pela internet de um arquivo em CSV no site do TSE: http://www.tse.jus.br/eleitor/biometria/biometria-atual-do-dia, mas o link é gerado a partir de um comando em javascrit, algo parecido com *javascript:window.location.href=apex.server.url({p_request: 'FLOW_EXCEL_OUTPUT_R69876138832319545_pt-br'},2); null;* Alguém já teve esse problema? Anteriormente o site não limitava o número de registros e conseguia fazer o web scraping via html , mas a partir desta semana o site foi atualizado e está limitando a 30 por página. Alguma dica? Muito Obrigado Olympio Neto

2 3

Vinculando R com OpenBLAS em sistemas GNU/Linux
by Pedro Rafael 03 Dez '19

03 Dez '19

O *ropenblas* (<https://prdm0.github.io/ropenblas/>) é um pacote projetado para facilitar a vinculação da biblioteca OpenBLAS (< https://www.openblas.net/>) com a linguagem R. O pacote, que funciona apenas para sistemas Linux, baixará automaticamente o código-fonte mais recente da biblioteca OpenBLAS e o compilará. O pacote ligará automaticamente o idioma R para usar a biblioteca OpenBLAS. Tudo será feito automaticamente, independentemente da distribuição do Linux que você estiver usando. Você também pode especificar versões mais antigas da biblioteca OpenBLAS. Automaticamente, se nenhuma versão for especificada, o pacote ropenblas considerará a versão mais recente da biblioteca OpenBLAS. A consideração de usar a biblioteca OpenBLAS em vez do BLAS (< http://www.netlib.org/blas/>) pode trazer otimizações extras para seu código e melhor desempenho computacional para suas simulações, pois o OpenBLAS é uma implementação otimizada da biblioteca BLAS. Algumas das razões pelas quais é conveniente vincular a linguagem R ao uso de alternativas otimizadas para BLAS podem ser encontradas aqui (< https://csantill.github.io/RPerformanceWBLAS/>). Vários outros benchmarks que apontam para um melhor desempenho computacional considerando a biblioteca OpenBLAS podem ser encontrados na internet. Algumas vantagens de usar a biblioteca *ropenblas*: [1] Tudo é feito dentro da linguagem R; [2] O procedimento será o mesmo para qualquer distribuição Linux; [3] A biblioteca OpenBLAS será compilada e você escolherá qual versão de compilação vincular ao R, independentemente da sua distribuição Linux; [4] Se a sua distribuição GNU / Linux não possui versões atualizadas do OpenBLAS, pouco importa. O pacote ropenblas busca a versão estável mais recente da conta de desenvolvimento da biblioteca OpenBLAS no GitHub; [5] Você não precisa conhecer bem o Linux. Em algumas distribuições, pode não ser tão simples para um usuário menos experiente compilar e vincular a biblioteca à biblioteca OpenBLAS com a linguagem R; [6] É muito mais fácil direcionar uma pessoa para vincular o OpenBLAS ao R dizendo "execute *ropenblas* () no R" do que pedir a essa pessoa que verifique se uma versão não otimizada do BLAS está instalada no sistema. Então você deve orientar a remoção da versão não otimizada do BLAS e guiá-lo para a instalação da biblioteca OpenBLAS através dos mais diversos procedimentos, dependendo da distribuição GNU / Linux usada; [7] Como afirmado anteriormente, o procedimento funciona para qualquer Linux e isso inclui o Android. Se o seu Android for capaz de executar comandos privilegiados (ROOT) e se você tiver o R instalado via Termux com as dependências necessárias, poderá compilar e vincular o OpenBLAS ao R usando *ropenblas*.

1 0

Validação cruzada
by Andre Oliveira 02 Dez '19

02 Dez '19

boa noite, estou com dificuldades fazer a validação cruzada de modelos gerados com as bibliotecas sommer e MCMCglmm. Alguém do grupo que tenha experiência poderia ma dar uma ajuda? Segue o CMR # Dados require(MCMCglmm)data(PlodiaPO) str(PlodiaPO) # Divisão dos dados library(caTools)divisao = sample.split(PlodiaPO$PO, SplitRatio = 0.75)conj_treinamento = subset(PlodiaPO,divisao==TRUE)conj_validacao = subset(PlodiaPO$PO,divisao==FALSE) # Modelo require(sommer)fit=mmer(PO ~ 1, random = ~ FSfamily, data=conj_treinamento)summary(fit) # Validação ???pred=predict(fit, new_data = conj_validacao, classify = "FSfamily")pred att,.André

2 1

Re: [R-br] Desdobramento de um fatorial triplo
by Walmes Zeviani 21 Nov '19

21 Nov '19

Respostas dentro da mensagem. À disposição. Walmes. On Mon, Nov 18, 2019 at 7:19 PM Maurício Lordêlo <mslordelo(a)gmail.com> wrote: > Obrigado Walmes! > Tenho feito bastante esta análise exploratória por achar fundamental > conhecer os dados. A questão que como dou muito suporte para trabalhos de > mestrado e doutorado, eles cobram os testes de normalidade e homogeneidade > de variâncias. > É uma cultura a ser mudada. Hábitos que de uma época com pouco recursos de visualização e do uso precedural da estatística (faça o texte X, se rejeitar H0, faça Y, caso contrário faça Z, etc). Também tenho contato com público assim, mas tô sempre argumentando, redijo a leitura dos gráficos, não reporto CV nas minhas análises (apenas quando o revisor não aceita as justificativas dadas para omissão). > Confesso que tenho muita dificuldade em explicar algumas funções e > resultados do R para os pesquisadores da área de Agrárias e Biológicas. > Muitas vezes é uma tarefa árdua pelos vícios que adquirem. > É uma cultura. Mudá-la leva tempo. > Tenho feito o uso do pacote ExpDes, devido ao fato de muitos já terem > utilizado o SisVar. > Rodei esta transformação proposta [(m_seca_raiz + 0.00001)^0.5] usando a > função a fat3dic deste pacote. > Aproveitando a oportunidade, como o coeficiente de variação deu bem alto > neste caso, qual a sua sugestão para convencer estes pesquisadores que esta > medida não é "a medida"? > Eles "idolatram" este coeficiente...rs rs rs. > Primeiro ponto é argumentar que não existe vínculo entre CV alto/baixo e não atendimento ou não dos pressupostos. Pode-se ter CV baixo com sérios afastamentos. Pode-se ter CV alto com total conformidade com as suposições. As pessoas criaram a regra no sentido inverso. Por exemplo, sempre viam que dados de contagem tem alto CV e geralmente não atendem os pressupostos. Aí entenderam que CV alto implica em falta de pressupostos, o que não é fundamentalmente verdade. Sempre que se usava uma transformação estabilizadora da variância, o CV era grande, mas porque a distribuição assumida para aplicações das funções estabilizadoras mais comuns é Poisson (contagem) ou proporção (Binomial), que não atendem os pressupostos. Portanto, a prática sem muita reflexão ou esclarecimento criou uma regra ao contrário, a de que CV alto é sinal de problema. O segundo ponto são as faixas para classificação do CV em baixo, médio e alto. Elas são arbitrárias. Aí os pesquisadores querem comparar um CV de 8% para produtividade de grãos em experimentos em vários locais com um CV de 38% de crescimento radicular em cultivo in vitro avaliando o efeito de doses homeopáticas de hormônio. Não tem como comparar isso. Uma coisa é feita em macro escala, grandes parcelas, variável resposta controlada por muitos genes e condições ambientais. O outro é micro escala. Um é produtividade de grãos (kg/ha) o outro é comprimento (mm de raíz). Como que uma estatística adquiriu tamanha e desproporcional importância? Se você pensar o que o CV está medindo, fica claro que ele não deveria ser usado em uma análise de experimentos. É o quociente entre desvio padrão residual e média amostral do experimento (CV 100 · DP/M). Mas se o experimento é feito com a premissa de existência de diferença entre médias (efeito dos tratamentos), que valor tem uma estatística que usa uma média global? Em alguns poucos contextos haverá justificativa (i.e. efeitos aleatórios). Outro ponto é que pressupõe uma relação 1:1 entre média e desvio padrão, justamente se contrapondo a suposição de ausência de relação entre média e variância. Existem muitos pontos frágeis sobre o CV. Eu poderia prosseguir aqui por muitas linhas. Mas eu geralmente me atenho ao primeiro: não existe relação entre valor do CV com atendimento dos pressupostos (a menos se for conhecida a distribuição dos dados). O CV estando alto ou não, se os pressupostos forem atendidos, a parte inferencial está assegurada (distribuições das estatísticas de teste conforme esperado, cobertura dos intervalos de confiança conforme esperado, níveis de significância conforme esperado, etc). Mas se ainda o sujeito quiser olhar pro CV alto e dizer que o experimento foi mal conduzido (é a coisa mais comum de ouvir), o que é lamentável, veja se consegue dialogar e desfazer a cultura aos poucos. > > Outro pedido é que gostaria de uma função diferente destas (do pacote > ExpDes) para obter estes resultados do desdobramento. Você tem alguma > sugestão cujas sintaxe e resultado > são mais ou menos fáceis de serem entendidas por este público? > Não tenho opções que sejam mais fáceis que o respostado pelo ExpDes ou SISVAR. Tenho apenas alternativas mais gerais para casos não gaussianos, balanceados, etc. O que posso recomendar para fazer o desdobramento no caso no glm(), survreg() e outros modelos paramétricos ou delineamentos mais completos é `emmeans` e `multcomp`. Veja alguns exemplos a seguir. GLM: http://leg.ufpr.br/~walmes/analises/WAALima/caiuae/caiuae.html Survreg: http://leg.ufpr.br/~walmes/analises/CCastellar/frutosmaduros.html LM: http://leg.ufpr.br/~walmes/analises/PSFLichtemberg/Colletotrichum/phy-gf.ht… Ancova: http://leg.ufpr.br/~walmes/mpaer/analise-de-covariancia.html > > Agradeço mais uma vez caso possa ajudar. > > Abraço, > > Maurício >

2 2

Erro no pacote RWeka
by Nilson Guiotoku 19 Nov '19

19 Nov '19

Pessoal, boa tarde, não estou conseguindo instalar o pacote RWeka. O que devo fazer? >* install.packages("RWeka")* WARNING: Rtools is required to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding: https://cran.rstudio.com/bin/windows/Rtools/ Installing package into ‘C:/Users/usuario/Documents/R/win-library/3.6’ (as ‘lib’ is unspecified) trying URL ' https://cran.rstudio.com/bin/windows/contrib/3.6/RWeka_0.4-41.zip' Content type 'application/zip' length 632650 bytes (617 KB) downloaded 617 KB package ‘RWeka’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\usuario\AppData\Local\Temp\Rtmp0w0FBo\downloaded_packages Então peguei a versão recomentada do Rtool35.exe ( https://cran.rstudio.com/bin/windows/Rtools/) e instalei dentro da pasta (C:/Users/usuario/Documents/R/win-library/3.6) do meu computador. E depois instalei o RWeka, usando o zip ( https://cran.rstudio.com/bin/windows/contrib/3.6/RWeka_0.4-41.zip) e também usando o comando como mencionei acima. As pastas Rtools e Rweka aparecem dentro de ( C:/Users/usuario/Documents/R/win-library/3.6) como se já tivesse instalado. Porém, quando chamo a biblioteca RWeka, aparece o seguinte erro. > *library(RWeka)* Error: package or namespace load failed for ‘RWeka’: .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry

1 0

Desdobramento de um fatorial triplo
by Maurício Lordêlo 18 Nov '19

18 Nov '19

Caros, Peço auxilio para um problema de um fatorial em que há um problema em atender aos pressupostos. Agradeço alguma sugestão. Maurício m_seca_raiz = c(0.000000000, 0.000000000, 0.000460000, 0.000000000, 0.000000000, 0.000140000, 0.000000000, 0.000300000, 0.000150000, 0.000000000, 0.000120000, 0.000500000, 0.000560000, 0.000300000, 0.000000000, 0.000160000, 0.000160000, 0.000600000, 0.000200000, 0.000200000, 0.000440000, 0.000000000, 0.000080000, 0.000240000, 0.000380000, 0.000080000, 0.000000000, 0.000250000, 0.000960000, 0.000800000, 0.000320000, 0.000760000, 0.000920000, 0.000366667, 0.000720000, 0.000540000, 0.000440000, 0.000620000, 0.000100000, 0.000600000, 0.000200000, 0.000280000, 0.000000000, 0.000520000, 0.000480000, 0.000000000, 0.000280000, 0.000280000, 0.000800000, 0.000280000, 0.000000000, 0.000680000, 0.000700000, 0.000380000, 0.000000000, 0.000360000, 0.000300000, 0.000340000, 0.000280000, 0.000360000, 0.000340000, 0.001780000, 0.000250000, 0.000260000, 0.000300000, 0.000280000, 0.000160000, 0.000120000, 0.000125000, 0.000080000, 0.000220000, 0.000320000) meio = factor((rep(c("m2","m3"),each=12, times=3))) sacarose = factor(rep(c("7,5","15,0","30,0"), each=24)) carvão = factor((rep(c("c1","c2"),each=6, times=6))) tratamento = meio:sacarose:carvão head(data.frame(m_seca_raiz, meio, sacarose, carvão, tratamento),18) hist(m_seca_raiz) boxplot(m_seca_raiz ~ tratamento,ylab="Massa seca da raiz", xlab="Tratamento", cex.axis=0.85 ) max(m_seca_raiz) modelo_m_seca_raiz = aov(m_seca_raiz ~ tratamento) summary(modelo_m_seca_raiz) hist(modelo_m_seca_raiz$residuals) shapiro.test(modelo_m_seca_raiz$residuals) require(car) leveneTest(modelo_m_seca_raiz$residuals ~ tratamento) anova(modelo_m_seca_raiz)$'Mean Sq'[2] res_pad_m_seca_raiz = modelo_m_seca_raiz$residuals/sqrt(anova(modelo_m_seca_raiz)$'Mean Sq'[2]) plot(modelo_m_seca_raiz$fitted.values, res_pad_m_seca_raiz, ylim = c(-5,5)) abline(h=c(-3,3),lty=3) #Elimando a observação 62 shapiro.test(modelo_m_seca_raiz$residuals[-62]) leveneTest(modelo_m_seca_raiz$residuals[-62] ~ tratamento[-62]) #Transformação Box-Cox require(forecast) #require(fpp) BoxCox.lambda(m_seca_raiz+0.5, method = "guerrero", lower=-3,upper = 3) BoxCox.lambda(m_seca_raiz+0.5, method = "loglik", lower=-3,upper = 3) #lambda = BoxCox.lambda(m_seca_raiz,lower=-3) #lambda m_seca_raizt = BoxCox(m_seca_raiz+0.5, lambda=-3) hist(m_seca_raizt) boxplot(m_seca_raizt ~ tratamento, ylab="m_seca_raiz (transformada)", xlab="Tratamentos") modelo_m_seca_raizt = aov(m_seca_raizt ~ tratamento) summary(modelo_m_seca_raizt) hist(modelo_m_seca_raizt$residuals) plot(modelo_m_seca_raizt$residuals) plot(modelo_m_seca_raizt$fitted.values, modelo_m_seca_raizt$residuals) shapiro.test(modelo_m_seca_raizt$residuals) leveneTest(modelo_m_seca_raizt$residuals ~tratamento) anova(modelo_m_seca_raizt)$'Mean Sq'[2] res_pad_m_seca_raizt = modelo_m_seca_raizt$residuals/sqrt(anova(modelo_m_seca_raizt)$'Mean Sq'[2]) plot(modelo_m_seca_raizt$fitted.values, res_pad_m_seca_raizt, ylim = c(-6,6)) abline(h=c(-3,3),lty=3)

2 1

Applied Linear Statistical Models
by Cleber Iack 13 Nov '19

13 Nov '19

Boa noite, eu estou escrevendo um artigo utilizando o R. Mas a minha formulação teórica para embasar a aplicação dos dados está faltando uma informação que necessito. Algum dos senhores teriam o Apêndice A do livro "Applied Linear Statistical Models", O pdf que eu tenho está com o Apêndice incompleto. Nao apresenta o inicio da Apendice A. Passa da pagina 1289 para a pagina 1306. Obrigado.. Ceber Iack

2 1

Organizar data.frame e sumarizar dados em transposição
by Marcelo Laia 13 Nov '19

13 Nov '19

Caro(a), Aqui na Universidade nós temos que sumarizar e organizar os dados dos discentes bolsistas por mês e valor recebido de bolsa. Os dados estão assim: CPF Nome Descrição Data Valor xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/02/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/03/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/04/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/07/2018 400,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Manutenção 01/08/2018 400,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Emergencial 01/06/2018 200,00 Nós precisamos da seguinte saída: CPF Nome Quantidade Parcelas Valor Mês Valor Ano xxx.xxx.xxx-xx Fulano de Tal 6 266,67 1600,00 Se puder ajudar, agradeço! -- Marcelo

6 10

Dados climáticos
by Ze Henrique 13 Nov '19

13 Nov '19

Bom doa Srs. alguém poderia me ajudar onde consigo achar a série histórica de dados climáticos do Brasil para alguns municípios entre 2000 e 2018, para uso no R. verifiquei que tem o pacote install "inmetr". Mas não to conseguindo baixar ele. Desde já agradeço. -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

2 1

Última Versão do R em ambiente Linux
by Eduardo Bearzoti 11 Nov '19

11 Nov '19

Prezados colegas: Utilizo o sistema operacional linux (Linux Mint 18.3), e gostaria (precisaria) muito de utilizar o R mais novo, versão 3.6. Embora use o Mint 18.3, tentei seguir as instruções deste link sugeridas para instalar o R no Mint 19, mas não deu certo. Após seguir as instruções, ao rodar: sudo apt-get install r-baseobtive: Os pacotes a seguir têm dependências desencontradas: r-base : Depende: r-base-core (>= 3.6.1-3bionic) mas não será instalado Depende: r-recommended (= 3.6.1-3bionic) mas não será instalado E: Impossível corrigir problemas, você manteve (hold) pacotes quebrados. E ao tentar rodar, por exemplo:sudo apt-get install r-base-coreobtive: Os pacotes a seguir têm dependências desencontradas: r-base-core : Depende: libc6 (>= 2.27) mas 2.23-0ubuntu11 está para ser instalado Depende: libcurl4 (>= 7.28.0) mas não é instalável Depende: libicu60 (>= 60.1-1~) mas não é instalável Depende: libreadline7 (>= 6.0) mas não é instalável E: Impossível corrigir problemas, você manteve (hold) pacotes quebrados. Parece-me assim que instalar a versão mais nova do R não é trivial, pois dependeria de bibliotecas que não estão disponíveis em minha versão do Mint. Como acredito que haja colegas aqui que também têm o Linux como sistema operacional, gostaria de consultá-los, se têm conseguido utilizar a versão mais nova do R, ou se isto realmente não é possível. Qualquer dica/informação será muito bem vinda. Abraços a todos,Eduardo. | | | | | | | | | | | How to Install R on Linux Mint 19 | Linux4one Table of ContentsPrerequisitesInstalling R on Linux MintRun Simple Hello World Program In RConclusion Install R ... | | | ******************************************** Prof. Eduardo Bearzoti Universidade Federal de Ouro Preto Departamento de Estatística - ICEB Campus Morro do Cruzeiro, s/n 35400-000 - Ouro Preto, MG, Brasil

3 2

densidade de uma variável
by Nilson Guiotoku 11 Nov '19

11 Nov '19

Pessoal, Estou com uma dúvida sobre a aplicação da função density, em um banco de dados de uma distribuição normal. Após fazer uma transformação encontrando o score (componente principal), ao aplicar a função density, eu esperava que os valores deveriam variar de 0 a 1. Plotei até o histograma da densidade. Por que a densidade não está no intervalo 0 a 1? Gerei um código reproduzível, porém, minha variável que estou aplicando, não necessariamente é exatamente uma variável normal. set.seed(23) X = cbind( 1:230, rnorm(230, 450, 20)) X[,2] <- (X[,2] - min(X[,2]))/(max(X[,2])-min(X[,2])) #plot(X) X[,1] <- X[,1]-mean(X[,1]) X[,2] <- X[,2]-mean(X[,2]) C <- cov(X) E <- eigen(C) score <- X %*% E$vectors densidade = density(score[,2], n=100) # os valores da densidade não estão dentro do intervalo [0,1] a = hist(score[,2], freq = F) lines(densidade)

3 2

Re: [R-br] Digest R-br, volume 107, assunto 7
by Vanessa Gomes 07 Nov '19

07 Nov '19

ah ótimo! Já deu pra entender. Muito obrigada. *Vanessa Matos Gomes* *Laboratório de Ecologia Evolutiva e Biodiversidade* Departamento de Ecologia, Genética e Evolução Instituto de Ciências Biológicas UFMG - Brasil *Ecosystem Restoration and Intervention Ecology * School of Biological Science University of Western Australia Em qui., 7 de nov. de 2019 às 22: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: Dúvida de 'warning message' (Cesar Rabak) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 6 Nov 2019 15:00:53 -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] Dúvida de 'warning message' > Message-ID: > <CAKrF98nhAu8hnzfXocqKyY5Jgnrvfz= > UG0hK96vnoMLaBfFx9A(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Você entende o mecanismo de criação de "postos" (*ranks* em inglês)? > > Quando dois ou mais deles têm a mesma posição devido ao cálculo como é > feito, você tem o *empate*. > > Caso você precise aprofundar, um livro texto de Estatística tem o material > para esclarecer esse assunto. > > HTH > -- > Cesar Rabak > > > On Wed, Nov 6, 2019 at 3:08 AM Vanessa Gomes por (R-br) < > r-br(a)listas.c3sl.ufpr.br> wrote: > > > Olá pessoal, > > > > após rodar correlação de Spearman no R, recebi a seguinte warning > message: In > > cor.test.default(EEG, DP, method = "spearman") : Impossível calcular o > > valor exato de p com empates > > > > Alguém sabe me dizer o que isso significa? > > > > Obrigada. > > > > *O que rodei foi:* > > > > > plot(EEG,DP) > > > cor(EEG,DP,method = "spearman") > > [1] 0.3688426 > > > cor.test(EEG,DP,method = "spearman") > > > > Spearman's rank correlation rho > > > > data: EEG and DP > > S = 6728.1, p-value = 0.01918 > > alternative hypothesis: true rho is not equal to 0 > > sample estimates: > > rho > > 0.3688426 > > > > Warning message: > > In cor.test.default(EEG, DP, method = "spearman") : > > Impossível calcular o valor exato de p com empates > > > > > > *Vanessa Matos Gomes* > > > > *Laboratório de Ecologia Evolutiva e Biodiversidade* > > Departamento de Ecologia, Genética e Evolução > > Instituto de Ciências Biológicas > > UFMG - Brasil > > *Ecosystem Restoration and Intervention Ecology * > > School of Biological Science > > University of Western Australia > > _______________________________________________ > > 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/20191106/808f721a/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 107, assunto 7 > ***************************************** >

1 0

Dúvida de 'warning message'
by Vanessa Gomes 06 Nov '19

06 Nov '19

Olá pessoal, após rodar correlação de Spearman no R, recebi a seguinte warning message: In cor.test.default(EEG, DP, method = "spearman") : Impossível calcular o valor exato de p com empates Alguém sabe me dizer o que isso significa? Obrigada. *O que rodei foi:* > plot(EEG,DP) > cor(EEG,DP,method = "spearman") [1] 0.3688426 > cor.test(EEG,DP,method = "spearman") Spearman's rank correlation rho data: EEG and DP S = 6728.1, p-value = 0.01918 alternative hypothesis: true rho is not equal to 0 sample estimates: rho 0.3688426 Warning message: In cor.test.default(EEG, DP, method = "spearman") : Impossível calcular o valor exato de p com empates *Vanessa Matos Gomes* *Laboratório de Ecologia Evolutiva e Biodiversidade* Departamento de Ecologia, Genética e Evolução Instituto de Ciências Biológicas UFMG - Brasil *Ecosystem Restoration and Intervention Ecology * School of Biological Science University of Western Australia

2 1

(sem assunto)
by Jersiton Tiago Pereira Matos 04 Nov '19

04 Nov '19

Bom dia, prezados. Tenho uma dúvida, é possível remover valores repetidos em uma coluna dado um critério de uma outra coluna do mesmo dataframe? att. Jersiton Mattos

2 3

Design of Experiments: replicatas no ponto central e intervalo de confiança
by Nei 03 Nov '19

03 Nov '19

Prezados e Prezadas, Estou estudando Planejamento Fatorial Completo de 2 níveis. Fiz um planejamento completo para 4 fatores em 2 níveis e para isso fiz 16 experimentos (2^4), sem replicatas. Agora pretendo fazer 3 replicatas no ponto central para estimar o desvio padrão e a partir daí calcular o intervalo de confiança e avaliar a significância dos parâmetros do modelo. Mas estou com dificuldades para encontrar um tutorial explicando como calcular o erro padrão e como definir o número de graus de liberdade para o valor de t para poder calcular o intervalo de confiança. Alguém saberia me dar uma dica? Obrigado, Nei

2 2

derivada da densidade de probabilidade
by Nilson Guiotoku 02 Nov '19

02 Nov '19

Pessoal, eu tenho uma variável que segue uma distribuição bimodal, e eu preciso encontrar o melhor ponto que separa essas bimodais, de forma lógica matematicamente. Então eu pensei em gerar a curva de densidade e encontrar a derivada dessa curva. O local em que a derivada for zero eu teria o valor critico e também o ponto de mínimo, portanto, o melhor valor para determinar o corte e encontrar os dois grupos distintos. Como encontrar a derivada primeira e segunda da densidade? # Código mínimo executável: X = c(rnorm(100, 10, 3), rnorm(50, 40, 5)) # Variável bimodal hist(X, freq = F) denfuncao = density(X) lines(denfuncao, lwd = 2) # como encontrar a derivada de uma curva se não conheço a função de # densidade? deriv(expr = denfuncao$x, namevec = denfuncao$x)

2 1

Analise conjunta R
by Jose Wellingthon Dos Santos 01 Nov '19

01 Nov '19

Prezado prof Walmes, bom dia. Em uma análise conjunta que fiz, a interação Cultivar x Local foi significativa. Precisa desdobrá-la e aplicar o teste Skott-Knott nas medias de cultivar dentro de cada Local. Como fazê-lo? Att -- José Wellingthon dos Santos Pesquisador - Métodos Quantitativos Embrapa Algodão Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Campina Grande, PB jose-wellingthon.santos(a)embrapa.br Telefone: +55 (83) 3182-4401 | Fax: +55 (83) 3182-4367 www.embrapa.br/algodao | https://twitter.com/embrapa_algodao Confira também: www.facebook.com/agrosustentavel ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

2 2

como exportar um raster categorico a partir do R para o SAGAGIS?
by Waldir De Carvalho Junior 31 Out '19

31 Out '19

Boa tarde a Todos estou tentando exportar um raster categorico elaborado por um modelo multinom para um SIG Mas os comandos sempre mostram erro por não ser um raster de números e sim de classes (classes de solos) Tanto o write.raster como write.asciigrid ou writeGDAL só aceitam se o raster for numérico... Se alguem puder me dar uma solução.... __________ Dr. Waldir de Carvalho Junior Engº Agroº D.Sc. Ciência do Solo - Mapeamento Digital Embrapa Solos Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Rio de Janeiro/RJ Waldir.Carvalho(a)embrapa.br Telefone: +55 (21) 2179-4614 | Fax: +55 (21) 2274-5291 www.embrapa.br/solos | | fb.com/embrapa | twitter.com/embrapa Call Send SMS Add to Skype You'll need Skype Credit Free via Skype ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

2 1

[OFF-TOPIC] 2º R Day · Curitiba/PR · 13/Dez/2019
by Walmes Zeviani 28 Out '19

28 Out '19

Saudações listeiRxs, Venho por meio deste divulgar que teremos o 2º R Day dia 13/Dez aqui na UFPR, Curitiba/PR. O evento visa congregar usuários, desenvolvedores, apoiadores e entusiastas para discutir o uso do R em ambiente acadêmico e profissional, além de meios para fortalecer e crescer a comunidade no Brasil e criar rede de colaboração. Detalhes sobre o evento estão em http://rday.leg.ufpr.br/. A inscrição é GRATUITA mas as vagas são LIMITADAS. Esperamos contar com a participação de vocês dessa lista que é um importante marco na história da comunidade R no Brasil. Sintam-se livres e estimulados para repassarem essa divulgação dos círculos dos quais participam. Pela atenção e certo da vossa solidariedade, agradeço. Walmes.

1 0

Como identificar colinearidade?
by Chiara Lubich 18 Out '19

18 Out '19

Boa tarde, Pessoal fiz uma análise de GLM, porém no resultado do modelo ficava aparecendo a seguinte mensagem: Coefficients: (2 not defined because of singularities). Buscando na internet, vi que se trata de colinearidade dos dados, porém não sei como identificar. Na mesma página de internet que encontrei a sugestão de identificar utilizando o comando pairs(), para tentar detectar colinearidade. Mas não sei como fazer a interpretação do gráfico gerado: Segue o gráfico com a variável Y1= CM: CM= Comprimento médio (contínua) D = distância do rio (contínua) R = rios (categórica - nome dos rios) TP = tempo de pesca (Dummy - 0= rios considerados tradicionais e 1= rios com pesca recente) NP= número de pescador por rio [image: image.png] Se puderem me ajuda na explicação de como identifico a colinearidade a partir do gráfico, agradeço. Ou até mesmo outra forma de encontrar... Abraço

2 1

Correlação de Pearson com escala Likert
by Angélica Ricarte 17 Out '19

17 Out '19

Prezados, boa noite. Por gentileza, poriam me ajudar? Preciso fazer uma análise de correlação de Pearson para resultados de uma pesquisa de satisfação para a qual foi utilizado um questionário e realizadas modificações em escala Likert, em que: [cid:image002.png@01D5846B.3C568D50] Alguém sabe me informar qual o pacote do R trabalha com tal análise e como fazê-la? Utilizo o “agricolae” para fazer correlação, mas não havia feito nada antes com escala Likert. Toda informação sobre esta escala me será útil. Grata, aguardo retorno. Abraços, Angélica.

2 1

Re: [R-br] RES: Digest R-br, volume 106, assunto 6
by sznelwar＠uol.com.br 10 Out '19

10 Out '19

2 1

RES: Digest R-br, volume 106, assunto 6
by Max 10 Out '19

10 Out '19

Marcelo Laia Preparei o script e tentei disponibilizar neste link: https://gist.github.com/Verginio/d671ff36b852c68b6ba7b7f59114a22e Pelo que entendi, daria para realizar esta tarefa com a função “group_by” seguidas de “summarise”. Em que a primeira define quais variáveis serão agrupadas e a segunda executa o cálculo para os grupos. library(dplyr) df_summ <- df_bolsa %>% group_by(CPF, Nome) %>% summarise( `Quantidade Parcelas` = n(), `Valor Mês` = mean(Valor), #média `Valor Ano` = sum(Valor) #somatório ) #Resultado head(df_summ) Abraços, Max Verginio De: r-br-request(a)listas.c3sl.ufpr.br Enviado:quinta-feira, 10 de outubro de 2019 12:00 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Digest R-br, volume 106, assunto 6 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. Organizar data.frame e sumarizar dados em transposição (Marcelo Laia) 2. Modelos lineares generalizados multivariados -literatura (Jose Wellingthon Dos Santos) 3. Re: Organizar data.frame e sumarizar dados em transposição (Cesar Rabak) 4. Re: Organizar data.frame e sumarizar dados em transposição (Marcelo Laia) 5. Re: Modelos lineares generalizados multivariados -literatura (Walmes Zeviani) ---------------------------------------------------------------------- Message: 1 Date: Wed, 9 Oct 2019 12:28:26 -0300 From: Marcelo Laia <marcelolaia(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Organizar data.frame e sumarizar dados em transposição Message-ID: <20191009152826.GC25795@localhost> Content-Type: text/plain; charset=iso-8859-1 Caro(a), Aqui na Universidade nós temos que sumarizar e organizar os dados dos discentes bolsistas por mês e valor recebido de bolsa. Os dados estão assim: CPF Nome Descrição Data Valor xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/02/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/03/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/04/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/07/2018 400,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Manutenção 01/08/2018 400,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Emergencial 01/06/2018 200,00 Nós precisamos da seguinte saída: CPF Nome Quantidade Parcelas Valor Mês Valor Ano xxx.xxx.xxx-xx Fulano de Tal 6 266,67 1600,00 Se puder ajudar, agradeço! -- Marcelo ------------------------------ Message: 2 Date: Wed, 9 Oct 2019 15:46:40 -0300 (BRT) From: Jose Wellingthon Dos Santos <jose-wellingthon.santos(a)embrapa.br> To: a lista Brasileira oficial de discussÃ£o do programa R. <r-br(a)listas.c3sl.ufpr.br> Cc: "walmes " <walmes(a)ufpr.br>, "walmeszeviani " <walmeszeviani(a)gmail.com> Subject: [R-br] Modelos lineares generalizados multivariados -literatura Message-ID: <1889187931.11565119.1570646800406.JavaMail.root(a)embrapa.br> Content-Type: text/plain; charset="iso-8859-1" Oi pessoal, Alguém conhece literatura para modelos lineares generalizados multivariados? Att -- José Wellingthon dos Santos Pesquisador - Métodos Quantitativos Embrapa Algodão Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Campina Grande, PB jose-wellingthon.santos(a)embrapa.br Telefone: +55 (83) 3182-4401 | Fax: +55 (83) 3182-4367 www.embrapa.br/algodao | https://twitter.com/embrapa_algodao Confira também: www.facebook.com/agrosustentavel ----- Mensagem original ----- De: "Marcelo Laia por (R-br)" <r-br(a)listas.c3sl.ufpr.br> Para: "a lista Brasileira oficial de discussão do programa R." <r-br(a)listas.c3sl.ufpr.br> Cc: "Marcelo Laia" <marcelolaia(a)gmail.com> Enviadas: Quarta-feira, 9 de outubro de 2019 12:28:26 Assunto: [R-br] Organizar data.frame e sumarizar dados em transposição Caro(a), Aqui na Universidade nós temos que sumarizar e organizar os dados dos discentes bolsistas por mês e valor recebido de bolsa. Os dados estão assim: CPF Nome Descrição Data Valor xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/02/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/03/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/04/2018 200,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/07/2018 400,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Manutenção 01/08/2018 400,00 xxx.xxx.xxx-xx Fulano de Tal Bolsa Emergencial 01/06/2018 200,00 Nós precisamos da seguinte saída: CPF Nome Quantidade Parcelas Valor Mês Valor Ano xxx.xxx.xxx-xx Fulano de Tal 6 266,67 1600,00 Se puder ajudar, agradeço! -- Marcelo _______________________________________________ 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. ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20191009/38dbdf38/atta…> ------------------------------ Message: 3 Date: Wed, 9 Oct 2019 16:32:38 -0300 From: Cesar Rabak <cesar.rabak(a)gmail.com> To: R-br <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Organizar data.frame e sumarizar dados em transposição Message-ID: <CAKrF98mdt6jM3Sne8+iiHxDMWEs0_HUmKehmg23Lq+ha-McyEA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Sem um CMR é difícil ser assertivo, mas o delineamento do que você deseja seria obter uma lista sem repetições da coluna dos CPF (*uniq*), e usando os itens dessa lista como parâmetro para *which* ou *subset* fazer a soma e contar o número de parcelas adicionando no novo dataframe esses resultados. HTH -- Cesar Rabak On Wed, Oct 9, 2019 at 12:28 PM Marcelo Laia por (R-br) < r-br(a)listas.c3sl.ufpr.br> wrote: > Caro(a), > > Aqui na Universidade nós temos que sumarizar e organizar os dados dos > discentes > bolsistas por mês e valor recebido de bolsa. > > Os dados estão assim: > > CPF Nome Descrição Data Valor > xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/02/2018 200,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/03/2018 200,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/04/2018 200,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/07/2018 400,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Manutenção 01/08/2018 400,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Emergencial 01/06/2018 200,00 > > > Nós precisamos da seguinte saída: > > CPF Nome Quantidade Parcelas Valor Mês Valor > Ano > xxx.xxx.xxx-xx Fulano de Tal 6 266,67 1600,00 > > Se puder ajudar, agradeço! > > -- > Marcelo > _______________________________________________ > 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/20191009/fa39d328/atta…> ------------------------------ Message: 4 Date: Wed, 9 Oct 2019 16:53:19 -0300 From: Marcelo Laia <marcelolaia(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] Organizar data.frame e sumarizar dados em transposição Message-ID: <20191009195319.GE20280@localhost> Content-Type: text/plain; charset=iso-8859-1 Então, Cesar. Eu tentei fazer um CMR mas não tive conhecimento para tal. Eu tentei com a família apply e com os pacotes dplyr, tidyr e plyr [1]. Não dei conta de fazer. Se puder ajudar, terá que dar mastigado. Isso não é um experimento. Precisamos fazer um relatório aqui na Universidade sobre bolsas dadas a discentes vulneráveis e preciso tratar os dados para inclui-los numa plataforma fixa do governo. Pra eu calcular isso na unha (calc) vai ser muito demorado e temos um prazo. Mas, eu sei que a lista não é pra esse tipo de coisa (dar mastigado). Abração Obrigado 1. http://www.rpubs.com/Kassio_Ferreira/tutorialdplyr Marcelo On 09/10/19 at 04:32, Cesar Rabak por (R-br) wrote: > Sem um CMR é difícil ser assertivo, mas o delineamento do que você > deseja seria obter uma lista sem repetições da coluna dos CPF (uniq), e > usando os itens dessa lista como parâmetro para which ou subset fazer a > soma e contar o número de parcelas adicionando no novo dataframe esses > resultados. > HTH > -- > Cesar Rabak > > On Wed, Oct 9, 2019 at 12:28 PM Marcelo Laia por (R-br) > <[1]r-br(a)listas.c3sl.ufpr.br> wrote: > > Caro(a), > Aqui na Universidade nós temos que sumarizar e organizar os dados > dos discentes > bolsistas por mês e valor recebido de bolsa. > Os dados estão assim: > CPF Nome Descrição Data Valor > xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/02/2018 200,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Permanência 01/03/2018 200,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/04/2018 200,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Integração 01/07/2018 400,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Manutenção 01/08/2018 400,00 > xxx.xxx.xxx-xx Fulano de Tal Bolsa Emergencial 01/06/2018 200,00 > Nós precisamos da seguinte saída: > CPF Nome Quantidade Parcelas Valor Mês > Valor Ano > xxx.xxx.xxx-xx Fulano de Tal 6 266,67 > 1600,00 > Se puder ajudar, agradeço! > -- > Marcelo > _______________________________________________ > R-br mailing list > [2]R-br(a)listas.c3sl.ufpr.br > [3]https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem ([4]http://www.leg.ufpr.br/r-br-guia) e > forneça código mínimo reproduzível. > > Referências > > 1. mailto:r-br@listas.c3sl.ufpr.br > 2. mailto:R-br@listas.c3sl.ufpr.br > 3. https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > 4. http://www.leg.ufpr.br/r-br-guia > _______________________________________________ > 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. -- Marcelo ------------------------------ Message: 5 Date: Wed, 9 Oct 2019 18:05:14 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: Jose Wellingthon Dos Santos <jose-wellingthon.santos(a)embrapa.br> Cc: a lista Brasileira oficial de discussÃ£o do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Modelos lineares generalizados multivariados -literatura Message-ID: <CAFU=EkbNPd5EA=_RRCuVSBOiS4RF4ACSTCiFxvcc0W=q+b+YKg(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" José Wellingthon O prof. Wagner e eu oferecemos um Curso da 62ª RBras que foi em Lavras sobre análise multivariada de dados experimentais http://cursos.leg.ufpr.br/mcglm4aed/ Dê uma olhada no material. Você tem uma necessidade específica de orientação para uma análise ou procura material geral? À disposição. Walmes. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20191009/546c75b5/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 106, assunto 6 *****************************************

1 0

MCMCglmm - Parameter Expansion
by Andre Oliveira 03 Out '19

03 Out '19

bom dia! Alguém da lista trabalhou com esta biblioteca? O que seria Parameter Expansion? Sei que é uma reparametrização de priori, mas, não consegui enxergar a distribuição de prob. que se prõe ao usar a Parameter Expansion. att,.André

1 0

Re: [R-br] Ajuda com análise conjunta de dois experimentos com parcelas subdividas
by sznelwar＠uol.com.br 02 Out '19

02 Out '19

1 0

fluxograma no Rmarkdown
by Luiz Leal 02 Out '19

02 Out '19

Prezados, bom dia.Alguém ja utilizou o DiagrammeR com o Rmarkdown?Há um espaço excessivo entre o fluxograma e o texto. Alguém pode me ajudar? (em anexo está o arquivo Rmd)Desde já agradeçoLuiz

2 2

Ajuda com análise conjunta de dois experimentos com parcelas subdividas
by José Lucas Safanelli 02 Out '19

02 Out '19

Caros colegas, Peço ajuda a vocês para conseguir entender e analisar o seguinte experimento: - são dois experimentos feitos lado a lado: um com plantio direto, outro com convencional. - em cada local, existe um experimento de parcelas subdividades, sendo o fator manejo na parcela e dose de nitrogenio na subparcela. Diante disso, gostaria de saber como expressar o modelo da ANOVA no R, a fim de fazer comparação de médias, e mais pra frente, análise multivariada via discriminante canonica. Sei que são necessárias duas coisas: - Dividir o erro dentro da abordagem da análise conjunta. - Dividir o erro na parcela subdividida. Tenho utilizado a seguinte expressão, mas não estou confiante quanto a expressão. Infelizmente nunca trabalhei com experimento mais complexos como este. aov_N <- aov(N_kg_ha ~ Manejo * Rotacao * Nitrogenio + Bloco + Error(Nitrogenio/Bloco/Manejo), data = dados) summary(aov_N) Tanto o croqui quanto os dados estão neste link: https://drive.google.com/drive/folders/1rjVv9KwtE5A7y-VPXTV9Y2GlilKLXK2F?us… Agradeço se alguém puder me ajudar com alguma orientação. José Lucas.

3 4

Ajuda: Escolha de distribuição de probabilidade (GLM)
by Chiara Lubich 30 Set '19

30 Set '19

Oi, Estou tendo dificuldade para achar uma distribuição certa para minha análise (GLM)... Mesmo dados não são normais e também não apresentam homoscedasticidade, mesmo após realizar transformação (log). Qual família devo escolher? Estava trabalhando com binomial negativa, é a correta? Obrigada e aguardo retorno -- *Chiara Lubich Cardoso Furtado* Mestranda em Biologia de Água Doce e Pesca Interior - INPA Engenheira de Pesca - UFAM Telefone: (92) 99110-6657 e-mail: lubichchiara(a)gmail.com Lattes: http://lattes.cnpq.br/0673892253719514

2 1

dic parcela subdividida no tempo
by ANDREA VITA REIS MENDONCA 20 Set '19

20 Set '19

Trata-se de um experimento deduzas temperaturas e seis períodos de avaliação, em DIC. Resolvi conforme mostrado a seguir. Gostaria de saber se está correto. As analises ficaram boas. pkg= c("lattice", "latticeExtra", "doBy", "multcomp", "reshape", "plyr", "nlme") sapply(pkg, require, character.only=TRUE) source("http://dl.dropboxusercontent.com/u/48140237/bandas.R") sessionInfo() trellis.device(color=FALSE) ang=read.csv("ang.csv") str(ang) dh=melt(subset(ang, select=1:8), id.vars=1:2) d=diff(as.Date(c("2017-01-01","2017-01-16", "2017-01-31","2017-02-15", "2017-03-02","2017-03-17"))); d d =cumsum(c(15, d)); d dh$dia=d[as.integer(dh$variable)] levels(dh$variable) xyplot(value~dia|Temperatura, groups=Temperatura, data=dh, type=c("p","a"), auto.key=TRUE) dh1=transform(dh, D=factor(dia), parc=interaction(Rep, Temperatura), subp=interaction(Rep, Temperatura)) str(dh1) m2=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="ML") m3=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="REML") -- Andrea Vita Reis Mendonça Profa. Associado I UFRB/CCAAB

2 1

Não carrega o biblioteca - lmer4
by Andre Oliveira 13 Set '19

13 Set '19

boa noite,de uma hora para outra não consigo carregar a biblioteca lme4.Alguém tem alguma sugestão de ajuda? (estou usando win10) require(lmer4)Carregando pacotes exigidos: lmer4Warning message:In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package called ‘lmer4’ att., André

5 5

gráficos com bandas de confiança no R - com sombra
by Nilson Guiotoku 12 Set '19

12 Set '19

Pessoal, gostaria de fazer um gráfico conforme o anexo, cujo o intervalo de confiança fique com formato de sombra dentro do gráfico. Conforme figura anexada: URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190904/bfde02f0/atta… > Como tentativa usei o *polygon*, mas não sei se existe um outro recurso melhor. Se alguém souber e puder contribuir, a lista do R agradece. Segue o código: x <- seq(1:30) Y <- x + rnorm(30, 0,3) model <- lm(Y ~ x) pred <- data.frame(x1 = x) pred$y1 <- predict(model, newdata = pred, interval = "confidence") plot(consumo, type = "n", main = "Consumo médio, kg", xlab = "Consumo estimado, kg ", ylab = "Consumo observado, kg" ) icfit = pred$y1 yy = c(as.numeric(as.vector(icfit[,2])), rev(as.numeric(as.vector(icfit[,3])))) xx = c(x, rev(x)) polygon(xx, yy, col = "lightsteelblue2", border = NA) points(consumo) abline(model, col="blue", lwd= 2) points(seq(2:15), seq(2:15), type = "l", col = "green", lwd = 2) legend("topright", legend = c("modelo linear base y = x", "modelo ajustado", "Bandas de confiança, 95%" ), lty = c(1,1,1), col = c("green", "blue", "lightsteelblue2"), lwd = c(2,2,NA), bty = "n", pch = c(NA,NA, 15))

4 3

Re: [R-br] gráficos com bandas de confiança no R - com sombra
by sznelwar＠uol.com.br 04 Set '19

04 Set '19

1 0

gráficos com bandas de confiança no R - com sombra
by Nilson Guiotoku 04 Set '19

04 Set '19

Pessoal, gostaria de fazer um gráfico conforme o anexo, cujo o intervalo de confiança fique com formato de sombra dentro do gráfico. Conforme figura anexada: [image: image.png] Como tentativa usei o *polygon*, mas não sei se existe um outro recurso melhor. Se alguém souber e puder contribuir, a lista do R agradece. Segue o código: x <- seq(1:30) Y <- x + rnorm(30, 0,3) model <- lm(Y ~ x) pred <- data.frame(x1 = x) pred$y1 <- predict(model, newdata = pred, interval = "confidence") plot(consumo, type = "n", main = "Consumo médio, kg", xlab = "Consumo estimado, kg ", ylab = "Consumo observado, kg" ) icfit = pred$y1 yy = c(as.numeric(as.vector(icfit[,2])), rev(as.numeric(as.vector(icfit[,3])))) xx = c(x, rev(x)) polygon(xx, yy, col = "lightsteelblue2", border = NA) points(consumo) abline(model, col="blue", lwd= 2) points(seq(2:15), seq(2:15), type = "l", col = "green", lwd = 2) legend("topright", legend = c("modelo linear base y = x", "modelo ajustado", "Bandas de confiança, 95%" ), lty = c(1,1,1), col = c("green", "blue", "lightsteelblue2"), lwd = c(2,2,NA), bty = "n", pch = c(NA,NA, 15))

1 0

Problema criar colunas: socorro!
by Elerson 02 Set '19

02 Set '19

Prezados: Trabalho um banco de dados com mais de 10000 linhas, representado pelo exemplo abaixo. *PROCESSO* *DATA_DISTRIBUICAO* *NOME_PARTE* *DATA_CRIME* *DATA_DENUNCIA* *ENQUADRAMENTO* 180022121 02/01/2018 A_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022121 02/01/2018 B_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.35 - 11343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.33 - 11343/06 180022097 02/01/2018 D_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022097 02/01/2018 E_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR.4º INC.II E IV - CPB 180022097 02/01/2018 F_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022501 03/01/2018 G_NONIMATO 02/01/2018 03/01/2018 ART.157 PAR.2 INC.I E II - CP 180022501 03/01/2018 G_NONIMATO 02/01/2018 03/01/2018 ART.14 - CP Hipóteses: a) mesmo PROCESSO, NOME_PARTE diferentes, ENQUADRAMENTO igual (linhas 1 e 2) b) mesmo PROCESSO, mesmo NOME_PARTE, ENQUADRAMENTO diferentes (linhas 3 e 4) c) mesmo PROCESSO, mesmo NOME_PARTE, ENQUADRAMENTO diferentes (linhas 8 e 9) Minha questão se refere à hipótese c: preciso criar um coluna (enquadramento2) para transferir o art. 14 transformando o caso em apenas uma linha… sem que se mexa na hipótese b. Em todo o banco sempre estará escrito ART.14 – CP. Para ficar assim: *PROCESSO* *DATA_DISTRIBUICAO* *NOME_PARTE* *DATA_CRIME* *DATA_DENUNCIA* *ENQUADRAMENTO* *ENQUADRAMENTO2* 180022121 02/01/2018 A_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022121 02/01/2018 B_NONIMATO 02/01/2018 02/01/2018 ART.33 PAR.CAPUT - 11.343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.35 - 11343/06 180022105 02/01/2018 C_NONIMATO 01/01/2018 02/01/2018 ART.33 - 11343/06 180022097 02/01/2018 D_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022097 02/01/2018 E_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR.4º INC.II E IV - CPB 180022097 02/01/2018 F_NONIMATO 01/01/2018 02/01/2018 ART.155 PAR. 4º INC. II E IV - CPB 180022501 03/01/2018 G_NONIMATO 02/01/2018 03/01/2018 ART.157 PAR.2 INC.I E II - CP ART.14 - CP Obrigado.

3 3

Estimativas de componentes de variância!
by Andre Oliveira 31 Ago '19

31 Ago '19

boa noite!Estou analisando as saídas destes dois modelos abaixo. Progenie=factor(Progenie)Row=factor(Row)Column=factor(Column) require(lmer4)fit1=lmer(Prod ~ (1|Progenie) + (1|Row) , data=dados)summary(fit1) Random effects: Groups Name Variance Std.Dev. Progenie (Intercept) 66.2 8.14 Row (Intercept) 47.6 6.90 Residual 100.1 10.00 Number of obs: 162, groups: Progenie, 81; Row, 19 require(sommer)fit=mmer(Prod ~ 1, random = ~ Progenie + Row, data=dados)summary(fit)Variance-Covariance components: VarComp VarCompSE Zratio ConstraintProgenie.Prod-Prod 62.11 20.60 3.016 PositiveRow.Prod-Prod 1.75 2.54 0.688 Positiveunits.Prod-Prod 110.69 17.49 6.328 Positive Alguém saberia dizer por qual motivo os componentes não são "iguais"? Seria algum erro na declaração do modelo ou seria devido o método de aproximação? att.André

1 0

Como reportar o valor do Odds Ratio
by Elias Carvalho 28 Ago '19

28 Ago '19

Fiz uma regressão logistica e a relação da variável "X" sobre o desfecho foi negativa: -0.53, porém o Odds ratio foi positivo: 0.559, por ser uma proporção entendo que OR deve ser sempre positivo. Então como se reporta isso (X tem uma influencia negativa de acordo com a RL e para cada unidade de alcool consumida há uma aumento dessa influencia de 0.559?) O que reportar na tabela de resultados? -- *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"*

4 3

Package Tibble
by Edmar Caldas 22 Ago '19

22 Ago '19

Boa noite! estou com um erro no código abaixo e não consigo descobrir o motivo. Já executei esse mesmo código em outra maquina. versao do R version.string R version 3.6.1 (2019-07-05) library(readxl)Salario <- read_excel("F:/FIAP/R/Excel/Salario.xlsx")View(Salario) attach(Salario)install.packages("dplyr")library(dplyr)agregar <-summarise(group_by(Salario, id), salario_m = mean(salario_m), data_pagamento_min=min(data_pagamento), data_pagamento_max=max(data_pagamento), contagem =n())agregar View(agregar) install.packages("rfm")install.packages("lubridate")library(rfm)library(lubridate) data_atual <- lubridate::as_date("2019-04-05", tz = "UTC-3")agregar$data_atual <-data_atualattach(agregar) data_compra <-lubridate::as_date(data_pagamento_max, tz="UTC")agregar$data_compra<-data_compra attach(agregar)View(agregar) rfm_result <- rfm_table_order(agregar, id, data_compra, salario_m, data_atual)rfm_result install.packages("tibble")library(tibble) rfm<-as.tibble(rfm_result)rfm erro Error: Tibble columns must have consistent lengths, only values of length one are recycled:* Length 0: Requested with `.rows` argument* Length 50: Column `rfm`Call `rlang::last_error()` to see a backtraceIn addition: Warning message:`as.tibble()` is deprecated, use `as_tibble()` (but mind the new semantics).This warning is displayed once per session. Edmar

3 2

Re: [R-br] Digest R-br, volume 104, assunto 3
by Edmar Caldas 20 Ago '19

20 Ago '19

Edmar Em terça-feira, 20 de agosto de 2019 12:00:28 BRT, r-br-request(a)listas.c3sl.ufpr.br <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. Package Tibble (Edmar Caldas) 2. Re: Package Tibble (Fernando Souza) ---------------------------------------------------------------------- Message: 1 Date: Mon, 19 Aug 2019 22:51:50 +0000 (UTC) From: Edmar Caldas <edmar_caldas(a)yahoo.com.br> To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Package Tibble Message-ID: <2090336582.86957.1566255110255(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Boa noite! estou com um erro no código abaixo e não consigo descobrir o motivo. Já executei esse mesmo código em outra maquina. versao do R version.string R version 3.6.1 (2019-07-05) library(readxl)Salario <- read_excel("F:/FIAP/R/Excel/Salario.xlsx")View(Salario) attach(Salario)install.packages("dplyr")library(dplyr)agregar <-summarise(group_by(Salario, id), salario_m = mean(salario_m), data_pagamento_min=min(data_pagamento), data_pagamento_max=max(data_pagamento), contagem =n())agregar View(agregar) install.packages("rfm")install.packages("lubridate")library(rfm)library(lubridate) data_atual <- lubridate::as_date("2019-04-05", tz = "UTC-3")agregar$data_atual <-data_atualattach(agregar) data_compra <-lubridate::as_date(data_pagamento_max, tz="UTC")agregar$data_compra<-data_compra attach(agregar)View(agregar) rfm_result <- rfm_table_order(agregar, id, data_compra, salario_m, data_atual)rfm_result install.packages("tibble")library(tibble) rfm<-as.tibble(rfm_result)rfm erro Error: Tibble columns must have consistent lengths, only values of length one are recycled:* Length 0: Requested with `.rows` argument* Length 50: Column `rfm`Call `rlang::last_error()` to see a backtraceIn addition: Warning message:`as.tibble()` is deprecated, use `as_tibble()` (but mind the new semantics).This warning is displayed once per session. Edmar -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190819/3d08833e/atta…> ------------------------------ Message: 2 Date: Mon, 19 Aug 2019 21:37:29 -0300 From: Fernando Souza <nandodesouza(a)gmail.com> To: "=?utf-8?Q?a_lista_Brasileira_oficial_de_discuss=C3=83=C2=A3o_do_programa_R.?=" <r-br(a)listas.c3sl.ufpr.br> Cc: Edmar Caldas <edmar_caldas(a)yahoo.com.br> Subject: Re: [R-br] Package Tibble Message-ID: <982FECB8-90E4-41AC-9120-62372264667C(a)getmailspring.com> Content-Type: text/plain; charset="utf-8" Seu exemplo não é reproduzível,, uma vez que não é possível acessar os dados. Att On Aug 19 2019, at 7:51 pm, Edmar Caldas por (R-br) <r-br(a)listas.c3sl.ufpr.br> wrote: > Boa noite! > > estou com um erro no código abaixo e não consigo descobrir o motivo. > Já executei esse mesmo código em outra maquina. > > versao do R > > version.string R version 3.6.1 (2019-07-05) > > > library(readxl) > Salario <- read_excel("F:/FIAP/R/Excel/Salario.xlsx") > View(Salario) > > attach(Salario) > install.packages("dplyr") > library(dplyr) > agregar <-summarise(group_by(Salario, id), > salario_m = mean(salario_m), > data_pagamento_min=min(data_pagamento), > data_pagamento_max=max(data_pagamento), > contagem =n()) > agregar > > View(agregar) > > > install.packages("rfm") > install.packages("lubridate") > library(rfm) > library(lubridate) > > data_atual <- lubridate::as_date("2019-04-05", tz = "UTC-3") > agregar$data_atual <-data_atual > attach(agregar) > > data_compra <-lubridate::as_date(data_pagamento_max, tz="UTC") > agregar$data_compra<-data_compra > > attach(agregar) > View(agregar) > > rfm_result <- rfm_table_order(agregar, id, data_compra, salario_m, data_atual) > rfm_result > > install.packages("tibble") > library(tibble) > > rfm<-as.tibble(rfm_result) > rfm > > > > erro > > Error: Tibble columns must have consistent lengths, only values of length one are recycled: > * Length 0: Requested with `.rows` argument > * Length 50: Column `rfm` > Call `rlang::last_error()` to see a backtrace > In addition: Warning message: > `as.tibble()` is deprecated, use `as_tibble()` (but mind the new semantics). > This warning is displayed once per session. > > > > Edmar > _______________________________________________ > 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/20190819/c2981dae/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 104, assunto 3 *****************************************

1 0

Re: [R-br] Digest R-br, volume 104, assunto 3
by Edmar Caldas 20 Ago '19

20 Ago '19

Edmar Em terça-feira, 20 de agosto de 2019 12:00:28 BRT, r-br-request(a)listas.c3sl.ufpr.br <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. Package Tibble (Edmar Caldas) 2. Re: Package Tibble (Fernando Souza) ---------------------------------------------------------------------- Message: 1 Date: Mon, 19 Aug 2019 22:51:50 +0000 (UTC) From: Edmar Caldas <edmar_caldas(a)yahoo.com.br> To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Package Tibble Message-ID: <2090336582.86957.1566255110255(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Boa noite! estou com um erro no código abaixo e não consigo descobrir o motivo. Já executei esse mesmo código em outra maquina. versao do R version.string R version 3.6.1 (2019-07-05) library(readxl)Salario <- read_excel("F:/FIAP/R/Excel/Salario.xlsx")View(Salario) attach(Salario)install.packages("dplyr")library(dplyr)agregar <-summarise(group_by(Salario, id), salario_m = mean(salario_m), data_pagamento_min=min(data_pagamento), data_pagamento_max=max(data_pagamento), contagem =n())agregar View(agregar) install.packages("rfm")install.packages("lubridate")library(rfm)library(lubridate) data_atual <- lubridate::as_date("2019-04-05", tz = "UTC-3")agregar$data_atual <-data_atualattach(agregar) data_compra <-lubridate::as_date(data_pagamento_max, tz="UTC")agregar$data_compra<-data_compra attach(agregar)View(agregar) rfm_result <- rfm_table_order(agregar, id, data_compra, salario_m, data_atual)rfm_result install.packages("tibble")library(tibble) rfm<-as.tibble(rfm_result)rfm erro Error: Tibble columns must have consistent lengths, only values of length one are recycled:* Length 0: Requested with `.rows` argument* Length 50: Column `rfm`Call `rlang::last_error()` to see a backtraceIn addition: Warning message:`as.tibble()` is deprecated, use `as_tibble()` (but mind the new semantics).This warning is displayed once per session. Edmar -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190819/3d08833e/atta…> ------------------------------ Message: 2 Date: Mon, 19 Aug 2019 21:37:29 -0300 From: Fernando Souza <nandodesouza(a)gmail.com> To: "=?utf-8?Q?a_lista_Brasileira_oficial_de_discuss=C3=83=C2=A3o_do_programa_R.?=" <r-br(a)listas.c3sl.ufpr.br> Cc: Edmar Caldas <edmar_caldas(a)yahoo.com.br> Subject: Re: [R-br] Package Tibble Message-ID: <982FECB8-90E4-41AC-9120-62372264667C(a)getmailspring.com> Content-Type: text/plain; charset="utf-8" Seu exemplo não é reproduzível,, uma vez que não é possível acessar os dados. Att On Aug 19 2019, at 7:51 pm, Edmar Caldas por (R-br) <r-br(a)listas.c3sl.ufpr.br> wrote: > Boa noite! > > estou com um erro no código abaixo e não consigo descobrir o motivo. > Já executei esse mesmo código em outra maquina. > > versao do R > > version.string R version 3.6.1 (2019-07-05) > > > library(readxl) > Salario <- read_excel("F:/FIAP/R/Excel/Salario.xlsx") > View(Salario) > > attach(Salario) > install.packages("dplyr") > library(dplyr) > agregar <-summarise(group_by(Salario, id), > salario_m = mean(salario_m), > data_pagamento_min=min(data_pagamento), > data_pagamento_max=max(data_pagamento), > contagem =n()) > agregar > > View(agregar) > > > install.packages("rfm") > install.packages("lubridate") > library(rfm) > library(lubridate) > > data_atual <- lubridate::as_date("2019-04-05", tz = "UTC-3") > agregar$data_atual <-data_atual > attach(agregar) > > data_compra <-lubridate::as_date(data_pagamento_max, tz="UTC") > agregar$data_compra<-data_compra > > attach(agregar) > View(agregar) > > rfm_result <- rfm_table_order(agregar, id, data_compra, salario_m, data_atual) > rfm_result > > install.packages("tibble") > library(tibble) > > rfm<-as.tibble(rfm_result) > rfm > > > > erro > > Error: Tibble columns must have consistent lengths, only values of length one are recycled: > * Length 0: Requested with `.rows` argument > * Length 50: Column `rfm` > Call `rlang::last_error()` to see a backtrace > In addition: Warning message: > `as.tibble()` is deprecated, use `as_tibble()` (but mind the new semantics). > This warning is displayed once per session. > > > > Edmar > _______________________________________________ > 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/20190819/c2981dae/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 104, assunto 3 *****************************************

1 0