search for: emmean

...hilesh Singh wrote: > I am a Professor of Statistics at Indira Gandhi Krishi Vishwavidyalaya, > Raipur, India. While teaching in class about analysis of variance using R, > I was doing a one-way analysis for the two data-sets given below in the > R-class. I got a typical error in "emmeans" package, please help: > > Data-set-1: > -------------- > Medley and Clements (1998) investigated the impact of zinc contamination > (and other heavy metals) on the diversity of diatom species in the USA > Rocky Mountains. The diversity of diatoms (number of species) and de...

...-Jan-2018 8:18 PM, "Sal Mangiafico" <salvatore.s.mangiafico at gmail.com> wrote: One way to avoid this error is to create the aov without using the with function, but instead use the data= option in the aov function. That is, medley2 = aov(diversity ~ zinc, data=medley.clementis) emmeans::emmeans(medley2, "zinc") You can see the difference in the calls: medley2$call medley.clementis.aov$call This works for the other data set as well, e.g. keough2 = aov(serpulid.ln ~ biofilm, data=keough.raimondi.ln) ~ Sal Mangiafico On 1/8/2018 4:44 PM, Rolf Turner wrote: On 07/0...

I am a Professor of Statistics at Indira Gandhi Krishi Vishwavidyalaya, Raipur, India. While teaching in class about analysis of variance using R, I was doing a one-way analysis for the two data-sets given below in the R-class. I got a typical error in "emmeans" package, please help: Data-set-1: -------------- Medley and Clements (1998) investigated the impact of zinc contamination (and other heavy metals) on the diversity of diatom species in the USA Rocky Mountains. The diversity of diatoms (number of species) and degree of zinc contamination (ca...

...also tried: (by (GROUP) cor.test (WK1FREQ,WK2FREQ)) (2) I also wanted to run the following analyses, written in SPSS syntax. I have searched the R archives and gather I may need to use lm instead of glm, but got lost. UNIANOVA avfreq BY sex WITH blik /METHOD = SSTYPE(3) /INTERCEPT = INCLUDE /EMMEANS = TABLES(OVERALL) /EMMEANS = TABLES(sex) /PRINT = DESCRIPTIVE ETASQ OPOWER HOMOGENEITY /CRITERIA = ALPHA(.05) /DESIGN = blik sex blik*sex . GLM avmth3fq avfreq BY SEX /WSFACTOR = timex 2 Polynomial /METHOD = SSTYPE (3) /PRINT = DESCRIPTIVE ETASQ TEST(MMATRIX) HOMOGENEITY /CRITER...

Dear R-Devel, I seem to no longer be able to access the bug-reporting system, so am doing this by e-mail. My report concerns models where variables are explicitly referenced (or is it "dereferenced"?), such as: cars.lm <- lm(mtcars[[1]] ~ factor(mtcars$cyl) + mtcars[["disp"]]) I have found that it is not possible to predict such models with new data. For example:

I read through the editorial. This is the one of the most mega-ultra-super-biased articles I've ever read. e.g. The authors encourage Baysian methods, and literally encourage subjective approaches. However, there's only one reference to robust methods and one reference to nonparametric methods, both of which are labelled as purely exploratory methods, which I regard as extremely