- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
大家好,我从酵母菌落板的图像中获得了一些强度值。我需要能够从强度值中找到峰值。下面是一个示例图像,显示了绘制时值的外观。
一些值的示例
5.75.38.216.534.258.8**75.4**7565.962.658.666.471.453.540.526.814.28.65.97.714.930.549.969.1**75.3**69.858.857.256.367.16945.127.613.485
这些值在 75.4 和 75.3 处显示了两个峰值,您可以看到值先增大后减小。变化并不总是一样的。
强度值图表
http://lh4.ggpht.com/_aEDyS6ECO8s/THKTLgDPhaI/AAAAAAAAAio/HQW7Ut-HBhA/s400/peaks.png来自 research我正在考虑做的一件事是将每个组(即山脉)存储在哈希中,然后查找组中的最大值。如果我看到的问题是如何确定每个组的边界。
这是我目前所拥有的代码的链接: http://paste-it.net/public/y485822/
这是一个完整数据集的链接: http://paste-it.net/public/ub121b4/
我正在用 Perl 编写代码。任何帮助将不胜感激。谢谢
最佳答案
您需要确定您希望山峰的局部性。此处的方法在数据的广泛区域内找到峰值和谷值。
use strict;
use warnings;
my @data = (
5.7, 5.3, 8.2, 16.5, 34.2, 58.8, 75.4, 75, 65.9, 62.6,
58.6, 66.4, 71.4, 53.5, 40.5, 26.8, 14.2, 8.6, 5.9, 7.7,
14.9, 30.5, 49.9, 69.1, 75.3, 69.8, 58.8, 57.2, 56.3, 67.1,
69, 45.1, 27.6, 13.4, 8, 5,
);
# Determine mean. Or use Statistics::Descriptive.
my $sum;
$sum += $_ for @data;
my $mean = $sum / @data;
# Make a pass over the data to find contiguous runs of values
# that are either less than or greater than the mean. Also
# keep track of the mins and maxes within those groups.
my $group = -1;
my $gt_mean_prev = '';
my @mins_maxs;
my $i = -1;
for my $d (@data){
$i ++;
my $gt_mean = $d > $mean ? 1 : 0;
unless ($gt_mean eq $gt_mean_prev){
$gt_mean_prev = $gt_mean;
$group ++;
$mins_maxs[$group] = $d;
}
if ($gt_mean){
$mins_maxs[$group] = $d if $d > $mins_maxs[$group];
}
else {
$mins_maxs[$group] = $d if $d < $mins_maxs[$group];
}
$d = {
i => $i,
val => $d,
group => $group,
gt_mean => $gt_mean,
};
}
# A fun picture.
for my $d (@data){
printf
"%6.1f %2d %1s %1d %3s %s\n",
$d->{val},
$d->{i},
$d->{gt_mean} ? '+' : '-',
$d->{group},
$d->{val} == $mins_maxs[$d->{group}] ? '==>' : '',
'.' x ($d->{val} / 2),
;
}
输出:
5.7 0 - 0 ..
5.3 1 - 0 ==> ..
8.2 2 - 0 ....
16.5 3 - 0 ........
34.2 4 - 0 .................
58.8 5 + 1 .............................
75.4 6 + 1 ==> .....................................
75.0 7 + 1 .....................................
65.9 8 + 1 ................................
62.6 9 + 1 ...............................
58.6 10 + 1 .............................
66.4 11 + 1 .................................
71.4 12 + 1 ...................................
53.5 13 + 1 ..........................
40.5 14 - 2 ....................
26.8 15 - 2 .............
14.2 16 - 2 .......
8.6 17 - 2 ....
5.9 18 - 2 ==> ..
7.7 19 - 2 ...
14.9 20 - 2 .......
30.5 21 - 2 ...............
49.9 22 + 3 ........................
69.1 23 + 3 ..................................
75.3 24 + 3 ==> .....................................
69.8 25 + 3 ..................................
58.8 26 + 3 .............................
57.2 27 + 3 ............................
56.3 28 + 3 ............................
67.1 29 + 3 .................................
69.0 30 + 3 ..................................
45.1 31 + 3 ......................
27.6 32 - 4 .............
13.4 33 - 4 ......
8.0 34 - 4 ....
5.0 35 - 4 ==> ..
关于perl - 在 Perl 中需要峰值信号检测方面的帮助,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/3549205/
我对构面有疑问,并根据构面进行了一些过滤。 我知道这是一个重复的问题,但我找不到答案。 我想知道如何在 flex 搜索中实现相同的功能。 假设我有一个有关汽车和某些方面的索引-例如模型和 颜色。 颜色
我正在尝试找到一种解决方案来为某些方面创建子方面列表。 我有一些产品的衣服尺码,它们存储在 solr 中 "Size_both":"W30L30","尺寸宽度":"W30","Size_length"
我正在尝试找到一种解决方案来为某些方面创建子方面列表。 我有一些产品的衣服尺码,它们存储在 solr 中 "Size_both":"W30L30","尺寸宽度":"W30","Size_length"
我对方面有疑问。他们不开火。我有小方面: @Aspect @Component public class SynchronizingAspect { @Pointcut("execution(
这是在 ruby 中启用散列自动生成的巧妙技巧(取自 facets): # File lib/core/facets/hash/autonew.rb, line 19 def self.a
这个问题在这里已经有了答案: 8年前关闭。 Possible Duplicate: Creating a facet_wrap plot with ggplot2 with different ann
XMLHttpRequest 能否从 http://mydomain.example/ 向 http://mydomain.example:81/ 发送请求? 最佳答案 要使两个文档被视为具有相同的来
我对 Elasticsearch 中的方面有一点问题。 我有一个表格视频,一个表格 channel ,一个 channel 有很多视频。 我只想在 X 个最新视频上显示每个 channel 的 %vi
假设我正在为 4 个人绘制数据图表:Alice、Bob、Chuck 和 Dana。我正在使用 ggplot2 制作一个多面图,每个人一个方面。我的磁盘上还有 4 张图像:Alice.png、Bob.p
我已经下载了收件箱,并且正在使用Pig和Hadoop处理电子邮件。我已经使用Pig和Wonderdog在ElasticSearch中为这些电子邮件编制了索引。 现在,我为收件箱中的每个电子邮件地址创建
我有一个模块如下: define([...], function(...){ function anothermethod() {...} function request() {....}
(defprotocol IAnimal "IAnimal" (report [o] (println (type o) " reporting.\n") (inner-repor
我有一个 Bean 需要向 InfluxDB 报告。数据库在表 INFLUX_DB_SERVER 中注册了 InfluxDB。如果你看一下代码,你会发现方法reportMemory做了很多工作,它构造
我的问题与分面有关。在下面的示例代码中,我查看了一些分面散点图,然后尝试在每个分面的基础上叠加信息(在本例中为平均线)。 tl;dr 版本是我的尝试失败了。要么我添加的平均线计算所有数据(不尊重方面变
假设我正在为 4 个人绘制数据图表:Alice、Bob、Chuck 和 Dana。我正在使用 ggplot2 制作一个多面图,每个人一个方面。我的磁盘上还有 4 张图像:Alice.png、Bob.p
尝试用两个方面包装服务类来获取此调用链: javanica..HystrixCommandAspect -> MyCustomAroundAspect -> MyService 遇到两个问题: Hys
我是 AspectJ 的初学者。我用它在我的网络驱动程序中截取屏幕截图。以下是我的包结构。 我想知道如何在 Browser 类中运行我的程序,以便它使用 Screenshots 类中定义的 Aspec
我在使用 spring aop 时遇到问题 (编辑:如果我的方法不是静态的,则代码可以正常工作) 我的包中有这个结构: aaa.bbb.ccc.Clase1.java aaa.bbb.ddd.Clas
我有一个通用存储库类,其中包含各种标记有 PostSharp 方面 (SecuredOperation) 的方法... public class Repository : IRepository, I
我有一个运行多线程的 Hibernate 事务方法“doImportImpl”。而某些记录需要依次导入,所以代码结构大致是这样的: public RecordResult doImportImpl(S
我是一名优秀的程序员,十分优秀!