search for: aovlist

I am having what I think is a strange problem with applying TukeyHSD to an aov fit with error strata. TukeyHSD is supposed to take "A fitted model object, usually an 'aov' fit." aov (with error strata) is supposed to generate an object of type aovlist, which is a list of objects of type aov. But I can't seem to feed components of my aovlist to TukeyHSD. I guess I wouldn't expect to be able to use the error strata, but I did expect to be able to use the final stratum. I have posted a complete example, which I hope explains why I am con...

I WROTE THIS FUNCTION BELOW test <- function(x, ...) UseMethod('test', x) test.data.frame = function(x, model, which, error, ...) { av <- aov(formula(model), data = x) res <- test.aovlist(av, which = which, error = error) return(res) } test.aovlist <- function(x, which, error, ...) { mm <- model.tables(x, "means") tabs <- mm$tables[-1] tabs <- tabs[which] res <- list(av=x, which=which) return(res) } GENERETING DATA var <- sort...

...and apparently ran correctly, but it failed "make check". Earlier versions of R failed "make check" as well, not always for the same reason. I just located the failure in the tests. It is in the file: R-2.6.0/tests/reg-tests-1.R and fails in the code: ## related checks on eff.aovlist example(eff.aovlist) # helmert contrasts eff1 <- eff.aovlist(fit) fit <- aov(Yield ~ A * B * C + Error(Block), data = aovdat) eff2 <- eff.aovlist(fit) stopifnot(all.equal(eff1, eff2)) # will have rounding-error differences options(contrasts = old) ## Were different in earlier versions The...

Hi friends, there is methods() function to see the all available methods for a particular function, for example: > head(methods("print")) [1] "print.acf" "print.anova" "print.aov" "print.aovlist" "print.ar" "print.Arima" In this list, there are some functions which are asterisked like print.acf(). How can I see the contents of those function? Thanks and regards, [[alternative HTML version deleted]]

...I use \Sexpr(rx) only the first value (in this case 21) is printed out. I have been trying to get xtable to work. > xtable(rx) Error in UseMethod("xtable") : no applicable method for "xtable" > methods(xtable) [1] xtable.anova* xtable.aov* xtable.aovlist* [4] xtable.coxph* xtable.data.frame* xtable.glm* [7] xtable.lm* xtable.matrix* xtable.prcomp* [10] xtable.summary.aov* xtable.summary.aovlist* xtable.summary.glm* [13] xtable.summary.lm* xtable.summary.prcomp* xtable.table* [16] xtable.ts*...

...action time (RT), there is a 3-level Condition factor (Cond; within-subjects), a number of subjects (Subj), and the dataframe is called WMU3C. The model is > RT.aov <- aov(RT~Cond + Error(Subj/Cond), WMU3C) I understand that TukeyHSD only works with an aov object, but that RT.aov is an aovlist object. > class(RT.aov) [1] "aovlist" "listof" I've tried to work around it using the "maiz" example in the MMC documentation of the HH package (a solution previously recommended), but I couldn't get it to work: My best shot was to calculate another ao...

...and 'data.frame': 1 obs. of 5 variables: .. ..$ Df : num 8 .. ..$ Sum Sq : num 36.4 .. ..$ Mean Sq: num 4.55 .. ..$ F value: num NA .. ..$ Pr(>F) : num NA ..- attr(*, "class")= chr [1:2] "summary.aov" "listof" - attr(*, "class")= chr "summary.aovlist" But for this dummy data > dput(BGBottles,"") structure(list(Machine = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L), .rk.invalid.fields = structure(list(), .Names = character(0)), .Label = c("1", "2", "3", "4"), class = &q...

...e. Can somebody let me know how to figure them? Value: An object of class 'c("aov", "lm")' or for multiple responses of class 'c("maov", "aov", "mlm", "lm")' or for multiple error strata of class '"aovlist"'. There are 'print' and 'summary' methods available for these. Regards, Peng

...a=MyData) My first naive question is: Is the Error(...) term correct for the design that I describe? I am not quite sure if I understand the syntax correctly. In any case, it actually seems to work fine, but I can't figure out how to do post hoc testing on this. TukeyHSD does not work for the aovlist which is returned, so I am kind of stuck. Is there a simple way to do the post hoc test based on the aovlist? I have been reading various questions on the list and I think that I have understood that I should be using lme from the nlme package and this is where I run into some problems. As I under...

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 Winsem...

