search for: psmirnov

HI, I've noticed what I think is an incorrect behavior of stats::psmirnov function and consequently of ks.test when run in an exact mode. For example: psmirnov(1, sizes=c(50, 50), z=1:100, two.sided = FALSE, lower.tail = F, exact=TRUE) produces 2.775558e-15 However, the exact value should be 1/combination(100, 50), which is 9.9e-30. While the absolute error is small,...

...I have two main questions. Could you provide me some more information? I download different versions of the r language source code through r language website (https://www.r-project.org/). By reading R-4.3.3/src/library/stats/R/ks.test.R, I encounter the following problem: before call the default psmirnov function (in two-sample case), z <- NULL. According to the T and F the TIES, determines the value of z assigned to w or not. However, when psmirnov is called, z=w is always used. I am curious whether the TIES parameter can be omitted. Compared with the previous ks.test, such as version 4.1.3, t...

...TIES <- TRUE line 74 if (TIES) line 75 z <- w But later this z variable is not used as a variable in the code. It looks to me that this TIES variable can be deleted without affecting anything else. What I gathered from the investigation is that probably now ties are handled by psmirnov() and for numeric x and y the computations are exact even with ties, however I am a bit puzzled about the warning for approximate values, when exact = FALSE is set anyway. So my question - is everything currently OK with the code and the documentation?

...kup by angle brackets. * The new experimental primitive function declare() is intended to eventually allow information about R code to be communicated to the interpreter, compiler, and code analysis tools. The syntax for declarations is still being developed. * Functions psmirnov(), qsmirnov() and rsmirnov() in package stats have had argument two.sided renamed to alternative, to take into account that the permutation distributions of the one-sided statistics can be different in the case of ties. Consequence of PR#18582. * sort() is now an impli...

...kup by angle brackets. * The new experimental primitive function declare() is intended to eventually allow information about R code to be communicated to the interpreter, compiler, and code analysis tools. The syntax for declarations is still being developed. * Functions psmirnov(), qsmirnov() and rsmirnov() in package stats have had argument two.sided renamed to alternative, to take into account that the permutation distributions of the one-sided statistics can be different in the case of ties. Consequence of PR#18582. * sort() is now an impli...