设为首页收藏本站

爱吱声

 找回密码
 注册
搜索
楼主: 晨枫
打印 上一主题 下一主题

[科研心得] 问题:如何从数据里估算普瓦松分布的均值?

[复制链接]
  • TA的每日心情
    慵懒
    6 小时前
  • 签到天数: 2251 天

    [LV.Master]无

    21#
    发表于 2019-2-4 13:21:29 | 只看该作者
    本帖最后由 松叶牡丹 于 2019-2-4 13:23 编辑 : C& E2 q, F# b; O$ C. I: A1 V
    晨枫 发表于 2019-2-4 13:181 h4 Z. {# I2 c+ d3 K
    正是:看数据点然后脑补。人工就是这么干的。但操作工要照管的事情太多,一会儿不看,可能就错过了,就可 ...

    4 _9 J( K& Z" g- J1 w" I7 v7 O. G. |" s6 `
    降低精度,把一系列正态分布离散化为有限几个histogram,每次测试点先第一步进行范围判断,然后归到不同的histogram里面。再反向把histogram变换为对应的大致正态分布,这样可行吗?
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    22#
     楼主| 发表于 2019-2-4 13:23:16 | 只看该作者
    tanis 发表于 2019-2-3 23:123 n) Y; M% Q/ r/ Q
    找温度第二,第三高的两个点,拟合个直线,取中间的那个位置
    ! n2 J$ J0 t% q! I2 y
    塔板是等间距分布的,这样我还不如直接取温度最高的那一块了
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    2023-4-1 00:01
  • 签到天数: 627 天

    [LV.9]渡劫

    23#
    发表于 2019-2-4 13:27:54 | 只看该作者
    晨枫 发表于 2019-2-4 13:23" o2 Q+ A5 Z  z0 A5 a
    塔板是等间距分布的,这样我还不如直接取温度最高的那一块了

      [2 W  a- z6 c. S2 S! R哦哦。。。 那就取最高和第二个的中间位置 有人工干预的时候就用人工,人忘了看就用这个。。。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    24#
     楼主| 发表于 2019-2-4 13:28:01 | 只看该作者
    松叶牡丹 发表于 2019-2-3 23:218 }( }; E$ b) ?* \2 W$ {# ?" _3 E! K
    降低精度,把一系列正态分布离散化为有限几个histogram,每次测试点先第一步进行范围判断,然后归到不同 ...
    ; V  A$ Q9 G  {. [2 a6 q

    6 x  y+ F1 s0 B9 C* ?+ F一把抓住:我的问题正是从histogram变换为某种分布,最想poison分布或者log normal distribution,不行的话正态分布也行。有什么现成的计算方法吗?
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    25#
     楼主| 发表于 2019-2-4 13:29:03 | 只看该作者
    tanis 发表于 2019-2-3 23:278 M6 J  r4 z- }! M% y: K
    哦哦。。。 那就取最高和第二个的中间位置 有人工干预的时候就用人工,人忘了看就用这个。。。 ...

    $ P" N% V1 s! I6 _* `6 T, K这也是个办法,总是有点“心有不甘”……
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    26#
     楼主| 发表于 2019-2-4 13:33:41 | 只看该作者
    松叶牡丹 发表于 2019-2-3 23:21( K; w- G/ Y1 ]* a; n' z8 ~
    降低精度,把一系列正态分布离散化为有限几个histogram,每次测试点先第一步进行范围判断,然后归到不同 ...

    2 c& f" B5 h/ q! Z唉,MATLAB里有histfit命令,干的正是我要的,可惜没法“偷”过来用啊
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    慵懒
    6 小时前
  • 签到天数: 2251 天

    [LV.Master]无

    27#
    发表于 2019-2-4 13:36:06 | 只看该作者
    晨枫 发表于 2019-2-4 13:33
    ! e/ {5 J5 |) B( ^8 a( X唉,MATLAB里有histfit命令,干的正是我要的,可惜没法“偷”过来用啊

    7 H3 t1 B; @7 k8 N" B4 S完全不会FORTRAN,这样啊,要不就得自己写算法了?
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    2023-4-1 00:01
  • 签到天数: 627 天

    [LV.9]渡劫

    28#
    发表于 2019-2-4 13:37:22 | 只看该作者
    晨枫 发表于 2019-2-4 13:33" ?4 A+ F2 Z2 |/ s) a1 d4 Z
    唉,MATLAB里有histfit命令,干的正是我要的,可惜没法“偷”过来用啊

    ! Y7 l4 i5 i. @* m5 r7 v+ b6 PMatlab 直接有poissfit啊。。。
    7 N1 e6 P5 _! o0 ]
    " E  L/ m, O) G7 h1 s我感觉人类脑补的水平估计和取中值也差不了多少了既然连忘了看也只是可能出问题,我脚的中值差不多 )(缺乏工匠精神的我。。。)
      u" H: B7 s  [* R
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    29#
     楼主| 发表于 2019-2-4 13:40:49 | 只看该作者
    tanis 发表于 2019-2-3 23:37
    0 R9 Q5 k, v( Z* X8 [/ UMatlab 直接有poissfit啊。。。/ m) M" z& S( Z% q

    " [4 V% U: [  M$ S! @, @我感觉人类脑补的水平估计和取中值也差不了多少了既然连忘了看也 ...
    5 v4 s, a2 @( M4 p# d% p+ L6 K
    问题是我没法把MATLAB装进DCS里用啊。嗯,如果不行的话,可能也只能内插取中值了。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    慵懒
    6 小时前
  • 签到天数: 2251 天

    [LV.Master]无

    30#
    发表于 2019-2-4 13:41:18 | 只看该作者
    您这个采样数据如果有之前很长时间的大批量数据的话,要不选一部分建模,建立比如10个distribution式子,用剩下数据validate,看看是不是符合六西格玛质量控制要求?如果不行就提高精度,设20,30或者更多式子?我这是笨办法。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    31#
     楼主| 发表于 2019-2-4 13:43:57 | 只看该作者
    松叶牡丹 发表于 2019-2-3 23:41, M  {7 J6 @0 I' j# V  F; ~
    您这个采样数据如果有之前很长时间的大批量数据的话,要不选一部分建模,建立比如10个distribution式子,用 ...
    & l" A6 v' Z7 ?) F8 ~9 C
    这个好像不行?我需要的是实时估计当前的峰值点,历史数据能帮什么忙呢?
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    2018-6-21 09:14
  • 签到天数: 2 天

    [LV.1]炼气

    32#
    发表于 2019-2-4 13:55:23 | 只看该作者
    如果问题是拟合若干点的光滑曲线 找到最大值,kernel density 就好,matlab里的histfit 也是基于这个算法。fortran里面library浩若烟海,查查肯定能找到。速度绝对够快,不需要POISSON 假设。
    6 S+ J- B1 Y3 o5 x# A* @* {+ |+ r8 B) k0 g7 Q6 h" A* U
    试着搜一下,这个fortran code 试试?: https://ideas.repec.org/c/boc/bocode/g00004.html

    评分

    参与人数 1爱元 +2 收起 理由
    tanis + 2 给力

    查看全部评分

    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    慵懒
    6 小时前
  • 签到天数: 2251 天

    [LV.Master]无

    33#
    发表于 2019-2-4 14:07:31 | 只看该作者
    晨枫 发表于 2019-2-4 13:43# r- T' _3 `& b. {7 O  s" N
    这个好像不行?我需要的是实时估计当前的峰值点,历史数据能帮什么忙呢? ...
      D: p) Z! w* J2 [
    历史数据建模自己写算法做验证用。不过如果可以像楼上那位说的可以找到现成library调用,kernel density estimation和 guassian mixture model应该都可以。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    34#
     楼主| 发表于 2019-2-4 14:16:53 | 只看该作者
    gnomegordon 发表于 2019-2-3 23:55
      p+ c6 _8 r* B+ z) h% @如果问题是拟合若干点的光滑曲线 找到最大值,kernel density 就好,matlab里的histfit 也是基于这个算法。 ...

    ( d! K3 P2 u* {/ S十分感谢!打开看到的是SRC文件,还在找FORTRAN……
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    35#
     楼主| 发表于 2019-2-4 14:19:36 | 只看该作者
    松叶牡丹 发表于 2019-2-4 00:07
    " g; D- W6 D6 \历史数据建模自己写算法做验证用。不过如果可以像楼上那位说的可以找到现成library调用,kernel density  ...
    " H1 `' q+ V- ^0 ?7 U
    多谢!这个kernel density以前没有玩过,正在研究中……
    回复 支持 1 反对 0

    使用道具 举报

  • TA的每日心情
    慵懒
    2020-7-26 05:11
  • 签到天数: 1017 天

    [LV.10]大乘

    36#
    发表于 2019-2-4 14:24:18 | 只看该作者
    本帖最后由 Dracula 于 2019-2-4 14:28 编辑
    5 i* N$ l3 h& ^& a- {! y& g4 X) t+ \( o
    说实话,你这个题目我没看懂。
    2 p% T& t' J5 m3 f# p
    - @& w% A% ~- d0 ~$ a4 p! Y但如果你你单纯就是想实现
    ( i6 y; }  M* }. R* I  Q' d% R" b! H4 q
    我就是要拟合一个“钟形曲线”,然后找峰值所在点。相当于上图中红线(或者lambda=2)里横轴1-2之间的位置。

    7 Q* }* p5 U1 a5 }0 V
    3 r3 Z5 u/ p+ u4 J其实应该挺容易。假设这个分布是log normal,如果你知道这个分布的参数, μ和 σ,你要想计算峰值(英语是mode),有现成公式,套公式就可以了(公式你可以上网查,比如Wikipedia上就有)。估计这个分布怎么做,我没看懂你的题目,不好说。但既然正态分布你知道能怎么做,log-normal没什么本质区别,就是公式不一样。简单的使用method of moments应该很容易就能推导出公式。实现就是简单的算术计算,不需要什么数值方法。' {& H5 ~- u' S7 ?$ U$ W4 D" |+ D

      I6 ~, \7 H! |  g2 o, ^4 c
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    37#
     楼主| 发表于 2019-2-4 14:34:17 | 只看该作者
    Dracula 发表于 2019-2-4 00:24
    : T9 L# j) ?' N' ~) W  |说实话,你这个题目我没看懂。
      O+ G2 O* U: [8 @+ e4 k# ^! {* n4 |9 r0 p" `: U8 @! F& v
    但如果你你单纯就是想实现
    $ _6 @& p8 f- W, V
    有了 μ和 σ想计算峰值就容易了,我的问题是如何从histogram计算log normal的 μ和 σ。看来这也是个办法,把数据log一下,再按normal做。我也想到过这个办法,就是不大确定。明天试试,上班的活儿用下班的时间试,可不划不来嘛。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    慵懒
    6 小时前
  • 签到天数: 2251 天

    [LV.Master]无

    38#
    发表于 2019-2-4 14:35:05 | 只看该作者
    Dracula 发表于 2019-2-4 14:24- u8 B2 j, P( R. ]9 m5 y
    说实话,你这个题目我没看懂。
    * I& [6 m: U" o& ^9 r9 E1 u6 m  c" f
    但如果你你单纯就是想实现
    * A/ i2 ?: h+ i+ t& n
    看题目意思应该是不知道参数的。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    慵懒
    6 小时前
  • 签到天数: 2251 天

    [LV.Master]无

    39#
    发表于 2019-2-4 14:36:11 | 只看该作者
    晨枫 发表于 2019-2-4 14:19* X# g5 I, k9 m6 s. c) B* d
    多谢!这个kernel density以前没有玩过,正在研究中……

    " T1 A8 ~1 r: t+ b/ C! {1 {晨大辛苦,您太客气了。祝新年快乐!
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    2018-6-21 09:14
  • 签到天数: 2 天

    [LV.1]炼气

    40#
    发表于 2019-2-4 14:39:24 | 只看该作者
    晨枫 发表于 2019-2-4 14:16
    1 ~, ?0 s! z9 c6 L8 x7 c十分感谢!打开看到的是SRC文件,还在找FORTRAN……
    2 d) [" c" C2 Z- d" h: m$ @
    apologize. 网上搜code太麻烦,还得验证。最好有本书可以翻翻 或者搜library
    回复 支持 反对

    使用道具 举报

    手机版|小黑屋|Archiver|网站错误报告|爱吱声   

    GMT+8, 2025-4-26 14:08 , Processed in 0.040535 second(s), 19 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

    快速回复 返回顶部 返回列表