05.05.2014 Views

Intervale de incredere, teste statistice

Intervale de incredere, teste statistice

Intervale de incredere, teste statistice

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Laborator biostatistica<br />

Intervalul <strong>de</strong> incre<strong>de</strong>re<br />

Principalii parametri ce <strong>de</strong>scriu o populatie statistica sunt media si abaterea standard. Acesti parametri<br />

se estimeazã pe baza <strong>de</strong>terminãrilor efectuate pe esantioane din respectiva populatie. Parametrii<br />

probelor extrase nefiind perfect i<strong>de</strong>ntici cu cei ai populatiei studiate, apare necesitatea <strong>de</strong> a calcula<br />

intervalul în care se pot încadra acesti parametri - interval <strong>de</strong> incre<strong>de</strong>re, cu o anumitã probabilitate,<br />

aleasa în functie <strong>de</strong> exactitatea dorita (<strong>de</strong> obicei <strong>de</strong> 95% sau 99%). Aceasta înseamna ca, pentru un<br />

numar mare <strong>de</strong> probe din aceeasi populatie, 95% (sau 99%) din ele vor avea parametrii care se<br />

incadreaza in intervalul <strong>de</strong> incre<strong>de</strong>re si doar 5% (sau 1%) din ele se vor gasi în afara acestui interval.<br />

Se <strong>de</strong>finesc:<br />

- Eroarea Standard (ES) - arata cat <strong>de</strong> precis aproximeaza media calculata din valorile unei serii, media<br />

populatiei din care a fost extras lotul pe care s-au facut masuratorile.<br />

ES =<br />

<br />

n<br />

un<strong>de</strong>:<br />

- abaterea standard<br />

n - numarul <strong>de</strong> valori din serie<br />

- Prag <strong>de</strong> semnificatie (q) - probabilitatea ca media sa fie in afara intervalului <strong>de</strong> incre<strong>de</strong>re<br />

q = 100 P<br />

un<strong>de</strong>: P - probabiliatea ca media sa fie in interiorul intervalului <strong>de</strong> incre<strong>de</strong>re<br />

- Intervalul <strong>de</strong> incre<strong>de</strong>re:<br />

(M - t x ES).....M......(M + t x ES)<br />

un<strong>de</strong>: M - media calculata pentru seria respectiva<br />

ES - eroarea standard calculata pentru seria respectiva<br />

t - testul <strong>de</strong> semnificatie (a se ve<strong>de</strong>a tabelul urmator)<br />

Exemplu <strong>de</strong> calcul:<br />

Glicemia masurata pe un esantion <strong>de</strong> 90 persoane are o medie <strong>de</strong> 110 mg/dl. Abaterea standard<br />

corespunzatoare este <strong>de</strong> 14 m/dl.<br />

Sa se calculeze intervalul <strong>de</strong> incre<strong>de</strong>re corespunzator unei probabilitati <strong>de</strong> 95%.<br />

Sa se calculeze intervalul <strong>de</strong> incre<strong>de</strong>re corespunzator unui prag <strong>de</strong> semnificatie <strong>de</strong> 1%.<br />

Rezolvare:<br />

1 - se calculeaza eroarea standard:<br />

Abaterea standard:<br />

14<br />

Media:<br />

M 110<br />

Numarul <strong>de</strong> elemente din serie:<br />

n 90<br />

Numarul <strong>de</strong> gra<strong>de</strong> <strong>de</strong> libertate (n-1):<br />

<br />

Eroarea standard: ES <br />

n<br />

G n 1<br />

ES 1.476<br />

G 89<br />

2 - se alege valoarea testului <strong>de</strong> semnificatie din tabel, pentru numarul <strong>de</strong> gra<strong>de</strong> <strong>de</strong> libertate si pragul <strong>de</strong><br />

semnificatie corespunzatoare.<br />

Pentru primul caz:<br />

P 95<br />

<strong>de</strong>ci<br />

q 100 P<br />

q 5<br />

se obtine din tabel:<br />

t 1 1.986


Pentru cel <strong>de</strong>-al doilea caz: q 1<br />

t 2 2.631<br />

se obtine din tabel:<br />

Tabel cu valorile testului <strong>de</strong> semnificatie "t"


3 - se calculeaza limitele intervalului <strong>de</strong> incre<strong>de</strong>re:<br />

