I alle eksemplerne ovenfor omkring to normalfordelte observationssæt er de forskellige tests lavet ved at bruge python som en lommeregner. I python laves disse test med funktionen ttest$\text{\textunderscore}$ind fra scipy.stats-modulet. Denne funktion laver imidlertid ikke det tilhørende konfidensinterval, og jeg har derfor selv kodet en funktion ttest2. Funktionen ligger i pytFunktioner.py omtalt i afsnit 1.6. Python har ikke en indbygget funktion til at lave $F$-testet for hypotesen om ens varianser i Statistisk Model 6.1.1 med to normalfordelte observationssæt. Jeg har derfor selv kodet en funktion med navnet vartest2. Funktionen ligger i pytFunktioner.py. Hvis data ligger i to vektorer x og y, og output placeres i fUD, bliver kaldet De forskellige dele af output er som følger: Det 95%-konfidensinterval, der angives i output, er for forholdet mellem de to variansparametre, det vil sige for $\sigma_1^2/\sigma_2^2$ (dette konfidensinterval har jeg ikke omtalt ovenfor). Gå nu tilbage til Eksempel 6.4.3, og find de beregnede værdier der i output fra et kald af vartest2.

6.5.1 Teste to varianser ens

For at teste at middelværdierne er ens i to normalfordelinger, skal man enten bruge $t$-testet, hvis de to varianser er ens, eller også bruge Welchs test, hvis de to varianser ikke er ens. Begge de to tests udregnes med funktionen ttest2. Hvis data ligger i to vektorer x1 og x2, bliver kaldet i tilfældet, hvor varianserne antages ens: Hvis de to varianser antages forskellige skal man ændre True til False. De forskellige dele af output er som følger: Nedre og øvre grænse er for et 95%-konfidensinterval for forskel i middelværdi mellem gruppe 1 og gruppe2, det vil sige for parameteren $\delta=\mu_1-\mu_2.$ Gå nu tilbage til Eksempel 6.2.2 og Eksempel 6.3.2 og gentag beregningerne ved hjælp af ttest2.

6.5.2 Teste to middelværdier ens

I eksempel 6.2.2 så vi, at mængden af ekstraheret cell free DNA var næsten dobbelt så stor ved Triton-metoden som ved Qiagen-metoden. Den empiriske spredning ved Triton-metoden er også noget højere end den empiriske spredning ved Qiagen-metoden, selvom et formelt test for hypotesen om ens varianser vil give en $p$-værdi over 0.05. At middelværdi og spredning "følges ad" er ikke helt atypisk, når data vedrører en positiv størrelse (her mængde). I sådanne situationer vil der ofte ske det, at hvis data logaritmetransformeres, vil der efterfølgende være større lighed mellem varianserne. Lad os betegne logaritmen til mængden af ekstraheret CFDNA med henholdsvis $\text{logQia}_{i}$ og $\text{logTri}_{i}$ for den $i$'te prøve i de to grupper. Vi betragter Statistisk Model 6.1.1, her skrevet som hvor $\nu_j$ er middelværdien af logaritmen til mængden. Man kan matematisk vise sammenhængen $\mu_j=\exp(\nu_j+\frac{1}{2}\tau_j^2),$ hvor $\mu_j$ er middelværdien af mængden. I kodevinduet nedenfor laves der qqplots for de logaritmetransformerede data, og disse giver ikke anledning til at forkaste modellen. Først undersøges hypotesen om samme varians i de to grupper for de logaritmetransformerde værdier. Beregningen er vist i kodevinduet nedenfor: $F$-teststørrelsen er 1.3427, og $p$-værdien (to gange sandsynlighed for værdi større end 1.3427) fra en $F(13,19)$-fordeling er 0.54. Da $p$-værdien er langt over 0.05, siger vi, at data ikke strider mod samme varians på logaritmeskalaen.I kodevinduet laves der også et 95%-konfidensinterval for forskel i middelværdi mellem gruppe 1 og gruppe2, $\delta=\nu_1-\nu_2,$ under antagelsen om samme varians. Konfidensintervallet er baseret på $t(32)$-fordelingen, og bliver $[-0.897,\, -0.586].$ Vi kan oversætte dette konfidensinterval til et konfidensinterval for forholdet mellem middelværdierne på den oprindelige skala. For data omkring mængden af CFDNA giver dette resultat 95%-konfidensintervallet $[e^{-0.897},\,e^{-0.586}]=[0.41,\, 0.56].$ Her står, at med 95% sikkerhed er middelværdien af mængden af CFDNA ved Qiagen-metoden mellem 41% og 56% af middelværdien ved brug af Triton-metoden.

Beregninger i python

ForegåendeNæste