具有正态分布错误的函数值

Question

我想说我想绘制一些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) 

Answer 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) )