rpois gives a large number repeatedly (PR#439)
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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._