Regressão probit (análise dose-resposta)
| Comando: | Estatísticas |
|---|
O procedimento de regressão probit ajusta uma curva dose-resposta sigmoide probit e calcula os valores (com IC de 95%) da variável de dose que correspondem a uma série de probabilidades. Por exemplo, a DE50 (dose efetiva mediana, ED50) ou a DL50 (dose letal mediana, LD50) são os valores que correspondem a uma probabilidade de 0,50, e o limite de detecção (CLSI, 2012) é o valor que corresponde a uma probabilidade de 0,95.
A equação de regressão probit tem a forma:
em que X é a variável de dose (eventualmente transformada em logaritmo) e probit(p) é o valor da função de distribuição cumulativa normal padrão inversa Φ−1 correspondente a uma probabilidade p:
Probit(p) pode ser transformado em uma probabilidade p com a função de distribuição cumulativa normal padrão Φ:
O MedCalc ajusta os coeficientes de regressão a e b com o método da máxima verossimilhança.
Como inserir os dados
Você pode inserir os dados em formato binário ou em formato agrupado.
Binário
No formato binário, há 2 variáveis: uma variável para a dose (concentração) e outra para a resposta binária.
Para cada medição individual, há uma linha com a dose e a resposta, que é codificada como 0 (sem resposta) e 1 (resposta).
Por exemplo:
| A | B | C | |
|---|---|---|---|
| Concentration | Positive | ||
| 1 | 0,1 | 0 | |
| 2 | 0,2 | 0 | |
| 3 | 0,2 | 1 | |
| 4 | 0,1 | 0 | |
| 5 | 0,5 | 1 | |
| 6 | 0,5 | 0 | |
| 7 | 0,5 | 0 | |
| 8 | 0,9 | 1 | |
| 9 | 0,9 | 1 |
Agrupado
No formato agrupado, há 3 variáveis: uma variável para a dose, outra para o número total de medições e outra para o número de medições com resposta.
Por exemplo:
| A | B | C | D | |
|---|---|---|---|---|
| Concentration | Subjects | Positive | ||
| 1 | 0,1 | 10 | 0 | |
| 2 | 0,2 | 10 | 1 | |
| 3 | 0,3 | 10 | 3 | |
| 4 | 0,4 | 10 | 5 | |
| 5 | 0,5 | 10 | 7 | |
| 6 | 0,6 | 10 | 8 | |
| 7 | 0,7 | 10 | 10 | |
| 8 | 0,8 | 10 | 10 | |
| 9 | 0,9 | 10 | 10 |
Dados necessários
Tipo de dados
Selecione a opção que corresponde à forma como você inseriu os dados: binários ou agrupados (veja acima).
Variável de dose
Selecione a variável de dose.
Variáveis no caso de dados binários
- Variável de resposta: a variável de resposta deve ser binária ou dicotômica e deve conter apenas dados codificados como 0 (sem resposta) ou 1 (resposta). Se os seus dados estiverem codificados de outra forma, você pode usar a ferramenta Definir estado para recodificá-los.
Variáveis no caso de dados agrupados
- Número total de casos: selecione a variável que contém o número de medições de cada dose.
- Número de respostas: selecione a variável que contém o número de respostas de cada dose.
Filtro
(Opcionalmente) insira um filtro de dados para incluir apenas um subgrupo selecionado de casos na análise.
Opções
- Transformação logarítmica: selecione esta opção se a variável de dose exigir uma transformação logarítmica. Quando a variável de dose contém valores 0, o MedCalc adiciona automaticamente um número pequeno aos dados para tornar possível a transformação logarítmica. Esse número pequeno é subtraído quando os resultados são reconvertidos para a apresentação.
- Diagrama dose-resposta: selecione esta opção para obter um diagrama dose-resposta.Marcadores: clique nesta opção para que os dados sejam representados no gráfico como marcadores. Observe que, quando você seleciona a transformação logarítmica e a variável de dose contém valores 0, esses valores não podem ser representados no gráfico como marcadores.
Resultados
Variável de dose | Concentration |
|---|---|
Variável para o número de casos | Subjects |
Variável para o número de casos positivos | Positive |
| Tamanho da amostra | 90 |
|---|---|
Casos positivos | 54 (60.00%) |
Casos negativos | 36 (40.00%) |
Ajuste geral do modelo
| −2 log-verossimilhança do modelo nulo | 121.142 |
|---|---|
−2 log-verossimilhança do modelo completo | 56.959 |
Qui-quadrado | 64.183 |
GL | 1 |
Nível de significância | P < 0.0001 |
Cox & Snell R2 | 0.5099 |
Nagelkerke R2 | 0.6893 |
Coeficientes e Erros Padrão
Variável | Coeficiente | Erro Padrão | Wald | P |
|---|---|---|---|---|
Concentration | 6.35651 | 1.16609 | 29.7147 | <0.0001 |
Constante | -2.60520 | 0.52209 | 24.8993 | <0.0001 |
Tabela dose-resposta
Probabilidade | Dose | Intervalo de confiança de 95% | |
|---|---|---|---|
| 0.010 | 0.043869 | -0.17191 | 0.15357 |
| 0.020 | 0.086754 | -0.10711 | 0.18674 |
| 0.025 | 0.10151 | -0.084896 | 0.19823 |
| 0.05 | 0.15108 | -0.010650 | 0.23723 |
| 0.10 | 0.20823 | 0.073895 | 0.28324 |
| 0.20 | 0.27744 | 0.17369 | 0.34155 |
| 0.25 | 0.30374 | 0.21041 | 0.36490 |
| 0.50 | 0.40985 | 0.34643 | 0.47127 |
| 0.75 | 0.51596 | 0.45600 | 0.60411 |
| 0.80 | 0.54225 | 0.47982 | 0.64035 |
| 0.90 | 0.61146 | 0.53893 | 0.73934 |
| 0.95 | 0.66861 | 0.58532 | 0.82352 |
| 0.975 | 0.71819 | 0.62452 | 0.89756 |
| 0.980 | 0.73294 | 0.63605 | 0.91973 |
| 0.990 | 0.77583 | 0.66933 | 0.98442 |
Tamanho da amostra e casos com desfecho negativo e positivo
Primeiro, o programa apresenta o tamanho da amostra e o número e a proporção de casos com e sem resposta.
Ajuste geral do modelo
O −2 log-verossimilhança do modelo nulo é dado por −2 * ln(L0), em que L0 é a verossimilhança de obter as observações no modelo “nulo”, um modelo sem a variável de dose.
O −2 log-verossimilhança do modelo completo é dado por −2 * ln(L), em que L é a verossimilhança de obter as observações com a variável de dose incorporada no modelo.
A diferença entre esses dois valores dá uma estatística qui-quadrado, que é uma medida de quanto a variável de dose afeta a variável de resposta.
O R2 de Cox & Snell e o R2 de Nagelkerke são outras medidas de qualidade do ajuste, conhecidas como pseudo R ao quadrado. Observe que o pseudo R ao quadrado de Cox & Snell tem um valor máximo que não é 1. O R2 de Nagelkerke ajusta o de Cox & Snell para que o intervalo de valores possíveis se estenda até 1.
Coeficientes de regressão
Os coeficientes de regressão são os coeficientes a (constante) e b (inclinação) da equação de regressão:
A estatística de Wald é o coeficiente de regressão dividido pelo seu erro padrão, ao quadrado: (b/EP)2.
Transformação logarítmica
Quando você seleciona a transformação logarítmica da variável de dose, a e b são, na verdade, os coeficientes da equação de regressão:
Uso da equação ajustada
A probabilidade prevista de uma resposta positiva pode ser calculada com a equação de regressão.
Quando a equação de regressão é, por exemplo:
para uma dose de 0,500, probit(p) é igual a 0,57. Probit(p) pode ser transformado em p com a função de planilha do MedCalc NORMSDIST(z) ou com a função equivalente do Excel.
Em alternativa, você pode usar a tabela a seguir.
| Probit(p) | |
|---|---|
| 2,326 | 0,99 |
| 1,645 | 0,95 |
| 1,282 | 0,90 |
| 0,842 | 0,80 |
| 0,000 | 0,50 |
| -0,842 | 0,20 |
| -1,282 | 0,10 |
| -1,645 | 0,05 |
| -2,326 | 0,01 |
No exemplo, com probit(p) igual a 0,57, p = 0,72.
Uma probabilidade p pode ser transformada em probit(p) com a tabela acima ou com a função de planilha do MedCalc NORMSINV(p) ou com a função equivalente do Excel. Para uma probabilidade p=0,5, você encontra na tabela que probit(p)=0. Quando a equação de regressão é
probit = −2,61 + 6,36 × Dose
então
Dose = (probit+2,61)/6,36
e, portanto, dose = 2,61/6,36 = 0,41.
Tabela dose-resposta
Esta tabela lista uma série de probabilidades com a dose correspondente, com um intervalo de confiança de 95% da dose (Finney, 1947).
Os valores em texto cinza-claro são valores de dose que ficam fora do intervalo observado da variável de dose.
Transformação logarítmica
Quando você seleciona a transformação logarítmica da variável de dose, o MedCalc reconverte os resultados e mostra a variável de dose na sua escala original na tabela dose-resposta.
Gráfico
Este gráfico mostra as probabilidades, que variam de 0 a 1, e a dose correspondente. Duas curvas adicionais representam o intervalo de confiança de 95% da dose.
A dose e o intervalo de confiança de 95% correspondentes a uma determinada probabilidade são obtidos a partir de uma linha horizontal nesse nível de probabilidade.
Literatura
- CLSI (2012) Evaluation of detection capability for clinical laboratory measurement procedures; Approved guideline - 2nd edition. CLSI document EP17-A2. Wayne, PA: Clinical and Laboratory Standards Institute.
- Finney DJ (1947) Probit Analysis. A statistical treatment of the sigmoid response curve. Cambridge: Cambridge University Press.