• No results found

Verdeling en onafhankelijkheid van kwadratensommen in de variantie-analyse

N/A
N/A
Protected

Academic year: 2021

Share "Verdeling en onafhankelijkheid van kwadratensommen in de variantie-analyse"

Copied!
12
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

Verdeling en onafhankelijkheid van kwadratensommen in de

variantie-analyse

Citation for published version (APA):

Doornbos, R. (1990). Verdeling en onafhankelijkheid van kwadratensommen in de variantie-analyse. (Memorandum COSOR; Vol. 9029). Technische Universiteit Eindhoven.

Document status and date: Gepubliceerd: 01/01/1990

Document Version:

Uitgevers PDF, ook bekend als Version of Record

Please check the document version of this publication:

• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.

• The final author version and the galley proof are versions of the publication after peer review.

• The final published version features the final layout of the paper including the volume, issue and page numbers.

Link to publication

General rights

Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain

• You may freely distribute the URL identifying the publication in the public portal.

If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:

www.tue.nl/taverne

Take down policy

If you believe that this document breaches copyright please contact us at:

openaccess@tue.nl

(2)

EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computing Science

Memorandum COSOR 90-29

Verdeling en onafhankelijkheid van kwadratensommen in de

variantie-analyse

R. Doornbos

Eindhoven University of Technology

Department of Mathematics and Computing Science P.O. Box 513

5600 MB Eindhoven The Netherlands

Eindhoven, August 1990 The Netherlands

(3)

Verdeling en onafhankelijkheid van

kwadratensommen in de variantie-analyse.

R. Doornbos

o.

Inleiding

In de meeste boeken over variantie-analyse en/of lineaire modellen vinden we een of meer inleidende hoofdstukken over matrixrekening en over de verdeling van kwadra-tische vormen. Dit wordt dan later bij de variantie-analyse alleen toegepast bij mo-dellen met factoren met vast nivo's. In deze notitie laten we zien hoe stochastische en gemengde modellen de voor de F-toets nodige eigenschappen eveneens rechtstreeks volgen uit de bekende stellingen over kwadratische vormen. Ook het geneste model kan op de zelfde wijze worden behandeld. Wij beperken ons, om de notatie eenvoudig te houden, tot het geval van 2 factoren met herhalingen. We zullen voor de gebruikte stellingen en notaties verwijzen naar de dictaten Lineaire Modellen (dictaatnr. 2.221, najaarssemester 1981) en Kansrekening en Statistiek 2, bestemd voor Bdk (2S750), herfsttrimester 1986), die we resp. met LM en K&S zullen aangeven.

1. Enkele stellingen

Stelling 1 1!.. '" Np(tt, V), V regulier

y"'Ay" '" X~(tt'Att) {:} AV idempotent en r = rCA) .

(L.M. Stelling 4.11).

Speciaal geval:

Stelling 1a

y.. '"

Np(tt, (721)

y'Py Ip

P orthogonale projector,

rep)

= r, dan - 2- '"

X~(6),

6'= tt / .