Pentru primul caz<br />

limita inferioara:<br />

L i M t 1 ES<br />

L i 107.069<br />

limita superioara:<br />

L s M t 1 ES<br />

L s 112.931<br />

latimea intervalului <strong>de</strong> incre<strong>de</strong>re:<br />

I c L s L i<br />

I c 5.862<br />

Pentru al doilea caz<br />

limita inferioara:<br />

L i M t 2 ES<br />

L i 106.117<br />

limita superioara:<br />

L s M t 2 ES<br />

L s 113.883<br />

latimea intervalului <strong>de</strong> incre<strong>de</strong>re:<br />

I c L s L i<br />

I c 7.765<br />

Exercitii <strong>de</strong> rezolvat:<br />

1. Sa se calculeze intervalul <strong>de</strong> incre<strong>de</strong>re 99% pentru masuratorile tensiunii sistolice, facute pe un<br />

esantion <strong>de</strong> 176 <strong>de</strong> pacienti, stiind ca media calculata este m = 125 mmHg iar abaterea standard este<br />

<strong>de</strong> 18 mmHg.<br />

2. Sa se calculeze intervalul <strong>de</strong> incre<strong>de</strong>re 95% pentru masuratorile continutului <strong>de</strong> calciu facute pe un<br />

esantion <strong>de</strong> 241 <strong>de</strong> pacienti, stiind ca media calculata este m = 120 mg/dl iar abaterea standard este <strong>de</strong><br />

16 mg/dl.<br />

Testul Grubbs<br />

Teste <strong>statistice</strong><br />

Permit verificarea valorii unei experiment în functie <strong>de</strong> probabilitatile aparitiei variatiei in natura. Pot<br />

preciza daca un anumit rezultat experimental poate fi consi<strong>de</strong>rat ca fiind datorat variabilitatii biologice<br />

normale si care nu va in<strong>de</strong>parta semnificativ rezultatul final al experimentului <strong>de</strong> “valoarea a<strong>de</strong>varata”<br />

sau este un rezultat aberant care se datoreaza unor factori acci<strong>de</strong>ntal aparuti in timpul experimentului.<br />

Testul consta in calcularea valorii Z si compararea acesteia cu valoarea admisa:<br />

Z =<br />

<br />

X <br />

<br />

Xi<br />

Conditia <strong>de</strong> verificare:<br />

Z Z admisibil<br />

un<strong>de</strong>:<br />

- abaterea standard<br />

<br />

X - valoarea medie calculata cu toate<br />

elementele din serie<br />

X i - valoarea suspecta<br />

n - numarul <strong>de</strong> elemente din serie<br />

Z admisibil - valoarea admisibila (vezi<br />

tabelul urmator)


Tabel cu valorile corespunzatoare testului Grubbs<br />

Exemplu <strong>de</strong> calcul:<br />

Sa se verifice valoarea <strong>de</strong> 160 kg din seria <strong>de</strong><br />

masuratori din tabelul alaturat.<br />

Rezolvare:<br />

1 - se introduc datele intr-o foaie <strong>de</strong> lucru Excel si<br />

se importa in LazStat<br />

2 - se calculeaza media = 73.93 kg si abaterea<br />

standard = 27.59 kg<br />

3 - se calculeaza Z si se alege Z admisibil din tabel<br />

Z <br />

73.93 160<br />

27.59<br />

Z admisibil 2.55


4 - se face verificarea:<br />

valoareaberanta "da" if Z Z admisibil<br />

"nu"<br />

otherwise<br />

valoareaberanta "da"<br />

Testul Stu<strong>de</strong>nt<br />

Testul Stu<strong>de</strong>nt este utilizat în analiza statisticã pentru compararea mediei unei caracteristici la doua<br />

populatii. Caracteristica studiata trebuie sã fie o caracteristicã cantitativa, masurabila.<br />

In cazul compararii a doua medii <strong>statistice</strong> se pot face doua ipoteze<br />

- Ipoteza H 0 : media µ1 a caracteristicii în populatia din care face parte grupul <strong>de</strong> test este egala cu<br />

media µ2 a caracteristicii în populatia din care face parte grupul martor (µ1 = µ2).<br />

- Ipoteza H 1 : media µ1 a caracteristicii în populatia din care face parte grupul <strong>de</strong> test este diferita <strong>de</strong><br />

media µ2 a caracteristicii în populatia din care face parte grupul martor ( 1 2 )”<br />

Formula <strong>de</strong> calcul a testului Stu<strong>de</strong>nt este:<br />

t =<br />

2<br />

s p =<br />

<br />

X Y<br />

2 2<br />

s p s p<br />

<br />

n 1 n 2<br />

n 1 1 1<br />

2<br />

n 2 1<br />

n 1 n 2 2<br />

un<strong>de</strong>:<br />

X - reprezintã media caracteristicii in grupul <strong>de</strong> test<br />

<br />

Y - reprezintã media caracteristicii in grupul martor<br />

n 1 - numarul <strong>de</strong> masuratori din grupul <strong>de</strong> test<br />

n 2 - numarul <strong>de</strong> masuratori din grupul martor<br />

1 - abaterea standard din grupul <strong>de</strong> test<br />

- abaterea standard din grupul martor<br />

Stabilirea concluziei testului Stu<strong>de</strong>nt:<br />

- Dacã probabilitatea P 0.05 , atunci se repinge ipoteza H 0 si se acceptã ipoteza H 1 , prin urmare<br />

exista o diferenta semnificativa intre mediile caracteristicii in cele doua populatii.<br />

- Daca probabilitatea P > 0.05 , atunci se accepta ipoteza H 0 , prin urmare nu exista o diferenta<br />

semnificativa intre mediile caracteristicii in cele doua populatii.<br />

Testul Stu<strong>de</strong>nt in LazStat<br />

Se activeaza prin comanda: Analyses - Comparisons - t-tests. Se <strong>de</strong>schi<strong>de</strong> urmatoarea fereastra:


Se pot face urmatoarele setari:<br />

"Data Entry By:" se poate opta pentru introducerea valorilor direct in fereastra sau din foia <strong>de</strong> calcul<br />

"Test Assumptions:" se poate opta pentru variabile in<strong>de</strong>pen<strong>de</strong>nte sau <strong>de</strong>pen<strong>de</strong>nte (se presupune ca<br />

populatiile au aceeasi dispersie).<br />

"Test Probability:" se poate opta pentru un test cu un capat (One-tailed) - in acest caz, in cazul<br />

respingerii ipotezei H 0 se indica media cea mai mare. La testul cu doua capete, in cazul respingerii<br />

ipotezei H 0 se consi<strong>de</strong>ra ca exista diferente intre mediile celor doua caracteristici, fara a se indica<br />

media cea mai mare. In general se utilizeaza testul cu doua capete.<br />

Exemplu <strong>de</strong> calcul:<br />

Intr-un studiu al efectului bumetami<strong>de</strong>i in secretia <strong>de</strong><br />

calciu in urina, 9 persoane alese aleator au primit<br />

fiecare cate o doza <strong>de</strong> 0,5 mg <strong>de</strong> medicament. S-a<br />

colectat în fiecare ora, timp <strong>de</strong> 6 ore, urina <strong>de</strong> la cele 9<br />

persoane. La fel s-a procedat cu alte 10 persoane care<br />

nu au primit medicamentul. Pentru fiecare persoana<br />

s-a calculat o medie (prin calculul mediei celor 6 valori<br />

citite). Datele obtinute sunt prezentate in tabelul<br />

alaturat.<br />

Sa se <strong>de</strong>termine daca secretia <strong>de</strong> calciu in urina difera la cele doua grupuri, <strong>de</strong>ci daca administrarea<br />

medicamentului are efect in cresterea secretiei <strong>de</strong> calciu. Se introduc datele intr-o foaie <strong>de</strong> calcul<br />

Excel, se importa in LazStat si se aplica testul Stu<strong>de</strong>nt <strong>de</strong> tipul 2 (Correlated scores) si cu un capat<br />

(One-tailed).<br />

Ipoteze:<br />

- H 0 : medicamentul nu are efect in cresterea secretiei <strong>de</strong> calciu.<br />

-H 1 : medicamentul are efect in cresterea secretiei <strong>de</strong> calciu.<br />

- P > 0,05 - se accepta ipoteza H 0, P < 0.05 - se accepta ipoteza H 1<br />

Se obtin urmatoarele date:<br />

***********************************************************************************************************************<br />

