Kruskal-Wallis para grupos com apenas 1 elemento
Prezados, preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso. Obrigada. -- Letícia Cavalari Pinheiro
Manda um CMR, assim facilita para os listeiros. CMR: Comando Mínimo Reproduzível [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro <leticia.ufmg@gmail.com> escreveu: Prezados, preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso. Obrigada. -- Letícia Cavalari Pinheiro _______________________________________________ R-br mailing list R-br@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.
O código que estou usando é o seguinte: dados<-data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12)) kruskal.test(dados$Doses50, dados$grupos, group=FALSE) Estou tentando reproduzir o resultado de um artigo, porém os resultados que obtenho dessa forma são completamente diferentes. Gostaria de saber se o teste de Kruskal-Wallis a ser feito nesse caso (de apenas 1 observação em cada grupo) seria dessa forma mesmo. []'s Em 30 de abril de 2014 13:10, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Manda um CMR, assim facilita para os listeiros. CMR: Comando Mínimo Reproduzível
[ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
Você quer comparar o que com o que? Ou você quer fazer um teste de ajustamento? Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 12:23, Letícia Pinheiro <leticia.ufmg@gmail.com> escreveu: O código que estou usando é o seguinte: dados<-data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12)) kruskal.test(dados$Doses50, dados$grupos, group=FALSE) Estou tentando reproduzir o resultado de um artigo, porém os resultados que obtenho dessa forma são completamente diferentes. Gostaria de saber se o teste de Kruskal-Wallis a ser feito nesse caso (de apenas 1 observação em cada grupo) seria dessa forma mesmo. []'s Em 30 de abril de 2014 13:10, Edson Lira <edinhoestat@yahoo.com.br> escreveu: Manda um CMR, assim facilita para os listeiros.
CMR: Comando Mínimo Reproduzível
[ ]'s.
Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro <leticia.ufmg@gmail.com> escreveu:
Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro _______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
Cada valor é referente a uma população. Quero saber se alguns dos valores são significativamente de outros, ou seja, se há diferenças entre os valores das populações examinadas. Em 30 de abril de 2014 13:32, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Você quer comparar o que com o que? Ou você quer fazer um teste de ajustamento?
Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 12:23, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: O código que estou usando é o seguinte:
dados<-data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12))
kruskal.test(dados$Doses50, dados$grupos, group=FALSE)
Estou tentando reproduzir o resultado de um artigo, porém os resultados que obtenho dessa forma são completamente diferentes. Gostaria de saber se o teste de Kruskal-Wallis a ser feito nesse caso (de apenas 1 observação em cada grupo) seria dessa forma mesmo.
[]'s
Em 30 de abril de 2014 13:10, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Manda um CMR, assim facilita para os listeiros. CMR: Comando Mínimo Reproduzível
[ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
-- Letícia Cavalari Pinheiro
Leticia, O teste permite ver isso sim (isso é, é possível calcular). Mas entenda que para qualquer conjunto de 12 observações tu vai ter o mesmo p-valor: dados <- data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12)) dados2 <- data.frame(Doses50 = rnorm(12), grupos = c(1:12)) kruskal.test(Doses50 ~ grupos, dados) # Teus dados kruskal.test(Doses50 ~ grupos, dados2) # Dados completamente aleatórios Isso porque o teste se baseia no ranking de cada observação, e não na observação em si. Um exemplo mais extremo: dados3 <- data.frame(Doses50 = c(1:11, 9999999999999999), grupos = c(1:12)) kruskal.test(Doses50 ~ grupos, dados3) Segundo o teste, as populações são iguais (p = 0.4433), mas temos, obviamente, um bem diferente dos demais. Em resumo: Com 1 observação de cada grupo, o teste não é poderoso o suficiente para detectar diferenças. 2014-04-30 13:36 GMT-03:00 Letícia Pinheiro <leticia.ufmg@gmail.com>:
Cada valor é referente a uma população. Quero saber se alguns dos valores são significativamente de outros, ou seja, se há diferenças entre os valores das populações examinadas.
Em 30 de abril de 2014 13:32, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Você quer comparar o que com o que? Ou você quer fazer um teste de
ajustamento?
Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 12:23, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: O código que estou usando é o seguinte:
dados<-data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12))
kruskal.test(dados$Doses50, dados$grupos, group=FALSE)
Estou tentando reproduzir o resultado de um artigo, porém os resultados que obtenho dessa forma são completamente diferentes. Gostaria de saber se o teste de Kruskal-Wallis a ser feito nesse caso (de apenas 1 observação em cada grupo) seria dessa forma mesmo.
[]'s
Em 30 de abril de 2014 13:10, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Manda um CMR, assim facilita para os listeiros. CMR: Comando Mínimo Reproduzível
[ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
Oi Rodrigo, obrigada pelo esclarecimento. Esse foi o resultado que encontrei mesmo. Porém o artigo que estou tentando reproduzir mostra um resultado bem diferente. Nesse caso, teria algum outro teste mais adequado para fazer a comparação entre as populações? []´s Em 30 de abril de 2014 14:00, Rodrigo Coster <rcoster@gmail.com> escreveu:
Leticia,
O teste permite ver isso sim (isso é, é possível calcular). Mas entenda que para qualquer conjunto de 12 observações tu vai ter o mesmo p-valor:
dados <- data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12)) dados2 <- data.frame(Doses50 = rnorm(12), grupos = c(1:12))
kruskal.test(Doses50 ~ grupos, dados) # Teus dados kruskal.test(Doses50 ~ grupos, dados2) # Dados completamente aleatórios
Isso porque o teste se baseia no ranking de cada observação, e não na observação em si. Um exemplo mais extremo:
dados3 <- data.frame(Doses50 = c(1:11, 9999999999999999), grupos = c(1:12)) kruskal.test(Doses50 ~ grupos, dados3)
Segundo o teste, as populações são iguais (p = 0.4433), mas temos, obviamente, um bem diferente dos demais.
Em resumo: Com 1 observação de cada grupo, o teste não é poderoso o suficiente para detectar diferenças.
2014-04-30 13:36 GMT-03:00 Letícia Pinheiro <leticia.ufmg@gmail.com>:
Cada valor é referente a uma população. Quero saber se alguns dos valores
são significativamente de outros, ou seja, se há diferenças entre os valores das populações examinadas.
Em 30 de abril de 2014 13:32, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Você quer comparar o que com o que? Ou você quer fazer um teste de
ajustamento?
Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 12:23, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: O código que estou usando é o seguinte:
dados<-data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12))
kruskal.test(dados$Doses50, dados$grupos, group=FALSE)
Estou tentando reproduzir o resultado de um artigo, porém os resultados que obtenho dessa forma são completamente diferentes. Gostaria de saber se o teste de Kruskal-Wallis a ser feito nesse caso (de apenas 1 observação em cada grupo) seria dessa forma mesmo.
[]'s
Em 30 de abril de 2014 13:10, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Manda um CMR, assim facilita para os listeiros. CMR: Comando Mínimo Reproduzível
[ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
Leticia, Provavelmente o artigo apenas resume um estudo mostrando as medidas por grupo mas no teste que usaram tinham mais alguma informação ou todos os dados para fazer o kruskal-wallis 2014-04-30 14:03 GMT-03:00 Letícia Pinheiro <leticia.ufmg@gmail.com>:
Oi Rodrigo,
obrigada pelo esclarecimento. Esse foi o resultado que encontrei mesmo. Porém o artigo que estou tentando reproduzir mostra um resultado bem diferente. Nesse caso, teria algum outro teste mais adequado para fazer a comparação entre as populações?
[]´s
Em 30 de abril de 2014 14:00, Rodrigo Coster <rcoster@gmail.com> escreveu:
Leticia,
O teste permite ver isso sim (isso é, é possível calcular). Mas entenda que para qualquer conjunto de 12 observações tu vai ter o mesmo p-valor:
dados <- data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12)) dados2 <- data.frame(Doses50 = rnorm(12), grupos = c(1:12))
kruskal.test(Doses50 ~ grupos, dados) # Teus dados kruskal.test(Doses50 ~ grupos, dados2) # Dados completamente aleatórios
Isso porque o teste se baseia no ranking de cada observação, e não na observação em si. Um exemplo mais extremo:
dados3 <- data.frame(Doses50 = c(1:11, 9999999999999999), grupos = c(1:12)) kruskal.test(Doses50 ~ grupos, dados3)
Segundo o teste, as populações são iguais (p = 0.4433), mas temos, obviamente, um bem diferente dos demais.
Em resumo: Com 1 observação de cada grupo, o teste não é poderoso o suficiente para detectar diferenças.
2014-04-30 13:36 GMT-03:00 Letícia Pinheiro <leticia.ufmg@gmail.com>:
Cada valor é referente a uma população. Quero saber se alguns dos valores
são significativamente de outros, ou seja, se há diferenças entre os valores das populações examinadas.
Em 30 de abril de 2014 13:32, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Você quer comparar o que com o que? Ou você quer fazer um teste de
ajustamento?
Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 12:23, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: O código que estou usando é o seguinte:
dados<-data.frame(Doses50 = c(0.58, 0.61, 0.64, 0.69, 0.70, 0.70, 0.70, 0.71, 0.75, 0.78, 0.85, 0.86), grupos = c(1:12))
kruskal.test(dados$Doses50, dados$grupos, group=FALSE)
Estou tentando reproduzir o resultado de um artigo, porém os resultados que obtenho dessa forma são completamente diferentes. Gostaria de saber se o teste de Kruskal-Wallis a ser feito nesse caso (de apenas 1 observação em cada grupo) seria dessa forma mesmo.
[]'s
Em 30 de abril de 2014 13:10, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Manda um CMR, assim facilita para os listeiros. CMR: Comando Mínimo Reproduzível
[ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 30 de Abril de 2014 11:57, Letícia Pinheiro < leticia.ufmg@gmail.com> escreveu: Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Letícia Cavalari Pinheiro
_______________________________________________ R-br mailing list R-br@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.
On 04/30/2014 12:56 PM, Letícia Pinheiro wrote:
Prezados,
preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso.
Obrigada.
-- Letícia Cavalari Pinheiro
Leticia, O teste em questão serve para compara se os postos de cada população são diferentes. Você acha que tendo apenas 1 observação em cada amostra é razoável a realização deste teste ? Na verdade, a DL50 corresponde a dose que seja suficiente para matar metade (50%) da população, desta forma seriam necessários - pelo menos - 2 observações por população. Com apenas 1 observação você vai obter varias letalidades de 0% [s] Tura
participantes (5)
-
Bernardo Rangel Tura -
Cesar Rabak -
Edson Lira -
Letícia Pinheiro -
Rodrigo Coster