Intervale de incredere, teste statistice
Intervale de incredere, teste statistice
Intervale de incredere, teste statistice
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