Come verificare un assunto proporzionali pericoli con SAS

February 22

modelli di pericolo aiutano statistici e altri ricercatori di valutare il rischio attraverso la comprensione delle variabili associate a rischio. Un tipo di modello di rischio, il rischio proporzionale di Cox modello, è una forma semplificata e specifico di questo modello che si basa su una forte assunzione: le variabili associate con pericolo per la sicurezza sono legati in modo moltiplicativo (si moltiplicano nel modello in contrasto con l'aggiunta di ). Prima di poter utilizzare con successo un modello di rischio proporzionale di Cox, è necessario verificare questa ipotesi. Utilizzando SAS, un potente pacchetto di software statistico, si può facilmente verificare questa ipotesi sui dati.

istruzione

1 Inserire i dati nel SAS. Fare clic su "File" e "Importa dati". Seguire i passaggi che appaiono e aprire il file di dati.

2 Digitare "proc phreg dati = yourdata;", dove "yourdata" è la matrice di dati che devono essere testati. Questo dice SAS quali dati si dovrà lavorare con.

3 Raddoppiare le variabili nei dati, moltiplicando ciascuno con il registro del tempo. Per ogni variabile, utilizzare il comando "newvar = oldvar

log (tempo);". Ad esempio, se si dispone di una variabile per stipendio chiamato "stipendio", immettere il comando "newsalary = log stipendio (tempo);". Fate questo per tutte le variabili nel set di dati.

4 Tipo "* censore tempo modello (0) = yourmodel", dove "yourmodel" è l'insieme di variabili che si desidera includere nel modello. Questo insieme di variabili dovrebbe includere tutte le variabili originali in dati, così come le nuove variabili che sono stati moltiplicati per il registro del tempo.

5 Preparare il test di proporzionalità per le nuove variabili. Tipo "proportionality_test: newvariables prova;", dove "newvariables" è l'insieme di tutte le nuove variabili create moltiplicando le vecchie variabili dal registro di tempo. Ad esempio, se si dispone di due variabili nel vostro studio, stipendio e l'età, si avranno due nuove variabili, newsalary e newage. Pertanto, si dovrebbe inserire "proportionality_test: prova newsalary newage;".

6 Eseguire il test. Tipo "correre;".

7 Osservare i valori di p in uscita. I valori di p si trovano nella colonna entitiled "Pr> ChiSq". Se si notano valori p sotto 0.05, il test ha esito negativo, il che significa l'assunzione di rischi proporzionali fallisce.