Your example code isn't consistent with the comments. You say
> ### still I find output values >1, even though yleft=1:
> range(y.out)
[1] 0.0000000 0.9816907
But the range is not bigger than 1. However, when I run the version from your mailing list posting, I get inconsistent results from run to run, so I'll take a look.

(In reply to comment #1)
> Your example code isn't consistent with the comments. You say
>
> > ### still I find output values >1, even though yleft=1:
> > range(y.out)
> [1] 0.0000000 0.9816907
>
> But the range is not bigger than 1. However, when I run the version from your
> mailing list posting, I get inconsistent results from run to run, so I'll take
> a look.
My mistake: I had found a very different results when running the script previously (as mentioned, the results are inconsitent): unfortunately I hadn't seen that this time the range is correct indeed... Sorry for that. In most cases, I get values above 1, though...
Thanks a lot,
Sam