Notação matricial para regressão logística

15

Na regressão linear (perda ao quadrado), usando matriz, temos uma notação muito concisa para o objetivo

minimizar  __UMAx-b__2

Onde é a matriz de dados, x são os coeficientes eb é a resposta.UMAxb

Existe uma notação matricial semelhante para o objetivo da regressão logística? Todas as anotações que eu vi não conseguem se livrar da soma em todos os pontos de dados (algo como dumatumaeulogístico(y,βTx) ).


Edição: obrigado por joceratops e grande resposta de AdamO. A resposta deles me ajudou a perceber que outra razão pela qual a regressão linear tem uma notação mais concisa é porque a definição da norma, que encapsula o quadrado e a soma ou ee . Mas na perda logística, não existe essa definição, o que torna a notação um pouco mais complicada.

Haitao Du
fonte

Respostas:

18

Na regressão linear, a solução Maximize Likelihood Estimation (MLE) para estimar x possui a seguinte solução de formulário fechado (supondo que A seja uma matriz com classificação completa da coluna):

x^lin=argminx__UMAx-b__22=(UMATUMA)-1UMATb

Isso é lido como "encontre o que minimiza a função objetivo, ". O bom de representar a função objetivo de regressão linear dessa maneira é que podemos manter tudo na notação matricial e resolver manualmente. Como Alex R. menciona, na prática geralmente não consideramos diretamente porque é computacionalmente ineficiente e geralmente não atende aos critérios de classificação completa. Em vez disso, nos voltamos para o pseudoinverso de Moore-Penrose . Os detalhes da resolução computacional para o pseudo-inverso podem envolver a decomposição de Cholesky ou a decomposição de valor singular.x__UMAx-b__22x lin ( A T A ) - 1 Umx^lin(UMATUMA)-1UMA

Alternativamente, a solução MLE para estimar os coeficientes na regressão logística é:

x^registro=argminxEu=1Ny(Eu)registro(1+e-xTuma(Eu))+(1-y(Eu))registro(1+exTuma(Eu))

onde (supondo que cada amostra de dados seja armazenada em linhas):

x é um vetor representa coeficientes de regressão

uma(Eu) i t h Um é um vetor que representa a / amostra / linha na matriz de dadosEuthUMA

y(Eu) { 0 , 1 } i t h i t h é um escalar em e o rótulo corresponde à amostra{0 0,1}EuthEuth

NA é o número de amostras de dados / número de linhas da matriz de dados de .UMA

Novamente, isso é lido como "encontre o que minimiza a função objetivo".x

Se você quisesse, poderia dar um passo adiante e representar na notação de matriz da seguinte maneira:x^registro

x^registro=argminx[1(1-y(1))1(1-y(N))][registro(1+e-xTuma(1))...registro(1+e-xTuma(N))registro(1+exTuma(1))...registro(1+exTuma(N))]

mas você não ganha nada fazendo isso. A regressão logística não possui uma solução de formulário fechado e não obtém os mesmos benefícios que a regressão linear, representando-a na notação matricial. Para resolver são utilizadas técnicas de estimativa como descida em gradiente e o método Newton-Raphson. Usando algumas dessas técnicas (por exemplo, Newton-Raphson), é aproximado e representado em notação matricial ( veja o link fornecido por Alex R. ).x^registro x logx^registro

