- R-br - DINF UFPR

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

05 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

Re: [R-br] Digest R-br, volume 104, assunto 2
by Maurício Sangiogo 08 Ago '19

08 Ago '19

Fernando Souza, quando não uso uma tabela com as medianas, uso um gráfico Box Plot, eu acho a apresentação desse gráfico interessante para esse tipo de dados. Att, _______________________________________________________ Eng. Agr. Maurício Sangiogo Doutorando em Fitossanidade - UFPel Analista Ambiental - FEPAM ________________________________ 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: quinta-feira, 8 de agosto de 2019 12:00 Para: r-br(a)listas.c3sl.ufpr.br <r-br(a)listas.c3sl.ufpr.br> Assunto: Digest R-br, volume 104, assunto 1 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. Contrastes - MLG (Andre Oliveira) 2. Graficos para não parametrica (Fernando Souza) 3. Re: Contrastes - MLG (Cesar Rabak) ---------------------------------------------------------------------- Message: 1 Date: Wed, 7 Aug 2019 19:47:49 +0000 (UTC) From: Andre Oliveira <andreolsouza(a)yahoo.com.br> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Contrastes - MLG Message-ID: <642002604.1904317.1565207269786(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Boa tarde! Tem como inverter a ordem do contraste que vem por padrão na função glht do pacote multcomp? Digo inverter para bater com o padrão do resultado do summary(fit) onde coloca t1-t2 e não t2 - t1 como fica no objeto p. dados <- read.table("http://www.leg.ufpr.br/~paulojus/dados/exemplo01.txt", head = T) attach(dados) boxplot(resp~trat, col=8) hist(resp, col=8) fit=glm(resp~trat, family = Gamma) summary(fit) require(multcomp) p <- glht(fit, linfct = mcp(trat = "Tukey")) plot(p, main= "", xlab = "") obg att.André -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190807/e5508270/atta…> ------------------------------ Message: 2 Date: Wed, 7 Aug 2019 17:20:06 -0300 From: Fernando Souza <nandodesouza(a)gmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Graficos para não parametrica Message-ID: <CAFrWFsm3a-MUoh1+rRRK4SnQikyfjt=H+BXfqNYbckSn5Tmq-A(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Pessoal Gostari de uma informação. Qual o melhor gráfico para utilizar com dados não paramétricos? Fiz uma análise utiizando a função art() do pacote ARTols do R o qual me permitiu testar interações não parametrica. Porém para mostrar os dados na forma gráfica e ilustrar uma interação significativa, há alguma forma? Desde já agradeço -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ========================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190807/655083be/atta…> ------------------------------ Message: 3 Date: Wed, 7 Aug 2019 18:57:00 -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] Contrastes - MLG Message-ID: <CAKrF98=k=uudnQ8aJECUSoj+j2-HGECsGSH=Sxh2LM3u3esWUA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" André, A função glht está usando a ordem "alfabética" dos seus fatores, como a função não deve ter meios de mudar isso (você já *olhou* no FM, não *olhou*?) ?, a única forma seria fazer uma cópia do objeto retornado 'p' neste teu exemplo, e inverter os sinais de *todos* os contrastes e colocar os *labels* manualmente no gráfico. HTH -- Cesar Rabak On Wed, Aug 7, 2019 at 4:50 PM Andre Oliveira por (R-br) < r-br(a)listas.c3sl.ufpr.br> wrote: > Boa tarde! > > Tem como inverter a ordem do contraste que vem por padrão na função glht > do pacote multcomp? Digo inverter para bater com o padrão do resultado do > summary(fit) onde coloca t1-t2 e não t2 - t1 como fica no objeto p. > > dados <- read.table("http://www.leg.ufpr.br/~paulojus/dados/exemplo01.txt > <https://l.facebook.com/l.php?u=http%3A%2F%2Fwww.leg.ufpr.br%2F%7Epaulojus%2…>", > head = T) > > attach(dados) > > boxplot(resp~trat, col=8) > > hist(resp, col=8) > > fit=glm(resp~trat, family = Gamma) > > summary(fit) > > require(multcomp) > > p <- glht(fit, linfct = mcp(trat = "Tukey")) > plot(p, main= "", xlab = "") > > obg > > att. > André > _______________________________________________ > 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/20190807/54b9507a/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 1 *****************************************

1 0

Contrastes - MLG
by Andre Oliveira 07 Ago '19

07 Ago '19

Boa tarde! Tem como inverter a ordem do contraste que vem por padrão na função glht do pacote multcomp? Digo inverter para bater com o padrão do resultado do summary(fit) onde coloca t1-t2 e não t2 - t1 como fica no objeto p. dados <- read.table("http://www.leg.ufpr.br/~paulojus/dados/exemplo01.txt", head = T) attach(dados) boxplot(resp~trat, col=8) hist(resp, col=8) fit=glm(resp~trat, family = Gamma) summary(fit) require(multcomp) p <- glht(fit, linfct = mcp(trat = "Tukey")) plot(p, main= "", xlab = "") obg att.André

2 1

Graficos para não parametrica
by Fernando Souza 07 Ago '19

07 Ago '19

Pessoal Gostari de uma informação. Qual o melhor gráfico para utilizar com dados não paramétricos? Fiz uma análise utiizando a função art() do pacote ARTols do R o qual me permitiu testar interações não parametrica. Porém para mostrar os dados na forma gráfica e ilustrar uma interação significativa, há alguma forma? Desde já agradeço -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

1 0