similar to: dbinom link

In principle a good idea, but I'm not sure the whereabouts of Catherine Loader are known at this point. Last peeps from her on the net seem to be about a decade old. .pd > On 18 May 2020, at 10:31 , Abby Spurdle <spurdle.a at gmail.com> wrote: > > This has come up before. > > Here's the last time: > https://stat.ethz.ch/pipermail/r-devel/2019-March/077478.html

What about using the Wayback Machine archive ? The web archive should be more stable than other links which also might disappear in the future. E.g. https://web.archive.org/web/20070610002602/http://www.herine.net/stat/software/dbinom.html , which also links to an archived copy of the PDF. Best regards, Hilmar On 18.05.20 10:57, peter dalgaard wrote: > In principle a good idea, but

This has come up before. Here's the last time: https://stat.ethz.ch/pipermail/r-devel/2019-March/077478.html I guess my answer to the following the question... Perhaps we should ask permission to nail the thing down somewhere on r-project.org? ...would be, to reproduce it somewhere. And then update the link in the binom help file. Given that the article was previously available

Hi Folks, The recent correspondence about "strange fisher.test result", and especially Peter Dalgaard's reply on Tue 03 April 2007 (which I want to investigate further) led me to take a close look at the code for binom.test(). I now have a query! The code for the two-sided case computes the p-value as follows: if (p == 0) (x == 0) else if (p == 1) (x == n)

Hi all, I am wondering whether there is a small bug in the binom.test function of the ctest library (I'm using R 1.6.0 on windows 2000, but Splus 2000 seems to have the same behaviour). Or perhaps I've misunderstood something. the command binom.test(11,100,p=0.1) and binom.test(9,100,p=0.1) give different p-values (see below). As 9 and 11 are equidistant from 10, the mean of the

I?ve installed R 3.7.0 on a new laptop running macos 10.14.5 and have managed to get most of my usual packages to compile from source with a ~/.R/Makevars file that looks like this: CC=/usr/local/clang8/bin/clang CXX=/usr/local/clang8/bin/clang++ LDFLAGS=-L/usr/local/clang8/lib CPPFLAGS=-I/usr/local/clang8/include -I/opt/X11/include/freetype2 FC=/usr/local/gfortran/bin/gfortran

Hi Folks, There has been past correspondence regarding Catherine Loader's Bell Labs (oops, Lucent) paper "Fast and Accurate Computation of Binomial Probabilities" which gives the algorithm on which R's dbinom() is based. The original URL given in the R documentation "?dbinom" is: http://cm.bell-labs.com/cm/ms/departments/sia/catherine/dbinom but this link is

I believe the binom.test procedure is producing one tailed p values rather than the two tailed value implied by the alternative hypothesis language. A textbook and SAS both show 2*9.94e-07 = 1.988e-06 as the two tailed value. As does the R summation syntax from R below. It looks to me like the alternative hypothesis language should be revised to something like " ... greater than or equal

Hi, With the following Sweave minimal file: ---<--------------------cut here---------------start------------------->--- \documentclass{article} \usepackage{Sweave} \begin{document} <<binom-sim>>= thetas <- seq(0, 1, by=0.001) prior <- rep(1, length(thetas)) / length(thetas) lik <- dbinom(1, 1, thetas) lik.p <- prior * lik post <- lik.p / sum(lik.p)

Full_Name: Uffe H?gsbro Thygesen Version: 2.2.0 OS: linux Submission from: (NULL) (130.226.135.250) Hello all. pbinom(q=0,size=0,prob=0.5) returns the value NaN. I had expected the result 1. In fact any value for q seems to give an NaN. Note that dbinom(x=0,size=0,prob=0.5) returns the value 1. Cheers, Uffe

On 28/06/2019 6:27 a.m., Koenker, Roger W wrote: > Thanks, I was just coming to that conclusion and beginning to look for a way to make a list of > exportable objects. After library(quantreg), ls("package:quantreg") will list all the names you currently export. So cat(ls("package:quantreg"), sep = ", ") will print the list in a form suitable for including

> ... reports that the fortran names of all the registered fortran functions are ?undocumented objects?... Those symbols should not be exported. In quantreg v5.33, NAMESPACE has 'exportPattern(".")'. Maybe it is a good time to remove that and export explicitly the symbols that are meant for export. Georgi Boshnakov ------------------------------ Message: 6 Date: Fri, 28

Hi, I trying to determine the best way to compute the power for a one-sample one-sided binomial test. Specifically I need to sample a population of individuals and ask whether a sample rate of 0% is compatable with a minimum threshold of 3% and how many samples are needed. I have made use of power.prop.test but I am not sure if a) that is the correct (or best) function to use and b) if the

On Mon, Mar 12, 2018 at 2:18 PM, Martin Maechler <maechler at stat.math.ethz.ch> wrote: > [...] > Is that so? Not according to my reading of the 'Writing R > Extensions' manual, nor according to what I have been doing in > all of my packages for ca. 2 years: > > The rule I have in my mind is > > 1) NAMESPACE Import(s|From) \ >

In a private response to Tony Plate's suggestion to replace version() output with sessionInfo() in R-help requests, > roger koenker wrote: >> Thanks for this, it would seem useful to have version numbers for >> the packages too? and Tony replied, > > Sounds sensible to me! If I were you I'd send a message to R-devel > suggesting this. AFAIK, some changes to

Thierry, I have this: if (require(MatrixModels) && require(Matrix)) { X <- model.Matrix(Terms, m, contrasts, sparse = TRUE) in my function rqss() I've tried variants of requireNamespace too without success. If I understand properly model.Matrix is from MatrixModels but it calls sparse.model.matrix which is part of Matrix, and it is the latter function that I'm not

My apologies, the data can now be found at: url <- "http://www.econ.uiuc.edu/~roger/research/ebayes/velo.d" x <- scan(url,skip = 1) If I could get each of the histograms to mimic what is produced by hist(x, 100, freq = FALSE) I?ve experimented with xlim, ylim, without success so far... url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu

I sent the message below to r-sig-mac yesterday, but having no reply I decided to explore a bit myself and found that editing: /Library/Frameworks/R.framework/Versions/2.3/Resources/share/make/ shlib.mk yzzy: diff shlib.mk shlib.mk~ 3c3 < include $(R_HOME)/etc/Makeconf --- > include $(R_HOME)/etc${R_ARCH}/Makeconf restored the functionality of R CMD INSTALL. Is this a known issue?

Dear All, Just when I thought I had the plague of gfortran-9 under control, I made the tactical error of allowing my mac mini to ?upgrade? to macOS 10.14.6 which apparently also upgraded Xcode to 10.3. In consequence I?m having difficulty building my packages. The current symptom is: /usr/local/clang7/bin/clang -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG

I am trying to simulate a series of ones and zeros (1 or 0) and I am using "rbinom" but realizing that the number of successes expected is not accurate. Any advice out there. This is the example: N<-500 status<-rbinom(N, 1, prob = 0.15) count<-sum(status) 15 percent of 500 should be 75 but what I obtain from the "count" variable is 77 that gives the probability of