- (16) 3376-2047
- [email protected]
- Portfólio de Serviços
- AT
O interesse dos pesquisadores em estudos cross-over é saber o quanto as formulações são semelhantes, diferentes ou equivalentes. A análise estatística proposta por Westlake (1972 e 1976), Metzler (1974) é a utilização de intervalos de confiança para a diferença e razão entre as médias das formulações. Nesta seção, trataremos sobre intervalos de confiança para a média e para a diferença de médias.
Seja o modelo (desprezando o efeito carry-over): $$ Y_{ijk}=\mu+S_{ik}+\pi_j+\tau_d{}_{[i,j]}+\varepsilon_{ijk}\left\{\begin{array}{cc}i=1,\cdots,s\quad\mbox{(Sequência)}\\j=1,\cdots,p\quad\mbox{(Período)}\\k=1,\cdots,n_i\quad\mbox{(Indivíduo)}\end{array}\right.$$
Considerando um planejamento cross-over, sem réplicas, para dois medicamentos (T = teste; R = referência), temos que cada indivíduo é aleatoriamente alocado para a sequência RT ou TR em dois períodos de dosagem. Isto é, indivíduos alocados na sequência RT (TR) recebem formulação R (T) no primeiro período de dosagem e formulação T (R) no segundo período de dosagem.
Seja:
Então: $$E(\bar{Y}_R)=\frac{1}{2}\left(2\mu+\pi_1+\pi_2+2F_R\right)$$
$$var(\bar{Y}_R)=\frac{(\sigma^2_S+\sigma^2)}{4}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)$$
$$E(\bar{Y}_T)=\frac{1}{2}\left(2\mu+\pi_1+\pi_2+2F_T\right)$$
$$var(\bar{Y}_T)=\frac{(\sigma^2_S+\sigma^2)}{4}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)$$
Considerando $\bar{Y}_l$ para $l=R$ ou $T$ a média de interesse, então, sabemos que $\quad\frac{\bar{Y}_l-E(\bar{Y}_l)}{var(\bar{Y}_l)}\sim N(0,1).$
Sendo $E(\bar{Y}_l)=\mu_l$ e $(\sigma^2_S+\sigma^2)$ estimado por $QM_{\mbox{intra}}$ ,então: $$T=\frac{\bar{Y}_l-\mu_l}{\sqrt{\frac{QM_{\mbox{intra}}}{4}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)}}\sim t_{n_1+n_2-2}.$$
Sabemos que $$\bar{Y}_T\sim N\left(2\mu+\pi_1+\pi_2+2F_T;\frac{1}{4}(\sigma^2_S+\sigma^2)\left(\frac{1}{n_1}+\frac{1}{n_2}\right)\right)$$ e
$$\bar{Y}_R\sim N\left(2\mu+\pi_1+\pi_2+2F_R;\frac{1}{4}(\sigma^2_S+\sigma^2)\left(\frac{1}{n_1}+\frac{1}{n_2}\right)\right)$$
Nesse sentido, $$\bar{Y}_T-\bar{Y}_R\sim N\left(F_T-F_R;\frac{1}{2}(\sigma^2_S+\sigma^2)\left(\frac{1}{n_1}+\frac{1}{n_2}\right)\right)$$
Temos então que, $$Z=\frac{\left(\bar{Y}_T-\bar{Y}_R\right)-\left(F_T-F_R\right)}{\sqrt{\frac{(\sigma^2_S+\sigma^2)}{2}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)}}\sim~N(0,1).$$ Vamos considerar que $(\sigma^2_S+\sigma^2)$ pode ser estimado por $QM_{\mbox{intra}}$. Com isso,
$$U=\frac{(n_1+n_2-2)QM_{\mbox{intra}}}{(\sigma^2_S+\sigma^2)}\sim\chi^2(n_1+n_2-2)$$
Sendo Z e U independentes, então $$T=\frac{Z}{\sqrt{\frac{U}{n_1+n_2-2}}}=\frac{\left(\bar{Y}_T-\bar{Y}_R\right)-\left(F_T-F_R\right)}{\sqrt{\frac{QM_{\mbox{intra}}}{2}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)}}\sim t(n_1+n_2-2)$$
A partir disso, fixando um valor $\alpha$ tal que, $0 \ \textless \ 1-\alpha \ \textless \ 1$, podemos encontrar um valor $t_{\frac{\alpha}{2}}$ tal que $P\left(-t_{\frac{\alpha}{2}} \ \textless \ T \ \textless \ t_{\frac{\alpha}{2}}\right)=1-\alpha$.
Assim, o intervalo de confiança para $(F_T-F_R)$, com nível de significância $\alpha$, é dado por:
$\left[\left(\bar{Y}_T-\bar{Y}_R\right)-t_{\frac{\alpha}{2}}\sqrt{\frac{QM_{\mbox{intra}}}{2}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)} \ \textless \ \left(F_T-F_R\right) \ \textless \ (\bar{Y}_T-\bar{Y}_R)+t_{\frac{\alpha}{2}}\sqrt{\frac{QM_{\mbox{intra}}}{2}\left(\frac{1}{n_1}+\frac{1}{n_2}\right)}\right]\label{ICdiff}.$
Voltando ao Exemplo 1.2.1.
clique aqui para efetuar o download dos dados utilizados nesse exemplo
Calculamos o intervalo de confiança para $\mu_l$.
Para $l = R, T$ temos que:
Sabemos que $\bar{y}_R=5,7748$,$n_1=12$,$n_2=12$,$QM_{\mbox{intra}}=0,005842685$ e para $\alpha=10\%$ temos que $t_{(\frac{\alpha}{2},22)}=1,7171$. Então, o intervalo é $$\left(5,7480;5,8016\right)$$
Sabemos que $\bar{y}_T=5,7873$,$n_1=12$,$n_2=12$,$QM_{\mbox{intra}}=0,005842685$ e para $\alpha=10\%$ temos que $t_{(\frac{\alpha}{2},22)}=1,7171$. Então, o intervalo é $$\left(5,7605;5,8141\right)$$
Agora calculamos o intervalo de confiança para $(F_T-F_R)$.
Sabemos que $\bar{y}_T=5,7873$,$\bar{y}_R=5,7748$,$n_1=12$,$n_2=12$,$QM_{\mbox{intra}}=0,005842685$ e para $\alpha=10\%$ temos que $t_{(\frac{\alpha}{2},22)}=1,7171$. Então, o intervalo é$$\left(-0,0254;0,0504\right)$$
O Portal Action é mantido pela Estatcamp - Consultoria Estatística e Qualidade, com o objetivo de disponibilizar uma ferramenta estatística em conjunto com uma fonte de informação útil aos profissionais interessados.