使用R去离子软件包进行反卷积的问题

问题描述

我在处理R中的反卷积问题上特别困难。

我有一套以磅为单位的重量，该重量被测量误差污染了。有72个数据点。我也有用于测量产品的比例的测量标准偏差。我想恢复对真实未污染分布的估计。

我假设测量误差呈正态分布，并且适用附加误差模型。

我正在尝试使用以下代码对数据进行反卷积：

Weight <- c(1.4760,1.4740,1.4700,1.4760,1.4720,1.4746,1.4710,1.4730,1.4705,1.4756,1.4770,1.4752,1.4735,1.4743,1.4747,1.4736,1.4732,1.4731,1.4709,1.4754,1.4707,1.4708,1.4753,1.4692,1.4764,1.4748,1.4749,1.4745,1.4671,1.4671)

plot(density(weight))

# Producing this image of the density

[Original Data Density Plot][1]


# Now Generate the noise vector

stddev <- 0.00000606
mn <- 0

noise <- rnorm(n = 72,mean = mn,sd = stddev)

# and finally perform the deconvolution

ST <- deamerSE(y = Weight,error = noise,grid.length = 100)

# Then view the plotted results

plot(density(Weight))
lines(ST)

[Original plus deconvoluted density][2]

 #####

即使测量的标准偏差很小，我也无法得到解卷积的密度看起来像原始密度。

任何帮助将不胜感激。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

deconvolution r r