[R-es] forest plot

[Jose Betancourt Bethencourt]

Gracias por su vliosa ayuda
saludos
José

El 20/6/22, Marcelino de la Cruz Rot <marcelino.delacruz using urjc.es> escribió:
> 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
>
> _______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>


-- 
Dr. Jose A. Betancourt Bethencourt
Universidad de Ciencias Medicas Carlos j. Finlay