> i = c ( 1, 1, 2, 3, 3, 3, 4, 4, 5, 5, 5, 5 ) # practice ID
> m = c ( 4, 4, 4, 3, 3, 4, 4, 3, 6, 4, 5, 3 ) # patients
> y = c ( 1, 3, 3, 2, 3, 3, 4, 3, 1, 4, 1, 0 ) # patients with response=1
> 
> vc = vcbd(i, m, y)
> vcbd.print(vc)
[1] Variance Components for Triply-Nested Binary Responses (Version 1.0)
[1] (c) Bahjat Qaqish (2005)
[1] Department of Biostatistics, CB 7420
[1] Chapel Hill,  NC 27599-7420
[1] 
[1] Number of practices:    5
[1] Number of physicians:   12
[1] Number of patients:     47
[1] Mean response           0.595744680851064
[1] Mean response(weighted) 0.710838468066934
[1] Estimates:
[1]     Theta:     0.595744680851064
[1]     Lambda_1:  0.133753985259339
[1]     Lambda_2:  0.325346871318211
[1] Variance components:
[1]     Between practices:                  13.3753985259339 %
[1]     Between providers within practice:  28.1830430687743 %
[1]     Between patients  within provider:  58.4415584052918 %
[1] Correlation:
[1]     Between providers within practice:  0.133753985259339
[1]     Between patients  within provider:  0.415584415947082
[1] Odds ratio:
[1]     Between providers within practice:  1.74013052360563
[1]     Between patients  within provider:  6.05242366240168
> 
> 
> 
> 
> m = c (4, 4, 4, 3, 3, 4, 4, 3, 6, 4, 5, 3)
> y = c (0, 4, 2, 0, 3, 0, 0, 3, 0, 4, 0, 3)
> vc = vcbd(i, m, y, print_flag = T)  # prints right away
[1] Variance Components for Triply-Nested Binary Responses (Version 1.0)
[1] (c) Bahjat Qaqish (2005)
[1] Department of Biostatistics, CB 7420
[1] Chapel Hill,  NC 27599-7420
[1] 
[1] Number of practices:    5
[1] Number of physicians:   12
[1] Number of patients:     47
[1] Mean response           0.404255319148936
[1] Mean response(weighted) 0.420536458905587
[1] Estimates:
[1]     Theta:     0.371986851603132
[1]     Lambda_1:  0
[1]     Lambda_2:  0.927034102886675
[1] Variance components:
[1]     Between practices:                  0 %
[1]     Between providers within practice:  92.7034102886675 %
[1]     Between patients  within provider:  7.29658971133249 %
[1] Correlation:
[1]     Between providers within practice:  0
[1]     Between patients  within provider:  0.927034102886675
[1] Odds ratio:
[1]     Between providers within practice:  1
[1]     Between patients  within provider:  746.348757915868