微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

R中两个密度混合的卷积

如何解决R中两个密度混合的卷积

我目前正在尝试使用软件包“ distr”在R中的两个对数正态树突的混合物之间进行卷积。混合物由下式给出:

library(distr)
M1 <- UnivarMixingdistribution(Lnorm(meanlog=8.1,sdlog=0.02),Lnorm(meanlog=8.15,sdlog=0.035),mixCoeff = c(0.45,0.55))

M2 <- UnivarMixingdistribution(Lnorm(meanlog=8.17,sdlog=0.03),mixCoeff = c(0.65,0.35))

这两个函数应该在X = [3000,4000]范围内提供正概率。但是,当我进行卷积时:

conv <- convpow(M1+M2,1)
M3  <- d(conv) 
M3(seq(3000,4000,by=1))

我得到的概率都是零,实际上PDF位于7000区域附近。有人可以帮我正确地进行卷积吗?

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。