<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />私は、およそ30から40の監視局にわたる約15年間のデータについて、糞便バクテリアの幾何平均を見つけようとしています。 私のスクリプトの目的は、2つのフィルター(1つはプロジェクトID列に、もう1つは特性名列に)を実行し、次に地理的手段をステーション、次に年ごとにグループ化し、最後の表にステーションID、年、特性を表示することです。名前、N、地理的意味、および中央値。
地理的意味を見つけるためのスクリプトをスクリプトの要約部分に含めることはできませんが、このスクリプトのどこに配置すればよいのかを把握できていません。 スクリプトを完全に再構築する必要があるのか、まったく別の方法でこれを実行する必要があるのかわかりません。
何かアドバイスや提案は大歓迎です!
myPath = "O:/計画/ BPB / IBP_Pilot_Files / JMPスクリプト/ Geomean /";
dt =データテーブル( "Data_dwq_20140502_133555_RegRes_For_Testing.jmp");
dt << Select Where(:プロジェクトID == "NCAMBNT" |:プロジェクトID == "NCISCO" |:プロジェクトID == "TPBA");
dt << Subset((Selected Rows)、出力テーブル名( "subset"));
dt =データテーブル(「サブセット」)。
dt << Select Where(:特性名== "糞便大腸菌" |:特性名== "合計大腸菌");
dt << Subset((Selected Rows)、出力テーブル名( "subset 2"));
dt =データテーブル(「サブセット2」)。
dtsum = dt <<概要(
グループ(:駅ID、:年、:特性名)、
N(:数値としての結果値)
Exp(Col Mean(対数(:結果の数値)))、
中央値(:数値としての結果値)
;
こんにちはメラニー。
1つの方法があります。データを要約する準備ができたら、目的の変数のログ列を作成します。次に、要約を使用してこのログデータのグループ平均を見つけます。次に、結果を指数化して新しい列に配置します。
これが一例です。
dt = open( "$ Sample_Data / Car Physical Data.jmp");
dt << new列( "log_weight"、式(log(:weight)));
dtSum =(dt << summary(平均(:log_weight)、グループ(:国、タイプ)));
dtsum << new列( "Geom_mean_weight"、式(exp(:name( "Mean(log_weight)")))););
乾杯、
ブレイディ