Codice SAS in JMP?Vorrei sapere se esiste un modo per eseguire un codice SAS (sotto) in JMP. O come posso emulare il codice SAS di seguito in JMP.
proc mixed
Class year week treatment rep;
model stdmisummax=year|week|treatment/ddfm=kr residual;
random rep(year) treatment*rep(year);
repeated week/type=ar(1) subject=treatment*rep(year);
run;
I miei fattori sono: trattamento, settimana (misura ripetuta), anno, ripetizione (replica).
Grazie in anticipo.
PS: ho provato a eseguire qualcosa di simile in JMP utilizzando l'aggiunta di misure ripetute, ma non funzionerà, perché non posso incrociare tre fattori per rendere l'ID soggetto.
Penso che la mia versione di JMP non sia pro o sia una sorta di versione pro per l'università.
Non ho il software SAS.
Re: codice SAS in JMP?Sì, l'effetto casuale, Year(Subject), o qualunque sia il primo, è la variazione stimata tra i soggetti .
Re: codice SAS in JMP?@jiancao Grazie per la tua spiegazione.
Ogni volta che includo più di un effetto casuale, ottengo "0" per Estimate, Std Error (vedi screen shot) per uno degli effetti casuali. Perché dovrebbe essere? È a causa della nidificazione?
Grazie per l'aiuto,
Isadora
Re: codice SAS in JMP?È possibile aggiungere più di un effetto casuale nell'analisi delle misure ripetute. Tuttavia riceverai una nota che dice " La colonna dell'oggetto ripetuto è stata nidificata in tutti gli effetti casuali .” Guarda lo screenshot qui sotto. Penso che sia per questo che Adam ha suggerito di abbandonare Rep(Anno) perché si anniderà Trattamento.
Hai provato a includere due effetti casuali? Cos'hai visto?
Re: codice SAS in JMP?CIAO @AdamMorris E @DonMcCormack
Vorrei sapere se puoi aiutarmi a capire un po' meglio perché hai usato un solo termine come effetto casuale. È perché lo sono è consentito utilizzare un solo termine casuale in Effetti casuali quando utilizzo AR1?
Grazie,
Isadora
Re: codice SAS in JMP?@DonMcCormack Questo è molto utile!
Grazie per averlo condiviso con me e per avermi dato la tua opinione su quale soluzione sia l'approccio migliore. Lo apprezzo davvero.
Isadora
Re: codice SAS in JMP?La mia soluzione è qui sotto:
Fit Model(
Y( :Response ),
Effects(
:Treatment,
:Week,
:Treatment * :Week,
:Year,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Rep[:Treatment, :Year] ),
NoBounds( 1 ),
Personality( "Mixed Model" ),
Subject( :Treatment, :Year, :Rep ),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Run( Variogram( AR1( 1 ) ) )
);
Ci sono solo due differenze. Per l'effetto casuale aggiuntivo, ho annidato Rep in Year. Ho usato Treatment, Year e Rep per definire implicitamente l'oggetto. Tuttavia, trovo l'approccio di Adam di definire esplicitamente il Soggetto più intuitivamente attraente.
Re: codice SAS in JMP?@AdamMorris , sei stato così utile!
Grazie mille per aver condiviso le tue conoscenze e tutte queste buone informazioni! Lo apprezzo.
Re: codice SAS in JMP?@DonMcCormack Molte grazie per la tua risposta.
Ti dispiacerebbe condividere anche la tua soluzione? Mi piacerebbe confrontarli, perché sono interessato a saperne di più sulla codifica in JMP e altri approcci per i miei dati.
Grazie!
Isadora
Re: codice SAS in JMP?Ciao Isadora,
Per condividere con altri lettori di questo thread, la soluzione proposta (dal supporto tecnico JMP) è creare una nuova colonna Oggetto che sia la concatenazione delle 3 variabili che identificano i soggetti (Trattamento, Anno e Rep).
È inoltre necessario creare una copia continua della colonna Settimana da utilizzare come variabile Ripetuta.
JMP Pro non può adattarsi esattamente allo stesso modello specificato dal codice SAS (per i modelli AR(1), JMP nidifica gli effetti casuali nella variabile soggetto).Tuttavia, JMP può adattare il modello omettendo il termine Rep(Year).Lo script seguente mostra come si potrebbe impostare questo.
Grazie,
Adamo
New Column("Subject", Character, Formula( :Treatment || ", " || Char( :Year ) || ", " || Char( :Rep ) ));
New Column("Week2", numeric, continuous, formula(:Week));
Fit Model(
Y( :response ),
Effects(
:Treatment,
:Week,
:Year,
:Treatment * :Week,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Treatment*:Year*:Rep ),
Subject(:Subject),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Personality( "Mixed Model" )
);Re: codice SAS in JMP?Isadora,
Ho parlato con Adam Morris della tua domanda. Indipendentemente, ho ottenuto gli stessi risultati di lui, anche se ho impostato il modello in modo leggermente diverso. Mi è piaciuta un po' di più la sua configurazione, quindi userei quella. Per favore fatemi sapere se avete altre domande.
Assistente
Re: codice SAS in JMP?
CIAO @DonMcCormack , ho allegato uno screenshot di come appare la tabella dei dati per ogni soggetto.
Ho 4 repliche, 6 trattamenti e 2 anni. I campioni sono stati raccolti in 3 giorni diversi dopo l'applicazione dei trattamenti. Chiamo questi giorni di campionamento del numero di campionamento nella tabella dei dati e si chiamano "settimana" nello script SAS.
Ho JMP Pro 13.
Grazie in anticipo,
Isadora
Re: codice SAS in JMP?Sarebbe utile sapere qualcosa in più sulla struttura della tabella dati. Come sarebbe per un soggetto? Le ripetizioni si verificano entro una settimana? Avrai bisogno di JMP Pro per utilizzare la struttura di covarianza AR(1).
Re: codice SAS in JMP?Non sono un esperto di modelli di effetti misti in generale o di questa piattaforma JMP Pro in particolare. Forse qualcuno che è un tale esperto si unirà alla discussione. In caso contrario, è possibile inviarlo all'assistenza tecnica di JMP per ulteriori suggerimenti.
Re: codice SAS in JMP?@Mark_Bailey Grazie per aver condiviso questa nota. L'ho letto prima.
Ho la piattaforma mista nella mia versione di JMP. Il problema è che non riesco a stabilire correttamente l'ID del soggetto in questa piattaforma, perché il mio soggetto è un'interazione di trattamento*rep(anno). Quella piattaforma non mi consente di incrociare e nidificare per aggiungere l'oggetto "trattamento*rap(anno)" nella casella ID soggetto. Ti permette di aggiungere solo un fattore.
Grazie
Re: codice SAS in JMP?Dispiace per la confusione. No, non è possibile eseguire programmi SAS in JMP.
JMP è in grado di eseguire alcune semplici analisi di misure ripetute. Guarda questo Nota .Situazioni più complesse (ovvero non semplici strutture di covarianza) richiedono la piattaforma Mixed Model in JMP Pro.
Re: codice SAS in JMP?@Mark_Bailey Grazie per la tua risposta. Il problema è che non voglio eseguire le analisi in SAS. Voglio sapere se è possibile eseguire questo codice in JMP in qualche modo o come posso emulare quel codice in uno script JMP o in Fit Model.
Grazie.
Re: codice SAS in JMP?Nel caso in cui ciò si applichi alla tua situazione, dovresti sapere che JMP è un client completo per SAS. Quindi, proprio come SAS viene eseguito come client, anche JMP può farlo. Richiede una connessione al server SAS. Quel server potrebbe essere eseguito sullo stesso computer, come nel caso del mio computer portatile da lavoro, oppure potrebbe essere remoto tramite una connessione di rete. Hai accesso a un server SAS remoto (azienda, università)? In tal caso, incolla il codice SAS in un editor di script JMP. Fare clic con il pulsante destro del mouse nell'editor e selezionare Invia SAS.
Selezionare File > Preferenze e configurare prima la connessione al server SAS.
Re: codice SAS in JMP?@uday_guntupalli Grazie per la tua risposta. Sì, lo so, ma non voglio entrare in SAS. Ho ricevuto questo codice da un collega perché non sono riuscito o non sono riuscito a scoprire come eseguire le stesse analisi con i fattori stabiliti in questo modo in JMP.
Grazie 
Re: codice SAS in JMP?@bordini ,
So che dici di non avere SAS, ma sapevi che puoi scaricare SAS University Edition gratuitamente? Se hai davvero bisogno di testare il codice SAS, puoi scaricare SAS university edition e provarlo.
https://www.sas.com/en_us/software/university-edition/download-software.html
Non so come eseguire il passaggio equivalente in JMP, quindi cercherò anche altre risposte
