1 love does not come by demanding from others, but it is a self initiation
DESCRIPTION
3 MODEL: i = 1,..., R j = 1,..., C k= 1,..., n In general, n observations per cell, R C cells. Y ijk = i j ij ijkTRANSCRIPT
![Page 1: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/1.jpg)
1
Love does not come by demanding from others, but it is a self initiation.
![Page 2: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/2.jpg)
2
Two Factor DesignsConsider studying the impact of two factors on the yield (response):
Here we have R = 3 rows (levels of the Row factor), C = 4 (levels of the column factor), and n = 2 replicates per cell
[nij for (i,j)th cell if not all equal]
NOTE: The “1”, “2”,etc...mean Level 1, Level 2,etc..., NOT metric values
1 2 3 417.9, 18.1 17.8, 17.8 18.1, 18.2 17.8, 17.9
18.0, 18.2 18.0, 18.3 18.4, 18.1 18.1, 18.5
18.0, 17.8 17.8, 18.0 18.1, 18.3 18.1, 17.9
BRAND
123
DEVICE
![Page 3: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/3.jpg)
3
MODEL:
i = 1, ..., Rj = 1, ..., Ck= 1, ..., n
In general, n observations per cell, R • C cells.
Yijk = ijijijk
![Page 4: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/4.jpg)
4
the grand meanithe difference between the ith
row mean and the grand meanj the difference between the jth
column mean and the grand meanij the interaction associated with
the i-th row and the j-th columnijij
![Page 5: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/5.jpg)
5
Where Y••• = Grand mean
Yi•• = Mean of row iY•j• = Mean of column jYij• = Mean of cell (i,j)
[All the terms are somewhat “intuitive”, except for (Yij• -Yi•• - Y•j• + Y•••)]
Yijk = Y•••+ (Yi•• - Y•••) + (Y•j• - Y•••)+ (Yij• - Yi•• - Y•j• + Y•••)+ (Yijk - Yij•)
![Page 6: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/6.jpg)
6
The term (Yij• -Yi•• - Y•j• + Y•••) is more intuitively written as:
how a cellmean differs
from grand mean
adjustmentfor “row
membership”
adjustmentfor “column
membership”
We can, without loss of generality, assume (for a moment) that there is no error (random part); why then might the above be non-zero?
(Yij• - Y•••) (Yi•• - Y•••) (Y•j• - Y•••)
![Page 7: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/7.jpg)
7
ANSWER:
Two basic ways to look at interaction:
BL BH
AL 5 8AH 10 ?
If AHBH = 13, no interactionIf AHBH > 13, + interactionIf AHBH < 13, - interaction
- When B goes from BLBH, yield goes up by 3 (58).- When A goes from AL AH, yield goes up by 5 (510).- When both changes of level occur, does yield go up by the sum, 3 + 5 = 8?
Interaction = degree of difference from sum of separate effects
1)
“INTERACTION”
![Page 8: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/8.jpg)
8
2)
- Holding BL, what happens as A goes from AL AH? +5
- Holding BH, what happens as A goes from AL AH? +9If the effect of one factor (i.e., the impact of changing its level) is DIFFERENT for different levels of another factor, then INTERACTION exists between the two factors.
BL BH
AL 5 8AH 10 17
NOTE:- Holding AL, BL BH has impact + 3- Holding AH, BL BH has impact + 7
(AB) = (BA) or (9-5) = (7-3).
![Page 9: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/9.jpg)
9
(Yijk - Y•••) = (Yi•• - Y•••) + (Y•j• - Y•••)
+ [(Yij• - Yi••) - (Y•j• - Y•••)]
+ (Yijk - Yij•)
Going back to the (model) equation on page 4, and bringing Y... to the other side of the equation, we get
If we then square both sides, triple sum both sides over i, j, and k, we get, (after noting that all cross-product terms cancel):
Effect of column j at row i. Effect of column j
![Page 10: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/10.jpg)
10
TSS = SSBRows + SSBCols + SSIR,C+ SSWError
and, in terms of degrees of freedom,
R.C.n-1 = (R-1) + (C-1) + (R-1)(C-1) + R.C.(n-1); DF of Interaction = (RC-1)-(R-1)-(C-1) = (R-1)(C-1).
OR,
(Yijk - Y•••)n.C.Yi•• - Y•••
i j k i
+ n.R.Y•j• - Y•••)2
+ n.Yij• - Yi•• - Y•j• +Y•••
i j
(Yijk - Yij•
i j k
j
![Page 11: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/11.jpg)
11
17.9, 18.1 17.8, 17.8 18.1, 18.2 17.8,17.9
18.1 17.8 18.15 17.85
18.2, 18.0 18.0, 18.3 18.4, 18.1 18.1, 18.5
18.1 18.15 18.25 18.3
18.0, 17.8 17.8, 18.0 18.1, 18.3 18.1, 17.9
17.9 17.9 18.2 18.0
1 2 3 4
18.00 17.95 18.20 18.05
1
2
3
In our example:
DEV ICE
17.95
18.20
18.00
18.05
BRAND
![Page 12: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/12.jpg)
12
SSBrows =2 4[(17.95-18.05) 2 + (18.20-18.05) 2 + (18.0-18.05) 2]
= 8 (.01 + .0225 + .0025) = .28
SSBcol =2•3[(18-18.05) 2+(17.95-18.05) 2+(18.2-18.05) 2+( 18.05-18.05) 2]
= 6 (.0025 + .001 + .0225 + 0) = .21
SSIR,C = 2(18-17.95-18+18.05)2 + (17.8-17.95-17.95+18.05)2 ....… + (18-18-18.05+18.05)2
[]
= 2 [.055] = .11
SSW = (17.9-18.0) 2 + (18.1-18.0) 2 + (17.8-17.8) 2 + (17.8-17.8) 2 + …....... (18.1-18.0) 2 + (17.9-18.0) 2
= .30
TSS = .28 + .21+ .11 + .30 = .90
•
•
![Page 13: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/13.jpg)
13
FTV (2, 12) = 3.89 Reject Ho
FTV (3, 12) = 3.49 Accept Ho
FTV (6, 12) = 3.00 Accept Ho
1) Ho: All Row Means EqualH1: Not all Row Means Equal
2) Ho: All Col. Means EqualH1: Not All Col. Means Equal
3) Ho: No Int’n between factorsH1: There is int’n between factors
ANOVA
.05
SOURCE SSQ df M.S. FcalcRows .28 2 .14 5.6COL .21 3 .07 2.8Int’n .11 6 .0183 .73Error .30 12 .025
![Page 14: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/14.jpg)
14
![Page 15: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/15.jpg)
15
Minitab: Stat >> Anova >> General Linear Model
General Linear Model: time versus brand, device
Factor Type Levels Valuesbrand fixed 4 1, 2, 3, 4device fixed 3 1, 2, 3
Analysis of Variance for time, using Adjusted SS for Tests
Source DF Seq SS Adj SS Adj MS F Pbrand 3 0.21000 0.21000 0.07000 2.80 0.085device 2 0.28000 0.28000 0.14000 5.60 0.019brand*device 6 0.11000 0.11000 0.01833 0.73 0.633Error 12 0.30000 0.30000 0.02500Total 23 0.90000
S = 0.158114 R-Sq = 66.67% R-Sq(adj) = 36.11%
![Page 16: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/16.jpg)
16
![Page 17: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/17.jpg)
Assumption:
ijk follows N(0, 2) for all i, j, k, and they are independent.
17
![Page 18: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/18.jpg)
18
Test for Normality
1. Restore residuals when doing Anova.2. Stat >> Basic Statistics >> Normality Test
Mean -4.44089E-16StDev 0.1142N 24AD 0.815P-Value 0.030
Not really normal but not too far from normal.
![Page 19: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/19.jpg)
19
Test for Equal Variances
Minitab: Stat >> Anova >> Test for Equal Variances
Test for Equal Variances: time versus device, brand
95% Bonferroni confidence intervals for standard deviations
device brand N Lower StDev Upper 1 1 2 0.0463363 0.141421 49.6486 1 2 2 * 0.000000 * … 3 4 2 0.0463363 0.141421 49.6486
Bartlett's Test (normal distribution)Test statistic = 2.33, p-value = 0.993
* NOTE * Levene's test cannot be computed for these data.
![Page 20: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/20.jpg)
Fixed Effect Model
20
![Page 21: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/21.jpg)
Random Effect ModelAdditional assumptions:i follows N(0, 2
) for all i, and they are independent.j follows N(0, 2
) for all j, and they are independent.ij follows N(0, 2
) for all I, j, and they are independent.•All these random components i j , ij ijk are (mutually) independent.
21
![Page 22: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/22.jpg)
Mixed Effect Model (fixed rows and random columns)
22
![Page 23: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/23.jpg)
23
Fixed Random Mixed MSRows
cnn
cn
ncn
MSColRnnRnn
Rn
MSRC
nnn
MSError
Another issue:Table 17.17 (O/L 6th ed., p. 1057)
MEANSQUARE EXPECTATIONS
col = randomrow= fixed
Reference: Design and Analysis of Experiments by D.C. Montgomery, 4th edition, Chapter 11.
![Page 24: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/24.jpg)
24
Fixed: Specific levels chosen by the experimenterRandom: Levels chosen randomly from a large
number of possibilities
Fixed: All Levels about which inferences are to be made are included in the experiment
Random: Levels are some of a large number possible
Fixed: A definite number of qualitatively distinguishable levels, and we plan to study them all, or a continuous set of quantitative settings, but we choose a suitable, definite subset in a limited region and confine inferences to that subset
Random: Levels are a random sample from an infinite ( or large) population
![Page 25: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/25.jpg)
25
“In a great number of cases the investigator may argue either way, depending on his mood and his handling of the subject matter. In other words, it is more a matter of assumption than of reality.”
Some authors say that if in doubt, assume fixed model. Others say things like “I think in most experimental situations the random model is applicable.” [The latter quote is from a person whose experiments are in the field of biology].
![Page 26: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/26.jpg)
26
My own feeling is that in most areas of management, a majority of experiments involve the fixed model [e.g., specific promotional campaigns, two specific ways of handling an issue on an income statement, etc.] . Many cases involve neither a “pure” fixed nor a “pure” random situation [e.g., selecting 3 prices from 6 “practical” possibilities].
Note that the issue sometimes becomes irrelevant in a practical sense when (certain) interactions are not present. Also note that each assumption may yield you the same “answer” in terms of practical application, in which case the distinction may not be an important one.
![Page 27: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/27.jpg)
How to Fit these Models in Minitab
27
• “Balanced ANOVA” can fit restricted and unrestricted version. By default, it shows unrestricted model.
• “General Linear Model” can only fit unrestricted model.
• There are no difference between restricted or unrestricted versions for fixed effect and random effect model. It only matters for the mixed effect model.
![Page 28: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/28.jpg)
More on Minitab
28
• The notation in EMS output under restricted model matches with ours but it under unrestricted model is different to ours.
• General suggestion: use “General Linear Model” to fit the models BUT use “Balanced ANOVA, option of restricted model” to find the EMS for fixed and random effect models.
![Page 29: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/29.jpg)
29
Two-Way ANOVA in Minitab
Stat>>Anova>>General Linear Model:
Model device brand device*brand
Random factors
Results
Factor plots
Graphs
device
Tick “Display expected mean squares and variance components”
Main effects plots & Interactions plots
Use standardized residuals for plots
![Page 30: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/30.jpg)
30
General Linear Model: time versus device, brand
Factor Type Levels Values device random 3 1 2 3brand fixed 4 1 2 3 4
Analysis of Variance for time, using Adjusted SS for Tests
Source DF Seq SS Adj SS Adj MS F Pdevice 2 0.28000 0.28000 0.14000 7.64 0.022brand 3 0.21000 0.21000 0.07000 3.82 0.076device*brand 6 0.11000 0.11000 0.01833 0.73 0.633Error 12 0.30000 0.30000 0.02500Total 23 0.90000
![Page 31: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/31.jpg)
31
Exercise: Lifetime of a Special-purpose Battery
It is important in battery testing to consider different temperatures and modes of use; a battery that is superior at one temperature and mode of use is not necessarily superior at other treatment combination. The batteries were being tested at 4 different temperatures for three modes of use (I for intermittent, C for continuous, S for sporadic). Analyze the data.
![Page 32: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/32.jpg)
32
Battery Lifetime (2 replicates)
Mode of use 1 2 3 4
I 12, 16 15, 19 31, 39 53, 55
C 15, 19 17, 17 30, 34 51, 49
S 11, 17 24, 22 33, 37 61, 67
Temperature
![Page 33: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/33.jpg)
33
M FInteresting Example:*
Frontiersman
April
50 peopleper cell
Mean Scores “Frontiersman” “April” “Frontiersman” “April”Dependent males males females femalesVariables (n=50) (n=50) (n=50) (n=50)
Intent-to-purchase 4.44 3.50 2.04 4.52
(*) Decision Sciences”, Vol. 9, p. 470, 1978
Brand Name Appeal for Men & Women:
![Page 34: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/34.jpg)
34
1 2
1 2
2
3
4
gender
brandM
ean
Interaction Plot - Data Means for y
12Y
![Page 35: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/35.jpg)
35
ANOVA Results
Dependent Source d.f. MS FVariable
Intent-to- Sex (A) 1 23.80 5.61* purchase Brand name (B) 1 29.64 6.99**(7 pt. scale) A x B 1 146.21 34.48***
Error 196 4.24
*p<.05**p<.01
***p<.001
![Page 36: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/36.jpg)
36
Two Factors with No Replication,
When there’s no replication, there is no “pure” way to estimate ERROR.Error is measured by considering more than one observation (i.e., replication) at the same “treatment combination” (i.e., experimental conditions).
1 2 31 7 3 42 10 6 83 6 2 54 9 5 7
A
B
![Page 37: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/37.jpg)
37
Our model for analysis is “technically”:
Yij = i j + Iij
i = 1, ..., R
j = 1, ..., C
We can write:
Yij = Y•• + (Yi• - Y••) + (Y•j - Y••)
+ (Yij - Yi• - Y•j+ Y••)
![Page 38: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/38.jpg)
38
After bringing Y•• to the other side of the equation, squaring both sides, and double summing over i and j,
We Find:Yij - Y••)2 = C • Yi•-Y••)2
+ R • Y•j - Y••)2
+ (Yij - Yi• - Y•j +
Y••)2
R
i = 1
C
j=1
R
i=1
C
j=1
R
i=1
C
j=1
![Page 39: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/39.jpg)
39
TSS = SSBROWS + SSBCol + SSIR, C
R•C - 1 = (R - 1) + (C - 1) + (R - 1) (C - 1)Degrees of Freedom :
We Know, E(MSInt.) = VInt.
If we assume VInt. = 0, E(MSInt.) = 2,
and we can call SSIR,C SSW
MSInt MSW
![Page 40: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/40.jpg)
40
And, our model may be rewritten:
Yij = + i + j + ij,
and the “labels” would become:
TSS = SSBROWS+ SSBCol + SSWError
In our problem: SSBrows = 28.67
SSBcol = 32
SSW = 1.33
![Page 41: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/41.jpg)
41
Source SSQ df MSQ Fcalc
rows
col
Error
28.67
32.00
1.33
9.55
16.00
00.22
3
2
6
43
72
TSS = 62 11
at = .01,
FTV (3,6)
= 9.78
FTV(2,6)
= 10.93
ANOVAand:
![Page 42: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/42.jpg)
42
What if we’re wrong about there being no interaction?
If we “think” our ratio is,
in Expectation, 2 + VROWS , (Say, for ROWS) 2
and it really is (because there’s interaction)
2 + VROWS,2 + Vint’n
being wrong can lead only to giving us an underestimated Fcalc.
![Page 43: 1 Love does not come by demanding from others, but it is a self initiation](https://reader036.vdocument.in/reader036/viewer/2022070615/5a4d1bd97f8b9ab0599dbea0/html5/thumbnails/43.jpg)
43
Thus, if we’ve REJECTED Ho, we can feel confident of our conclusion, even if there’s
interaction
If we’ve ACCEPTED Ho, only then could the no interaction assumption be CRITICAL.