modelli di pericolo aiutano gli statistici e gli altri ricercatori di valutare il rischio attraverso la comprensione delle variabili associate a rischio . Un tipo di modello di rischio , il modello dei rischi proporzionali di Cox , è una forma semplificata e specifica di questo modello che si basa su una forte assunzione : le variabili associate con il rischio di pericolo sono legati in modo moltiplicativo ( si moltiplicano nel modello invece di aggiungere ) . Prima di poter utilizzare con successo un modello di rischio proporzionale di Cox , è necessario verificare questa ipotesi . Utilizzando SAS , un potente pacchetto software per le statistiche , si può facilmente verificare questa ipotesi sui dati. Istruzioni
1
Inserire i dati in SAS . Fare clic su "File" e " Importa dati " . Seguire le istruzioni visualizzate e aprire il file di dati
2
Tipo " proc phreg data = yourdata ; " . , Dove " yourdata " è la matrice di dati che deve essere testato . Questo dice SAS i dati che si dovrà lavorare con .
3
Raddoppia le variabili nei dati , moltiplicando ciascuno con il registro di tempo . Per ogni variabile , utilizzare il comando " nuovavar = vecchiavar * log ( tempo ) ; " . Ad esempio , se si dispone di una variabile di salario chiamato " stipendio " , immettere il comando " newsalary = stipendio * log ( tempo ) ; " . Fate questo per tutte le variabili nel dataset .
4
tipo "tempo di modello * censurare ( 0 ) = yourmodel " , dove " yourmodel " è l'insieme di variabili che si desidera includere nel modello . Questo insieme di variabili dovrebbe includere tutte le variabili originali nei 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 per il registro del tempo . Ad esempio , se si dispone di due variabili nel vostro studio , lo stipendio e l'età , si avranno due nuove variabili , newsalary e newage . Pertanto , si deve inserire " proportionality_test : prova newsalary newage ; "
6
Eseguire il test. . Tipo "run ";
7
Osservare i valori di p in uscita . . I valori di p si trovano nella colonna entitiled "Pr > ChiSq " . Se si notano i valori- p sotto 0.05 , il test ha esito negativo , il che significa l'assunzione di rischi proporzionali non riesce.