如何解决使用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 举报,一经查实,本站将立刻删除。