dies ist mein erster Eintrag und ich hoffe, dass ich an alle benötigten Infos gedacht habe
Ich arbeite zur Zeit an einer Meta-Analyse und mit dem metafor package. Als vereinfachte Datenstruktur habe ich zwei Variablen:
yi für die Effektstärken
vi für die Varianz der Effektstärken
Hier findet ihr mein reproduzierbares Minimalbeispiel:
Code: Alles auswählen
data <- structure(list(yi = structure(c(-3.97957983193277, -3.78316455696203,
-3.44876946564885, -3.19228965517241, -2.97759120879121), format.spss = "F18.15", display_width = 16L),
vi = structure(c(1.46494705176188, 0.418507803236661, 0.279065509002972,
0.367668965517241, 0.170291317473735), format.spss = "F18.15", display_width = 16L)), row.names = c(NA,
-5L), class = c("tbl_df", "tbl", "data.frame"))
Ich möchte jetzt gern ein Funnel Plot mit der Trim and Fill Methode erstellen. Beim Vorgehen habe ich mich an dem Beispiel auf der Seite http://www.metafor-project.org/doku.php ... ]=trimfill orientiert. Das random-effects model zu erstellen ist kein Problem, das habe ich für andere Analysen auch schon gemacht und hatte damit nie Probleme. Sobald ich aber die Funktion trimfill ausführe, bekomme ich immer die folgende Fehlermeldung:
Fehler in if (!is.null(mods) && (nrow(mods) != k)) stop(mstyle$stop("Number of rows of the model matrix does not match length of the outcome vector.")) : Fehlender Wert, wo TRUE/FALSE nötig ist
Durch debug & traceback glaube ich herausgefunden zu haben, dass das Problem an folgender Stelle in der Funktion liegt:
if (is.vector(mods))
mods <- cbind(mods)
Die erstellte Variable "mods" ist nämlich numerisch und kein Vector...ein einfaches Überschreiben in "if (is.numeric(mods))" hilft an dieser Stelle aber leider nicht. Wenn ich das oben genannte Beispiel mit der Datei dat.hackshaw1998 ausführe, bekomme ich keine Fehlermeldung.
Ich hoffe sehr, dass mir jemand helfen kann!
Ganz liebe Grüße,
Laura