library(pwr) pwr.t.test(d = {effectd}, sig.level = {alpha}, power = {power}, type = "two.sample")Gebruikershandleiding
Wanneer te gebruiken. Bepaal de minimale steekproefomvang voordat u met een onderzoek begint, rekening houdend met uw beoogde power en verwachte effectgrootte.
Aannames. Je hebt een geplande alfa (meestal 0.05), een streefvermogen (meestal 0.80) en een schatting van de effectgrootte op basis van pilotgegevens of literatuur.
Veelgemaakte fouten. Bereken de power niet achteraf om een niet-significant resultaat te rechtvaardigen (post-hoc power geeft geen informatie).
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Cohen J. Statistische poweranalyse voor de gedragswetenschappen. 2e editie. 1988.
Verwacht: 2 kolommen: Waarden van Groep 1 | Waarden van Groep 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden van Groep 1 | Waarden van Groep 2
g1 <- c({g1}) g2 <- c({g2}) t.test(g1, g2, var.equal = {eqvar})Gebruikershandleiding
Wanneer te gebruiken. Om de gemiddelden van twee onafhankelijke groepen te vergelijken op basis van een continue uitkomstvariabele.
Aannames: Onafhankelijkheid van de waarnemingen; bij benadering normale verdelingen in elke groep; gelijke varianties (of gebruik de Welch-correctie).
Veelgemaakte fouten. Niet gebruiken voor gepaarde gegevens: gebruik in plaats daarvan de gepaarde t-toets. Controleer eerst de normaliteit bij kleine steekproeven.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Student. De waarschijnlijke fout van een gemiddelde. Biometrika. 1908;6(1):1-25.
- Welch BL. De generalisatie van het probleem van Student. Biometrika. 1947;34:28-35.
Verwacht: 2 kolommen: Voor / Voorwaarde A | Na / Voorwaarde B
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Voor / Voorwaarde A | Na / Voorwaarde B
pre <- c({pre}) post <- c({post}) t.test(pre, post, paired = TRUE)Gebruikershandleiding
Wanneer te gebruiken. Vergelijk twee metingen die bij dezelfde proefpersonen zijn uitgevoerd (voor/na, gepaarde metingen).
Aannames. De verschillen tussen paren zijn bij benadering normaal verdeeld.
Veelgemaakte fouten. Alleen geldig wanneer de paren daadwerkelijk overeenkomen; niet gebruiken voor twee onafhankelijke steekproeven.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Student. De waarschijnlijke fout van een gemiddelde. Biometrika. 1908;6(1):1-25.
{p_r} # z: 2 * pnorm(-abs(z)) t: 2 * pt(-abs(t), df) # chi-kwadraat: pchisq(x, df, lower.tail = FALSE) F: pf(x, df1, df2, lower.tail = FALSE)Gebruikershandleiding
Wanneer te gebruiken. Zet een toetsingsstatistiek (z, t, χ² of F) om in een p-waarde onder de nulhypothese.
Aannames. Je hebt al een correct berekende statistiek en het juiste aantal vrijheidsgraden.
Veelgemaakte fouten. Een p-waarde is niet de waarschijnlijkheid dat de nulhypothese waar is, noch de grootte van een effect.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Wasserstein RL, Lazar NA. De ASA-verklaring over p-waarden. Am Stat. 2016;70:129-133.
Verwachte opmaak: Waargenomen aantallen (rijen op nieuwe regels, gescheiden door komma's)
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Waargenomen aantallen (rijen op nieuwe regels, gescheiden door komma's)
tbl <- matrix(c({obs_flat}), nrow = {nrows}, byrow = TRUE) chisq.test(tbl)Gebruikershandleiding
Wanneer te gebruiken. Test of twee categorische variabelen met elkaar samenhangen, met behulp van een contingentietabel van tellingen.
Aannames. Het verwachte aantal in elke cel moet minstens 5 zijn. De waarnemingen zijn onafhankelijk.
Veelgemaakte fouten. Als de verwachte aantallen lager zijn dan 5 (kleine steekproeven), gebruik dan de exacte toets van Fisher.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Pearson K. Over het criterium dat een gegeven systeem van afwijkingen… Philos Mag. 1900;50:157-175.
Verwacht formaat: Voorbeeldwaarden
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht formaat: Voorbeeldwaarden
x <- c({data}) t.test(x, conf.level = {conf})$conf.intGebruikershandleiding
Wanneer te gebruiken. Om de plausibele bandbreedte voor een populatiegemiddelde te schatten op basis van een steekproef.
Aannames. De steekproef is willekeurig; voor kleine n is de onderliggende verdeling bij benadering normaal.
Veelgemaakte fouten. Een 95% betrouwbaarheidsinterval betekent niet dat er een kans van 95% is dat de werkelijke waarde binnen dit ene interval ligt.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Neyman J. Schets van een theorie van statistische schatting. Philos Trans R Soc A. 1937;236:333-380.
Verwacht: 2 kolommen: Waarden van Groep 1 | Waarden van Groep 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden van Groep 1 | Waarden van Groep 2
g1 <- c({g1}) g2 <- c({g2}) library(effsize) cohen.d(g1, g2) cohen.d(g1, g2, hedges.correction = TRUE) # Hedges' gGebruikershandleiding
Wanneer te gebruiken. Om de gestandaardiseerde omvang van het verschil tussen twee groepen te kwantificeren.
Aannames. Continue data; de gepoolde standaarddeviatie is relevant (vergelijkbare varianties).
Veelgemaakte fouten. Effectgrootte meet de omvang, niet de statistische significantie of klinische relevantie.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Cohen J. Statistische poweranalyse voor de gedragswetenschappen. 1988.
Verwacht: één rij per record, Label: v1, v2 …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, Label: v1, v2 …
df <- stack(list({anova_groups})) summary(aov(values ~ ind, data = df))Gebruikershandleiding
Wanneer te gebruiken. Om de gemiddelden van drie of meer onafhankelijke groepen tegelijk te vergelijken.
Aannames: Onafhankelijkheid; normaliteit binnen groepen; homogeniteit van variantie.
Veelgemaakte fouten. Voer geen meerdere paarsgewijze t-toetsen uit: dit verhoogt de kans op een type I-fout.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Fisher RA. Statistische methoden voor onderzoekers. 1925.
Verwacht: 2 kolommen: Variabele X | Variabele Y
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Variabele X | Variabele Y
x <- c({x}) y <- c({y}) cor.test(x, y, method = "{method}"{spearman_exact})Gebruikershandleiding
Wanneer te gebruiken. Om de sterkte en richting van een lineair (Pearson) of monotoon (Spearman) verband tussen twee variabelen te meten.
Aannames. Pearson: lineariteit, normaliteit, geen grote uitschieters. Spearman: alleen een monotone relatie.
Veelgemaakte fouten. Correlatie impliceert geen causaliteit. Controleer de spreidingsgrafiek op niet-lineariteit.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Pearson K. Notitie over regressie en erfelijkheid. Proc R Soc Lond. 1895;58:240-242.
- Spearman C. Het bewijs en de meting van associatie. Am J Psychol. 1904;15:72-101.
tp <- {tp}; fp <- {fp}; fn <- {fn}; tn <- {tn} sens <- tp/(tp+fn); spec <- tn/(tn+fp) c(sensitivity = sens, specificity = spec)Gebruikershandleiding
Wanneer te gebruiken? Evalueer de prestaties van de diagnostische test aan de hand van een vergelijking van het testresultaat met de werkelijke ziektestatus (2x2).
Aannames. Een geldige referentiestandaard (gouden standaard) classificeert de werkelijke ziektestatus.
Veelgemaakte fouten. De positief voorspellende waarde (PPV) en negatief voorspellende waarde (NPV) zijn sterk afhankelijk van de prevalentie van de ziekte: rapporteer deze.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Altman DG, Bland JM. Diagnostische tests 1: sensitiviteit en specificiteit. BMJ. 1994;308:1552.
Verwachte opmaak: Score, Label: één per regel (1=positief, 0=negatief)
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Score, Label: één per regel (1=positief, 0=negatief)
library(pROC) scores <- c({roc_scores}) labels <- c({roc_labels}) roc(labels, scores, direction = "{roc_dir}") # {roc_dir_note}Gebruikershandleiding
Wanneer te gebruiken. Evalueer het onderscheidend vermogen van een continue test of score.
Aannames. Binaire werkelijke uitkomst (0/1); continue voorspellende variabele of score.
Veelgemaakte fouten. Een AUC-waarde dicht bij 0.5 betekent geen onderscheidingsvermogen; interpreteer een enkele drempelwaarde niet te sterk.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Hanley JA, McNeil BJ. De betekenis en het gebruik van het gebied onder een ROC-curve. Radiologie. 1982;143:29-36.
tbl <- matrix(c({a},{b},{c},{d}), nrow = 2, byrow = TRUE) fisher.test(tbl) # of epitools::oddsratio(tbl)Gebruikershandleiding
Wanneer te gebruiken. Kwantificeer de associatie in een 2×2-tabel; OR voor case-controlstudies, RR voor cohortstudies.
Aannames: De tellingen zijn onafhankelijk; de cellen zijn niet nul (anders moet een continuïteitscorrectie worden toegepast).
Veelgemaakte fouten. Interpreteer OR niet als RR wanneer de uitkomst vaak voorkomt (>10%).
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Woolf B. Over het schatten van de relatie tussen bloedgroep en ziekte. Ann Hum Genet. 1955;19:251-253.
cer <- {cer}; eer <- {eer} arr <- cer - eer # > 0 voordeel, < 0 nadeel als (arr > 0) c(ARR = arr, NNT_voordeel = plafond(afronden(1/arr, 6))) anders c(ARR = arr, NNH_nadeel = plafond(afronden(1/(eer - cer), 6)))Gebruikershandleiding
Wanneer te gebruiken. Druk het behandelingsvoordeel uit als het aantal patiënten dat behandeld moet worden om één incident te voorkomen.
Aannames. De incidentiecijfers zijn afkomstig van vergelijkbare groepen over een vastgestelde tijdsperiode.
Veelgemaakte fouten. Vermeld altijd de tijdshorizon; NNT is zonder die informatie betekenisloos.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Laupacis A, Sackett DL, Roberts RS. Een beoordeling van klinisch bruikbare meetinstrumenten. N Engl J Med. 1988;318:1728-1733.
Verwacht: één rij per record, TP, FP, FN, TN
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, TP, FP, FN, TN
bibliotheek(mada) dat <- read.csv(text = "TP,FP,FN,TN {dm_csv}") madad(dat); reitsma(dat)Gebruikershandleiding
Wanneer te gebruiken. Combineer de diagnostische nauwkeurigheid (sensitiviteit en specificiteit) van verschillende studies met behulp van een random-effects model.
Aannames. Elk onderzoek rapporteert een 2×2-vergelijking van testresultaten versus referentiestandaard.
Veelgemaakte fouten. Het afzonderlijk samenvoegen van sensitiviteit en specificiteit negeert hun correlatie; een bivariate/HSROC-model heeft de voorkeur wanneer dit mogelijk is.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Reitsma JB, et al. Bivariate analyse van sensitiviteit en specificiteit. J Clin Epidemiol. 2005;58:982-990.
Verwachte opmaak: Onderwerpen: één per regel, aantal beoordelaars per categorie
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Onderwerpen: één per regel, aantal beoordelaars per categorie
library(irr) # rijen = onderwerpen, kolommen = categorieën (aantal beoordelaars) kappam.fleiss(ratings)
Gebruikershandleiding
Wanneer te gebruiken. Om de overeenstemming te meten tussen drie of meer beoordelaars die onderwerpen in categorieën indelen (Fleiss' kappa).
Aannames. Elk onderwerp wordt beoordeeld door hetzelfde aantal beoordelaars; categorieën sluiten elkaar uit.
Veelgemaakte fouten. Kappa is prevalentiegevoelig; rapporteer de categoriefrequenties erbij.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Fleiss JL. Het meten van overeenstemming op nominale schalen tussen meerdere beoordelaars. Psychol Bull. 1971;76:378-382.
Verwacht: 2 kolommen: Metingen met methode 1 | Metingen met methode 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Metingen met methode 1 | Metingen met methode 2
m1 <- c({m1}) m2 <- c({m2}) d <- m1 - m2 c(bias = mean(d), lower = mean(d)-1.96*sd(d), upper = mean(d)+1.96*sd(d)) cor.test(m1, m2) coef(lm(m2 ~ m1)) # gewone kleinste-kwadratenlijnGebruikershandleiding
Wanneer te gebruiken? Vergelijk twee meetmethoden met behulp van bias, 95%-betrouwbaarheidsgrenzen, correlatie en de gewone kleinste-kwadratenmethode.
Aannames. Gepaarde metingen van dezelfde grootheid met beide methoden.
Veelgemaakte fouten. Beoordeel overeenstemming niet alleen op basis van correlatie: geef de grenzen van de overeenstemming aan.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Bland JM, Altman DG. Statistische methoden voor het beoordelen van overeenstemming tussen twee methoden. Lancet. 1986;1:307-310.
Verwacht: 2 kolommen: X-waarden | Y-waarden
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: X-waarden | Y-waarden
x <- c({x}) y <- c({y}) summary(lm(y ~ x))Gebruikershandleiding
Wanneer te gebruiken. Om een continue uitkomst te modelleren en te voorspellen op basis van één continue voorspellende variabele.
Aannames: lineariteit, onafhankelijkheid, constante variantie (homoscedasticiteit), normale residuen.
Veelgemaakte fouten. Extrapoleer voorspellingen niet buiten het waargenomen bereik van x.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Galton F. Regressie naar middelmatigheid in erfelijke lengte. J Anthropol Inst. 1886;15:246-263.
Verwacht: één rij per record, y, x1, x2, …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, y, x1, x2, …
df <- read.csv(text = "y,x1,x2 {mlr_csv}") summary(lm(y ~ ., data = df))Gebruikershandleiding
Wanneer te gebruiken. Modelleer een continue uitkomstvariabele op basis van meerdere voorspellende variabelen tegelijk; lees elke coëfficiënt als het aangepaste effect.
Aannames: Lineariteit, onafhankelijke fouten, constante variantie, normale residuen, geen ernstige multicollineariteit.
Veelgemaakte fouten. Vermijd het toevoegen van meer voorspellende variabelen dan de steekproef aankan; controleer de variantie-inflatiecoëfficiënten.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Draper NR, Smith H. Toegepaste regressieanalyse. 3e editie. 1998.
Verwacht: één rij per record, y(0/1), x1, x2, …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, y(0/1), x1, x2, …
df <- read.csv(text = "y,x1 {logit_csv}") m <- glm(y ~ ., data = df, family = binomial) summary(m); exp(cbind(OR = coef(m), confint(m)))Gebruikershandleiding
Wanneer te gebruiken. Modelleer een binaire uitkomst en bereken de aangepaste odds ratio's voor elke voorspellende variabele.
Aannames. De uitkomst is binair; de log-odds zijn lineair in de voorspellende variabelen; de waarnemingen zijn onafhankelijk.
Veelgemaakte fouten. Vuistregel: minimaal ~10 gebeurtenissen per voorspeller om instabiele schattingen te voorkomen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Hosmer DW, Lemeshow S, Sturdivant RX. Toegepaste logistische regressie. 3e editie. 2013.
Verwacht: één rij per record, aantal, x1, x2, …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, aantal, x1, x2, …
df <- read.csv(text = "count,x1 {pois_csv}") m <- glm(count ~ ., data = df, family = poisson) summary(m); exp(coef(m))Gebruikershandleiding
Wanneer te gebruiken. Modelleer het aantal uitkomsten en bereken de relatieve risico's voor elke voorspellende variabele.
Aannames. De tellingen volgen een Poisson-verdeling; het gemiddelde is gelijk aan de variantie (controleer op overdispersie).
Veelgemaakte fouten. Als de variantie veel groter is dan het gemiddelde, gebruik dan een negatief binomiaal model of een quasi-Poisson model.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Cameron AC, Trivedi PK. Regressieanalyse van telgegevens. 2e editie. 2013.
Verwacht: 2 kolommen: Waarden van Groep 1 | Waarden van Groep 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden van Groep 1 | Waarden van Groep 2
g1 <- c({g1}) g2 <- c({g2}) wilcox.test(g1, g2)Gebruikershandleiding
Wanneer te gebruiken? Niet-parametrische vergelijking van twee onafhankelijke groepen wanneer normaliteit twijfelachtig is.
Aannames. Onafhankelijke waarnemingen; ordinale of continue gegevens.
Veelgemaakte fouten. Het vergelijkt rangverdelingen, niet gemiddelden: beschrijf de resultaten als een verschuiving in de verdeling.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Mann HB, Whitney DR. Over een toets om te bepalen of een van twee willekeurige variabelen stochastisch groter is. Ann Math Stat. 1947;18:50-60.
Verwacht: 2 kolommen: Waarden vóór | Waarden na
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden vóór | Waarden na
pre <- c({pre}) post <- c({post}) wilcox.test(pre, post, paired = TRUE)Gebruikershandleiding
Wanneer te gebruiken. Niet-parametrische gepaarde toets voor gepaarde gegevens wanneer de verschillen niet normaal verdeeld zijn.
Aannames. Verschillen zijn symmetrisch rond de mediaan; gepaard ontwerp.
Veelgemaakte fouten. Nulverschillen worden weggelaten; gelijke scores verminderen de statistiek.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Wilcoxon F. Individuele vergelijkingen door middel van rangschikkingsmethoden. Biometrics Bull. 1945;1:80-83.
Verwacht: één rij per record, Label: v1, v2 …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, Label: v1, v2 …
df <- stack(list({anova_groups})) kruskal.test(values ~ ind, data = df)Gebruikershandleiding
Wanneer te gebruiken? Niet-parametrische vergelijking tussen drie of meer groepen.
Aannames. Onafhankelijke groepen; ordinale of continue gegevens; gelijkvormige verdelingen.
Veelgemaakte fouten. Een significante H-waarde geeft aan welke groepen verschillen, niet welke: voer daarom de Dunn-test uit.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Kruskal WH, Wallis WA. Gebruik van rangen in een variantieanalyse met één criterium. J Am Stat Assoc. 1952;47:583-621.
tbl <- matrix(c({a},{b},{c},{d}), nrow = 2, byrow = TRUE) fisher.test(tbl)Gebruikershandleiding
Wanneer te gebruiken. Nauwkeurige associatietoets in een 2×2-tabel, ideaal voor kleine steekproeven.
Aannames. Vaste marginale totalen; onafhankelijke waarnemingen.
Veelgemaakte fouten. Voor grote steekproeven is de chi-kwadraat toets rekenkundig eenvoudiger en equivalent.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Fisher RA. De logica van inductieve gevolgtrekking. JR Stat Soc. 1935;98:39-82.
Verwacht: één rij per record, tijd, gebeurtenis (1=gebeurtenis, 0=gecensureerd)[, groep]
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, tijd, gebeurtenis (1=gebeurtenis, 0=gecensureerd)[, groep]
library(survival) df <- read.csv(text = "time,event,group {km_csv}") fit <- survfit(Surv(time, event) ~ group, data = df) summary(fit); survdiff(Surv(time, event) ~ group, data = df)Gebruikershandleiding
Wanneer te gebruiken. Schat de overlevingskans in de loop van de tijd op basis van gecensureerde gegevens over de tijd tot de gebeurtenis, en vergelijk groepen met de log-ranktoets.
Aannames. Censurering is niet-informatief; overlevingstijden worden gemeten vanaf een gemeenschappelijk beginpunt.
Veelgemaakte fouten. Beschouw gecensureerde onderwerpen aan het eind niet als gebeurtenisvrij: de methode houdt er al rekening mee.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Kaplan EL, Meier P. Niet-parametrische schatting op basis van onvolledige waarnemingen. J Am Stat Assoc. 1958;53:457-481.
Verwacht: één rij per record, tijd, gebeurtenis (1/0), x1, x2, …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, tijd, gebeurtenis (1/0), x1, x2, …
library(survival) df <- read.csv(text = "time,event,x1 {cox_csv}") summary(coxph(Surv(time, event) ~ ., data = df))Gebruikershandleiding
Wanneer te gebruiken. Schat de hazard ratio's voor covariaten op de tijd tot een gebeurtenis, zonder de basishazard te specificeren.
Aannames. Proportionele hazards (covariaateffecten constant in de tijd); niet-informatieve censurering.
Veelgemaakte fouten. Test de proportional-hazards-aanname (bijvoorbeeld Schoenfeld-residuen) voordat u de HR's vertrouwt.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Cox DR. Regressiemodellen en sterftetabellen. JR Stat Soc B. 1972;34:187-220.
Verwachte opmaak: Studies: Label, Effect, LowerCI, UpperCI (één per regel)
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Studies: Label, Effect, LowerCI, UpperCI (één per regel)
library(meta) studlab <- c({meta_studlab}) TE <- c({meta_te}) lower <- c({meta_lower}) upper <- c({meta_upper}) seTE <- (upper - lower) / (2 * 1.96) # SE van het 95% betrouwbaarheidsinterval, zoals de calculator doet m <- metagen(TE, seTE, studlab = studlab, sm = "SMD", common = {meta_common}, random = {meta_random}, method.tau = "DL") summary(m) # gepoolde effect, Cochran's Q, I2, tau2Gebruikershandleiding
Wanneer te gebruiken? Combineer de effectgroottes over verschillende studies en kwantificeer de heterogeniteit (I-kwadraat).
Aannames. De onderzoeken zijn klinisch en methodologisch voldoende vergelijkbaar om te kunnen worden gecombineerd.
Veelgemaakte fouten. Een hoge I-kwadraatwaarde (>75%) duidt op aanzienlijke heterogeniteit: interpreteer de gepoolde schatting met de nodige voorzichtigheid.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- DerSimonian R, Laird N. Meta-analyse in klinische onderzoeken. Control Clin Trials. 1986;7:177-188.
Verwacht: één rij per record, label, effect, standaardfout
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, label, effect, standaardfout
library(meta) studlab <- c({funnel_studlab}) TE <- c({funnel_te}) seTE <- c({funnel_se}) m <- metagen(TE, seTE, studlab = studlab, common = TRUE, random = FALSE) summary(m) # gepoolde schatting met vaste effecten en 95% betrouwbaarheidsinterval, zoals de calculator rapporteert funnel(m) metabias(m, method.bias = "linreg", k.min = 3) # Egger's regressietest: intercept, t, pGebruikershandleiding
Wanneer te gebruiken. Detecteer asymmetrie in trechterdiagrammen (een indicator van mogelijke publicatiebias) via de regressietest van Egger.
Aannames. Minimaal ~10 studies voor een goed onderbouwde Egger-test; effectgroottes en standaardfouten per studie.
Veelgemaakte fouten. Asymmetrie kan ook wijzen op echte heterogeniteit of effecten van kleine studies, en niet alleen op vertekening.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Egger M, Davey Smith G, Schneider M, Minder C. Vertekening in meta-analyse gedetecteerd door een eenvoudige, grafische test. BMJ. 1997;315:629-634.
library(PRISMA2020) # lever tellingen aan PRISMA_flowdiagram()
Gebruikershandleiding
Wanneer te gebruiken. Stel de aantallen records samen voor een PRISMA 2020-stroomschema voor studiekeuze en controleer de rekenkundige balans.
Aannames. De aantallen zijn afkomstig uit uw screeningslogboek bij elke fase.
Veelgemaakte fouten. De cijfers moeten kloppen: gescreende teksten min uitgesloten teksten moet gelijk zijn aan volledig beoordeelde teksten.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Page MJ, et al. De PRISMA 2020-verklaring. BMJ. 2021;372:n71.
# GRADE is een gestructureerd oordeel, geen berekening. # Start: RCT = Hoog (4), Observationeel = Laag (2). # Verlaging (-1/-2) bij risico op bias, inconsistentie, # indirectheid, onnauwkeurigheid, publicatiebias. # Verhoging (+1/+2) bij groot effect, dosis-responsrelatie.
Gebruikershandleiding
Wanneer te gebruiken. Beoordeel de zekerheid van een reeks bewijzen (Hoog / Matig / Laag / Zeer laag) met behulp van het GRADE-raamwerk.
Aannames. U heeft elk GRADE-domein beoordeeld op de gewenste uitkomst.
Veelgemaakte fouten. Beoordeel de zekerheid per uitkomst, niet per onderzoek; upgrades gelden normaal gesproken alleen voor observationeel bewijs.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Guyatt GH, et al. GRADE: een opkomende consensus over het beoordelen van de kwaliteit van bewijsmateriaal. BMJ. 2008;336:924-926.
Verwachte opmaak: Overeenkomstmatrix: rijen op nieuwe regels, gescheiden door komma's.
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Overeenkomstmatrix: rijen op nieuwe regels, gescheiden door komma's.
library(irr) tbl <- matrix(c({mat_flat}), nrow = {nrows}, byrow = TRUE) r1 <- rep(row(tbl), tbl); r2 <- rep(col(tbl), tbl) # één paar beoordelingen per item kappa2(cbind(r1, r2))Gebruikershandleiding
Wanneer te gebruiken. Om de categorische overeenstemming tussen beoordelaars te meten, voorbij wat op basis van toeval verwacht kan worden.
Aannames. Twee beoordelaars; dezelfde categorieën; onafhankelijke beoordelingen.
Veelgemaakte fouten. Kappa is gevoelig voor prevalentie: een hoge mate van overeenstemming kan een lage kappa opleveren.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Cohen J. Een overeenkomstcoëfficiënt voor nominale schalen. Educ Psychol Meas. 1960;20:37-46.
Verwachte opmaak: Beoordelaars (één per regel): Naam beoordelaar: v1, v2 …
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Beoordelaars (één per regel): Naam beoordelaar: v1, v2 …
library(irr) ratings <- {icc_ratings} # rijen = onderwerpen, kolommen = beoordelaars icc(ratings, model = "oneway", unit = "single")Gebruikershandleiding
Wanneer te gebruiken. Kwantificeer de betrouwbaarheid van continue metingen door meerdere beoordelaars.
Aannames. Continue beoordelingen; consistente beoordelaars voor alle onderwerpen.
Veelgemaakte fouten. Specificeer de ICC-vorm (1,1 / 2,1 / 3,1): deze zijn niet uitwisselbaar.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Shrout PE, Fleiss JL. Intraclasscorrelaties. Psychol Bull. 1979;86:420-428.
Verwachte opmaak: Respondenten × items: één per regel, items gescheiden door komma's.
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwachte opmaak: Respondenten × items: één per regel, items gescheiden door komma's.
library(psych) # rijen = respondenten, kolommen = items alpha(data_matrix)
Gebruikershandleiding
Wanneer te gebruiken. Beoordeel de interne consistentie (betrouwbaarheid) van een schaal met meerdere items.
Aannames. De items meten één onderliggend construct; scoren op intervalniveau.
Veelgemaakte fouten. De alfawaarde neemt toe met het aantal items: interpreteer hoge waarden op lange schalen niet te sterk.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Cronbach LJ. Coëfficiënt alfa en de interne structuur van tests. Psychometrika. 1951;16:297-334.
te <- {te}; vooroordeel <- {vooroordeel}; cv <- {cv} sigma <- (te - abs(bias))/cv sigmaGebruikershandleiding
Wanneer te gebruiken. Bereken de sigma-waarde voor een laboratoriumtest en selecteer de juiste kwaliteitscontrole-regels.
Aannames. TEa, bias en CV zijn uitgedrukt in dezelfde procentuele eenheden.
Veelgemaakte fouten. Een Sigma lager dan 3 vereist intensieve kwaliteitscontrole met meerdere regels; werk niet met afzonderlijke regels.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Westgard JO. Basisprincipes van kwaliteitscontrole. 4e editie. 2016.
Verwacht: 2 kolommen: Waarden van methode 1 | Waarden van methode 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden van methode 1 | Waarden van methode 2
m1 <- c({m1}) m2 <- c({m2}) d <- m1 - m2 c(bias = gemiddelde(d), onder = gemiddelde(d)-1.96*sd(d), boven = gemiddelde(d)+1.96*sd(d))Gebruikershandleiding
Wanneer te gebruiken? Beoordeel de overeenstemming tussen twee meetmethoden aan de hand van de gemiddelde bias en de betrouwbaarheidsgrenzen.
Aannames. De verschillen zijn bij benadering normaal verdeeld en constant over het meetbereik.
Veelgemaakte fouten. Gebruik nooit de correlatiecoëfficiënt om de overeenstemming tussen methoden te beoordelen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Bland JM, Altman DG. Statistische methoden voor het beoordelen van overeenstemming. Lancet. 1986;1:307-310.
Verwacht: 2 kolommen: Waarden van methode 1 | Waarden van methode 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden van methode 1 | Waarden van methode 2
library(mcr) m1 <- c({m1}) m2 <- c({m2}) mcreg(m1, m2, method.reg = "PaBa")Gebruikershandleiding
Wanneer te gebruiken. Niet-parametrische methodevergelijking die robuust is tegen uitschieters en meetfouten in beide methoden.
Aannames. Continue gepaarde metingen; geen aanname over de foutverdeling.
Veelgemaakte fouten. Een redelijk meetbereik is nodig om de helling betrouwbaar te kunnen inschatten.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Passing H, Bablok W. Een nieuwe biometrische procedure voor het testen van de gelijkheid van metingen. J Clin Chem Clin Biochem. 1983;21:709-720.
events <- {events}; pt <- {pt} poisson.test(events, pt) # exact betrouwbaarheidsinterval per eenheid persoon-tijd poisson.test(events, pt)$conf.int * {mult_n} # 95% betrouwbaarheidsinterval per {mult_n}Gebruikershandleiding
Wanneer te gebruiken. Bereken het aantal nieuwe incidenten per eenheid persoonsjaren met risico.
Aannames: De persoons-tijd wordt correct geaccumuleerd; de gebeurtenissen zijn onafhankelijk.
Veelgemaakte fouten. Verwar incidentie (nieuwe gevallen in de loop van de tijd) niet met prevalentie (bestaande gevallen).
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Rothman KJ, Greenland S, Lash TL. Moderne epidemiologie. 3e editie. 2008.
cases <- {cases}; pop <- {pop} prop.test(cases, pop)Gebruikershandleiding
Wanneer te gebruiken. Om het percentage van een populatie met een bepaalde aandoening op een bepaald moment of in een bepaalde periode te schatten.
Aannames. Het aantal gevallen en de populatie worden over dezelfde periode gemeten.
Veelgemaakte fouten. Prevalentie weerspiegelt zowel de incidentie als de duur, niet het aantal nieuwe gevallen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Rothman KJ, Greenland S, Lash TL. Moderne epidemiologie. 3e editie. 2008.
a <- {a}; b <- {b}; c <- {c}; d <- {d} Re <- a/(a+b); Ru <- c/(c+d) c(AR = Re-Ru, ARP = (Re-Ru)/Re, RR = Re/Ru)Gebruikershandleiding
Wanneer te gebruiken. Kwantificeer het extra risico dat aan een blootstelling kan worden toegeschreven (toewijsbaar risico, toewijsbare fractie, populatie-AR).
Aannames. Een causale interpretatie vereist dat het verband tussen blootstelling en uitkomst niet verstoord wordt door confounders.
Veelgemaakte fouten. Toewijsbare fracties veronderstellen causaliteit: ze zijn niet van toepassing op puur associatieve data.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Rothman KJ, Greenland S, Lash TL. Moderne epidemiologie. 3e editie. 2008.
ea <- {ea}; na <- {na}; ec <- {ec}; nc <- {nc} fp <- function(x1, x0) fisher.test(matrix(c(x1, na - x1, x0, nc - x0), 2, byrow = TRUE))$p.value p0 <- fp(ea, ec); x1 <- ea; x0 <- ec; fi <- 0 als (p0 < 0.05) terwijl (fp(x1, x0) < 0.05) { als (ea <= ec) x1 <- x1 + 1 anders x0 <- x0 + 1; fi <- fi + 1 } c(fisher_p = p0, fragility_index = if (p0 < 0.05) fi anders NA)Gebruikershandleiding
Wanneer te gebruiken? Beoordeel hoe robuust een statistisch significant resultaat van een gerandomiseerde gecontroleerde studie is.
Aannames. Binaire uitkomst 2×2 uit een gerandomiseerd onderzoek.
Veelgemaakte fouten. Een lage fragiliteitsindex betekent op zich niet dat een onderzoek ongeldig is.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Walsh M, et al. De statistische significantie van de resultaten van gerandomiseerde, gecontroleerde onderzoeken is vaak fragiel. J Clin Epidemiol. 2014;67:622-628.
Verwacht formaat: Waarden
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht formaat: Waarden
x <- c({data}) shapiro.test(x)Gebruikershandleiding
Wanneer te gebruiken. Beoordeel of de gegevens aannemelijk afkomstig zijn van een normale verdeling voordat u een parametrische toets kiest.
Aannames. Continue data; n tussen 3 en 5000 voor Shapiro-Wilk.
Veelgemaakte fouten. Een p-waarde groter dan 0.05 bewijst geen normaliteit: bekijk altijd ook het histogram en de QQ-plot.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Shapiro SS, Wilk MB. Een variantieanalyse-test voor normaliteit. Biometrika. 1965;52:591-611.
Verwacht formaat: Waarden
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht formaat: Waarden
x <- c({data}) summary(x); sd(x)Gebruikershandleiding
Wanneer te gebruiken. Om een continue variabele samen te vatten voor een tabel (gemiddelde, standaarddeviatie, mediaan, interkwartielbereik, bereik).
Aannames. Geen: uitsluitend beschrijvend.
Veelgemaakte fouten. Geef het gemiddelde ± standaarddeviatie alleen weer bij nagenoeg symmetrische gegevens; gebruik de mediaan (interkwartielbereik) bij scheve gegevens.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Amerikaanse Psychologische Vereniging. Publicatiehandleiding. 7e editie. 2020.
x <- {val}; mu <- {mu}; sigma <- {sig} z <- (x - mu)/sigma c(z = z, percentiel = pnorm(z))Gebruikershandleiding
Wanneer te gebruiken. Standaardiseer een waarde en bepaal het bijbehorende percentiel in een normale verdeling.
Aannames: Het populatiegemiddelde en de standaarddeviatie zijn bekend; de verdeling is normaal.
Veelgemaakte fouten. Gebruik in plaats daarvan een t-score bij het schatten van parameters op basis van een kleine steekproef.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Moore DS, McCabe GP, Craig BA. Inleiding tot de praktijk van de statistiek. 9e editie. New York: WH Freeman; 2017.
Verwacht formaat: Waarden
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht formaat: Waarden
x <- c({data}) hist(x); summary(x)Gebruikershandleiding
Wanneer te gebruiken? Om de vorm, het centrum en de spreiding van een continue variabele te visualiseren.
Aannames. Continue data; voldoende steekproefomvang voor zinvolle indelingen.
Veelgemaakte fouten. De breedte van de bak heeft een grote invloed op het uiterlijk: probeer er een paar uit voordat u een definitieve keuze maakt.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Sturges HA. De keuze van een klasse-interval. J Am Stat Assoc. 1926;21:65-66.
Verwacht: 2 kolommen: Waarden van methode A | Waarden van methode B
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Waarden van methode A | Waarden van methode B
x <- c({x}) y <- c({y}) n <- length(x) sx <- mean((x-mean(x))^2); sy <- mean((y-mean(y))^2); sxy <- mean((x-mean(x))*(y-mean(y))) ccc <- 2*sxy/(sx+sy+(mean(x)-mean(y))^2) ccc cor(x, y) # precisie; Cb = ccc / r (nauwkeurigheid)Gebruikershandleiding
Wanneer te gebruiken. Om de overeenstemming (niet alleen de correlatie) tussen twee methoden die dezelfde continue grootheid meten te kwantificeren.
Aannames: Gepaarde metingen bij dezelfde proefpersonen; een bij benadering bivariate normale verdeling voor het Fisher-z betrouwbaarheidsinterval; een redelijk breed bereik aan waarden.
Veelgemaakte fouten. Gebruik de Pearson-correlatiecoëfficiënt (r) niet als enige maatstaf voor overeenstemming. Bij zeer kleine steekproeven is het betrouwbaarheidsinterval onbetrouwbaar; voeg een Bland-Altman-grafiek toe.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Lin LI. Een concordantiecorrelatiecoëfficiënt om de reproduceerbaarheid te evalueren. Biometrics. 1989;45(1):255-268.
- McBride GB. Een voorstel voor criteria voor de mate van overeenstemming voor Lin's concordantiecorrelatiecoëfficiënt. NIWA Client Report HAM2005-062. 2005.
Verwacht: 2 kolommen: Referentie / methode X | Testmethode Y
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Referentie / methode X | Testmethode Y
x <- c({x}) y <- c({y}) lambda <- {delta} sxx <- var(x); syy <- var(y); sxy <- cov(x, y) b <- (syy - lambda*sxx + sqrt((syy - lambda*sxx)^2 + 4*lambda*sxy^2)) / (2*sxy) a <- gemiddelde(y) - b*mean(x) c(snijpunt = a, helling = b)Gebruikershandleiding
Wanneer te gebruiken. Om twee meetmethoden te vergelijken wanneer beide een meetfout bevatten (in tegenstelling tot de gewone kleinste kwadratenmethode, die ervan uitgaat dat X exact is).
Aannames. De verhouding van de foutvarianties (λ) van de twee methoden is bekend of kan worden geschat op basis van herhaalde metingen; de fouten zijn ruwweg constant over het bereik.
Veelgemaakte fouten. Als de fout evenredig is met de concentratie, gebruik dan gewogen Deming of Passing-Bablok. Gebruik geen OLS als beide methoden ruis bevatten.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Deming WE. Statistische aanpassing van gegevens. Wiley; 1943.
- Linnet K. Evaluatie van regressieprocedures voor methodenvergelijkingsstudies. Clin Chem. 1993;39(3):424-432.
Verwacht: één rij per record, naam, waarde, type
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: één rij per record, naam, waarde, type
u <- c(0.12, 0.20, 0.05/sqrt(3), 0.10/sqrt(6)) # standaardonzekerheden uc <- sqrt(som(u^2)) k <- {k} c(gecombineerd = uc, uitgebreid = k*uc)Gebruikershandleiding
Wanneer te gebruiken. Stel een GUM-onzekerheidsbudget op: converteer elk onderdeel naar een standaardonzekerheid en combineer ze door middel van de wortel van de som der kwadraten.
Aannames. De componenten zijn onafhankelijk (geen covariantietermen) en de meetgrootheid is een som van de bijdragen met gevoeligheidscoëfficiënten van 1.
Veelgemaakte fouten. Meng geen uitgebreide en standaard onzekerheden. Voeg gecorreleerde of niet-eenheidsgevoelige termen handmatig toe voordat u ze invoert.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- JCGM 100:2008. Evaluatie van meetgegevens — Richtlijn voor de uitdrukking van meetonzekerheid (GUM).
- Eurachem/CITAC-handleiding CG 4. Kwantificering van onzekerheid in analytische metingen. 3e editie. 2012.
Verwacht: 4 kolommen: Blanco herhalingen | Herhalingen van monsters met lage concentratie | Kalibratieconcentraties | Instrumentrespons
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 4 kolommen: Blanco herhalingen | Herhalingen van monsters met lage concentratie | Kalibratieconcentraties | Instrumentrespons
blank <- c({blank}) low <- c({low}) LoB <- mean(blank) + 1.645 * sd(blank) cp <- 1.645 / (1 - 1 / (4 * (length(low) - 1))) LoD <- LoB + cp * sd(low) c(LoB = LoB, LoD = LoD) # ICH Q2 kalibratiecurve alternatief: conc <- c({x}); resp <- c({y}) fit <- lm(resp ~ conc) c(LoD = 3.3 * summary(fit)$sigma / coef(fit)[2], LoQ = 10 * summary(fit)$sigma / coef(fit)[2])Gebruikershandleiding
Wanneer te gebruiken. Schat de blanco-limiet en de detectielimiet op basis van blanco- en laaggedoseerde herhalingen met behulp van de parametrische formules (CLSI EP17-A2 beschrijft het onderzoeksontwerp), of LoD/LoQ op basis van een kalibratiecurve (ICH Q2(R1)).
Aannames. EP17: de resultaten zijn bij benadering Gaussisch, gebaseerd op ten minste 60 blanco metingen en 60 metingen op laag niveau (20 van elk ter verificatie), verspreid over meerdere dagen en verschillende reagentia-batches. ICH: lineaire respons nabij de limiet.
Veelgemaakte fouten. De rekenmachine accepteert 10 of meer herhalingen; de EP17-richtlijnen gebruiken ongeveer 60 blanco- en 60 laag-niveau-resultaten om een limiet vast te stellen en minstens 20 om een geclaimde limiet te verifiëren. Deze parametrische methode gaat uit van één enkel laag-niveau-monster; bij meerdere monsters of niet-Gaussiaanse blanco's dient de niet-parametrische procedure te worden gebruikt. Bevestig de limieten experimenteel.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- CLSI. Evaluatie van de detectiemogelijkheden voor klinische laboratoriummeetprocedures; EP17-A2. 2012.
- ICH Q2(R1). Validatie van analytische procedures: tekst en methodologie. 2005.
- Armbruster DA, Pry T. Grenswaarde van blanco, detectielimiet en kwantificatielimiet. Clin Biochem Rev. 2008;29(Suppl 1):S49-S52.
Verwacht formaat: Referentiepopulatiewaarden
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht formaat: Referentiepopulatiewaarden
x <- c({data}) quantile(x, c(0.025, 0.975), type = 6) # niet-parametrisch (CLSI C28-A3) mean(x) + c(-1, 1) * qnorm(0.975) * sd(x) # parametrischGebruikershandleiding
Wanneer te gebruiken. Stel een referentie-interval vast (de centrale 95% van een gezonde referentiepopulatie) voor een laboratoriumanalyt.
Aannames. De waarden zijn afkomstig van ten minste 120 goed gedefinieerde referentiepersonen (CLSI C28-A3); uitschieters zijn beoordeeld en niet blindelings verwijderd.
Veelgemaakte fouten. Bij minder dan 120 waarden moeten de grenzen als benadering worden beschouwd. Maak een indeling op leeftijd of geslacht wanneer de verdelingen duidelijk verschillen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- CLSI. Het definiëren, vaststellen en verifiëren van referentie-intervallen in het klinisch laboratorium; Goedgekeurde richtlijn C28-A3. 2008.
- Horn PS, Pesce AJ. Referentie-intervallen: een update. Clin Chim Acta. 2003;334:5-23.
Verwacht: 2 kolommen: Verwachte concentraties | Gemeten resultaten
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Verwachte concentraties | Gemeten resultaten
x <- c({x}) y <- c({y}) summary(lm(y ~ x + I(x^2)))$coëfficiënten # rij 3 = niet-lineariteitstest summary(lm(y ~ x))$r.squaredGebruikershandleiding
Wanneer te gebruiken. Controleer of een test lineair reageert over het gehele meetbereik door te testen op een kwadratische component (de polynomiale benadering van de eerdere CLSI EP6-A).
Aannames. Minstens vijf even informatieve niveaus; herhalingen gemiddeld per niveau; constante variantie over het bereik.
Veelgemaakte fouten. Een hoge R²-waarde alleen bewijst geen lineariteit. Niet significant bewijst geen lineariteit bij zeer weinig significantieniveaus: controleer de residuen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- CLSI. Evaluatie van de lineariteit van kwantitatieve meetprocedures; EP06. 2e editie. 2020.
Verwacht: 3 kolommen: Variabele X | Variabele Y | Controlevariabele Z
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 3 kolommen: Variabele X | Variabele Y | Controlevariabele Z
x <- c({x}) y <- c({y}) z <- c({z}) rx <- resid(lm(x ~ z)); ry <- resid(lm(y ~ z)) r <- cor(rx, ry) t <- r * sqrt((length(x) - 3) / (1 - r^2)) c(partial_r = r, t = t, p = 2 * pt(-abs(t), length(x) - 3))Gebruikershandleiding
Wanneer te gebruiken. Om de samenhang tussen X en Y te schatten na het verwijderen van het lineaire effect van een derde variabele Z (een verstorende factor).
Aannames: Lineaire verbanden; variabelen bij benadering normaal verdeeld; Z wordt zonder fouten gemeten.
Veelgemaakte fouten. Partiële correlatie bewijst niet dat Z een verstorende factor is of dat X Y veroorzaakt. Gebruik voor meerdere covariaten meervoudige regressie.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Kim S. ppcor: een R-pakket voor een snelle berekening van semi-partiële correlatiecoëfficiënten. Commun Stat Appl Methods. 2015;22(6):665-674.
e1 <- {e1}; n1 <- {n1}; e0 <- {e0}; n0 <- {n0} margin <- {margin}/100 p1 <- e1/n1; p0 <- e0/n0 se <- sqrt(p1*(1-p1)/n1 + p0*(1-p0)/n0) z <- (p1 - p0 + margin)/se # non-inferioriteit (hoger is beter) c(diff = p1 - p0, z = z, p_one_sided = 1 - pnorm(z)) # Equivalentie (TOST): z2 <- (margin - (p1 - p0))/se max(1 - pnorm(z), 1 - pnorm(z2)) # TOST p-waardeGebruikershandleiding
Wanneer te gebruiken. Aantonen dat een nieuwe behandeling niet meer dan een vooraf vastgestelde marge slechter is dan de controlegroep (non-inferioriteit), of dat twee behandelingen in beide richtingen minder dan die marge van elkaar verschillen (equivalentie, TOST), voor een binaire uitkomst.
Aannames. De marge werd in het protocol vastgelegd voordat de gegevens beschikbaar waren en is klinisch gerechtvaardigd; voldoende grote steekproeven voor de Wald-benadering.
Veelgemaakte fouten. Bepaal de marge niet pas na het bekijken van de gegevens. Gebruik bij kleine cellen of percentages dicht bij 0 of 100% een scoremethode (Farrington-Manning). Niet-significantie in een superioriteitstest is geen bewijs van equivalentie.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Blackwelder WC. Het bewijzen van de nulhypothese in klinische onderzoeken. Control Clin Trials. 1982;3:345-353.
- Schuirmann DJ. Een vergelijking van de twee eenzijdige testprocedures en de power-benadering voor het beoordelen van de equivalentie van de gemiddelde biologische beschikbaarheid. J Pharmacokinet Biopharm. 1987;15:657-680.
- Piaggio G, et al. Rapportage van gerandomiseerde onderzoeken naar non-inferioriteit en equivalentie (CONSORT-uitbreiding). JAMA. 2012;308(24):2594-2604.
Verwacht: 4 kolommen: Reeks AB, Periode 1 | Reeks AB, Periode 2 | Reeks BA, Periode 1 | Reeks BA, Periode 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 4 kolommen: Reeks AB, Periode 1 | Reeks AB, Periode 2 | Reeks BA, Periode 1 | Reeks BA, Periode 2
ab1 <- c({ab1}); ab2 <- c({ab2}) # volgorde AB: periode 1 (A), periode 2 (B) ba1 <- c({ba1}); ba2 <- c({ba2}) # volgorde BA: periode 1 (B), periode 2 (A) dA <- ab1 - ab2; dB <- ba1 - ba2 t.test(dA, dB, var.equal = TRUE) # behandelingseffect A - B = (gemiddelde(dA) - gemiddelde(dB)) / 2 t.test(dA, -dB, var.equal = TRUE) # periode-effect t.test(ab1 + ab2, ba1 + ba2, var.equal = TRUE) # carryover (verschil in subjecttotalen)Gebruikershandleiding
Wanneer te gebruiken. Analyseer een crossover-onderzoek met twee perioden en twee behandelingen (AB/BA) met een continue uitkomst: schat het behandelingseffect en test op periode- en carry-over-effecten (Hills-Armitage / Grizzle).
Aannames: Deelnemers willekeurig toegewezen aan de twee behandelingsreeksen; voldoende uitwasperiode; nagenoeg normale verschillen binnen de proefpersonen; geen interactie tussen behandeling en periode.
Veelgemaakte fouten. De carryover-test heeft een lage power en mag niet als poortwachter voor de analyse worden gebruikt (Senn); bepaal de washout in het onderzoeksontwerp. Niet gebruiken voor aandoeningen die vanzelf verdwijnen of behandelingen met blijvende effecten.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Senn S. Cross-overonderzoek in klinisch onderzoek. 2e editie. Wiley; 2002.
- Hills M, Armitage P. Het tweeperioden cross-over klinisch onderzoek. Br J Clin Pharmacol. 1979;8:7-20.
- Grizzle JE. Het tweefasen-overgangsontwerp en het gebruik ervan in klinische onderzoeken. Biometrics. 1965;21:467-480.
power.t.test(delta = {delta}, sd = sqrt(2) * {sigma}, sig.level = {alpha}, power = {power}, type = "paired")$n # totaal NGebruikershandleiding
Wanneer te gebruiken. Plan een crossover-onderzoek met twee perioden en twee volgordes (AB/BA) waarbij twee behandelingen worden vergeleken op basis van een continue uitkomstmaat.
Aannames: Geen carry-over-effect of interactie tussen periode en behandeling; gelijke aantallen in de twee sequenties; binnen-subject standaarddeviatie uit een eerder cross-over onderzoek.
Veelgemaakte fouten. Niet gebruiken voor parallelgroepstudies of wanneer overdrachtseffecten waarschijnlijk zijn (langdurige behandelingen, geen uitwasperiode). Bio-equivalentie vereist een TOST-gebaseerde berekening.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Senn S. Cross-overonderzoek in klinisch onderzoek. 2e editie. Wiley; 2002.
- Chow SC, Shao J, Wang H. Steekproefgrootteberekeningen in klinisch onderzoek. 2e editie. 2008.
# Klassieke (vaste vorm) groepssequentiële grenzen, gelijkmatig verdeelde looks K <- {looks}; alpha <- {alpha}; type <- "{interim_type}" # "OF" = O'Brien-Fleming, "Pocock" w <- if (type == "Pocock") rep(1, K) else sqrt(K / seq_len(K)) surv <- function(C, NG = 240) { simp <- c(1, rep(c(4, 2), length.out = NG - 1), 1) for (k in seq_len(K)) { b <- C * w[k] * sqrt(k); hn <- 2 * b / NG; ng <- -b + (0:NG) * hn nd <- als (k == 1) dnorm(ng) anders sapply(ng, functie(s) som(simp * dens * dnorm(s - grid)) * h / 3) s <- som(simp * nd) * hn / 3; grid <- ng; dens <- nd; h <- hn } s } C <- uniroot(functie(C) 1 - surv(C) - alpha, c(1, 6), tol = 1e-10)$root round(C * w, 3) # kritieke |z| bij elke lookGebruikershandleiding
Wanneer te gebruiken. Stel effectiviteitsgrenzen in voor een onderzoek met geplande tussentijdse controles, zodat de algehele type I-fout op α blijft.
Aannames. Gelijkmatig verdeelde observaties (gelijke informatie per fase); tweezijdige toets van een normaal verdeelde statistiek; grenzen vastgesteld vóór de ontblinding.
Veelgemaakte fouten. Voeg geen ongeplande looks toe. Gebruik voor ongelijke spatiëring de alfa-spending-methode (Lan-DeMets) in gespecialiseerde software.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- O'Brien PC, Fleming TR. Een procedure voor meervoudige testen voor klinische onderzoeken. Biometrics. 1979;35(3):549-556.
- Pocock SJ. Groepssequentiële methoden bij het ontwerp en de analyse van klinische onderzoeken. Biometrika. 1977;64(2):191-199.
- Jennison C, Turnbull BW. Groepssequentiële methoden met toepassingen in klinische onderzoeken. CRC; 2000.
set.seed({seed}) arms <- trimws(strsplit("{arms}", ",")[[1]]); blocks <- c({blocks}); n <- {n} strata <- trimws(strsplit("{strata}", ",")[[1]]); scheme <- "{scheme}" gen <- function(n) { if (grepl("Simple", scheme)) return(arms[sample.int(length(arms), n, replace = TRUE)]) out <- character(0) while (length(out) < n) { b <- blocks[sample.int(length(blocks), 1)] out <- c(out, sample(rep(arms, length.out = b))) } out[1:n] } if (grepl("Stratified", scheme)) lapply(setNames(strata, strata), function(s) gen(n)) else gen(n)Gebruikershandleiding
Wanneer te gebruiken. Genereer een reproduceerbare toewijzingslijst (eenvoudig, gepermuteerd blok of gestratificeerd) met een vastgelegde seed. Het R-commando reproduceert de exacte lijst.
Aannames. De lijst wordt beheerd door iemand die onafhankelijk is van de werving (toewijzingsverhulling); de startwaarde en de lijst worden gearchiveerd.
Veelgemaakte fouten. Gebruik geen vaste blokgrootte in niet-geblindeerde onderzoeken (de laatste toewijzingen worden dan voorspelbaar). Eenvoudige randomisatie kan leiden tot onevenwichtige resultaten bij kleine onderzoeken. Dynamische methoden zoals minimalisatie vereisen speciale software.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Schulz KF, Grimes DA. Het genereren van toewijzingssequenties in gerandomiseerde onderzoeken. Lancet. 2002;359:515-519.
- Kim J, Shin W. Hoe voer je willekeurige toewijzing (randomisatie) uit? Clin Orthop Surg. 2014;6(1):103-109.
e1 <- {e1}; n1 <- {n1}; e0 <- {e0}; n0 <- {n0} pt <- e1/n1; pc <- e0/n0 arr <- pc - pt se <- sqrt(pt*(1-pt)/n1 + pc*(1-pc)/n0) ci <- arr + c(-1, 1) * qnorm(0.975) * se c(ARR = arr, NNT = 1/arr) 1/rev(ci) # NNT 95% BI (geldig als de ARR CI 0 uitsluit)Gebruikershandleiding
Wanneer te gebruiken. Rapporteer de behandelingseffecten als absolute risicoreductie, relatief risico, odds ratio en number needed to treat (NNT), elk met een 95% betrouwbaarheidsinterval, plus een Cates-grafiek die laat zien wat er gebeurt met 100 behandelde personen.
Aannames. Binaire uitkomst over een vaste follow-up; redelijk grote aantallen voor de Wald-intervallen.
Veelgemaakte fouten. Wanneer het ARR-interval nul omvat, loopt het NNT-interval door tot in het oneindige (NNTB tot NNTH): geef het niet weer als een eenvoudig bereik. NNT is specifiek voor de follow-uptijd.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Altman DG. Betrouwbaarheidsintervallen voor het aantal patiënten dat behandeld moet worden om één gunstig resultaat te bereiken. BMJ. 1998;317:1309-1312.
- Laupacis A, Sackett DL, Roberts RS. Een beoordeling van klinisch bruikbare meetinstrumenten voor de gevolgen van een behandeling. N Engl J Med. 1988;318:1728-1733.
library(effectsize) d <- {value} # als er wordt geconverteerd van d d_to_r(d) # gelijke groepsgroottes d_to_oddsratio(d) J <- 1 - 3/(4*({n1}+{n2}-2) - 1); d * J # Hedges gGebruikershandleiding
Wanneer te gebruiken. Omrekenen tussen Cohen's d, Hedges' g, de correlatiecoëfficiënt r en de odds ratio, zodat studies in een meta-analyse kunnen worden samengevoegd.
Aannames. De onderliggende continue uitkomst is in beide groepen bij benadering normaal verdeeld (de logistische conversie naar OR gaat hiervan uit).
Veelgemaakte fouten. Converteer niet wanneer de oorspronkelijke gegevens echt binair zijn: gebruik de gerapporteerde odds ratio direct. Conversies tussen d en r gaan uit van een ontwerp met twee groepen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Borenstein M, Hedges LV, Higgins JPT, Rothstein HR. Inleiding tot meta-analyse. Wiley; 2009. Hoofdstuk 7.
- Chinn S. Een eenvoudige methode om een odds ratio om te zetten naar een effectgrootte. Stat Med. 2000;19:3127-3131.
tab <- matrix(c({a}, {c}, {b}, {d}), nrow = 2) # rijen: voor +/-, kolommen: na +/- mcnemar.test(tab) binom.test({b}, {b} + {c}) # exacte p-waardeGebruikershandleiding
Wanneer te gebruiken. Vergelijk gepaarde binaire uitkomsten (dezelfde proefpersonen voor/na, of overeenkomende case-controlparen).
Aannames. Paren zijn onafhankelijk van elkaar; de uitkomst is dichotoom op beide tijdstippen; alleen discordante paren bevatten informatie.
Veelgemaakte fouten. Niet gebruiken voor onafhankelijke groepen (gebruik de chi-kwadraat- of Fisher-toets). Bij minder dan ongeveer 25 afwijkende paren is de exacte p-waarde de voorkeur.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- McNemar Q. Opmerking over de steekproeffout van het verschil tussen gecorreleerde proporties. Psychometrika. 1947;12(2):153-157.
- Fagerland MW, Lydersen S, Laake P. De McNemar-test voor binaire matched-pairs-gegevens. BMC Med Res-methode. 2013;13:91.
Verwacht: 3 kolommen: Werkelijke status, één per onderwerp | Scores beoordelaar 1 | Scores beoordelaar 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 3 kolommen: Werkelijke status, één per onderwerp | Scores beoordelaar 1 | Scores beoordelaar 2
status <- c({status}) s1 <- c({score1}) s2 <- c({score2}) if (grepl("^Lower", "{dir1}")) s1 <- -s1 # maak "hoger = positiever" if (grepl("^Lower", "{dir2}")) s2 <- -s2 alpha <- {alpha} # DeLong (1988) structurele componenten, alleen basis R psi <- function(x, y) (x > y) + 0.5 * (x == y) comp <- function(s) { M <- outer(s[status == 1], s[status == 0], psi) # positieven x negatieven list(V10 = rowMeans(M), V01 = colMeans(M)) } c1 <- comp(s1); c2 <- comp(s2) m <- sum(status == 1); n <- sum(status == 0) auc <- c(mean(c1$V10), mean(c2$V10)) S10 <- cov(cbind(c1$V10, c2$V10)) S01 <- cov(cbind(c1$V01, c2$V01)) S <- S10 / m + S01 / n # covariantiematrix van (AUC1, AUC2) d <- auc[1] - auc[2] se <- sqrt(S[1, 1] + S[2, 2] - 2 * S[1, 2]) z <- d / se p <- 2 * pnorm(-abs(z)) q <- qnorm(1 - alpha / 2) ci <- function(a, v) pmin(1, pmax(0, a + c(-1, 1) * q * sqrt(v))) ci1 <- ci(auc[1], S[1, 1]); ci2 <- ci(auc[2], S[2, 2]) cid <- d + c(-1, 1) * q * se rho <- S[1, 2] / sqrt(S[1, 1] * S[2, 2]) cat(sprintf("n positieven = %d, n negatieven = %d\n", m, n)) cat(sprintf("AUC1 = %.4f [%.4f, %.4f]\n", auc[1], ci1[1], ci1[2])) cat(sprintf("AUC2 = %.4f [%.4f, %.4f]\n", auc[2], ci2[1], ci2[2])) cat(sprintf("Verschil = %.4f [%.4f, %.4f], z = %.4f, p = %.4f\n", d, cid[1], cid[2], z, p)) cat(sprintf("Correlatie van AUC's = %.4f\n", rho)) # Optionele kruiscontrole als pROC is geïnstalleerd: # pROC::roc.test(pROC::roc(status, s1, direction = "<"), # pROC::roc(status, s2, direction = "<"), method = "delong")Gebruikershandleiding
Wanneer te gebruiken. Test of twee diagnostische tests of biomarkers, gemeten bij dezelfde proefpersonen, verschillen in onderscheidend vermogen (oppervlakte onder de ROC-curve).
Aannames. Elk subject heeft zowel scores als een werkelijke 0/1-status; de subjecten zijn onafhankelijk van elkaar; de twee AUC's zijn gecorreleerd omdat ze dezelfde subjecten delen (DeLong houdt hier rekening mee). Gelijkstanden worden afgehandeld met de midrank-conventie (0.5).
Veelgemaakte fouten. Niet gebruiken voor markers die bij verschillende proefpersonen zijn gemeten (gebruik een ongepaarde vergelijking). Bij zeer weinig positieve of negatieve resultaten is de normale benadering onbetrouwbaar; een niet-significant resultaat is geen bewijs van gelijke prestaties. Vergelijk partiële AUC's of klinisch relevante drempelwaarden wanneer de volledige curve niet relevant is.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- DeLong ER, DeLong DM, Clarke-Pearson DL. Het vergelijken van de oppervlakten onder twee of meer gecorreleerde receiver operating characteristic-curven: een niet-parametrische benadering. Biometrics. 1988;44(3):837-845.
- Hanley JA, McNeil BJ. Een methode om de oppervlakten onder receiver operating characteristic-curven te vergelijken die zijn afgeleid van dezelfde gevallen. Radiologie. 1983;148(3):839-843.
- Robin X, Turck N, Hainard A, et al. pROC: een open-sourcepakket voor R en S+ om ROC-curven te analyseren en te vergelijken. BMC Bioinformatics. 2011;12:77.
Verwacht: 2 kolommen: Ziektestatus per proefpersoon | Markerwaarde per proefpersoon
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 2 kolommen: Ziektestatus per proefpersoon | Markerwaarde per proefpersoon
status <- c({status}) score <- c({score}) higher <- "{dir}" == "Hoger = positiever" method <- "{method}" target <- {target} / 100 cr <- {costratio} prev <- suppressWarnings(as.numeric("{prev}")) / 100 # leeg -> NA -> prevalentie van de steekproef P <- sum(status == 1); N <- sum(status == 0) if (is.na(prev)) prev <- P / (P + N) u <- sort(unique(score)) cut <- (head(u, -1) + tail(u, -1)) / 2 # kandidaten: middelpunten tussen verschillende waarden posf <- function(s, c) if (higher) s >= c else s <= c sens <- sapply(cut, function(c) mean(posf(score[status == 1], c))) spec <- sapply(cut, function(c) mean(!posf(score[status == 0], c))) obj <- if (grepl("Youden", method)) sens + spec - 1 else if (grepl("top-left", method)) -((1 - sens)^2 + (1 - spec)^2) else if (grepl("sensitivity", method)) ifelse(sens >= target - 1e-12, spec, -Inf) anders als (grepl("specificiteit", methode)) als anders (spec >= target - 1e-12, sens, -Inf) anders -(cr * prev * (1 - sens) + (1 - prev) * (1 - spec)) stop als niet (is.finite(max(obj))) # vaste methoden: doel moet haalbaar zijn i <- which(obj >= max(obj) - 1e-12)[1] # gelijkstanden -> laagste afkapwaarde c(cutoff = cut[i], sensitivity = sens[i], specificity = spec[i], J = sens[i] + spec[i] - 1) cp <- function(x, n) c(als (x == 0) 0 anders qbeta(.025, x, n - x + 1), als (x == n) 1 anders qbeta(.975, x + 1, n - x)) cp(round(sens[i] * P), P); cp(round(spec[i] * N), N) # exact (Clopper-Pearson) 95% CI # AUC met DeLong CI: library(pROC); r <- roc(status, score, direction = ifelse(higher, "<", ">")); auc(r); ci.auc(r, method = "delong")Gebruikershandleiding
Wanneer te gebruiken. Kies een enkele beslissingsdrempel voor een continue diagnostische marker en bekijk de sensitiviteit, specificiteit, voorspellende waarde, likelihood ratio en AUC bij die drempel. Vijf criteria worden naast elkaar vergeleken: Youden-index, de waarde die het dichtst bij de linkerbovenhoek ligt, vaste sensitiviteit, vaste specificiteit en de verwachte kosten gegeven de prevalentie en de relatieve kosten van een vals negatief resultaat.
Aannames. Eén markerwaarde per proefpersoon met een echte 0/1 referentiestandaardstatus; onafhankelijke proefpersonen. De drempelwaarden zijn de middelpunten tussen opeenvolgende, afzonderlijke, gesorteerde markerwaarden (dezelfde drempelwaarden die pROC gebruikt); een proefpersoon is testpositief wanneer de markerwaarde ≥ de drempelwaarde is (of ≤ als lagere waarden positiever zijn). Voorspellende waarden bij een andere prevalentie zijn alleen geldig als de sensitiviteit en specificiteit ook van toepassing zijn op die populatie.
Veelgemaakte fouten. De grenswaarde en de nauwkeurigheid ervan worden gekozen en geschat op basis van dezelfde gegevens en zijn optimistisch; valideer met nieuwe gegevens of door middel van bootstrapping, vooral bij minder dan ongeveer 30 proefpersonen per klasse. Verschillende criteria leiden terecht tot verschillende grenswaarden (Perkins & Schisterman). Gebruik de steekproefprevalentie niet voor PPV/NPV wanneer gevallen en controles afzonderlijk zijn bemonsterd; voer in plaats daarvan de werkelijke prevalentie in. Het AUC-betrouwbaarheidsinterval is een normale benadering (DeLong) afgekapt tot [0, 1] en is onbetrouwbaar in de buurt van 1 of bij zeer weinig proefpersonen.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Youden WJ. Index voor het beoordelen van diagnostische tests. Kanker. 1950;3(1):32-35.
- Clopper CJ, Pearson ES. Het gebruik van betrouwbaarheids- of fiduciaire grenzen geïllustreerd aan de hand van de binomiale verdeling. Biometrika. 1934;26(4):404-413.
- DeLong ER, DeLong DM, Clarke-Pearson DL. Het vergelijken van de oppervlakten onder twee of meer gecorreleerde receiver operating characteristic-curven: een niet-parametrische benadering. Biometrics. 1988;44(3):837-845.
- Perkins NJ, Schisterman EF. De inconsistentie van “optimale” grenswaarden verkregen met behulp van twee criteria gebaseerd op de receiver operating characteristic curve. Am J Epidemiol. 2006;163(7):670-675.
- Robin X, Turck N, Hainard A, et al. pROC: een open-sourcepakket voor R en S+ om ROC-curven te analyseren en te vergelijken. BMC Bioinformatics. 2011;12:77.
Verwacht: 3 kolommen: Waargenomen uitkomst, één waarde per proefpersoon | Voorspelde waarschijnlijkheid, model/marker 1 | Voorspelde waarschijnlijkheid, model 2
Sleep hier een CSV-bestand naartoe of klik om te kiezen. Het bestand wordt in uw browser verwerkt: er wordt niets geüpload.
Verwacht: 3 kolommen: Waargenomen uitkomst, één waarde per proefpersoon | Voorspelde waarschijnlijkheid, model/marker 1 | Voorspelde waarschijnlijkheid, model 2
y <- c({status}) # waargenomen uitkomst, 1 = gebeurtenis p1 <- c({p1}) # voorspeld risico, model 1 p2 <- c({p2}) # voorspeld risico, model 2 (laat c() leeg voor één model) thr <- round(seq({tmin}, {tmax}, by = {tstep}), 10) nb <- function(p, y, pt) { # netto voordeel; Behandelen wanneer voorspeld risico >= drempelwaarde behandelen <- p >= pt - 1e-12 som(behandelen & y == 1) / lengte(y) - som(behandelen & y == 0) / lengte(y) * pt / (1 - pt) } vorige <- gemiddelde(y) dca <- data.frame(drempelwaarde = thr, alles behandelen = vorige - (1 - vorige) * thr / (1 - thr), niets behandelen = 0, model1 = sapply(thr, functie(t) nb(p1, y, t))) als (lengte(p2)) dca$model2 <- sapply(thr, functie(t) nb(p2, y, t)) print(dca[dca$drempelwaarde %in% c(0.05, 0.10, 0.20, 0.30, 0.40, 0.50), ]) matplot(dca$threshold, dca[, -1], type = "l", lty = c(2, 1, 1, 1), lwd = 2, col = c("grey50", "black", "#1a5276", "#c0504d"), ylim = c(-0.05, prev), xlab = "Drempelwaarschijnlijkheid", ylab = "Netto voordeel") legend("topright", c("Behandel alles", "Behandel niets", "{label1}", if (length(p2)) "{label2}"), lty = c(2, 1, 1, 1), col = c("grey50", "black", "#1a5276", "#c0504d"), lwd = 2, bty = "n")Gebruikershandleiding
Wanneer te gebruiken? Beoordeel of een voorspellingsmodel of biomarker klinische beslissingen kan verbeteren: het netto voordeel weegt de ware positieven af tegen de valse positieven bij de risicodrempel waar een patiënt of arts zou handelen, en vergelijkt het model met het behandelen van iedereen of niemand.
Aannames. De voorspelde waarden zijn gekalibreerde waarschijnlijkheden tussen 0 en 1 voor een binaire uitkomst die bij elke proefpersoon is waargenomen (geen censurering); dezelfde proefpersonen worden gebruikt voor alle modellen; het drempelbereik weerspiegelt de risico's waarbij een beslissing plausibel zou veranderen. De curves op basis van de gegevens die zijn gebruikt om het model te fitten, zijn optimistisch.
Veelgemaakte fouten. Gebruik geen ruwe scores die geen waarschijnlijkheden zijn, lees geen curves buiten een klinisch relevant drempelbereik en claim geen voordeel op basis van kleine verschillen zonder validatie: deze tool toont geen betrouwbaarheidsintervallen en corrigeert niet voor overfitting. Netto voordeel is geen nauwkeurigheidsmaatstaf zoals AUC.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Vickers AJ, Elkin EB. Beslissingscurve-analyse: een nieuwe methode voor het evalueren van voorspellingsmodellen. Med Decis Making. 2006;26(6):565-574. doi:10.1177/0272989X06295361
- Vickers AJ, van Calster B, Steyerberg EW. Netto-batenbenaderingen voor de evaluatie van voorspellingsmodellen, moleculaire markers en diagnostische tests. BMJ. 2016;352:i6. doi:10.1136/bmj.i6
- Vickers AJ, van Calster B, Steyerberg EW. Een eenvoudige, stapsgewijze handleiding voor het interpreteren van beslissingscurve-analyse. Diagn Progn Res. 2019;3:18. doi:10.1186/s41512-019-0064-7
AUC: Var(A) = [A(1−A) + (n1−1)(Q1A2) + (n0−1)(Q2A2)] / (N1n0), Q1 = A/(2−A), Q2 = 2A2/(1+A). Kleinste n1 (n0 = plafond(verhouding·n)1)) met A−A0 ≥ z1−α/2√Var0 + vanenergie√Var1; precisie-only: z1−α/2√Var(A) ≤ d. Gekoppeld: Var(A1A2) = Var(A1) + Var(A2) − 2r√(Var(A1)Var(A2)).
# Steekproefomvang voor diagnostische nauwkeurigheidsstudies (basis R, geen pakketten) doel <- "{doel}"; alfa <- {alfa} za <- qnorm(1 - alfa / 2) als (grepl("Buderer", doel)) { # Buderer 1996 / Flahault 2005: n = z^2 P (1 - P) / d^2, dan delen door prevalentie of 1 - prevalentie se <- {se} / 100; sp <- {sp} / 100; prev <- {prev} / 100; d <- {d} / 100 n_dis <- za^2 * se * (1 - se) / d^2 # zieke personen nodig voor sensitiviteit (niet afgerond) n_non <- za^2 * sp * (1 - sp) / d^2 # niet-zieke personen nodig voor specificiteit (niet afgerond) tot_se <- ceiling(n_dis / prev); tot_sp <- ceiling(n_non / (1 - prev)) c(diseased = ceiling(n_dis), nondiseased = ceiling(n_non), total_sensitivity = tot_se, total_specificity = tot_sp, required_total = max(tot_se, tot_sp)) } else { # Hanley & McNeil (1982) variantie van een AUC; n0 = ceiling(ratio * n1) A <- {auc}; A0 <- {auc0}; A2 <- {auc2}; r <- {rho}; ratio <- {ratio} pw <- {power} / 100; dd <- {d} / 100; zb <- qnorm(pw) hmvar <- function(A, n1, n0) { Q1 <- A / (2 - A); Q2 <- 2 * A^2 / (1 + A) (A * (1 - A) + (n1 - 1) * (Q1 - A^2) + (n0 - 1) * (Q2 - A^2)) / (n1 * n0) } # Hanley & McNeil (1983): Var(A1 - A2) = Var(A1) + Var(A2) - 2 r SE1 SE2 vd <- functie(A1, A2, n1, n0) { v1 <- hmvar(A1, n1, n0); v2 <-hmvar(A2, n1, n0); v1 + v2 - 2 * r * sqrt(v1 * v2) } if (grepl("gecorreleerd", doel)) { Ab <- (A + A2) / 2; D <- abs(A - A2) ok <- functie(n1) { n0 <- plafond(verhouding * n1); D >= za * sqrt(vd(Ab, Ab, n1, n0)) + zb * sqrt(vd(A, A2, n1, n0)) } okp <- functie(n1) { n0 <- plafond(verhouding * n1); za * sqrt(vd(A, A2, n1, n0)) <= dd } } else { ok <- functie(n1) { n0 <- plafond(verhouding * n1); A - A0 >= za * sqrt(hmvar(A0, n1, n0)) + zb * sqrt(hmvar(A, n1, n0)) } okp <- functie(n1) { n0 <- plafond(verhouding * n1); za * sqrt(hmvar(A, n1, n0)) <= dd } } n1 <- 2; while (!ok(n1)) n1 <- n1 + 1 p1 <- 2; while (!okp(p1)) p1 <- p1 + 1 c(positives = n1, negatives = ceiling(ratio * n1), total = n1 + ceiling(ratio * n1), precision_only_total = p1 + ceiling(ratio * p1)) }Gebruikershandleiding
Wanneer te gebruiken. Bepaal hoeveel deelnemers een onderzoek naar diagnostische nauwkeurigheid nodig heeft: om de sensitiviteit en specificiteit te schatten met een gekozen precisie (Buderer), om aan te tonen dat een AUC een nulwaarde overschrijdt, of om de AUC's van twee tests die bij dezelfde proefpersonen zijn afgenomen te vergelijken.
Aannames. Buderer: Wald-intervallen (normale benadering), dus d moet klein genoeg zijn zodat het interval binnen 0 tot 100% blijft; prevalentie is de prevalentie in de gerekruteerde steekproef. AUC: de Hanley-McNeil-variantie is een vloeiende benadering die uitgaat van ruwweg binormale scoreverdelingen; de resultaten zijn voor de empirische (Mann-Whitney) AUC en een normale toetsingsstatistiek. Gepaarde vergelijking: r is de correlatie tussen de twee AUC-schattingen, zoals ingevoerd, en deze moet afkomstig zijn uit pilotgegevens of de literatuur.
Veelgemaakte fouten. Lees het Buderer-resultaat niet als een powerberekening voor een hypothesetoets (het is alleen gericht op de breedte van het betrouwbaarheidsinterval). Gebruik de AUC-resultaten niet voor zeer kleine steekproeven (een paar gevallen van de ziekte), voor gedeeltelijke AUC of voor geclusterde/meerdere beoordelaarsgegevens. Gebruik voor een exacte of zeer hoge/lage sensitiviteit een exact-binomiale planning of een specifiek softwarepakket.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Buderer NMF. Statistische methodologie: I. Het meenemen van de prevalentie van de ziekte in de berekening van de steekproefomvang voor sensitiviteit en specificiteit. Acad Emerg Med. 1996;3(9):895-900.
- Flahault A, Cadilhac M, Thomas G. Een steekproefomvangberekening is noodzakelijk voor een nauwkeurig onderzoeksontwerp bij diagnostische teststudies. J Clin Epidemiol. 2005;58(8):859-862.
- Hanley JA, McNeil BJ. De betekenis en het gebruik van het gebied onder een receiver operating characteristic (ROC)-curve. Radiologie. 1982;143(1):29-36.
- Hanley JA, McNeil BJ. Een methode om de oppervlakten onder receiver operating characteristic-curven te vergelijken die zijn afgeleid van dezelfde gevallen. Radiologie. 1983;148(3):839-843.
- Obuchowski NA. Berekeningen van de steekproefomvang in studies naar de nauwkeurigheid van tests. Stat Methods Med Res. 1998;7(4):371-392.
tp <- {tp}; fp <- {fp}; fn <- {fn}; tn <- {tn}; pretest <- {pretest} # pre-testkans, % if (any(c(tp, fp, fn, tn) == 0)) { tp <- tp + .5; fp <- fp + .5; fn <- fn + .5; tn <- tn + .5 } # correctie voor nulcellen sens <- tp/(tp+fn); spec <- tn/(tn+fp) # Vanuit sens/spec in plaats daarvan: sens <- {sens}/100; spec <- {spec}/100 lrp <- sens/(1-spec); lrn <- (1-sens)/spec # Vanuit LRs in plaats daarvan: lrp <- {lrp}; lrn <- {lrn} se_p <- sqrt(1/tp - 1/(tp+fn) + 1/fp - 1/(fp+tn)); se_n <- sqrt(1/fn - 1/(tp+fn) + 1/tn - 1/(fp+tn)) ci <- function(l, se) exp(log(l) + c(-1, 1) * qnorm(0.975) * se) odds <- (pretest/100) / (1 - pretest/100) post <- function(l) { o <- odds * l; 100 * o / (1 + o) } list(LRplus = c(lrp, ci(lrp, se_p)), LRminus = c(lrn, ci(lrn, se_n)), post_positive = post(lrp), post_negative = post(lrn))Gebruikershandleiding
Wanneer te gebruiken. Zet een testresultaat om in een post-testkans: voer de sensitiviteit/specificiteit, een 2×2-tabel of gepubliceerde likelihood ratio's in, stel de pre-testkans in en lees het antwoord af uit de tabel of het Fagan-nomogram.
Aannames. Gevoeligheid en specificiteit zijn overdraagbaar naar uw situatie; de pre-testkans is realistisch (een prevalentie of een klinische schatting); het resultaat is dichotoom. Bij 2×2 tellingen gebruikt het betrouwbaarheidsinterval de logaritmische methode (Simel 1991), die uitgaat van voldoende celtellingen. Als een cel nul is, wordt 0.5 bij elke cel opgeteld voordat de likelihood ratio (LR) en het betrouwbaarheidsinterval worden berekend (continuïteitscorrectie), en dat staat ook op de pagina vermeld.
Veelgemaakte fouten. Combineer likelihood ratio's van gecorreleerde tests niet alsof ze onafhankelijk zijn. De likelihood ratio (LR) van sensitiviteit en specificiteit, ingevoerd als percentages, heeft geen betrouwbaarheidsinterval (aantallen zijn nodig). Een specificiteit van 100% geeft LR+ = oneindig; in de tellingmodus geeft de correctie van 0.5 een eindige waarde, maar met een zeer breed interval.
Uitgewerkt voorbeeld
Klik hierboven op 'Voorbeeldgegevens laden' en vervolgens op 'Analyse uitvoeren' om een volledig uitgewerkt resultaat te zien met grafiek, formule en interpretatie voor deze test.
Referenties
- Fagan TJ. Nomogram voor de stelling van Bayes. N Engl J Med. 1975;293(5):257.
- Simel DL, Samsa GP, Matchar DB. Likelihood ratio's met betrouwbaarheid: schatting van de steekproefomvang voor diagnostische teststudies. J Clin Epidemiol. 1991;44(8):763-770.
- Deeks JJ, Altman DG. Diagnostische tests 4: waarschijnlijkheidsverhoudingen. BMJ. 2004;329(7458):168-169.
Alle gratis onderzoekstools
Alle onderstaande tools werken volledig in uw browser — geen account nodig, geen uploads vereist en er worden geen gegevens van uw apparaat verzonden. Elke tool heeft een eigen pagina met de formule, de onderliggende aannames, een uitgewerkt voorbeeld, primaire bronnen en de reproduceerbare R-opdracht.
Kernstatistieken
Klinisch en diagnostisch
Overleving en regressie
Niet-parametrische tests
Overeenstemming en betrouwbaarheid
Laboratorium- en methodevalidatie
Klinisch proefontwerp
Meta-analyse en bewijssynthese
Stroomschema's voor rapportagerichtlijnen
Een test kiezen
Verder lezen
Handleidingen en door vakgenoten beoordeelde bronnen om u te helpen bij het kiezen, uitvoeren en rapporteren van de juiste analyse.
- Welke statistische toets moet ik gebruiken? — interactieve selectie
- Calculator voor steekproefomvang en poweranalyse
- Inzicht in p-waarden en statistische significantie
- Meta-analyse: bosplots en heterogeniteit
- Alle 44 gratis statistische rekenmachines
- Greenland S, et al. Statistische toetsen, p-waarden, betrouwbaarheidsintervallen en power: een gids voor misinterpretaties. Eur J Epidemiol. 2016.PMID 27209009 · gratis volledige tekst
- Franc JM. Analyse van gepaarde gegevens. Prehosp Rampenmedicatie. 2025.PMID 40213959 · de fout bij gekoppelde gegevens
- Altman N, Krzywinski M. Inzicht in p-waarden en significantie. Lab Anim. 2024.PMID 39315628 · effectgrootte & meervoudige testen
Vraag om statistische hulp.
Vertel ons over uw data en de analyse die u nodig heeft. Een biostatisticus met een doctoraat reageert met het juiste pakket en een duidelijk plan – meestal binnen één werkdag.
