similar to: symmetric dotplot? (Wilkinson 1999, Am Stat 53 (3) 276-281

Given this example #start code a<-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940, 760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430) b<-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90, 3220,490,20790,290,740,5350,940,3910,0,640,850,260) wilcox.test(a, b, paired=FALSE) #sum of rank for first sample sum.rank.a <-

Hi all, I have considered neural network to classify the health status of the cow. I found a very neatly written R codes for classification method in here<http://home.strw.leidenuniv.nl/~jarle/IAC/RRoutines/classification-example.R> . It would be very helpful if you can answer some of the questions, that I am struggling with, I have set of time series data from different animals, I use

>>>>> Gustavo Fernandez Bayon <gbayon at gmail.com> >>>>> on Thu, 24 Aug 2017 16:42:36 +0200 writes: > Hello, > While doing some enrichment tests using chisq.test() with simulated > p-values, I noticed some strange behaviour. The computed p-value was > extremely small, so I decided to dig a little deeper and debug >

Hi, My name is Moise; I am trying to create a jittered stripchart of the average profit per employee. I enter the following code > stripchart(AvgProfitPerEmployee, method = "jitter", pch = 1, main="Average Profit in dollars per person", xlab = "Number of Employees", ylab = "Profit per Employee") > when I press enter , no error is shown and now chart

> On 25 Aug 2017, at 11:23 , Jari Oksanen <jari.oksanen at oulu.fi> wrote: > > It is not about "really arge total number of observations", but: > > set.seed(4711);tabs <- r2dtable(1e6, c(2, 2), c(2, 2)); A11 <- vapply(tabs, function(x) x[1, 1], numeric(1));table(A11) > > A11 > 0 1 2 > 166483 666853 166664 > > There are

Hi, I understand that Simon may be discontinuing his OpenSSH work. Does anyone know if someone plans to maintain the patch? Thank you, -- ******************************************************* Quellyn L. Snead UNIX Effort Team ( unixeffort at lanl.gov ) CCN-2 Enterprise Software Management Team Los Alamos National Laboratory (505) 667-4185 Schedule B

It is not about "really arge total number of observations", but: set.seed(4711);tabs <- r2dtable(1e6, c(2, 2), c(2, 2)); A11 <- vapply(tabs, function(x) x[1, 1], numeric(1));table(A11) A11 0 1 2 166483 666853 166664 There are three possible matrices, and these come out in proportions 1:4:1, the one with all cells filled with ones being most common. Cheers, Jari

Hi STATISTICA has a function called "Two-way joining" (see http://www.statsoft.com/TEXTBOOK/stcluan.html#twotwo) and the reference material states that this is based on the method as published by Hartigan (found this paper: http://www.jstor.org/pss/2284710 through wikipedia). What is the relationship (if any) between the "heatmap" function in R and this technique? Is there an

Hi all, Is there any code to run fixed effects Tobit models in the style of Honore (1992) in R? (The original Honore article is here: http://www.jstor.org/sici?sici=0012-9682%28199205%2960%3A3%3C533%3ATLALSE%3E2.0.CO%3B2-2) Cheers David [[alternative HTML version deleted]]

I'm not sure if this is the right mailing list. My apologies if I'm sending the question to the wrong mailing list. (The r-packages list moderator suggested using the R-help list.) My question is about the "chickwts" data in the datasets package. According to help(chickwts), the source of the data is ?Anonymous (1948) Biometrika, 35, 214. But this reference seems irrelevant.

Greetings, I'm working on some analyses where I need to calculate wilcox tests for paired samples. In my current literature search I've found a few papers on sample size determination for the wilcox test notably: Sample Size Determination for Some Common Nonparametric Tests Gottfried E. Noether Journal of the American Statistical Association

Dear all, I would like to run a simple regression model y~x1+x2+x3+... The problem is that I have a lot of independent variables (xi) -- around one hundred -- and that some of them are categorical with a lot of categories (like, for example, ZIP code). One straightforward way would be to (a) transform all categorical variables into 1/0 dummies and (b) enter all the variables into an lm model.

I've noticed that the implementation of Kendall's Tau in R is O(N^2). The following reference describes how it can be done in O(N log N): A Computer Method for Calculating Kendall's Tau with Ungrouped Data William R. Knight Journal of the American Statistical Association, Vol. 61, No. 314, Part 1 (Jun., 1966), pp. 436-439 http://www.jstor.org/pss/2282833 I'm interested in

I'd like to use stripplot for some plots because I want to use the jitter parameter. On the other hand, I'd like to use dotplot because I'd like to have the horizontal lines that it includes. dotplot doesn't have a jitter option and I'm not having any success with getting panel.grid(h=-1) with stripplot. Can anyone show me how to make dotplot-like lines on a stripplot? Or

For each level of the factor in dotplot, I have time points I'd like to connect with a line. In the example below, 'x' represents a starting time and 'd' a duration, and I wish to connect 'x' to 'x+d'. Ordinarily I would use Dotplot from hmisc for this, but I have not been able to find a time class that Dotplot will allow. I can get lattice dotplot to put

Hi All. Markus has commited the long-awaited GSSAPI patch to OpenBSD's ssh. There are patches. The first [1] is a straightforward port of the OpenBSD code to Portable. The second [2] contains the parts I've stolen from Simon Wilkinson's portable GSSAPI patch in an attempt to make it build. It is incomplete and doesn't currently work. The PAM support is not there and

Will Simon Wilkinson's GSSAPI Key Exchange patch ever be incorporated into the OpenSSH source? http://www.sxw.org.uk/computing/patches/openssh.html I'm sure I'm not the only one that uses it and would like to see it become part of the OpenSSH source. Is there something missing or is there some technical/philosophical reason for not including it?

Hello, I have trouble adding an abline to a dotplot() from "lattice" package. For example, I would like to draw a line at x=3: > library(lattice) > x<-1:5 > names(x) <- c("a", "b", "c", "d", "e") > dotplot(sample(x)) > panel.abline(v=3) Produces a line on the left hand side of the origin. What would be a correct

Dear list, is it possible to change the background color of dotplot's points? I tried in many ways but unsuccessfully Thanks in advance Gianandrea require(lattice) dotplot(variety ~ yield | site, data = barley, groups = year, pch=21) dotplot(variety ~ yield | site, data = barley, groups = year, pch=21, bg=c("2","3")) ??!!! -- View this message in context:

I found this paper on ANOVA on unequal error variance. Has this be incorporated to any R package? Is there any textbook that discuss the problem of ANOVA on unequal error variance in general? http://www.jstor.org/stable/2532947?cookieSet=1