similar to: Faster version of rgeom()

Dear all, I''m working with the geometric distribution for the time being, and I''m confused. This may have more to do with statistics than R itself, but since I''m getting results from R I find counterintuitive (well, yeah, my statistical intuition has not been properly sharpened), I feel like asking. The point first: If I do > rgeom(1,prob=1) I get: [1] NaN Warning

Dear all, Using rmultinom() in a stochastic model, I found this function returns an error message 'NA in probability' for an infinite probability. Maybe, a more precise message will be helpful when debugging. > rmultinom(1, 3:5, c(1/2, 1/3, Inf)) Error in rmultinom(1, 3:5, c(1/2, 1/3, Inf)) : NA in probability vector > rmultinom(1, 3:5, c(1/2, 1/3, NA)) Error in rmultinom(1,

>>>>> Christophe Dutang >>>>> on Sat, 8 Apr 2023 14:21:53 +0200 writes: > Dear all, > Using rmultinom() in a stochastic model, I found this function returns an error message 'NA in probability' for an infinite probability. > Maybe, a more precise message will be helpful when debugging. >> rmultinom(1, 3:5, c(1/2, 1/3,

On 08/04/2023 5:53 p.m., Martin Maechler wrote: >>>>>> Christophe Dutang >>>>>> on Sat, 8 Apr 2023 14:21:53 +0200 writes: > > > Dear all, > > > Using rmultinom() in a stochastic model, I found this function returns an error message 'NA in probability' for an infinite probability. > > > Maybe, a more

>>>>> Benjamin Tyner >>>>> on Mon, 20 Jan 2020 08:10:49 -0500 writes: > On 1/20/20 4:26 AM, Martin Maechler wrote: >> Coming late here -- after enjoying a proper weekend ;-) -- >> I have been agreeing (with Spencer, IIUC) on this for a long >> time (~ 3 yrs, or more?), namely that I've come to see it as a >>

>>>>> Ben Bolker >>>>> on Mon, 20 Jan 2020 12:54:52 -0500 writes: > Ugh, sounds like competing priorities. indeed. > * maintain type consistency > * minimize storage (= current version, since 3.0.0) > * maximize utility for large lambda (= proposed change) > * keep user interface, and code, simple (e.g., it would be easy

On 2020-01-22 02:54, Martin Maechler wrote: >>>>>> Martin Maechler >>>>>> on Tue, 21 Jan 2020 09:25:19 +0100 writes: >>>>>> Ben Bolker >>>>>> on Mon, 20 Jan 2020 12:54:52 -0500 writes: > >> Ugh, sounds like competing priorities. > > > indeed. > > >> * maintain type

On 1/20/20 4:26 AM, Martin Maechler wrote: > Coming late here -- after enjoying a proper weekend ;-) -- > I have been agreeing (with Spencer, IIUC) on this for a long > time (~ 3 yrs, or more?), namely that I've come to see it as a > "design bug" that rpois() {and similar} must return return typeof() "integer". > > More strongly, I'm actually pretty

Ugh, sounds like competing priorities. * maintain type consistency * minimize storage (= current version, since 3.0.0) * maximize utility for large lambda (= proposed change) * keep user interface, and code, simple (e.g., it would be easy enough to add a switch that provided user control of int vs double return value) * backward compatibility On 2020-01-20 12:33 p.m., Martin Maechler

>>>>> Martin Maechler >>>>> on Tue, 21 Jan 2020 09:25:19 +0100 writes: >>>>> Ben Bolker >>>>> on Mon, 20 Jan 2020 12:54:52 -0500 writes: >> Ugh, sounds like competing priorities. > indeed. >> * maintain type consistency >> * minimize storage (= current version, since 3.0.0) >> *

Hi, I am trying to use the Fisher scoring method with a geometric distribution, with p = .07, 100 observations from the geom distrib, and 10 iterations. I cannot quite get the code to work. Can anyone see the mistake? n <- 100 p <- 0.07 x <- rgeom(n, p) s <- sum(x) f <- function(x, p) p*(1-p)^x L <- function(p) p^n*(1-p)^s logL <- function(p)

On 1/20/20 12:33 PM, Martin Maechler wrote: > > It's really something that should be discussed (possibly not > here, .. but then I've started it here ...). > > The NEWS for R 3.0.0 contain (in NEW FEATURES) : > > * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), > rsignrank() and rwilcox() now return integer (not double) > vectors.

Hi there, I am trying to do multiple pairwise Wilcoxon signed rank tests in a manner similar to: a <- c(runif(1000, min=1,max=50), rnorm(1000, 50), rnorm(1000, 49.9, 0.5), rgeom(1000, 0.5)) b <- c(rep("group_a", 1000), rep("group_b", 1000), rep("group_c", 1000), rep("group_d", 1000)) pairwise.wilcox.test(a, b, alternative="two.sided",

Hi! I have just downloaded the R-0.60.1 sources and have problems compiling R on a Sun Ultra 1 running Solaris 2.6 and gcc 2.7.2.1. I have not been able to find to find any compiling hints in the documentation or the FAQ. After ./configure I use make and get the output below. Any hints are welcome. I am not on the list, so please answer me directly too. Best regards Jens --- Jens Lund

On a recent FreeBSD 8.0-CURRENT (i386) building R (any version) breaks with the following messages: ---------------------------------------------------------------------- [...snip...] gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -g -O2 -c wilcox.c -o wilcox.o gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include

Full_Name: Madeleine Yeh Version: 1.9.1 OS: AIX 5.2 Submission from: (NULL) (151.121.225.1) After compiling R-1.9.1 on AIX 5.2 using the IBM cc compiler, I ran the checks. One of them failed. Here is the output from running the check solo. root@svweb:/fsapps/test/build/R/1.9.1/R-1.9.1/tests/Examples: ># ../../bin/R --vanilla < stats-Ex.R R : Copyright 2004, The R

Full_Name: Kjetil Kjernsmo Version: Version 1.0.0 OS: osf1 Submission from: (NULL) (129.240.28.227) Hello again! This is a follow-up to my message on R-help about a problem with match.arg() I have a little more on the topic, but not much really. First, it was not entirely true what I wrote in r-help that I call match.arg() in my ramp(), because ramp() just calls qamp() with runif() as

Dear helpers, thanks. Installation of R and packages (like integrate) successful done. The LD_LIBRARY_PATH wasn't set in the correct order: system libs ahead of the compiler lib. Cross checks provided by Brian Ripley, Oxford, helped a lot. Yours Bodo gordon.harrington at uni.edu wrote: > Without tracing the loader system calls and knowing where the functions are > located I can

Hi again ! Thanks for the info on updating the config.site file which I have done. I have also added -lm in the Makeconf manually because this is needed explicitly for DEC cc. However, there are still a few problems when linking some of the files as you can see from the enclosed log. Ciao, Andreas ------------------------------------------------------- R-0.63.3>make make[1]: Entering

Hi, Just downloaded the R-0.62.4 of R and tried to install it. With the standard procedure : ./configure make At the end of the compilation I got the following message : ld: ../lib/libunix.a(system.o): main: multiply defined fort: Severe: Failed while trying to link. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. I attach the printout after the ./configure and make