similar to: proposed pbirthday fix

Full_Name: Andy Lynch Version: 1.9.1 OS: Windows Submission from: (NULL) (131.111.86.211) As I understand it, pbirthday(n,c,k) gives the approximate probability that we see a class with k coicident people in it when n people are sorted into c classes. so the command > pbirthday(4,classes=3,coincident=4) should give the approximate probability that when four people fall into three

Hi, pbirthday(, coincident = 2) starts to issue warnings (see (*) below) for larger number of classes (R 4.0.0, R-devel ./src/library/stats/R/birthday.R:47). The default coincident = 2 is computed as 1 - prod((c:(c - n + 1))/rep(c, n)) where c = classes. Using exp(log(...)), one can derive the return value if(n > 0) 1 - exp(sum(log1p(-(0:(n-1))/c))) else 0. Simplifying this a bit further one

I'm looking at algorithms for determining coincidences. In educational testing, it is interesting to look at cheating via the birthday problem where I can assess the probability of n students having the same test score in a class of size k. I was writing my own code for the b-day problem until I ran into the qbirthday() function, which has solutions for the overflow problems I kept running

I've rolled up R-2.13.1.tar.gz a short while ago. This is a maintenance release to consolidate various minor fixes to 2.13.0. The Licencing change to GPL-2 | GPL-3 may be of interest to some. You can get it from http://cran.r-project.org/src/base/R-2/R-2.13.1.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course.

I've rolled up R-2.13.1.tar.gz a short while ago. This is a maintenance release to consolidate various minor fixes to 2.13.0. The Licencing change to GPL-2 | GPL-3 may be of interest to some. You can get it from http://cran.r-project.org/src/base/R-2/R-2.13.1.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course.

Hi, I tried to calculate the formula for the birthday problem (the probability that at least two people out of a group of n people share the same birthday) But the factorial-function allows me only to calculate factorials up to 170. So is there a way to push that limit? to solve this formula: (factorial(365) / factorial((365-23))) / (365^23) (n=23)

Does anyone know of another version of the Mcnemar test that provides: 1. Odds Ratios 2. 95% Confidence intervals of the Odds Ratios 3. Sample probability 4. 95% Confidence intervals of the sample probability Obviously the Odds Ratios and Sample probabilities are easy to calculate from the contingency table, but I would appreciate any help on how to calculate the confidence

> De : Bill Venables <venables at acland.qld.cmis.csiro.au> > Objet : Re: [R] paste ? > Date?: mardi 25 avril 2000 08:45 > (...) > Secondly, I'm curious about the history of this kind of plot. > I've only heard it called a "Tukey mean difference" plot, (and > Trellis graphics has a function, tmd(), that does it, but no one > knows about it...).

Dear R-list members, I have a problem with translating a mathematica script into R. The whole script is at the end of the email (with initial values for easy reproduction) and can be pasted directly into R. The problematic part (which is included below of course) is <--- Original Mathematica ---> (* p_svbar *) UiA = Ni (Dsi - 2Di A + A^2)/2; UiiA = Nii (Dsii - 2Dii A + A^2)/2; psvbar =

ear ALL, The last part of my thesis analysis is the cross validation. Right now I am having difficulty using the cross validation of gstat. Below are my commands with the tsport_ace as the variable: nfold <- 3 part <- sample(1:nfold, 69, replace = TRUE) sel <- (part != 1) m.model <- x2[sel, ] m.valid <- x2[-sel, ] t<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0)) cv69

Ralf I don't doubt this is expected with the current implementation, I doubt the implementation is desirable. Suggesting to turn this to pbirthday(1e6, classes = 2^53) ## [1] 5.550956e-05 (which is still non-zero, but much less likely to cause confusion.) Best regards Kirill On 26.02.19 10:18, Ralf Stubner wrote: > Kirill, > > I think some level of collision is actually

Hello, I have placed a package on CRAN about two surveys of mortality in Iraq that were published in the Lancet. http://cran.at.r-project.org/src/contrib/Descriptions/lancet.iraqmortality.html > install.packages("lancet.iraqmortality") ... > library(lancet.iraqmortality) Loading required package: foreign > ?lancet.iraqmortality > vignette("mortality") This is

Hello, I have placed a package on CRAN about two surveys of mortality in Iraq that were published in the Lancet. http://cran.at.r-project.org/src/contrib/Descriptions/lancet.iraqmortality.html > install.packages("lancet.iraqmortality") ... > library(lancet.iraqmortality) Loading required package: foreign > ?lancet.iraqmortality > vignette("mortality") This is

I came across a post by Karl Knoblick regarding the modeling of longitudinal data (see https://stat.ethz.ch/pipermail/r-help/2007-May/132137.html). I am often asked by physicians to perform what Karl refers to in his post as option 1: to perform paired t-tests against baseline at each follow up time point (30 days, 90 days, 6 months, etc.). Unlike Karl's example, however, many of the trials

I get ties in output from runif() when I generate as few as 10^5 variates and get quite a lot when I generate 10^6. Is this expected?? I haven't seen any duplication with rnorm(10^6), but see varying amounts of duplication using rexp(), rbeta() and rgamma(). I would have thought that there'd be enough precision that one wouldn't get ties until generating samples larger than this..

I am running 1.3.1 on a Windows (NT 4.0) machine. I've fit a nonlinear model intended to predict crop yield from nutrient information, and want to use the profile function. If I type say, profile(simparj.fm) I get the following error message: "Error in prof$getProfile(): number of iterations exceeded maximum of 5.25515e-308" I used the profiler function to profile simparj,fm step

I am trying to implement a simple r-svm example using the iris data (only two of the classes are taken and data is within the code). I am running into some errors. I am not an expert on svm's. If any one has used it, I would appreciate their help. I am appending the code below. Thanks../Murli ####################################################### ### R-code for R-SVM ### use leave-one-out

Gabe As mentioned on Twitter, I think the following behavior should be fixed as part of the upcoming changes: R.version.string ## [1] "R Under development (unstable) (2019-02-25 r76160)" .Machine$double.digits ## [1] 53 set.seed(123) RNGkind() ## [1] "Mersenne-Twister" "Inversion"??????? "Rejection" length(table(runif(1e6))) ## [1] 999863 I don't

I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance release containing mainly bugfixes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz (give it some time to arrive there) or wait for it to be mirrored at a CRAN site nearer to you. If you're *really* impatient, www.biostat.ku.dk/~pd/R-release should work

I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance release containing mainly bugfixes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz (give it some time to arrive there) or wait for it to be mirrored at a CRAN site nearer to you. If you're *really* impatient, www.biostat.ku.dk/~pd/R-release should work