yuen does not work when there unused levels in the factors given to it. > yuen(GoalsGame ~ League, tr=0.2, data = SpainGer) Call: yuen(formula = GoalsGame ~ League, data = SpainGer, tr = 0.2)
Test statistic: NaN (df = NA), p-value = NA Trimmed mean difference: NaN 95 percent confidence interval: NaN NaN > yuen(GoalsGame ~ League, tr=0.2, data = droplevels(SpainGer)) Call: yuen(formula = GoalsGame ~ League, data = droplevels(SpainGer), tr = 0.2) Test statistic: 0.8394 (df = 16.17), p-value = 0.4135 Trimmed mean difference: -0.11494 95 percent confidence interval: -0.405 0.1751 Complain to WSR2's maintainer. > maintainer("WRS2") [1] "Patrick Mair <m...@fas.harvard.edu>" Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, Aug 30, 2016 at 2:07 PM, James Henson <jfhens...@gmail.com> wrote: > Greetings R community, > Here is a small but perplexing problem with the ‘yuen’ function in the > ‘WRS2’ package. > I begin with the ‘eurosoccer’ data frame from the ‘WRS2’ package. > Then make a subset that contains only two Leagues Spain and Germany > (subset data frame is ‘SpainGer’). The ‘yuen’ function cannot read > the data in subset data frame ‘SpainGer’. My code is below. > Thanks you for your help. > James F. Henson > > > # Examples from 'Robust Statistical Methods on R Using the WRS2 Package’ > # Robust t-test, and ANOVA (pages5-13) > library(WRS2) > data("eurosoccer") > class(eurosoccer) > print(eurosoccer) > library("digest") > library("DT") > datatable(eurosoccer) > str(eurosoccer) > # make a subset with only Spain and Germany > SpainGer <- subset (eurosoccer, subset = League =="Spain" | League == > "Germany") > print(SpainGer) > class(SpainGer) > str(SpainGer) > # The 'yuen' function can not read the data in the subset data.frame > "SpainGer" > yuen(GoalsGame ~ League, tr=0.2, data = SpainGer) > # the 'yuen' function works on the orginioal data.frame > yuen(GoalsGame ~ League, tr=0.2, data = eurosoccer) > # the 'aov' function reads the data in the subset data.frame "SpainGer" > Goals.fit <- aov(GoalsGame ~ League, data = SpainGer) > summary(Goals.fit) > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/ > posting-guide.html > and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.