मॉडल में यादृच्छिक शब्दों को शामिल करना ग्रेड के लिए कुछ सहसंयोजक संरचना को प्रेरित करने का एक तरीका है। स्कूल के लिए यादृच्छिक कारक एक ही स्कूल के विभिन्न छात्रों के बीच एक गैर शून्य कोवरियन को प्रेरित करता है, जबकि स्कूल अलग होने पर यह होता है।0
चलो अपने मॉडल लिखने
जहां रों अनुक्रमित स्कूल और मैं अनुक्रमित छात्रों (प्रत्येक विद्यालय में)। शर्तों स्कूल रों स्वतंत्र यादृच्छिक एक में तैयार चर हैं एन ( 0 , τ ) । ई एस , मैं स्वतंत्र यादृच्छिक एक में तैयार चर हैं एन ( 0 , 2
Ys,i=α+hourss,iβ+schools+es,i
sischoolsN(0,τ)es,iN(0,σ2) ।
इस सदिश का मूल्य अपेक्षित है
[α+hourss,iβ]s,i
जो काम किया घंटे की संख्या से निर्धारित होता है।
के बीच सहप्रसरण और वाई एस ' , मैं ' है 0 जब रों ≠ रों 'Ys,iYs′,i′0s≠s′ , जिसका अर्थ है कि उम्मीद मूल्यों से ग्रेड के प्रस्थान स्वतंत्र जब छात्रों को एक ही स्कूल में नहीं हैं।
के बीच सहप्रसरण और वाई एस , मैं ' है τ जब मैं ≠ मैं ' , और का विचरण वाई एस , मैं है τ + σYs,iYs,i′τi≠i′Ys,iτ+σ2 : एक ही स्कूल से छात्रों के ग्रेड उनकी उम्मीद मूल्यों से प्रस्थान सहसंबद्ध होगा ।
उदाहरण और सिम्युलेटेड डेटा
यहां पांच स्कूलों से पचास छात्रों के लिए एक छोटी अनुसंधान अनुकरण है (यहाँ मैं ले ); चर के नाम स्व दस्तावेज हैं: σ2=τ=1
set.seed(1)
school <- rep(1:5, each=10)
school_effect <- rnorm(5)
school_effect_by_ind <- rep(school_effect, each=10)
individual_effect <- rnorm(50)
हम प्रत्येक छात्र के लिए उम्मीद ग्रेड, है कि मामले से प्रस्थान साजिश , एक साथ (बिंदीदार रेखा) प्रत्येक स्कूल के लिए मतलब प्रस्थान के साथ:schools+es,i
plot(individual_effect + school_effect_by_ind, col=school, pch=19,
xlab="student", ylab="grades departure from expected value")
segments(seq(1,length=5,by=10), school_effect, seq(10,length=5,by=10), col=1:5, lty=3)
schoolsα+hoursβ
इस उदाहरण के लिए विचरण मैट्रिक्स
schoolses,i
⎡⎣⎢⎢⎢⎢⎢⎢A00000A00000A00000A00000A⎤⎦⎥⎥⎥⎥⎥⎥
10×10AA=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢2111111111121111111111211111111112111111111121111111111211111111112111111111121111111111211111111112⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥.