similar to: New package GAD (General ANOVA Design)

Jose - Before implementing SNK and Duncan's, you may want to be aware of some criticisms of these methods: >From Hsu (1996), "Newman-Keuls multiple range test is not a confident inequalities method and cannot be recommended." "Duncan's multiple range test is not a confident inequalities method and cannot be recommended either. In the words of Tukey (1991),

How can I perform a Student-Newman-Keuls-Test for multiple comparision of means in R? (I did not manage to find any specific function in the libraries) Jörg Peter Baresel Technische Universität München Institut für Ackerbau und INformatik im Pflanzenbau D-85354 Freising Lange Point 51 [[alternative HTML version deleted]]

Here is a little patch against 2.9p1 that performs the SNK (also known as TIS authserv) challenge-response automaticly instead of asking the user. hope you find it useful. --larry -------------- next part -------------- diff -NuBw openssh-2.9p1/Makefile.in openssh/Makefile.in --- openssh-2.9p1/Makefile.in Thu Apr 26 20:31:08 2001 +++ openssh/Makefile.in Wed Jun 6 16:15:56 2001 @@ -43,9 +43,9

Dear R-list, Could anybody tell me (or give me a tip) of how to implement the Duncan distribution in R? I've been trying to make a new and more flexible function for multiple comparison of means: Tukey, SNK and Duncan, from 'aov' objects, like TukeyHSD function. For while, it is running nice (Tukey and SNK), for simple design (completely randomized, randomized block and Latin

Hello Everyone,   I've been learning to use R in my spare time over the past several months. I've read about 7-8 books on the subject. Lately I've been testing what I've learned by trying to replicate the analyses from some of my SAS books. This helps me make sure I know how to use R properly and also helps me to understand how the two programs are similar and different.   Below is

I have Samba PDC and use tdbsam backend. I use the version 3.0.11, but recently I have exchanged it for the version 3.0.23c. Now I have a number of problem. Here is one of then. I wish to move one user "gad" from Domain Users (RID 513) to Domain Guests (RID 514). Now Primary group of user "gad" is "Domain Users": pdbedit -Lv gad -d0 Unix username: gad NT

Queridos R-users: Tengo una duda que hace mucho tiempo que estoy intentando resolver, os explico a modo de ejemplo: Tengo estos efectos: Año(5 niveles),Localidad (10 niveles) y genotipo (3 niveles), año y localidad son aleatorios y genotipo es fijo (los he escogido yo). Me gustaría hacer obtener una tabla parecida a la Tabla Anova donde aparezca cada factor y sus interacciones y

Hi all, I just upgraded more servers, and doing some tests I found that my setup for kerberos/ldap authentication against Active Directory is no more working. I don't know why... I followed some times ago scott Lowe blog for this setup : http://blog.scottlowe.org/2007/01/15/linux-ad-integration-version-4/ And it was working correctly until the upgrade. What is curious is that id command

On 8/3/07 22:12, "Gad Abraham" <g.abraham at ms.unimelb.edu.au> wrote: > Laura Hill wrote: >> >> >> On 7/3/07 00:15, "Gad Abraham" <g.abraham at ms.unimelb.edu.au> wrote: >> >>>> On 6 Mar 2007, at 08:54, Laura Hill wrote: >>>> >>>>> Hi, >>>>> >>>>> My name is

Hi all, I have a Dell server, CentOS 5.6 64 bits, on which I configured the CR repository. I just tried to update to 5.7, and had this error : # yum update .... http://mirror.centos.org/centos/5/cr/x86_64/repodata/filelists.sqlite.bz2: [Errno 14] HTTP Error 404: Not Found Trying other mirror. Error: failure: repodata/filelists.sqlite.bz2 from cr: [Errno 256] No more mirrors to try. You could

Hi all How can I change the limits (xlim or ylim) in a plot that has been already created? For example, consider this naive example curve(dbeta(x,2,4)) curve(dbeta(x,8,13),add=T,col=2) When adding the second curve, it goes off the original limits computed by R for the first graph, which are roughly, c(0,2.1) I know two obvious solutions for this, which are: 1) passing a sufficiently large

hello, I need to groupe some means and I wonder how to do to get critical range table for tukey or for Neuman-Keuls I think it's possible to do that with the qtukey() function? thanks. _____________________________________________________________________________ [[alternative HTML version deleted]]

Hello, I think that to apply the Neuman-Keuls test under R it's possible with the TukeyHSD function and I know that to do a Neuman-Keuls test I have to sort the means but I can't use the TukeyHSD function I don't understand how to do to apply this function to my data. In the R help they use the TukeyHSD function to a aov() do you think that I have to do teh same? thanks.

Hello guys, Is there any way to match a render :nothing? I coudn''t find any way to do this so i''ve just changed my controllers to do a "head :ok", but it would be nice to know if there is any other way :) -- Maur?cio Linhares http://alinhavado.wordpress.com/ (pt-br) | http://blog.codevader.com/ (en) Jo?o Pessoa, PB, +55 83 8867-7208

Hello everyone, Like a lot of people, I have been looking for functions in R doing ANOVA (ok) and performing multicomparisons (like Student-Newman-Keuls, etc.). As I have been a little bit disappointed, I have bee looking through the net for such "open source" softwares. I found one in: http://www.statpages.org/miller/openstat/OS4.html I have begun to use it, and it seems good and

Hi, Can anyone explain how the standard error in arima() is calculated? Also, how can I extract it from the Arima object? I don't see it in there. > x <- rnorm(1000) > a <- arima(x, order = c(4, 0, 0)) > a Call: arima(x = x, order = c(4, 0, 0)) Coefficients: ar1 ar2 ar3 ar4 intercept -0.0451 0.0448 0.0139 -0.0688 0.0010 s.e.

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn

Hi, I'm using survreg() from the survival package for parametric survival regression (modelling inter-arrival times of patients to a waiting list as exponentially distributed, with various regressors such as queue size and season). Does anyone know which algorithm survreg() uses for this? Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics The University of Melbourne

hello, I have programmed this function to calculate the Neuman-Keuls test but I have a problem the function return an empty list and I don't know why. summary(fm1) E <- sqrt((summary(fm1)[[1]]["Residuals","Mean Sq"])/length(LR)) lst <- list() lst1 <- list() lst2 <- list() NK <- function (x) { if (length(x) == 2) { Tstudent <- t.test(subset(exple,