Thanks, The 'droplevels' statement works. Best regards, James Henson On Tue, Aug 30, 2016 at 5:34 PM, William Dunlap <wdun...@tibco.com> wrote: > 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. > >
______________________________________________ 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.