O procedimento básico de correspondência de pontuação de propensão funciona com dados de seção transversal (ou seja, coletados em um determinado momento). O comando psmatch2 popular usa uma variável dummy indicando que uma observação pertence ao grupo de tratamento ou controle.
No entanto, no meu conjunto de dados, essa função indicadora varia no tempo. Os dados têm a seguinte aparência: identifico indivíduos e os acompanho por um período de 12 anos. Em algum momento durante esse período, os indivíduos podem ser tratados e, portanto, "transferem do controle para o grupo de tratamento". Observe que isso pode ocorrer a qualquer momento (sic!) Para cada indivíduo.
Acima disso, minha análise subsequente seria dependente do tempo: a partir do momento em que os indivíduos estão sendo tratados, quero estimar o TCA, mas apenas para o ano seguinte (sic!). Assim, não me importo mais com os indivíduos tratados 1 ano após o início do tratamento.
psmatch2
parece ser muito restritivo quando se trata disso. Você vê alguma possibilidade de como fazer isso? Talvez também haja métodos de correspondência que me permitam fazer isso.
PS: Há uma pergunta semelhante aqui no fórum, no entanto, isso não me ajudou a resolver esse problema.
fonte
Respostas:
Talvez o artigo a seguir seja relevante para o seu caso: Lu B. Pontuação de propensão correspondente a covariáveis dependentes do tempo. Biometrics 2005; 61, 721-728 .
Na situação considerada no artigo, os sujeitos podem iniciar o tratamento a qualquer momento durante um período de observação. Um indivíduo que fica exposto no momento é correspondido a vários controles selecionados no conjunto de riscos correspondente, ou seja, de todos os indivíduos que ainda correm o risco de ficar expostos no tempo .tt t
A correspondência refere-se a um escore de propensão dependente do tempo, definido como o risco de ficar exposto no tempo calculado a partir de um modelo de riscos proporcionais de Cox: onde é um vetor de preditores potencialmente variáveis no tempo do status do tratamento. Em cada conjunto de riscos, a correspondência é efetivamente realizada na escala preditora linear de acordo com a métrica h ( t ) = h 0 ( t ) exp ( β ' x ( t ) ) x ( t ) d ( x i ( t ) , x j ( t ) ) = ( β ' x i ( t ) - β ′ x j ( t ) ) 2 .t
fonte
O Stata 13 possui um estimador de efeitos de tratamento com vários valores . Pode ser possível reformular o seu problema como tratamento de múltiplos valores, em que o tratamento é indexado por tempo (tratado no ano 1, tratado no ano 2, ..., tratado no ano 12) em vez de binário. O resultado será medido no ano após o tratamento.
fonte