问题描述
我目前正在尝试使用软件包“ 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 (将#修改为@)