joceratops
fonte
Ótimo. Obrigado. Acho que a razão pela qual não temos algo como resolver é a razão pela qual não damos mais esse passo para fazer a notação da matriz e evitar o símbolo da soma. UMAUMAx=UMAb
Haitao Du
Nós temos alguma vantagem de dar um passo adiante, transformá-lo em multiplicação de matrizes para simplificar o código, e em muitas plataformas como o matlab, para loop com soma em todos os dados, é muito mais lento que as operações de matriz.
Haitao Du
5
@ hxd1011: Apenas um pequeno comentário: reduzir as equações da matriz nem sempre é sábio. No caso de , você não deve tentar procurar inversa de matriz , mas sim algo como uma decomposição de Cholesky que será muito mais rápida e estável numericamente. Para a regressão logística, existem vários esquemas de iteração diferentes que de fato usam cálculos matriciais. Para uma grande revisão ver aqui: research.microsoft.com/en-us/um/people/minka/papers/logreg/...A T AUMATUMAx=UMATbUMATUMA
Alex R.
1
@AlexR. Muito obrigado. Aprendi que o uso da equação normal tornará o número condicional da matriz ao quadrado. E QR ou Cholesky seria muito melhor. Seu link é ótimo, essa revisão com métodos numéricos é sempre o que eu queria.
Haitao Du
15

A resposta do @joceratops se concentra no problema de otimização da máxima probabilidade de estimativa. Esta é realmente uma abordagem flexível que é passível de muitos tipos de problemas. Para estimar a maioria dos modelos, incluindo modelos de regressão linear e logística, existe outra abordagem geral baseada no método de estimativa de momentos.

O estimador de regressão linear também pode ser formulado como a raiz da equação de estimativa:

0 0=XT(Y-Xβ)

Nesse sentido, é visto como o valor que recupera um resíduo médio de 0. Ele não precisa contar com nenhum modelo de probabilidade subjacente para ter essa interpretação. É, no entanto, interessante deduzir as equações de pontuação para uma probabilidade normal; você verá que elas assumem exatamente a forma exibida acima. Maximizar a probabilidade de família exponencial regular para um modelo linear (por exemplo, regressão linear ou logística) é equivalente a obter soluções para suas equações de pontuação.β

0 0=Eu=1nSEu(α,β)=βregistroeu(β,α,X,Y)=XT(Y-g(Xβ))

Onde possui o valor esperado . Na estimativa de GLM, é considerado o inverso de uma função de link. Nas equações de probabilidade normal, é a função de identidade e na regressão logística é a função de logit. Uma abordagem mais geral seria exigir que permita a especificação incorreta do modelo.YEug(XEuβ)gg-1g-10 0=Eu=1nY-g(XEuβ)

Além disso, é interessante notar que, para famílias exponenciais regulares, que é chamado de relação de variância média. De fato, para a regressão logística, a relação de variação média é tal que a média está relacionada à variação por . Isso sugere uma interpretação de um modelo GLM mal especificado como sendo aquele que fornece um resíduo médio de 0 Pearson. Isso sugere ainda uma generalização para permitir derivadas médias funcionais não proporcionais e relações de variação média.g(Xβ)β=V(g(Xβ))p=g(Xβ)var(YEu)=pEu(1-pEu)

Uma abordagem de equação de estimativa generalizada especificaria modelos lineares da seguinte maneira:

0 0=g(Xβ)βV-1(Y-g(Xβ))

Com uma matriz de variações baseada no valor ajustado (média) dado por . Essa abordagem de estimativa permite escolher uma função de link e uma relação de variação média como nos GLMs.Vg(Xβ)

Na regressão logística, seria o logit inverso e seria dado por . As soluções para esta equação de estimativa, obtidas por Newton-Raphson, produzirão o obtido a partir da regressão logística. No entanto, uma classe um pouco mais ampla de modelos é estimada sob uma estrutura semelhante. Por exemplo, a função de link pode ser considerada o log do preditor linear, de modo que os coeficientes de regressão são riscos relativos e não razões de chances . O que - dadas as armadilhas bem documentadas de interpretar as ORs como RRs - me cabe perguntar por que alguém se encaixa mais nos modelos de regressão logística.gVEuEug(XEuβ)(1-g(Xβ))β

AdamO
fonte
1
+1 ótima resposta. formulá-lo como uma descoberta raiz sobre derivado é realmente novo para mim. e a segunda equação é realmente concisa.
Haitao Du