Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteSto cercando di fare un'analisi di sopravvivenza, in cui sto combinando alcune analisi in un copione. Quindi quello che sto cercando di fare è inserire le colonne nella piattaforma di sopravvivenza, quindi prendo quelle colonne e le uso nella piattaforma dei rischi proporzionali. Il mio script funziona alla grande tranne che per la risposta Y, non troverà la variabile di risposta Y. Se cambio Y( As Column( Ycolname ) ), a qualcosa come Y("Time to Failure"), funziona alla grande.Ma quando provo a usare la variabile, non funziona.
</CODE></P><P>Ho due domande: perché il mio metodo funziona con Censor e Xs, ma non con Ys.</P><P>Seconda domanda, se non inseriscono una colonna per o frequenza, come faccio a dire al modello di adattamento che non c'è una colonna da inserire. Se uso As Column (variabile) dice che non riesce a trovare la mia colonna, il che è giusto.</P><P> </P><P> </P><P> </P><PRE><CODE class=" language-jsl">Names Default To Here( 1 );
Clear Log();
dt = Current Data Table();
If( Not( Is Scriptable( dt ) ),
dtsel = Pick File( "Open File", {"JMP|jmp", "All files|*"} );
If( dtsel == "",
Stop(),
Try( dt = Open( dtsel ), Stop() )
);
);
dtname = Data Table( dt ) << Get Name;
validated = 0;
New Window( "Survival/Reliability",
<<MODAL></MODAL>
AnaliticaMac OSScriptingRi: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteGrazie per la tua guida, sei stato utilissimo!
Il mio problema era: se nella prima schermata non c'era nessuna colonna messa in frequenza,
Freqname = Arg (sopravvivenza[ 4 ], 1 );
Questo ritornerebbe come vuoto.
Quindi, se ho provato a eseguire lo script:
Piattaforma = dt << Modello adatto (
frequenza( Colonna ( Nomefrequenza ) ),
Di( Colonna ( Per nome ) ),
Y( Colonna ( Ycolname ) ),
Censurare( Colonna (Nome Censore)),
Effetti( Colonna ( Nome effetto ) ),
Personalità( "Rischio proporzionale" ),
Codice censore( "1" ),
Esegui modello( Test del rapporto di verosimiglianza( 1 ), Intervalli di confidenza di probabilità( 1 ), rapporti di rischio( 1 ) )
)
Direbbe che non esiste una colonna per la frequenza, il che era vero.
Quindi quello che ho fatto è stato inserire 8 cicli if nidificati cercando di vedere quali colonne erano vuote.
Forse non è esattamente il modo migliore, ma ha funzionato.
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteNon sono sicuro di aver capito il problema. Ecco cosa ho provato e visto.
Per prima cosa modifico Rats.jmp per simulare un caso in modo da poter testare i tuoi codici:
Ho aggiunto una colonna Per e ho diviso equamente i dati in due gruppi.
Ho modificato l'ultimo pezzo dei tuoi codici nel seguente.
if (isempty(byname),
Platform = dt << Fit Model(
// Freq(),
// By(),
Y( Column( Ycolname ) ),
Censor( Column( Censorname ) ),
Effects( Column( Effectname ) ),
Personality( "Proportional Hazard" ),
Censor Code( "1" ),
Run Model( Likelihood Ratio Tests( 1 ), Likelihood Confidence Intervals( 1 ), Risk Ratios( 1 ) )
),
Platform = dt << Fit Model(
// Freq(),
By( column(byname) ),
Y( Column( Ycolname ) ),
Censor( Column( Censorname ) ),
Effects( Column( Effectname ) ),
Personality( "Proportional Hazard" ),
Censor Code( "1" ),
Run Model( Likelihood Ratio Tests( 1 ), Likelihood Confidence Intervals( 1 ), Risk Ratios( 1 ) )
)
);
Ho provato due percorsi per eseguire i codici:
Innanzitutto, esegui il codice e inserisci i ruoli come segue:
Nessun "da". e il comportamento dei codici è come prima.
In secondo luogo, esegui il codice e inserisci i ruoli come segue:
Con "By", e vedo un messaggio di errore alla fine:
Immagino che i risultati intermedi siano diversi e che i tuoi codici non abbiano gestito correttamente la situazione.
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteHai qualche suggerimento sulla mia seconda domanda?
se non inseriscono una colonna per o frequenza, come faccio a dire al modello di adattamento che non c'è una colonna da inserire. Se uso Colonna (variabile) dice che non riesce a trovare la mia colonna, il che è giusto.
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteSo di aver provato questo e non ha funzionato per me, ma quando ho incollato il tuo codice ha funzionato perfettamente. Non ho idea del perché, ma funziona quindi GRAZIE!
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteNon riesco a riprodurre il problema se cambio As Column in Column. Ecco i codici completi che eseguo. Ho appena cambiato gli ultimi tre As Column in Column.
Names Default To Here( 1 );
Clear Log();
dt = Current Data Table();
If( Not( Is Scriptable( dt ) ),
dtsel = Pick File( "Open File", {"JMP|jmp", "All files|*"} );
If( dtsel == "",
Stop(),
Try( dt = Open( dtsel ), Stop() )
);
);
dtname = Data Table( dt ) << Get Name;
validated = 0;
New Window( "Survival/Reliability",
<<MODAL></MODAL>
Ecco i miei passi e quello che ho visto. Ho controllato JMP14, JMP15 e JMP16. Credo che la tua versione non sia più indietro a causa del modo in cui i tuoi codici usano i nuovi nomi letterali.
- Aprire Rats.jmp
- Esegui lo script
- Specificare la finestra di dialogo quando viene visualizzata, quindi fare clic su OK
- Vengono visualizzati alcuni rapporti. E questa è la finestra di dialogo di avvio del modello di adattamento simile a:
Se non riesci a vedere quello che sto vedendo qui, posso solo sospettare che ci sia qualcosa di cui nessuno di noi due è a conoscenza. In tal caso, potresti voler contattare support@jmp.com per avere un po' di aiuto.
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteLo so, ho ottenuto la stessa cosa in cui il modello di adattamento non è venuto fuori, se ho spostato Eval (PH) in basso, motivo per cui l'ho spostato dove l'ho fatto. Ho anche riscontrato lo stesso errore con Rats, ma non questo errore con i miei file, non so perché. Se estraggo l'espressione e la lascio funzionare in questo modo, l'errore scompare e ottengo lo stesso problema di prima, senza Expr(PH) è incollato sotto l'immagine dello schermo, le altre colonne vengono visualizzate, ma non la Y . Questo script è corretto e il mio più recente. Se avete qualche idea, per favore fatemelo sapere! Questo è molto frustrante per me e non ho idea del perché!
Names Default To Here( 1 );
Clear Log();
dt = Current Data Table();
If( Not( Is Scriptable( dt ) ),
dtsel = Pick File( "Open File", {"JMP|jmp", "All files|*"} );
If( dtsel == "",
Stop(),
Try( dt = Open( dtsel ), Stop() )
);
);
dtname = Data Table( dt ) << Get Name;
validated = 0;
New Window( "Survival/Reliability",
<<MODAL></MODAL>Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancantePotresti per favore ricontrollare il tuo script? Non vedo come arrivare alla finestra di dialogo che mostri qui.
Ho ricevuto un messaggio di errore:
Se sposto Eval(PH) in fondo allo script e cambio As Column in Column, non vedo affatto la finestra di dialogo Fit Model.
A proposito, io uso Rats.jmp per testare. "days" su Y, "Censor" su Censor, "Group" su Grouping.
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteColumn fa la stessa cosa di As Column, inserisce le informazioni ad eccezione della Y:
Ri: Rischi proporzionali del modello di adattamento variabile della risposta Y mancanteSe usi Column, invece di As Column, ottieni quello che vuoi? In realtà non vedo in quale situazione, si deve usare As Column, piuttosto che Column. Forse, uno dei miei colleghi che ha più conoscenze su questa funzione può commentare.