...it being non-visible. I ran its help file example, which printed visible output. I did not notice how it differed from other functions, like summary.data.frame that is not marked non-invisible. What does the non-visible mean? Thanks, Bob > methods(summary) [1] summary.aov summary.aovlist [3] summary.connection summary.data.frame [5] summary.Date summary.default [7] summary.ecdf* summary.factor [9] summary.glm summary.infl [11] summary.lm summary.loess* [13] summary.manova summary.m...

...Sq Mean Sq F value Pr(>F) trat 4 1.580 0.3950 0.821 0.541 Residuals 10 4.811 0.4811 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Residuals 30 15.78 0.5259 Althougth the Anova Model object it isn't a object of class "aov" but a object of class"aovlist": #> attr(AnovaModel.2,"class") [1] "aovlist" "listof" This class isn't suitable for the glht function neither TukeyHSD function. I can extract a object of class "aov" through: #> ls.str(pat="AnovaModel.2") AnovaModel.2 : List of...

...a=MyData) My first naive question is: Is the Error(...) term correct for the design that I describe? I am not quite sure if I understand the syntax correctly. In any case, it actually seems to work fine, but I can't figure out how to do post hoc testing on this. TukeyHSD does not work for the aovlist which is returned, so I am kind of stuck. Is there a simple way to do the post hoc test based on the aovlist? I have been reading various questions on the list and I think that I have understood that I should be using lme from the nlme package and this is where I run into some problems. As I under...

It doesn't print anything: the summary.aov (or summary.aovlist) print method does. ?summary.aov tells you the structure of the objects they return. On Tue, 15 Feb 2005, RenE J.V. Bertin wrote: > I'd like to annotate a plot with the output of summary(aov(model)), > ideally just with the significant effects. I don't find a means to > redir...

Hello, I am trying to print out ANOVA summaries --- as returned by summary (aov(...))--- in tab-delimited format, either to a file or to the console, but so far I haven't been able to figure out a solution. My particular ANOVA has repeated measures, so it is an instance of class aovlist. Can anyone offer any suggestions? I am relatively new to R, so please forgive me if I am missing something obvious. Regards, -Chris Chris R. Sims Doctoral Student, Department of Cognitive Science Rensselaer Polytechnic Institute Troy, NY 12180 [[alternative HTML version deleted]]

Hi all, I search the archive for finding a simple solution for using TukeyHSD with a multistratum aov result (a repeated emasure anova). The Question have been asked but I've found no clear answer. res<-aov(y~Fa*Fb+Error(Subject/(Fa*Fb)) ) I think that the problem is that res is an aovlist object instead of the "aov" object required by TukeyHSD. Is there an easy solution to this problem ? Sylvain Cl?ment Sylvain CLEMENT (MCF) JE Neuropsychologie et Cognition Auditive UFR de Psychologie BP 60149, Universit? Ch. de Gaulle Lille 3 Domaine universitaire "Pont de Bois&...

Hello: I s there a way to use xtable with objects from the psych package, particularly principal()? Is there a difference between princomp and principal? xtable seems to play better with princomp. Thank you. Yours, Simon Kiss ********************************* Simon J. Kiss, PhD Assistant Professor, Wilfrid Laurier University 73 George Street Brantford, Ontario, Canada N3T 2C9

...ta_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen)+Error(Rep/Year/code),data=Rotationdata_R) summary(model) LSD.test(y=model,trt="code",console=TRUE) the error; Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : cannot coerce class "c("aovlist", "listof")" to a data.frame What is wrong with this. Ahmed Attia, Ph.D. Agronomist & Soil Scientist

> 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 missed my point that LSD.test is not in the packages loaded by...

I'm using the current version of xtable for 1.9.0 and I have an interesting error: Error in match.names(clabs, names(xi)) : names don't match previous names: F value, Pr(>F) In addition: Warning message: longer object length is not a multiple of shorter object length in: clabs == nmi This is produced in the following manner: >data.trans <-