<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Sto cercando di trovare la media geometrica dei batteri fecali per circa 15 anni di dati su circa 30-40 stazioni di monitoraggio. L'obiettivo del mio script è di condurre due filtri (uno sulla colonna ID progetto e un altro sulla colonna Nome caratteristica), quindi raggruppare il geo per stazione, poi per anno, e nella tabella finale mostrare ID stazione, Anno, Caratteristica Nome, N, geo media e mediana.
Ora so che non posso includere la sceneggiatura per trovare la geo mean nella parte di riepilogo della sceneggiatura, ma non sono stato in grado di capire come dovrei collocarla in questo script. Non sono sicuro se ho bisogno di ristrutturare completamente la sceneggiatura o di farlo in un modo completamente diverso.
Qualsiasi consiglio o suggerimento sarebbe molto apprezzato!
myPath = "O: / Planning / BPB / IBP_Pilot_Files / JMPScrips / Geomean /";
dt = Tabella dati ("Data_dwq_20140502_133555_RegRes_For_Testing.jmp");
dt << Selezionare Dove (: ID progetto == "NCAMBNT" |: ID progetto == "NCISCO" |: ID progetto == "TPBA");
dt << Sottoinsieme ((Righe selezionate), Nome tabella di output ("sottoinsieme"));
dt = data table ("sottoinsieme");
dt << Selezionare Dove (: Nome caratteristico == "Coliformi fecali" |: Nome caratteristico == "Coliformi totali");
dt << Sottoinsieme ((Righe selezionate), Nome tabella di output ("sottoinsieme2"));
dt = data table ("subset2");
dtsum = dt << Sommario (
Gruppo (: identificativo della stazione,: anno,: nome caratteristico),
N (: valore risultante come numero),
Exp (Col Mean (Log (: Valore risultato come numero))),
Mediana (: valore risultato come numero)
);
Ciao Melanie.
Ecco un modo: una volta pronti i dati per riepilogare, creare una colonna di registro per la variabile desiderata. Quindi utilizzare Sommario per trovare il mezzo di gruppo di questi dati di registro. Quindi esponenzia i risultati e inseriscili in una nuova colonna.
Ecco un esempio:
dt = open ("$ Sample_Data / Car Physical Data.jmp");
dt << nuova colonna ("peso_log", formula (log (: peso)));
dtSum = (dt << summary (mean (: peso_log), gruppo (: Paese,: tipo)));
dtsum << nuova colonna ("Geom_mean_weight", formula (exp (: name ("Mean (log_weight)"))));
Saluti,
Brady