[18616] in s-news-athena

home help back first fref pref prev next nref lref last post

[S] QQplot for two two unequal samples

daemon@ATHENA.MIT.EDU (ziv shkedy)
Tue Jul 13 09:35:08 1999

Date: Tue, 13 Jul 1999 15:32:20 +0200 (MET DST)
Message-Id: <199907131332.PAA00690@alpha.luc.ac.be>
Mime-Version: 1.0
Content-Type: text/plain; charset="us-ascii"
To: s-news@wubios.wustl.edu
From: ziv shkedy <ziv.shkedy@luc.ac.be>


Dear fellow S+ users,
I'm using S+ 4.5 for windows and my question is about QQplot for two samples
with unequal size.
I'm using the S+ fuel.frame dataset and I want to plot the QQplot for the
"Sporty" and the "Compact" car.
I used the following code and QQplot that I got contains 9 points. There are
9 observations in the "Sporty" group and 15 in the "Compact" group.
Does anybody know how the qq() work in this case? Which observation were
deleted for the "Compact" group" ?
You can find below the S+ code with sample size and sorted data.
Thanks in advance, Ziv.

> attach(fuel.frame)

> qq(Type~Weight,data=fuel.frame,aspect=1,                        # this
gives QQplot with 9 points
+           subset = (Type == "Sporty")|(Type == "Compact"))




######## SAMPLE SIZE  #############

> length(Weight[Type == "Compact"])
[1] 15

> length(Weight[ Type == "Sporty"])
[1] 9



#########The SORTED DATA###########

> sort(Weight[ Type == "Sporty"])
 Honda Civic CRX Si 4 Subaru XT 4 Ford Probe Honda Prelude Si 4WS 4 Nissan
240SX 4 Plymouth Laser Dodge Daytona Ford Mustang V8 
                 2170        2485       2695                   2710
2775           2840          2885            3310
 Chevrolet Camaro V8 
                3320
> sort(Weight[ Type == "Compact"])
 Peugeot 405 4 Buick Skylark 4 Oldsmobile Calais 4 Chevrolet Beretta 4 Audi
80 4 Mitsubishi Galant 4 Ford Tempo 4 Mazda 626 4 Honda Accord 4 
          2575            2640                2645                2655
2670                2745         2750        2780           2920
 Nissan Stanza 4 Toyota Camry 4 Subaru Legacy 4 Volvo 240 4 Chrysler Le
Baron V6 Mitsubishi Sigma V6 
            2920           2920            2935        2985
3065                3110
=========================================================
Ziv Shkedy
Biostatistics                   
Center for Statistics           
Limburgs Universitair Centrum
Universitaire Campus, department WNI
B-3590 Diepenbeek, Belgium
Tel: +32-(0)11-26.82.57
e-mail: ziv.shkedy@luc.ac.be 
=========================================================

-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu.  To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message:  unsubscribe s-news

home help back first fref pref prev next nref lref last post