weighted.bagplot: Compare two bagplots

In pbiecek/PISAtools: A set of tools useful when dealing with PISA data.

Description

This function compares two bagplots.

Usage

  weighted.bagplot(x, y, weights, resolution = 10) 

  compareBagPlots(bp.selected, bp.reference, col.selected="#55555555", col.reference="#ffffff", ...)

Author(s)

Przemyslaw Biecek

Examples

  library(PISA2009lite)
  sel <- student2009[student2009$OECD == "OECD",]
  rPV <- removeGroupMeans(sel$PV1MATH, factor(paste(sel$CNT, sel$SCHOOLID)), sel$W_FSTUWT)
  rES <- removeGroupMeans(sel$ESCS, factor(paste(sel$CNT, sel$SCHOOLID)), sel$W_FSTUWT)
  rBagPlot <- weighted.bagplot(rPV, rES, sel$W_FSTUWT)

  sel <- student2009[student2009$CNT == "Poland",]
  pPV <- removeGroupMeans(sel$PV1MATH, factor(paste(sel$CNT, sel$SCHOOLID)), sel$W_FSTUWT)
  pES <- removeGroupMeans(sel$ESCS, factor(paste(sel$CNT, sel$SCHOOLID)), sel$W_FSTUWT)
  pBagPlot <- weighted.bagplot(pPV, pES, sel$W_FSTUWT)

  compareBagPlots(rBagPlot, pBagPlot, 
        xlim=c(-100,100), ylim=c(-1,1), xlab="PV1MATH diff", ylab="ESCS diff")

pbiecek/PISAtools documentation built on May 24, 2019, 10:35 p.m.