data ARL; input DMU $ 1-22 X1 X2 X3 Y1 Y2 Y3; label DMU='도서관' X1='총장서數' X2='노동력' X3='운영비용' Y1='총대출건수' Y2='도서관간 자료협조건수' Y3='자료관련 상담건수'; cards; Boston U 2129423 289 1842734 456571 35172 66185 Brigham Young 2500849 304 1821800 733720 45880 398718 Brown 2978970 241 3621081 272192 30727 63341 Case Western Reserve 1996479 172 1252339 179170 130145 83223 Chicago 6271045 336 3757486 724930 53253 69004 Columbia 7018408 516 3194047 798503 52160 150475 Cornell 6260779 562 4052863 1105680 49962 193650 Duke 4764033 359 3896238 522942 82546 137037 Emory 2442504 297 4560948 386109 53166 108293 George Washington 1848636 246 2148307 388245 72248 179404 Georgetown 2363799 262 2979320 530208 24942 147644 Howard 2372112 188 5229215 216351 11530 95900 Johns Hopkins 3275082 316 4001586 634903 44489 132058 MIT 2532175 217 1828445 564119 37099 79651 Miami 2165040 221 1615182 184067 58822 190782 Northwestern 3954204 351 2029635 366614 63729 123857 Notre Dame 2644486 237 2405703 375760 31819 56516 Pennsylvania 4672777 393 7544145 497286 56736 394240 Rochester 2992304 214 2223800 387131 43216 62975 Southern California 3480853 339 5844566 643330 24713 162706 Stanford 6997003 601 9023677 1111611 32426 114247 Tulane 2148660 203 2079857 357693 25513 57867 Vanderbilt 2512072 291 2152493 345503 28582 92528 Washington-St. Louis 3296358 287 5258647 412624 50201 123062 Yale 10108371 611 6827000 741583 34753 325216 ; run; proc iml; use ARL; read all var {dmu} into dmuname; read all var {X1 X2 X3} into X; read all var {Y1 Y2 Y3} into Y; reset storage=mlib.dea; call ccrie(Theta, Lambda, X, Y); ATheta = Theta[:]; /* 평균계산 */ print dmuname Theta[f=10.4] ATheta; n = nrow(X); file print; do k = 1 to n; dmunamek = dmuname[k]; put k 3. ' ' dmunamek 22. @; do kk = 1 to n; if Lambda[k,kk] > 0.0001 then do; Lambdakk = Lambda[k,kk]; put kk 3. '(' Lambdakk 5.2 ')' @; end; end; put; end;