• A cabbages adathalmaz a MASS-ból
  • 60 egyed, két faktor és két mennyiségi változó
  • szimultán dobozábrák
In [64]:
library(MASS)
#help(cabbages)
#cabbages
In [65]:
by(cabbages$HeadWt,cabbages$Cult,summary)
boxplot(cabbages$HeadWt~cabbages$Cult)

cabbages$Cult: c39
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.600   2.500   2.800   2.907   3.450   4.300 
------------------------------------------------------------ 
cabbages$Cult: c52
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.000   1.625   2.100   2.280   2.800   4.200
In [66]:
by(cabbages$HeadWt,cabbages$Date,summary)
boxplot(cabbages$HeadWt~cabbages$Date)

cabbages$Date: d16
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.700   2.150   2.450   2.720   3.125   4.300 
------------------------------------------------------------ 
cabbages$Date: d20
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  2.000   2.600   2.800   2.955   3.200   4.200 
------------------------------------------------------------ 
cabbages$Date: d21
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.000   1.500   1.650   2.105   2.450   4.200
In [71]:
by(cabbages$VitC,cabbages$Date,summary)
boxplot(cabbages$VitC~cabbages$Date)

cabbages$Date: d16
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  42.00   50.75   55.50   56.40   61.50   72.00 
------------------------------------------------------------ 
cabbages$Date: d20
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  41.00   46.50   53.00   54.15   60.25   72.00 
------------------------------------------------------------ 
cabbages$Date: d21
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  45.00   54.00   66.00   63.30   71.25   84.00
In [72]:
by(cabbages$VitC,cabbages$Cult,summary)
boxplot(cabbages$VitC~cabbages$Cult)

cabbages$Cult: c39
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  41.00   46.00   51.00   51.50   54.75   68.00 
------------------------------------------------------------ 
cabbages$Cult: c52
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  47.00   58.00   64.50   64.40   70.75   84.00
In [78]:
qqnorm(cabbages$VitC)
qqline(cabbages$VitC)
In [84]:
sel=cabbages$Cult=="c52"
data=cabbages$VitC[sel]
qqnorm(data)
qqline(data)