Forestplot : deux scripts similaires avec un qui tourne et l'autre pas

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:

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:

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 :calim2:

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