[R-br] Kruskal-Wallis para grupos com apenas 1 elemento

Cesar Rabak cesar.rabak em gmail.com
Quinta Maio 1 19:20:44 BRT 2014


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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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 em 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/20140501/dd7b8776/attachment.html>


Mais detalhes sobre a lista de discussão R-br