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

[Rodrigo Coster]

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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140430/da76a37e/attachment.html>