具有正态分布误差的函数值答案

【问题标题】：Function values with normally distributed errors具有正态分布误差的函数值
【发布时间】：2014-03-23 17:02:10
【问题描述】：

假设我想绘制 Michaelis-Menten（或任何其他）函数的一些值。然而，不完全是这些值，而是具有正态分布误差的值。

我想出了以下代码。还有更直接的吗？

f<-function(x,a,b){a*x/(1+b*x)}
datax<-seq(1:10)
datay<-f(datax,2,3)
errors<-rnorm(10,0,0.01)
datayerrors<-datay+errors
plot(datax,datayerrors)

【问题讨论】：

【解决方案1】：

您可以将错误构建到函数中，并且真的不应该使用 seq(1:10) ... seq(1,10) 或 1:10 .... 我曾经遇到过一个令人讨厌的错误。

 f<-function(x,a,b){a*x/(1+b*x) + rnorm(length(x), 0, .01) }
 plot( 1:10, f(1:10) )

【讨论】：