Bonjour

Je réalise une méta-analyse via R et j'utilise le forestplot. J'ai deux jeux de données et j'applique exactement le même script pour chacun des jeux. Pour le jeu de données "white", le script tourne et crée le forestplot. Pour le jeu de données "MME", j'ai le message d'erreur suivant :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
#programme MME# 
#values
 mean_figure_M=c(gMME, ybarMME)
lower_figure_M=c(lower_CI95_gMME, lowerCI95_ybarMME)
upper_figure_M=c(upper_CI95_gMME, lowerCI95_ybarMME)
weight_figure_M=c(wMME, 0.2)
 
 #corrections 
 study_namesMME=c("Summary","Agban","Kenealy", "Lee", "Levya", Plockinger", "Weinstock", "West"," ")
 blanksMME=rep(" ", length(study_namesMME))
 textMME=cbind(study_namesMME, blanksMME)
 
 #Forest 
forestplot (labeltext=textMME, mean=mean_figure_M,lower=lower_figure_M,upper=upper_figure_M,align=rep(1,5),is.summary=subtotal,
+ zero=0,xlog=F,xticks=c(-0.2, 0, 0.2, 0.4, 0.6, 0.8,1),col=meta.colors(lines="black", box="gray", zero="gray", summary="blue"))[/COLOR]
[COLOR="Blue"]Erreur dans if ((convertX(unit(OR, "native") + unit(0.5 * size, "lines"),  : 
  valeur manquante là où TRUE / FALSE est requis
L'autre programme jumeau, qui lui tourne sans soucis, est le suivant :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
> #values
 mean_figure=c(gWhite,ybarWhite)
 lower_figure=c(lower_CI95_gWhite, lowerCI95_ybarWhite) 
 upper_figure=c(upper_CI95_gWhite, upperCI95_ybarWhite)
 weight_figure=c(wWhite, 0.2)
 
 #corrections 
study_names=c("Summary","Agban","Kenealy", "Lee", "Levya", "Plockinger", "Weinstock", "West"," ")
blanks=rep(" ", length(study_names))
text=cbind(study_names, blanks)
 
#Forest
forestplot (labeltext=text, mean=mean_figure, lower=lower_figure, upper=upper_figure, align=rep(1,5), is.summary=subtotal, zero=0, xlog=F, xticks=c(-0.2, 0, 0.2, 0.4, 0.6, 0.8, 1),  col=meta.colors(lines="black", box="gray", zero="gray", summary="blue"))
Peut-être que l'erreur est évident mais là, je cale sérieusement

Si l'une ou l'un d'entre vous a une super idée, je le remercie d'avance !!!

Merciiiiii

PS: je ne programme en R que depuis le 10/09/12