<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />我试图在大约30到40个监测站中找到大约15年数据的粪便细菌的几何平均值。 我的脚本的目标是进行两个过滤器(一个在项目ID列上,另一个在特征名称列上),然后按地点分组地理位置然后按年份分组,并在最终表格中显示工作站ID,年份,特征名称,N,地理平均值和中位数。
我现在知道我不能在脚本的摘要部分中包含用于查找地理平均值的脚本,但是我无法确定如何在此脚本中放置它。 我不确定我是否需要完全重构脚本或者需要以完全不同的方式解决这个问题。
任何建议或建议将不胜感激!
myPath =“O:/ Planning / BPB / IBP_Pilot_Files / JMPScrips / Geomean /”;
dt =数据表(“Data_dwq_20140502_133555_RegRes_For_Testing.jmp”);
dt <<选择位置(:项目ID ==“NCAMBNT”|:项目ID ==“NCISCO”|:项目ID ==“TPBA”);
dt <<子集((选定行),输出表名(“子集”));
dt =数据表(“子集”);
dt <<选择位置(:特征名称==“粪大肠菌群”|:特征名称==“总大肠菌群”);
dt << Subset((Selected Rows),输出表名(“subset2”));
dt =数据表(“subset2”);
dtsum = dt <<摘要(
组(:站号,:年份,:特征名称),
N(:结果值为数字),
Exp(Col Mean(Log(:结果值为数字))),
中位数(:结果值为数字)
);
嗨,梅兰妮。
这是一种方法:一旦准备好汇总数据,就为所需变量创建一个日志列。然后使用Summary查找此日志数据的组平均值。然后对结果取幂并将它们放在一个新列中。
这是一个例子:
dt = open(“$ Sample_Data / Car Physical Data.jmp”);
dt << new column(“log_weight”,formula(log(:weight)));
dtSum =(dt << summary(mean(:log_weight),group(:Country,:type)));
dtsum << new column(“Geom_mean_weight”,formula(exp(:name(“Mean(log_weight)”))));
干杯,
布雷迪