[R] QQ plots and boxcox

Erin Hodgess erinm.hodgess at gmail.com
Tue Dec 2 22:28:05 CET 2008


Dear R People:

In the DASL library, there is a story about hot dogs.

Here are the data:
Beef	186	495
Beef	181	477
Beef	176	425
Beef	149	322
Beef	184	482
Beef	190	587
Beef	158	370
Beef	139	322
Beef	175	479
Beef	148	375
Beef	152	330
Beef	111	300
Beef	141	386
Beef	153	401
Beef	190	645
Beef	157	440
Beef	131	317
Beef	149	319
Beef	135	298
Beef	132	253
Meat	173	458
Meat	191	506
Meat	182	473
Meat	190	545
Meat	172	496
Meat	147	360
Meat	146	387
Meat	139	386
Meat	175	507
Meat	136	393
Meat	179	405
Meat	153	372
Meat	107	144
Meat	195	511
Meat	135	405
Meat	140	428
Meat	138	339
Poultry	129	430
Poultry	132	375
Poultry	102	396
Poultry	106	383
Poultry	94	387
Poultry	102	542
Poultry	87	359
Poultry	99	357
Poultry	107	528
Poultry	113	513
Poultry	135	426
Poultry	142	513
Poultry	86	358
Poultry	143	581
Poultry	152	588
Poultry	146	522
Poultry	144	545


Here is my work:
> dog1.df <- read.table(file="dog1.dat",as.is=F,header=F,
+ col.names=c("type","cal","sodium"))
> dog1.aov <- aov(cal~type,data=dog1.df)
> summary(dog1.aov)
            Df  Sum Sq Mean Sq F value    Pr(>F)
type         2 17692.2  8846.1  16.074 3.862e-06 ***
Residuals   51 28067.1   550.3
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> TukeyHSD(dog1.aov)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = cal ~ type, data = dog1.df)

$type
                   diff       lwr       upr     p adj
Meat-Beef      1.855882 -16.82550  20.53726 0.9688129
Poultry-Beef -38.085294 -56.76667 -19.40391 0.0000277
Poultry-Meat -39.941176 -59.36515 -20.51720 0.0000239

> plot(dog1.aov)

When I look at the QQ plot, it's quite "snaky".
I thought that doing a transformation (based on boxcox) would help
with the snakiness.

However, when I re-run with the transformation, it's still snaky.

Any suggestions, please?

Thanks,
Sincerely,
Erin


-- 
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: erinm.hodgess at gmail.com



More information about the R-help mailing list