On Sat, 25 Apr 2020 19:02:59 +0000 (UTC) varin sacha via R-help <r-help@r-project.org> wrote:
> plot(function(d) > k*exp(-0.5*(d-m)^2/sd^2),col=2,add=TRUE,xlim=range(d)) Remove the `add=TRUE` or calculate the values manually and you'll see that the function value is zero because the values in `d` are too big compared to `m` (where the function peak is). Compare your result to the following: plot( function(d) k*exp(-0.5*(d-m)^2/sd^2), xlim = m + 3*sd * c(-1,+1) ) Your parameters don't quite fit the observed values. -- Best regards, Ivan ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.