COMPARISON OF TWO MEANS<br />

Variable Mean Variance Std.Dev. S.E.Mean N<br />

VAR 1 7.50 36.56 6.05 2.02 9<br />

VAR 2 6.17 3.13 1.77 0.59 9<br />

Assuming <strong>de</strong>pen<strong>de</strong>nt samples, t = 0.748 with probability = 0.4761 and 8 <strong>de</strong>grees of freedom<br />

Correlation between VAR 1 and VAR 2 = 0.517<br />

Difference = 1.33 and Standard Error of difference = 1.78<br />

Confi<strong>de</strong>nce interval = ( -2.78, 5.45)<br />

t for test of equal variances = 4.836 with probability = 0.0019<br />

**********************************************************************************************************************<br />

Rezultate: P = 0,4761 - se accepta ipoteza H 0 , <strong>de</strong>ci medicamentul nu are efect in cresterea secretiei <strong>de</strong><br />

calciu.


Exercitii <strong>de</strong> rezolvat<br />

Exercitiul 1<br />

Se efectueaza un studiu al nivelului <strong>de</strong> digoxin ser, dupa<br />

efectuarea rapida a unei injectii intravenoase cu acest<br />

medicament. Sa se stabileasca daca nivelul <strong>de</strong> digoxin<br />

ser la 4 ore dupa injectare difera semnificativ <strong>de</strong> nivelul<br />

<strong>de</strong> la 8 ore dupa injectare. Datele obtinute în urma<br />

studiului pe 10 subiecti sunt prezentate in tabelul alaturat.<br />

Pentru a obtine rezultatul studiului, aplicati testul Stu<strong>de</strong>nt<br />

cu doua capete si pentru grupuri <strong>de</strong>pen<strong>de</strong>nte. Daca<br />

probabilitatea P obtinuta este mai mica <strong>de</strong>cât 0,05 atunci<br />

exista diferente semnificative.<br />

Exercitiul 2<br />

Concentratia hemoglobinei în g/100 ml sânge, la un numar<br />

<strong>de</strong> 12 persoane cu anemie<br />

feripriva, a crescut dupa tratament conform tabelului<br />

alaturat.<br />

Se poate afirma ca tratamentul este eficace? Pentru a<br />

putea raspun<strong>de</strong> la aceasta intrebare utilizati testul<br />

STUDENT cu doua capete si pentru grupuri <strong>de</strong>pen<strong>de</strong>nte.<br />

Exercitiul 3<br />

S-a masurat glicemia la un lot <strong>de</strong> 5 persoane sanatoase,<br />

alese aleator. Apoi s-a<br />

masurat glicemia la un lot <strong>de</strong> 8 persoane alese <strong>de</strong> asemenea<br />

aleator, dar bolnave <strong>de</strong> diabet<br />

zaharat. Rezultatele obtinute sunt prezentate în tabelul<br />

alaturat<br />

Sa se stabileasca daca mediile celor doua loturi difera<br />

semnificativ, cu un risc <strong>de</strong> 0,05.<br />

Se va utiliza testul Stu<strong>de</strong>nt cu doua capete, variabile<br />

in<strong>de</strong>pen<strong>de</strong>nte.


Exercitiul 4<br />

S-a masurat urimia la doua loturi <strong>de</strong> câte 10 bolnavi <strong>de</strong> guta,<br />

dintre care unii au fost<br />

tratati si altii nu si s-au obtinut rezultatele prezentate alaturat.<br />

Sa se aprecieze cu un risc <strong>de</strong> 0,05 daca medicamentul a<br />

avut efect. Se va aplica testul STUDENT cu doua capete,<br />

variabile in<strong>de</strong>pen<strong>de</strong>nte.<br />

Exercitiul 5<br />

S-a masurat valoarea sistolica la grup <strong>de</strong> pacienti<br />

dignosticati cu stenoza si la un grup <strong>de</strong> pacienti<br />

asimptomatici. Sa se <strong>de</strong>termine daca media celor doua serii<br />

<strong>de</strong> date difera<br />

semnificativ. Se va utiliza testul STUDENT cu doua capete,<br />

pentru grupruri in<strong>de</strong>pen<strong>de</strong>nte).<br />

Rezultatele masuratorilor sunt exemplificate în tabelul<br />

alaturat

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!