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)
James: To be clear, I know nothing about this package, but the following may be helpful for those who do: 1. What exacty was the error message that you received? Or was it a warning? "cannot read" may be too vague to allow others to respond. 2. ?maintainer to find the package maintainer. It may be necessary to contact them directly with the problem. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Aug 30, 2016 at 2:07 PM, James Henson <jfhenson1 at 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 at 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.
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 <mair at fas.harvard.edu>" Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, Aug 30, 2016 at 2:07 PM, James Henson <jfhenson1 at 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 at 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]]
Thanks, The 'droplevels' statement works. Best regards, James Henson On Tue, Aug 30, 2016 at 5:34 PM, William Dunlap <wdunlap at 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 <mair at fas.harvard.edu>" > > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Tue, Aug 30, 2016 at 2:07 PM, James Henson <jfhenson1 at 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 at 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. > >