J'écris actuellement une méta-analyse en utilisant une méta-analyse d'effets aléatoires par paires pour comparer les taux de complication de 5 modalités de traitement, avec un comme étalon-or. J'ai réussi à faire sortir RR de l'analyse mais je n'ai pas pu en extraire les valeurs p. Comment puis-je obtenir les p-values de ce modèle? J'ai essayé d'utiliser pval.random
mais cela n'a pas fonctionné. Et je n'ai trouvé aucun autre code sur CRAN qui puisse m'aider. Quelqu'un peut-il m'aider avec le code R?P-valeur d'une méta-analyse de réseau d'effets aléatoires par paires
drf <- read.csv("drf zonder moroni.csv", sep = ";", header = TRUE, as.is = TRUE)
##
drf <- drf[, 1:5]
names(drf) <- c("study", "type", "treat", "events", "n")
compl <- subset(drf, type == "Complications")
library(netmeta)
p.compl <- pairwise(treat = treat, event = events, n = n,
studlab = study, data = compl)
n.compl <- netmeta(p.compl, reference = "PC", comb.random=TRUE)
n.compl
netgraph(n.compl, iterate = TRUE, number = TRUE)
partie du jeu de données
Study| Event Type| Treatment| Number of Events (n)| N| n/N|
Kumaravel| Complications| EF| 3| 23| 0,1304348|
Franck| Complications| EF| 2| 20| 0,1|
Schonnemann| Complications| EF| 8| 30| 0,2666667|
Aita| Complications| EF| 1| 16| 0,0625|
Hove| Complications| EF| 31| 39| 0,7948718|
Andersen| Complications| EF| 26| 75| 0,3466667|
Krughaug| Complications| EF| 22| 75| 0,2933333|
Moroni| Complications| EF| 0| 20| 0|
Plate| Complications| IMN| 3| 30| 0,1|
Chappuis| Complications| IMN| 4| 16| 0,25|
Gradl| Complications| IMN| 12| 66| 0,1818182|
Schonnemann| Complications| IMN| 6| 31| 0,1935484|
Aita| Complications| IMN| 1| 16| 0,0625|
Dremstrop| Complications| IMN| 17| 44| 0,3863636|
Wong| Complications| PC| 1| 30| 0,0333333|
Kumaravel| Complications| PC| 4| 25| 0,16|