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 :
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
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
Peut-être que l'erreur est évident mais là, je cale sérieusement
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"))
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
Partager