(7 (7

Stelling 2

y.. '"

Np(tt, V)

y'"Ay"en y'"By"onafhankelijk {:} AVB = 0 . (L.M. Stelling 4.13*).

Speciaal geval:

Stelling 2a1!.. '" Np(tt, (721)

(4)

In het bijzonder: A = P1 , B = P2, orthogonale projectoren op onderling loodrechte

deelruimten.

Stelling 3

EUAu. = sp(AV)

+

JL'AJ.l (L.M. Eigenschap 2.36).

De toepassing van deze stellingen is zeer eenvoudig, zodra bij de verschillende modellen de variantie-covariantiematrixV is bepaald.

2. De tweevoudige variant ie-analyse

We beschouwen 4 gevallen:

a. Vaste nivo's

b. Stochastische nivo's

c. Gemengd model

d. Geneste proefopzet.

In aIle gevallen geldt de splitsing van kwadratensommen (vgl. K&S)

waarbij in het geneste geval, met B : A geldt:

De corresponderende splitsing van de vrijheidsgraden is

abn - 1= (a - 1)

+

(b - 1)

+

(a - 1)(b - 1)

+

ab( n - 1)

en dit komt meer overeen met de notatie in bruto kwadratensommen:

(5)

Zoals in L.M. stelt P:=

.!.uu'

de orthogonale projector voor op

u

:= {l, ... , l}', zodat n

Py =

{y, ..., y}'.

PA is de projector die middelt per nivo van de factorA, PB middelt per nivo van Ben

PAB per nivocombinatie van A en B. Als we voor de waarnemingvector Y als volgorde afspreken

YIll, ••• , Ylln, Y121, ••• , Y12n, ... , YIbl, ... , YIbn ,

Y2ll, ... , Y21n, ... , Yabl, ... , Yabn ,

B

A 1 2 b 1

-2

--

-a dan is

1---1~et WL-~

- --.

(~X~J

1---1

'/

1 PA=-nb

PB bestaat uit a x a blokken van de gedaante

(6)

1 na 1- - (

'_.1

h

·ld

1 PAB = -n

1- - .1

f - - 1

1- - -1

{- -1 \.

R

~

In andere woorden:

PA is de projectie op de deelruimte opgespannen dooravectoren, die elke constant zijn

(bv

=

1) op een nivo van A en 0 op de (a - 1) andere.

PB is de projectie op de deelruimte opgespannen door b vectoren, elk constant (i 0) op een van de bnivo's van B en = 0 op de andere. PAB is de projectie op de deelruimte opgespannen doorab vectoren, constant

(i

0) op een van de abnivocombinaties van A

en B en

=

0 op de andere combinaties.

De producten van P, PA, PB en PAB volgen het schema:

P PA PB PAB

P P P P P

PA P PA P PA

PB P P PB PB

PAB P PA PB PAB

(7)

Het is gemakkelijk te verifieren dat de onderlinge producten van de 4 leden van het rechterlid gelijk zijn aan O. Het zijn 4 projecties op onderling loodrechte deelruimten. Verder is

SSA = y'(PA - P)y SSB = y'(PB - P)y

SSR = y'(I - PAB)y

De modellen met de gebruikelijke veronderstellingen voor de 4 gevallen zijn te vinden in K&5.

Voor de verwachtingen van de gemiddelde kwadratensommen is het voor ons doel beter onderscheid te blijven maken tussen vaste en stochastische nivo's en

bV

-1-1:al niet

a - I

te vervangen door (7~.

Deze verwachtingen zijn dan:

vast stochastisch gemengd*J

GJ(A --1:a~a - Inb + 0'2 nb0'2 + n(72

+

(72 - - G'inb 1: 2

+

n(72AB + (72

l A AB a - I GJ(B - -na 1:/32. +(72 na(72 + (72 + (72 na(7~

+

0'2 b - 1 3 B AB GJ(AB n )2 2 nO'~B + (72 n(72

+

(72 (a_l)(b_l)1:(a/3 ij+(7 AB GJ(R (72 (72 (72 *) A vast, B stochastisch genest GJ(A a - Inb 1:ai2+ n(7B:A2 + (72 GJ(B:A n(7B:A2 + (72 GJ(R (72

Al deze verwachtingen zijn zonder moeite uit stelling 3 af te lei den. Hierbij levert J-t' AJ-t steeds de uidtrukkingen in 1:a~, 1:/3; en 1:(a/3)~j op, terwijl sp AV de variantie-componenten geeft.

a. Bij het model met vaste nivo's voIgt uit stelling Ia dat de vier kwadratensommen, gedeeld door (72, een x2-verdeling hebben.

(8)

De restkwadratensom heeft altijd een centrale x2-verdeling, de andere drie een niet- centrale, behalve onder de drie te toetsen nulhypothesen.

De onderlinge onafhankelijkheid voIgt direct uit stelling 2a.

b. Bij het model met stochastische nivo's is de variantie-covariantiematrix

De verdelingen van de kwadratensommen volgen nu uit stelling 1. Neem by.

Dus SSA/(nbCT~

+

nCT~B

+

CT2) '" X~-l'

Rier zijn dus alle vier de kwadratensommen, gedeeld door de verwachting van de gemiddelde kwadratensom, centraal X2 verdeeld.

De onafhankelijkheid van bnPA en PAB voIgt uit stelling 2 omdat

= (PA - P)(PAB - PA - PB

+

P) V

=

0 , want alle vermenigvuldigingen zijn commutatief.

SSA (a-l)(b-l) nCT~B+CT2 .

Dus F = - S S . ( ) . b 2 2 2 heeft een centrale F-verdelmg

-AB a - 1 n CTA

+

CTAB

+

CT

en derhalve isGKA / GKAB onder de hypotheseHo : CT~

=

0 eveneens centraal F-verdeeld.

c. Ret gemengde model

a-l

met als voorwaarde ~i(a[!)ij= 0 en var(af!)ij = --CT~B kan (zie K&S) ook

a

geschreven worden als

met var 2j = CTb en var(

al)ij

= CT~C en alle stochastische variabelen onafhanke-lijk. Ret verband tussen beide modellen is

(9)

_ ....2 - vAC

v

Met stelling 1 vinden we weer de verdelingen van de kwadratensommen, uit stel-ling 2 voIgt hun onderline onafhankelijkheid.

d. Inhet genest model

Y"1e = J.L

+

Qi

+

f3 '.'

+

~3'1e

-'3 -3·'

is de variantie-covariantiematrix

en

Verdelingen en onafhankelijkheid volgen opnieuw onmiddellijk uit de stellingen 1 en 2.

Eindhoven, augustus 1990

(10)

3:INDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computing Science

PROBABILITY THEORY, STATISTICS, OPERATIONS RESEARCH AND SYSTEMS

fHEORY

P.O. Box 513

5600MB Eindhoven - The Netherlands Secretariate: Dommelbuilding 0.03 felephone: 040 - 47 3130

~istof COSOR-memoranda - 1990

~umber Month Author Title

\190-01 January I.J.B.F. Adan Analysis of the asymmetric shortest queue problem J. Wessels Part 1: Theoretical analysis

W.H.M.Zijm

\190-02 January D.A. Overdijk Meetkundige aspecten van de productie van kroonwielen

\1190-03 February I.J.B.F.Adan Analysis of the assymmetric shortest queue problem 1. Wessels Part II: Numerical analysis

W.H.M.Zijm

\1190-04 March P. van der Laan Statistical selection procedures for selecting the best variety L.R. Verdooren

\190-05 March W.H.M.Zijm Scheduling a flexible machining centre E.H.L.B. Nelissen

\1190-06 March O.Schuller The design of mechanizations: reliability, efficiency and flexibility W.H.M.Zijm

\1190-07 March W.H.M.Zijm Capacity analysis of automatic transport systems in an assembly fac-tory

\190-08 March OJ.v. Houtum Computational procedures for stochastic multi-echelon production W.H.M.Zijm systems

(11)

Number Month Author Title

M 90-09 March P.J.M. van Production preparation and numerical control in PCB assembly Laarhoven

W.H.M.Zijm

M90-10 March F.A.W. Wester A hierarchical planning system versus a schedule oriented planning J. Wijngaard system

W.H.M.Zijm

M90-11 April A. Dekkers Local Area Networks

M90-12 April P. v.d. Laan On subset selection from Logistic populations

M90-13 April P. v.d. Laan De Van Dantzig Prijs

M 90-14 June P. v.d. Laan Beslissen met statistische selectiemethoden

M 90-15 June F.W. Steutel Some recent characterizations of the exponential and geometric distributions

M 90-16 June J. van Geldrop Existence of general equilibria in infinite horizon economies with C. Withagen exhaustible resources. (the continuous time case)

M 90-17 June

p.e.

Schuur Simulated annealing as a tool to obtain new results in plane geometry M 90-18 July F.W. Steutel Applications of probability in analysis

M 90-19 July !.J.B.F. Adan Analysis of the symmetric shortest queue problem J. Wessels

W.H.M.Zijm

M90-20 July I.J.B.F. Adan Analysis of the asymmetric shortest queue problem with threshold J. Wessels jockeying

W.H.M.Zijm

M 90-21 July K.vanHarn On a characterization of the exponential distribution F.W. Steutel

M90-22 July A. Dekkers Performance analysis of a volume shadowing model J. van der Wal

(12)

Number Month M 90-23 July M 90-24 July M 90-25 July M 90-26 July M 90-27 August M 90-28 August M 90-29 August Author A Dekkers 1. van der Wal

D.A Overdijk J. van Oorschot A Dekkers 1. van Oorschot A Dekkers D.A Overdijk AW.J. Kolen J.K. Lenstra R. Doornbos Title

Mean value analysis of priority stations without preemption

Benadering van de kroonwielflank met behulp van regeloppervlakken

inkroonwieloverbrengingen met grote overbrengvertlOuding

Cake, a concurrent Make CASE tool

Measuring and Simulating an 802.3 CSMA/CD LAN

Skew-symmetric matrices and the Euler equations of rotational motion for rigid systems

CombinatoricsinOperations Research

Verdeling en onafhankelijkheid van kwadratensommen in de variantie-analyse

Referenties

GERELATEERDE DOCUMENTEN

Deze dalingen worden natuurlijk niet door één factor veroorzaakt, andere factoren, waaronder het toeval, spelen hierbij zeker ook een rol · Toch kan een

Kwelmilieus komen voor waar grondwater uittreedt in het rivier- bed langs hoger gelegen gronden langs de Maas en IJssel of in de overgang van de gestuwde Utrechtse Heuvelrug naar

Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers) Please check the document version of this publication:.. • A submitted manuscript is

The composition of copolymers of vinylacetate (VAc) and vinylbutyrate (VB) can be determined rapidly and satisfactorily by dynamic thermogravimetric analysis (TGA), in

Naar aanleiding van dit commentaar hebben we besloten om toch een werkblad (werkblad 5.3) met behuip van VU- Grafiek te maken. Op deze manier leren de laptopleerlingen het

Omdat in dit geval de kwadratische afwij- kingen binnen de steekproeven relatief een groter deel bijdragen aan de totale kwadratische afwijking, zou men de nulhypothese

The breast cancer data set was analyzed by employing an extended Markov renewal illness-death model (surgery, tumor recurrence and death). Prediction results be- tween the

The communication between these components flows as follow: In order to be able to get notifications, the client application requests the notification component to subscribe