الإصدار: 0.1.2 |التاريخ: 2016-10-22
Christopher T. Moore [aut, cre] “moor0554@umn.edu”
2016-10-30 22:16:14
أسامة عبدالعزيز الوسيدي “alwusaidio@rcyci.edu.sa”
عبدالرحمن النفيعي “aanofei@uqu.edu.sa”
library(gtheory)
# Loading Brennan.3.2 Data
data(Brennan.3.2)
head(Brennan.3.2)## Task Person Rater Score
## 1 1 1 1 5
## 2 1 2 1 9
## 3 1 3 1 3
## 4 1 4 1 7
## 5 1 5 1 9
## 6 1 6 1 3sapply(Brennan.3.2, class)## Task Person Rater Score
## "factor" "factor" "factor" "integer"# Loading Rajaratnam.2 Data
data(Rajaratnam.2)
head(Rajaratnam.2)## Person Subtest Item Score
## 1 1 1 1 4
## 2 2 1 1 2
## 3 3 1 1 2
## 4 4 1 1 1
## 5 5 1 1 3
## 6 6 1 1 1sapply(Rajaratnam.2, class)## Person Subtest Item Score
## "factor" "factor" "factor" "integer"#A univariate D study.
#Compare to results on page 116 of Brennan (2001).
data(Brennan.3.2)
formula.Brennan.3.2 <- "Score ~ (1 | Person) + (1 | Task) + (1 | Rater:Task) +
(1 | Person:Task)"
gstudy.out <- gstudy(data = Brennan.3.2, formula = formula.Brennan.3.2)
dstudy(gstudy.out, colname.objects = "Person", data = Brennan.3.2, colname.scores = "Score")## $components
## source var percent n
## 1 Person:Task 0.18648698 18.3 3
## 2 Rater:Task 0.05396961 5.3 12
## 3 Person 0.47297277 46.4 1
## 4 Task 0.10844755 10.6 3
## 5 Residual 0.19835839 19.4 12
##
## $var.universe
## [1] 0.4729728
##
## $generalizability
## [1] 0.5513672
##
## $var.error.rel
## [1] 0.3848454
##
## $sem.rel
## [1] 0.6203591
##
## $see.rel
## [1] 0.4606421
##
## $dependability
## [1] 0.4635919
##
## $var.error.abs
## [1] 0.5472625
##
## $sem.abs
## [1] 0.7397719
##
## $see.abs
## [1] 0.5036928
##
## attr(,"class")
## [1] "dstudy" "list"#A multivariate D study.
#Compare to results on pages 270-272 of Brennan (2001).
data(Rajaratnam.2)
formula.Rajaratnam.2 <- "Score ~ (1 | Person) + (1 | Item)"
gstudy.out <- gstudy(data = Rajaratnam.2, formula = formula.Rajaratnam.2,
colname.strata = "Subtest", colname.objects = "Person")
dstudy(gstudy.out, colname.objects = "Person", data = Rajaratnam.2, colname.scores = "Score",
colname.strata = "Subtest", weights = c(0.25, 0.5, 0.25))## $within
## $within$`1`
## $within$`1`$components
## source var percent n
## 1 Person 1.5714367 75.9 1
## 2 Item 0.2142874 10.3 2
## 3 Residual 0.2857130 13.8 2
##
##
## $within$`2`
## $within$`2`$components
## source var percent n
## 1 Person 2.78572919 89.8 1
## 2 Item 0.04985121 1.6 4
## 3 Residual 0.26785674 8.6 4
##
##
## $within$`3`
## $within$`3`$components
## source var percent n
## 1 Person 1.8571501 80.9 1
## 2 Item 0.2232052 9.7 2
## 3 Residual 0.2142860 9.3 2
##
##
##
## $between
## $between$var.obs
## 1 2 3
## 1 1.857143 1.4821429 0.5000000
## 2 1.482143 3.0535714 0.9464286
## 3 0.500000 0.9464286 2.0714286
## attr(,"colname.objects")
## [1] "Person"
##
## $between$var.universe
## 1 2 3
## 1 1.571437 1.4821429 0.5000000
## 2 1.482143 2.7857292 0.9464286
## 3 0.500000 0.9464286 1.8571501
## attr(,"colname.objects")
## [1] "Person"
##
## $between$generalizability
## 1 2 3
## 1 0.8461551 0.0000000 0.0000000
## 2 0.0000000 0.9122813 0.0000000
## 3 0.0000000 0.0000000 0.8965519
##
## $between$var.error.rel
## 1 2 3
## 1 0.285713 0.0000000 0.000000
## 2 0.000000 0.2678567 0.000000
## 3 0.000000 0.0000000 0.214286
##
## $between$sem.rel
## 1 2 3
## 1 0.5345213 0.0000000 0.0000000
## 2 0.0000000 0.5175488 0.0000000
## 3 0.0000000 0.0000000 0.4629104
##
## $between$see.rel
## 1 2 3
## 1 0.4916884 0.0000000 0.0000000
## 2 0.0000000 0.4943285 0.0000000
## 3 0.0000000 0.0000000 0.4383133
##
## $between$dependability
## 1 2 3
## 1 0.7586215 0.0000000 0.0000000
## 2 0.0000000 0.8976271 0.0000000
## 3 0.0000000 0.0000000 0.8093422
##
## $between$var.error.abs
## 1 2 3
## 1 0.5000004 0.0000000 0.0000000
## 2 0.0000000 0.3177079 0.0000000
## 3 0.0000000 0.0000000 0.4374912
##
## $between$sem.abs
## 1 2 3
## 1 0.7071071 0.0000000 0.0000000
## 2 0.0000000 0.5636559 0.0000000
## 3 0.0000000 0.0000000 0.6614312
##
## $between$see.abs
## 1 2 3
## 1 0.6158823 0.0000000 0.0000000
## 2 0.0000000 0.5340255 0.0000000
## 3 0.0000000 0.0000000 0.5950463
##
##
## $composite
## $composite$var.universe
## [,1]
## [1,] 1.580362
##
## $composite$generalizability
## [,1]
## [1,] 0.9414896
##
## $composite$var.error.rel
## [1] 0.09821412
##
## $composite$sem.rel
## [1] 0.3133913
##
## $composite$see.rel
## [,1]
## [1,] 0.3040848
##
## $composite$dependability
## [,1]
## [1,] 0.9196801
##
## $composite$var.error.abs
## [1] 0.1380202
##
## $composite$sem.abs
## [1] 0.3715107
##
## $composite$see.abs
## [,1]
## [1,] 0.3562786
##
##
## attr(,"class")
## [1] "dstudy" "list"#Conduct a univariate G study.
#Compare to results on page 116 of Brennan (2001).
data(Brennan.3.2)
formula.Brennan.3.2 <- "Score ~ (1 | Person) + (1 | Task) +
(1 | Rater:Task) + (1 | Person:Task)"
gstudy(data = Brennan.3.2, formula = formula.Brennan.3.2)## $components
## source var percent n
## 1 Person:Task 0.5594609 12.8 1
## 2 Rater:Task 0.6476353 14.8 1
## 3 Person 0.4729728 10.8 1
## 4 Task 0.3253427 7.4 1
## 5 Residual 2.3803007 54.3 1
##
## attr(,"class")
## [1] "gstudy" "list"#Conduct a multivariate G study.
#Compare to results on page 270 of Brennan (2001).
data(Rajaratnam.2)
formula.Rajaratnam.2 <- "Score ~ (1 | Person) + (1 | Item)"
gstudy(data = Rajaratnam.2, formula = formula.Rajaratnam.2, colname.strata = "Subtest",
colname.objects = "Person")## $within
## $within$`1`
## $within$`1`$components
## source var percent n
## 1 Person 1.5714367 61.1 1
## 2 Item 0.4285749 16.7 1
## 3 Residual 0.5714260 22.2 1
##
##
## $within$`2`
## $within$`2`$components
## source var percent n
## 1 Person 2.7857292 68.7 1
## 2 Item 0.1994048 4.9 1
## 3 Residual 1.0714269 26.4 1
##
##
## $within$`3`
## $within$`3`$components
## source var percent n
## 1 Person 1.8571501 68.0 1
## 2 Item 0.4464104 16.3 1
## 3 Residual 0.4285721 15.7 1
##
##
##
## $between
## $between$var.obs
## 1 2 3
## 1 1.857143 1.4821429 0.5000000
## 2 1.482143 3.0535714 0.9464286
## 3 0.500000 0.9464286 2.0714286
## attr(,"colname.objects")
## [1] "Person"
##
##
## attr(,"class")
## [1] "gstudy" "list"