similar to: ks.test (PR#1004)

There is a fairly new bug in wilcox.test in R-2.2.1 (stable). It wasn't there when I last taught nonparametrics in fall 2003. Line 86 of wilcox.test.R achieved.alpha<-2*psignrank(trunc(qu),n) It should be achieved.alpha<-2*psignrank(trunc(qu)-1,n) If you don't see why, decode the cookbook instructions p. 56 in Hollander and Wolfe (2nd ed.) or see

Full_Name: Thomas Waterhouse Version: 2.9.1 OS: OS X 10.5.7 Submission from: (NULL) (216.239.45.4) ks.test uses a biased approximation to the p-value in the case of the two-sample test with ties in the pooled data. This has been justified in R in the past by the argument that the KS test assumes a continuous distribution. But the two-sample test can be extended to arbitrary distributions by a

The point estimates produced by wilcox.test are perverse (not wrong, just brain damaged). The Hodges-Lehmann estimator that goes with the signed rank test is the median of the Walsh averages. The Hodges-Lehmann estimator that goes with the rank sum test is the median of the pairwise differences. wilcox.test agrees except that it uses the following very peculiar definition of "sample

I have to start by saying that I am new to R, so I might miss something crucial here. It seems to me that the results of friedman.test and ks.test are "wrong". Now, obviously, the first thing which crossed my mind was "it can't be, this is a package used by so many, someone should have observed", but I can't figure out what it might be. Problem: let's start with

I have a suse 8.0 server running samba 2.2.3, its ip adress is in 134.84.86.x and I'm trying to get it so that it works for both its subnet and 160.94.148.x. I have it configured how I want it for its own subnet, but I can't get it to work in the 160.94.148 subnet. This is what my smb.conf file looks like: [global] netbios name = muskrat workgroup = Statistics server string

I'm working with some old code that worked under R 2.15.3 but gives an error under 3.0.0 and 3.0.1 : > 1.0 - .C("psmirnov2x", p = as.double(10.0), as.integer(1000), as.integer(1000))$p Error in .C("psmirnov2x", p = as.double(10), as.integer(1000), as.integer(1000)) : C symbol name "psmirnov2x" not in load table I don't understand why that function would

Hello, I am working on a simple non-parametric (Theil) regression function and and am following Hollander and Wolfe 1999 text. I would like some help making my function faster. I have compared with pre-packaged version from "MBLM", which isnt very fast either, but it appears mine is faster with N = 1000 (see results below). I plan on running this function repeatedly, and I generally

Full_Name: Charles J. Geyer Version: 1.8.0 OS: i686-pc-linux-gnu (Suse 8.2) Submission from: (NULL) (134.84.86.22) Two bugs (perhaps related, perhaps independent) revealed by the same Poisson regression with offset mydata <- read.table(url("http://www.stat.umn.edu/geyer/5931/mle/seeds.txt")) out.fubar <- glm(seedlings ~ burn01 + vegtype * burn02 + offset(log(totalseeds)),

Here's a simpler example that should reproduce that error for you: ans <- utils::select.list(c("hello", "world", "again"), graphics=TRUE) Does it? FYI, I installed R 3.6.3 from source on CentOS 7 a few hours ago, and for me the above works just fine. For your immediate needs of selecting a CRAN mirror, you can set: options(menu.graphics = FALSE) as a

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. --1132542651-1468968864-1260896436=:8788 Content-Type: TEXT/PLAIN; format=flowed; charset=US-ASCII I am using R version 2.6.2 (2008-02-08) on a Ubuntu Linux system. I seemed to be finding occasional errors in the p-values produced by ks.test(a,b)

Is it possible to implement the Jonckheere-Terpstra test for ordered alternatives using the coin package: Conditional Inference Procedures in a Permutation Test Framework? I found jonckheere.test{clinfun}, but it uses a normal approximation when ties are present in the data. To make this concrete, I've include a small dataset. Thanks. --Dale Hollander and Wolfe, 1999 Table 6.6, pg 205

Just built 3.6.3 from source and tcl doesn't work. Worked fine with the same laptop in 3.6.2. Here's the exact error. blurfle$ R --vanilla R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to

I knew I could work around. But this shouldn't happen. And yes. Same problem with your example. blurfle$ R --vanilla R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain

OK. I give up. I'll ask a stupid question. How do I get the $!#@*$ signif stars line printed by summaries to not look extremely bizarre in the latex produced by Sweave? For example, see p. 7 of http://www.stat.umn.edu/geyer/aster/library/aster/doc/tutor.pdf I can see what the problem is. R emits non-ascii characters (as it is supposed to do), Sweave puts them in the tex file, and

>>>>> Charles Geyer writes: > In that case, perhaps the question could be changed to could CC0 be > added to the list of R licences. Right now the only CC licence that > is in the R licenses is CC-BY-SA-4.0. Hmm, I see Name: CC0 FSF: free_and_GPLv3_compatible (https://www.gnu.org/licenses/license-list.html#CC0) OSI: NA (https://opensource.org/faq#cc-zero) URL:

The subject line says it all. How can I find what Error: cannot set length of non-vector means? RTFS is no help. I can find out of course that it comes from "lengthgets", but who called that? Not me! The situation is as follows. I am trying to get a package ready for CRAN. Every time I run "R CMD check" I get this error when building the package vignette. But (!)

Can anyone tell me if the following C code (which proved very useful when using the snow package -- use it to nice slaves) compiles and dyn.loads under Windoze or Mac? It is (apparently) POSIX, so I suppose it is fine in OS X, but does Windoze's advertized POSIX compliance mean anything here? If it doesn't work on some platform, how do I deal with that? I want to submit as contributed

This is just to alert the list about a new technical report http://www.stat.umn.edu/geyer/gdor/ that shows how to determine, using exact infinite-precision rational arithmetic when the MLE exists in a GLM (or for that matter in any exponential family, such as aster models -- contributed package aster) and what to do when it does not. This suggests many changes to glm, glm.fit, predict.glm,

Hello group, Can someone please put me through how to estimate the nuisance parameters (Saturation parameter and radius) for the Geyer saturation process? It seems quite confusing to me how these parameters are achieved, but they are very important in determining the interaction between points in a point process. I am using Spatstat package. Thanks. neba [[alternative HTML version deleted]]

Does anyone know if any R package has a function that will conduct Hollander's test of bivariate symmetry? (Hollander, Biometrika, 1971) Either the exact test or an asymptotic version would be sufficient for my purposes. Thanks. Joe Boyer Statistical Sciences Renaissance Bldg 510, 3233-D Mail Stop RN0320 8-275-3661 cell: (610) 209-8531 [[alternative HTML version deleted]]