similar to: Do I need a special package to run LSD tests?

Hola a todos Estoy trabajando con el package agricolae para realizar pruebas de medias y obtengo un error que no entiendo al trabajar las interacciones. Os explico: Poseo un Diseño Completamente al Azar para el arreglo de tratamientos factorial de la forma 2*2*4. El modelo consta de 6228 observaciones, el esquema para el anova es: A= 2 niveles={2,5} B= 4 niveles={F1,F2,F3,F4} C= 2

> On Dec 25, 2017, at 2:09 PM, Ahmed Attia <ahmedatia80 at gmail.com> wrote: > > The model should be class aov or lm and my model class is aovlist. > tried tidy from broom library but did not work. To make it class aov, > I had to remove the error term; > > model <- aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen),data=Rotationdata_R) You seemed to have

The model should be class aov or lm and my model class is aovlist. tried tidy from broom library but did not work. To make it class aov, I had to remove the error term; model <- aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen),data=Rotationdata_R) Ahmed Attia, Ph.D. Agronomist & Soil Scientist On Mon, Dec 25, 2017 at 7:38 PM, David Winsemius <dwinsemius at

Hi, I did the analysis of variance of a split-plot and the effect of treatment was significant. I would like compare treatment means using Tukey. I can't extract the mean square to apply HSD.test to use in agricolae package. anava = aov(ganhos ~ Blocos + Trat*Supl + Error(Blocos/Trat)) names(anava) summary(anava) require(agricolae) HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05)

Hi, I'm designing a experiment in order to compare the growing of several clones of a tree specie. It will be a complete randomized block design. How can I decide what model of mean comparision to choose? LSD, HSD,TukeyHSD, Duncan,...? Thanks in advance

Hi, I have a matrix of results of multiple 2x2 chi^2 tests, non- significant tests are marked as TRUE. Is there a function for grouping those variables in a similar way LSD.test from agricolae library does? I reviewed LSD.test's source but it's not helpful for me. This is my matrix: [,1] [,2] [,3] [,4] [,5] [,6] [,7] 1 TRUE FALSE TRUE TRUE TRUE FALSE TRUE 2 FALSE TRUE

Hello, My boss wants me to do a Duncan's test, which is under the agricolae package. Unfortunately I am not versed enough in R to run my data. I have 7 subspecies of deer mouse for which I have 23 measurements which are my variables of interest. I have run an ANOVA for each of the set of subspecies and variables, my data look like this: subspecies WMF 1 rowleyi 2.50 2 rowleyi 2.30

> On Dec 25, 2017, at 10:21 AM, Ahmed Attia <ahmedatia80 at gmail.com> wrote: > > LSD.test ?LSD.test No documentation for ?LSD.test? in specified packages and libraries: you could try ???LSD.test? -- David Winsemius Alameda, CA, USA 'Any technology distinguishable from magic is insufficiently advanced.' -Gehm's Corollary to Clarke's Third Law

LSD-test produces error for this code; code <- as.factor(Rotationdata_R$`Rot/code`) #factor in the main Nitrogen <- as.factor(Rotationdata_R$Nitrogen) #factor in the sub Rep <- as.factor(Rotationdata_R$REP) #blocks Year <- as.factor(Rotationdata_R$YEAR) #years model <- aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen)+Error(Rep/Year/code),data=Rotationdata_R)

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

Please help. I really like R and I have been looking at how to do LSD multiple comparison test with data that has more than one factor. So far, I am unsuccessful. Please help! Me

Hi there, Is there a function that can do a Fisher's LSD multiple comparisons in a two-way ANOVA? I hope to get a result similar with TukeyHSD(). Especially, I hope to know the significance of comparisons between the interactions of two factors. In the following example: x <- c(76, 84, 78, 80, 82, 70, 62, 72, 71, 69, 72, 74, 66, 74, 68, 66, 69, 72, 72, 78, 74, 71, 73, 67, 86, 67, 72,

Hi there, ============================= x=c(145,40,40,120,180, 140,155,90,160,95, 195,150,205,110,160, 45,40,195,65,145, 195,230,115,235,225, 120,55,50,80,45 ) y2=c( rep(as.character(1),5), rep(as.character(2),5), rep(as.character(3),5), rep(as.character(4),5), rep(as.character(5),5), rep(as.character(6),5) ) crd2=data.frame(x,y2) model1=aov(x~y2,data=crd2) TukeyHSD(model1)

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20111102/25eb7284/attachment.pl>

I am returning this to the R-help list. Please keep followups on the list. Yes, it can be done. It is not currently easy because multcomp doesn't have the syntax yet. Making this easy is on Torsten's to-do list for the multcomp package. See the MMC.WoodEnergy example in the HH package. The current version on CRAN is HH_1.17. Please see the discussion of this example in R-help:

If i insert \u2264 inside the text, like this lsd_label <- "LSD (P \u2264 0.05) = " the result is "LSD (P = 0.05) = " instead "LSD (P ? 0.05) = " how can i solve this problem? Thanks. -- View this message in context: http://r.789695.n4.nabble.com/less-than-or-equal-to-glyph-tp3073557p3073557.html Sent from the R help mailing list archive at Nabble.com.

I have no problem setting up my mixed model, or performing anova or lsmeans on my model?s outputs. However, performing lsd mean separation is giving me fits. So I do not have a problem when using two-way anova model. When using the code: fit.yield.add <- lm(data = ryzup, Yield ~ Rep + Nitrogen + Treatment) LSD.test(fit.yield.add, trt = "Nitrogen", alpha = 0.1, console = TRUE)

Hi all. I'm working on an assignment for a psychology class and I am not sure how to adjust the y-axis so it displays the range: 0, 5, 10, 15 The code below is almost ideal: lsd = c(3, 5, 13) mar = c(1, 2, 3) g_range <- range(0, lsd, mar) plot(lsd, type="o", col="blue", ylim=g_range, axes=FALSE, ann=FALSE) axis(1, at=1:3,

On 16/07/13 6:31 PM, Martijn van Beurden wrote: > On 16-07-13 09:07, Leigh Dyer wrote: >> Hi, >> >> On a particular input file, FLAC (testing with current git) greatly >> inflates its output if I encode at level 7, which enables >> --exhaustive-model-search. The source is a 24-bit WAV file of about >> 60MB; flac -6 encodes this to a 43MB FLAC file, but flac

On 16/07/13 8:10 PM, Erik de Castro Lopo wrote: > Leigh Dyer wrote: > >> Certainly -- I've uploaded the analysis files for both the -6 and -7 >> encodes, in case you wanted to compare: >> >> http://wootangent.net/~lsd/blah/6.ana >> http://wootangent.net/~lsd/blah/7.ana >> >> The encode seems to proceed normally until 59% of the way through the