similar to: Help with a more flexible funtion for multiple comparision of means

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),

Dear R list, Please, could someone help me with SSQ decomposition and contrasts. Below my data, graphic, ANOVAs and my doubt: # Data a = paste('a', gl(3, 8), sep='') b = paste('b', gl(2, 4, 24), sep='') tra = sort(paste('t', rep(1:6, 4), sep='')) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6,

# Dear R list, # # A have a doubt about SSQ decomposition and contrasts with ANOVAs. # So, I would like a tip from more advanced R users. # Below my data, the basic script and my doubts: # Data r = paste('r', gl(3, 8), sep='') e = paste('e', rep(gl(2, 4), 3), sep='') tra = sort(paste('t', rep(1:6, 4), sep='')) y = c(26.2, 26.0, 25.0, 25.4,

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, I'm trying to calculate Mahalanobis distances for 'Species' of 'iris' data as obtained below: Squared Distance to Species From Species: Setosa Versicolor Virginica Setosa 0 89.86419 179.38471 Versicolor 89.86419 0 17.20107 Virginica 179.38471 17.20107 0 These distances were obtained with proc 'CANDISC'

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've been using R under WindowsXP. So, where the R stores the graphic archives (don't saved) if I use the option options(graphics.record=T) inside of Rprofile.site file? Regards, -- /////\\\\\/////\\\\\/////\\\\\/////\\\\\ Jose Claudio Faria Brasil/Bahia/UESC/DCET Estatistica Experimental/Prof. Titular joseclaudio.faria em terra.com.br joseclaudio.faria em

Dear R list, I would like to illustrate the origin of the Student t distribution using R. So, if (sample.mean - pop.mean) / standard.error(sample.mean) has t distribution with (sample.size - 1) degree free, what is wrong with the simulation below? I think that the theoretical curve should agree with the relative frequencies of the t values calculated: #== begin options===== # parameters

Dear R users, The solution is probably simple but I need someone to point me to it. How can I to generate a matrix from a numeric sequence of 1:10 like 'A' or 'B' below: A) |--------------------| | 1 2 3 4 5 | |--------------------| | 1 | 0 | | 2 | 1 0 | | 3 | 2 5 0 | | 4 | 3 6 8 0 | | 5 | 4 7 9 10 0 | |--------------------|

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi. I have been using SAS for some time, and now I have discovered R. I am very happy with it, but I have not found out how to perform some of the multiple comparisons I was used to do in SAS. With the SAS/STAT, I generally used the MEANS (for comparison of arithmetic means) and the LSMEANS (for adjusted means) statements of the GLM procedure (I

Hi all, I would like to suggest that all R functions/etc like: codes-deprecated grid-internal ns-alt ns-dblcolon ns-hooks ns-internals ns-lowlev ns-reflect.Rd tools-internal ts-defunct utils-deprecated utils-internal ... and another i.e, function/word separate for '-' were all substituted by codes_deprecated grid_internal ns_alt

# Dear R list, # # I'm needing help with lattice, regression and respective lines. # My data is below: bra = gl(2, 24, label = c('c', 's')) em = rep(gl(3, 8, label = c('po', 'pov', 'ce')), 2) tem = rep(c(0, 0, 30, 30, 60, 60, 90, 90), 6) tem2 = tem^2 r = rep(1:2, 24) y = c(40.58, 44.85, 32.55, 35.68, 64.86, 51.95, 42.52, 52.21,

Dear list, For educational purposes I have been working with the script below. I have a observation: line 31 #CScal[i] = (amo^2) # IT IS WRONG, I KNOW, BUT IT MAKE R TO CRASHES! I'm thinking this is a possible bug in the R! And I have a couple of doubts: 1) line 51 #curve(dchisq(x, n-1), add = T, col = 'red' I think that it is correct, but the function no

Dear list, I'm trying to do a set of generic functions do make contingency tables from data.frames. It is just running "nice" (I'm learning R), but I think it can be better. I would like to filter the data.frame, i.e, eliminate all not numeric variables. And I don't know how to make it: please, help me. Below one of the my functions ('er' is a mention to EasieR,

Dears friends, I have been translating and extended, with educational porposes only, the available material in the www (http://zoonek2.free.fr/UNIX/48_R/all.html) "Statistcs avec R" for the Portuguese of Brazil. I consider the material very good and would like to share this work with anothers users R. Already I translated and I extended two chapters: Introduction to the R and analysis

Dear list, I'm needing submit values (V1 = 8,6,4,3,1,2,9) (Id = 2:8) of a data.frame (DF), like below Id V1 V2 ... 1 0 1 ... 2 8 10 ... 3 6 2 ... 4 4 4 ... 5 3 7 ... 6 1 8 ... 7 2 6 ... 8 9 7 ... 9 6 1 ... 10 5 4 ... to selection (>=2 and <8) for remanescents like below: Id V1 V2 ... 1 0 1 ... 2 .

Dears members of R list, I would like that a more experienced R user help me to complete this analysis: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4, 22.8, 21.3) df =

Hi! I try to explain the efffect of (1) forest where i took samples's soils (* Lugar*: categorical variable with three levels), (2) nitrogen addition treatments (*Tra*: categorical variable with two levels) on total carbon concentration's soil samples (*C: *continue* *variable) during four months of sampling (*Time:* categorical and ordered variable with four levels). I fitted the

Hi! First I would like to thank you for the help with my first question. However,I am afraid I am still into trouble. I think I have properly used the getGropus commad but seems to be a fatal error when calculating the fixed effects (see bellow) and I cannot understand why. Any inputs? Thanks in advance, /Jose hongos1=read.table("c:/temp/hongos1.txt", header=T) > attach(hongos1)

Hi, everyone. GAD package analyses complex ANOVA models with any combination of orthogonal/nested and fixed/random factors, as described by Underwood (1997). There are two restrictions: (i) data must be balanced; (ii) fixed nested factors are not allowed. Homogeneity of variances is checked using Cochran's C test and "a posteriori" comparisons of means are done using