13 Gentagne målinger
Eksempel: Calcium uge 10.
13.1 Spaghettiplot: Individuelle profiler
Spaghettiplot laves ved brug af menuen Graph/Chart Builder/Line og valg af ikonet til højre (Multiple Line). Herefter klikkes der på fanen Groups/Point ID og der skal sættes flueben i Columns panel variable, hvorved der kommer et felt i øverste venstre hjørne med navnet Panel?. Sæt nu grp
i Panel?, bmd
på Y-aksen, visit
på X-aksen og girl
i Set color-boksen (kræver, at girl
gøres Nominal - dette gøres i Variable View-vinduet under kolonnen Measure).
NB: I dette plot er det tilsyneladende ikke muligt at markere nogle profiler i andre farver (f.eks. de, som ikke har fuldstændig follow up).
13.2 Figur af gennemsnitskurver
Benyt menuen Graph/Chart Builder/Line og vælg ikonet til højre (Multiple Line). Sæt nu bmd
på Y-aksen, visit
på X-aksen og grp
i Set color-boksen. Ovre i boksen Element Properties til højre i vinduret, vælges under Statistics Mean. Akserne bliver muligvis tossede, da de er afpassede efter selve observationerne, men dette kan afhjælpes som beskrevet i Lenes slides s. 100.
13.3 Mixed Model
Benyt menuen Analyze/Mixed Models/Linear sæt girl
over i Subjects og klik Continue.
Herefter sættes bmd
i Dependent Variable-boksen og såvel grp
visit
og girl
i Factor(s). Under fanen Fixed markeres både grp
og visit
, der vælges Factorial (defualt) og klikkes Add. Under fanen Random skal vi vælge Build nested terms i Random effects-boksen. Hernæst markeres først girl
og der trykkes pil ned. Derefter klikkes Within, man markerer grp
, trykker pil ned og klikker Add/Continue. Således at effekten fremstår som girl(grp). Under Statistics afkrydses Parameter estimates.
13.3.1 Prædikterede værdier
Når mixed modellen bygges benyttes også Save-knappen, hvor
man under Fixed Predicted Values sætter flueben ved
Predicted Values Derved fås en ny variabel i datasættet med navn FIXPRED_1
.
Brug herefter sædvanligt Lineplot for variablen FIXPRED_1
vs.
visit, med forskellige farver afhængigt af gruppe.
13.3.2 Modelkontrol i mixed model
Ved brug af Save-fanen afkrydes Predicted Values under Fixed Predicted Values . Under Predicted Values & Residuals sættes flueben ved Predicted Values og Residuals (som får navnene hhv PRED_1
og RESID_1
). Derved får man de små residualer (=individuelle/betingede).
Store residualer må vi selv danne: Observerede minus prædikterede værdier. Disse danner vi ved hjælp af menuen Transform/Compute, hvor vi i Target Variable skriver det nye variabelnavn, her fixed_resid
, og i feltet Numeric Expression sættes: bmd-Fixed Predicted Values
(Fixed Predicted Values er FXPRED_1
-variablen).
Herefter benyttes menuen Graph/Chart Builder/ til de forskellige figurer (residualplot og QQ-plot).
13.4 Specifikation af CS-struktur
Tager udgangspunkt i model fra kapitel 13.3.
Benyt menuen Analyze/Mixed Models/Linear, sæt girl
over i Subjects, visit
over i Repeated og skift Repeated Covariance Type til Compound Symmetry. Vælg Continue. Gå herefter ind i fanen Random og sørg for, at der ikke står noget her (fjern evt. girl fra denne). Nu sættes (eller bibeholdes) bmd
i Dependent Variable-boksen og såvel grp
, visit
og girl
i Factor(s). I fanen Fixed markeres både grp
og visit
, der vælges Factorial (defualt) og klikkes Add. Under Statistics afkrydses Parameter estimates og evt. Covariances of residuals.
13.5 Specifikation af kovariansstrukturer
Typen er her Unstructured. Fremgangsmåden/opsætningen er som ovenfor ved specifikation af CS-struktur, men i Repeated skal Repeated Covariance Type nu sættes til Unstructured.
Typen er her AR1.
Fremgangsmåden/opsætningen er som ovenfor ved specifikation af CS-struktur, men i Repeated skal Repeated Covariance Type nu sættes til AR(1). Man kan godt have en random effect af girl
samtidig med en AR(1)-korrelationsstruktur.
I visse versioner af SPSS kan man også vælge strukturen SP_POWER, som er den autoregressive struktur for irregulære tidspunkter.
13.6 Random regression, med ny tid
Benyt menuen Analyze/Mixed Models/Linear sæt girl
over i Subjects og klik Continue. Nu sættes bmd
i Dependent Variable-boksen og såvel grp
og girl
i Factor(s), men years
i Covariate(s). I fanen Fixed markeres både grp
og years
, der vælges Factorial (default) og klikkes Add. I fanen Random sættes girl
fra Subjects over i Combinations-boksen nederst i vinduet og years
sættes over i Model. Desuden sættes flueben ved Include intercept, og (meget
vigtigt): Covariance Type skiftes til Unstructured. Under Statistics afkrydses vi Parameter estimates.