[Rd] rpois gives a large number repeatedly (PR#439)

[Kjetil Kjernsmo]

  This message is in MIME format.  The first part should be readable text,
  while the remaining parts are likely unreadable without MIME-aware tools.
  Send mail to mime@docserver.cac.washington.edu for more info.

--0-759262972-950616086=:2636
Content-Type: TEXT/PLAIN; charset=US-ASCII

On 15 Feb 2000, Peter Dalgaard BSA wrote:

>kjetikj@astro.uio.no writes:
>
>> Full_Name: Kjetil Kjernsmo
>> Version: 0.65.1
>           !!!!!!

Yep, it's old.... When v1.0.0 comes, I guess I'll get the UNIX people to
install it.

>I can't seem to generate anything unusual with the latest R snapshot.

It is really strange, this. I can do 
> max(rpois(20000000,3))
[1] 16

>Perhaps it would help if you told us what parameters to call
>photonsincidenttodetectorbin with?

Whoops, sorry, I use 
x <- cbind(-50:50)
apply(x, 1, photonsincidenttodetectorbin, 100000, 3)

Anyway, coming back from lunch, I was unable to reproduce it with this
call, after trying about 50 times... This call made the problem show up
almost every time half an hour ago. I bet there is someone sitting inside
the computer having a good time right now... :-) However, I was able to
reproduce it once in about 10 trails by calling 
plot(lensedprofile(100000,3,1), type='l')
which is a function defined by the rest of my code, that I have attached. 

Best,

Kjetil
-- 
Kjetil Kjernsmo
Graduate astronomy-student                    Problems worthy of attack
University of Oslo, Norway            Prove their worth by hitting back
E-mail: kjetikj@astro.uio.no                                - Piet Hein
Homepage <URL:http://www.astro.uio.no/~kjetikj/>
Webmaster@skepsis.no 

--0-759262972-950616086=:2636
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="first.R"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.OSF.4.05.10002151301260.2636@rukbat.uio.no>
Content-Description: Full code
Content-Disposition: attachment; filename="first.R"

bmNsb3Vkc2JpbiA8LSBmdW5jdGlvbihiaW5ubywgbnRvdGFsY2xvdWRzLCBu
dW1iZXJvZmJpbnMgPSAxMDAsIGxpbmV3aWR0aCA9IG51bWJlcm9mYmlucyAv
IDYpDQogIHJldHVybihudG90YWxjbG91ZHMgKiBkbm9ybShiaW5ubywgc2Qg
PSBsaW5ld2lkdGgpKQ0KDQpwaG90b25zaW5jaWRlbnR0b2RldGVjdG9yYmlu
IDwtIGZ1bmN0aW9uKGJpbm5vLCBudG90YWxjbG91ZHMsIGF2Z3Bob3RvbmZy
b21jbG91ZCwgYW1wbW9kZSA9IDEsIGFtcG1lYW4gPSAyLCBudW1iZXJvZmJp
bnMgPSAxMDAsIGxpbmV3aWR0aCA9IG51bWJlcm9mYmlucyAvIDYpDQp7DQog
IHNjIDwtIGFtcG1lYW4gLSBhbXBtb2RlDQogIGlmKHNjIDw9IDApIHN0b3Ao
Ik1lYW4gQW1wbGlmaWNhdGlvbiBtdXN0IGJlIGdyZWF0ZXIgdGhhbiBNb2Rl
IEFtcGxpZmljYXRpb24iKSANCiAgbmNiIDwtIG5jbG91ZHNiaW4oYmlubm8s
IG50b3RhbGNsb3VkcywgbnVtYmVyb2ZiaW5zLCBsaW5ld2lkdGgpDQogIHJn
IDwtIHJnYW1tYShuY2IsIGFtcG1lYW4gLyBzYywgc2MpDQogIHJwIDwtIHJw
b2lzKG5jYiwgYXZncGhvdG9uZnJvbWNsb3VkKQ0KICBkciA8LSByZyAqIHJw
DQogIGlmKHN1bShkcikgPiAyNTAwMCkgcHJpbnQoY2JpbmQocmcsIHJwLCBk
ciwgYXZncGhvdG9uZnJvbWNsb3VkLCBuY2IpKQ0KICByZXR1cm4oc3VtKGRy
KSkNCn0NCg0KDQoNCnBob3RvbnNkZXRlY3RlZGJpbiA8LSBmdW5jdGlvbihi
aW5ubywgbnRvdGFsY2xvdWRzLCBhdmdwaG90b25mcm9tY2xvdWQsIHNlbnNp
dGl2aXR5LCBhbXBtb2RlID0gMSwgYW1wbWVhbiA9IDIsIG51bWJlcm9mYmlu
cyA9IDEwMCwgbGluZXdpZHRoID0gbnVtYmVyb2ZiaW5zIC8gNikNCiAgew0K
ICAgIHBkYiA8LSBzZW5zaXRpdml0eSAqIHBob3RvbnNpbmNpZGVudHRvZGV0
ZWN0b3JiaW4oYmlubm8sIG50b3RhbGNsb3VkcywgYXZncGhvdG9uZnJvbWNs
b3VkLCBhbXBtb2RlLCBhbXBtZWFuLCBudW1iZXJvZmJpbnMsIGxpbmV3aWR0
aCkNCiAgIHByaW50KHBkYikNCiAgICAjIFBhdGNoZXMgYnVnIGluIFIgdmVy
c2lvbiA8IDEuMC4wDQogICAgcmV0dXJuKGlmZWxzZShwZGIgPD0gMCwgMCwg
cnBvaXMoMSwgcGRiKSkpDQoNCiAgfQ0KbGVuc2VkcHJvZmlsZSA8LSBmdW5j
dGlvbihudG90YWxjbG91ZHMsIGF2Z3Bob3RvbmZyb21jbG91ZCwgc2Vuc2l0
aXZpdHksIGFtcG1vZGUgPSAxLCBhbXBtZWFuID0gMiwgbnVtYmVyb2ZiaW5z
ID0gMTAwLCBsaW5ld2lkdGggPSBudW1iZXJvZmJpbnMgLyA2KQ0KICByZXR1
cm4oYXBwbHkoY2JpbmQoLWFzLmludGVnZXIobnVtYmVyb2ZiaW5zLzIpOmFz
LmludGVnZXIobnVtYmVyb2ZiaW5zLzIpKSwgMSwgcGhvdG9uc2RldGVjdGVk
YmluLCBudG90YWxjbG91ZHMsIGF2Z3Bob3RvbmZyb21jbG91ZCwgc2Vuc2l0
aXZpdHksIGFtcG1vZGUsIGFtcG1lYW4sIG51bWJlcm9mYmlucywgbGluZXdp
ZHRoKSkNCg==
--0-759262972-950616086=:2636--
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._