[R-es] forest plot

Marcelino de la Cruz Rot m@rce||no@de|@cruz @end|ng |rom urjc@e@
Lun Jun 20 09:54:07 CEST 2022


Hola Jose:

Si miras la ayuda de metacont() verás que antes del argumento "studlab" 
requiere otros seis argumentos adicionales 
(n.e,mean.e,sd.e,n.c,mean.c,sd.c).

Tú le estás suministrando sólo 4 y de ahí que se queje de que no le 
suministras el 5º, "mean.c".


Un saludo,
Marcelino

El 20/06/2022 a las 0:43, Jose Betancourt Bethencourt escribió:
> Estimados
> Intentaba hacer el gráfico de forest plot imitando los códigos de un
> libro, lo del libro sale bien (parte de arriba, lo mio (parte de abajo
> no sale y no veo porque, apreciaria su ayuda.
> Adjunto ambas bases de datos
>
>
>
>
>
> LA DEL LIBRO QUE FUNCIONA
>
> rm(list = ls())
>
> library(readr)
> data1 <- read_csv("000 EN DISEñO/00 Meta/LIBRO/dataset01.csv")
> View(data1)
>
>
> wd=("D:/Public/Documents/000 EN DISEñO/00 Meta/"LIBRO")
> # 1. Read in the data
> data1 <- read.csv("dataset01.csv", as.is=TRUE)
>
>   2. Print structure of R object data1
> str(data1)
> # 3. To view an R object, just type its name:
> data1
>
>   # 1. Load add-on package meta
>     library(meta)
>
>   # 2. Do meta-analysis
>     m <- metacont(Ne, Me, Se, Nc, Mc, Sc,
>                     studlab=paste(author, year),
>                     data=data1)
>
>       # 3. Produce forest plot
>     forest(m, xlab="Maximum % fall in FEV1")
>
>
> #=====================================================================
> HAGO ESTO CON UNA BASE DE DATOS MIA EN EXCEL Y NO ME FUNCIONA
>
> rm(list = ls())
>
> library(readxl)
> data1 <- read_excel("000 EN DISEñO/00 BIA/guia1.xlsx")
> View(data1)
>
> 2. Print structure of R object data1
> str(data1)
> # 3. To view an R object, just type its name:
> data1
> attach(data1)
> # 1. Load add-on package meta
> library(meta)
>
>
> # 2. Do meta-analysis
> m <- metacont(Fat,SD1, Sist, SD2,
>                studlab=paste(author, year),
>                data=data1)
> ##AQUI ME DICE
> Error: Provide either argument 'mean.c' or 'median.c'
>
>
> # 3. Produce forest plot
> forest(m, xlab="Maximum % fall in FEV1")
>
> _______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


-- 
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España



Más información sobre la lista de distribución R-help-es