Probabilidade de Cobertura dos Intervalos de Confiança ... - Uem

76
77 proc p r i n t t o ; run ;
78
79 data n u l l ;
80 s e t emvs ( keep=parameter estimate s t a n d a r d e r r o r ) ;
81 nda = &i n d i c e ;
82 f i l e ” c : \ emv b &n&p . . dat ” mod ;
83 put nda ” ; ” parameter ” ; ” estimate ” ; ” s t a n d a r d e r r o r ;
84 run ;
85
86 data n u l l ;
87 s e t adic ( keep=l a b e l estimate s t a n d a r d e r r o r ) ;
88 nda = &i n d i c e ;
89 f i l e ” c : \ a d i c b &n&p . . dat ” mod ;
90 put nda ” ; ” l a b e l ” ; ” estimate ” ; ” s t a n d a r d e r r o r ;
91 run ;
92 %mend bootstrap ;
93 %macro geracao (mu=, beta =,n=, seed =,p=);
94 data o r i g i n a l ( drop=i U) ;
95 do i = 1 to &n ;
96 U = uniform(&seed ) ;
97 time = &mu∗(− l o g (1−U))∗∗(1/& beta ) ;
98 d e l t a = ranbin(&seed ,1 ,&p ) ;
99 output ;
100 end ;
101 %mend geracao ;
102 %macro cobertura (mu=, beta =,p=,n=,nboot1=,nboot2 =);
103 %d e l e t a (n=&n , p=&p ) ;
104 %do j = 1 %to &nboot1 ;
105 %geracao (mu=&mu, beta=&beta , n=&n , seed=&j ∗1000 ,p=&p ) ;
106 %estimacao dados ( dados=o r i g i n a l ,mu=&mu, beta=&beta , i n d i c e=&j , n=&n , p=&p ) ;
107 %bootstrap ( inpt=o r i g i n a l , seed=&j , s a m p l e s i z e=&n , r e p l i c a c a o=&nboot2 ,mu=&mu,
108 beta=&beta , i n d i c e=&j , n=&n , p=&p ) ;
109 %end ;
110 %mend cobertura ;
Para a contagem dos intervalos de confiança que continham o verdadeiro valor do parâmetro
de interesse, foi utilizado a proc iml do SAS.
92