重复分析<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />大家好,
我想根据附带的示例处理表中包含的数据。该表包含100列+一列作为控制措施,并包含30到80行。
目的是以迭代的方式应用统计测试,从对照与列51,对照与列52,对照与列53等开始,以便检测何时出现显着差异,在连续列上存在多长时间以及当意义丧失。
例如,可以从列63到列78检测差异,然后从c79到c86丢失,然后再次出现在c87到c98上。
我不知道如何编写一个脚本,允许自动分析这种类型的表,给出一个报告差异发生的时刻。
有什么建议吗?
Re:重复分析这是交换建议的建设性方式!谢谢你,凯文。有趣的是,我们可以与R和JMP来回走动。我一定会仔细研究一下变更点程序。
Re:重复分析嗨,保罗!
没有更多的信息,我所做的不仅仅是在黑暗中拍摄。
但是,如果我怀疑您正在搜索顺序收集的数据以试图辨别数据生成过程的变化,那么变换点方法可能更合理。
有许多方法可以检测变化点,并且有许多参考资料可以追溯到很多年。
CRAN有一个变点在R中实现几个最近研究过的方法的包.JMP有一些简洁的R集成,你可以在JMP数据集上执行R代码并将结果返回到JMP中。 尝试R.Killick,P.Fearnhead&IAEckley(2012)Optimal Detection of Changepoints with Linear Computational Cost,Journal of the American Statistical Association,107:500,1590-1598中引用的Pruned Exact Linear Time(PELT)方法。
Re:重复分析谢谢凯文,任何建议表示赞赏。
根据我的理解,保罗没有选择任何测试,但更多的是要求如何在JMP中运行一系列双样本测试(一个样本总是相同),给出他的原始数据集。我的观点是在这种情况下不需要编写脚本。
再次感谢您的贡献。
Re:重复分析谢谢你的留言Kevin,
我不想打扰杰出的统计学家,无论他们是活着还是死了。
一世将会满意至知道你的建议是什么。
Re:重复分析有人应该提一下,你们所采取的方法可能会导致相当数量的死去的杰出统治者像他们各自的坟墓中的车床一样旋转。
我无法准确地看出你想从你的描述中做出什么决定,但这听起来更像是对我的变换点检测。 如果是这样,从任何双样本测试中分析许多顺序聚合的p值可能被认为是最好的形式。
Re:重复分析嗨保罗,
谢谢,这帮了很多忙。
似乎脚本在这里不是强制性的。您可以执行以下操作:
1。表>堆栈:获取一个新数据表,其中包含来自控制列的数据和堆叠的所有测量列。这将创建一个包含列名称的列标签,以及包含所有数据的第二列。
附件显示我的结果。
2。分析>适合Y乘X,X =标签,Y =数据>热点>比较均值>使用控制,Dunnett:选择您的控制列,您将获得所有成对测试。您可以右键单击LSD阈值矩阵下的表,例如对其进行排序或将其转换为数据表(可能用于在Graph Builder中绘制p值图形)。
对于Dunnett的测试也见比较手段 。
希望有所帮助,
沃尔克
Re:重复分析嗨Volker,谢谢你照顾我的问题,
这是一个已经分析过的数据表,其中Ctrl / test比较显示出显着差异(在测试中评估,但Wilcoxon更合适)。
下面的每一行给出了Ctrl数据(控制列)和给定测试列(编号从100到200)之间的比较
Col nbr P值
...
160 0.01
161 0,0022
162 ns
163 ns
164 0,03
165 ns
166 0,008
167 0,04
168 0,04
169 ns
170 0,004
171 ns
172 ns
173 ns
174 0,04
175 ns
176 <0,0001
177 0,008
178 ns
179 0,004
180 0,02
181 0,003
182 0,004
183 0,009
184 0,04
185 0,005
186 ns
187 0,006
188 0,03
189 0,02
190 0,0002
191 ns
192 0,001
等等...
Re:重复分析嗨保罗,
这是我从聊天中学到的东西:
差异的定义:Wilcoxon测试将Ctrl列中的值与其他100列中的值进行比较(重复比较)。
示例表仅用于说明数据结构。上述行为不适用于此表。
目的是对每种分析表的差异的起始和偏移时间进行排序。
如果您可以提供另一个表格以及一些分析结果,您希望看到它们会有很大帮助。
BR,Volker