Fabry

Author: GNtem2

Fabry/Fabry.Rmd

@@ -0,0 +1,53 @@
+---
+title: "Fabry"
+author: "GNtem2"
+date: "10/06/2021"
+output: html_document
+---
+
+```{r setup, include=FALSE}
+knitr::opts_chunk$set(echo = TRUE)
+```
+
+## R Markdown
+
+This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see <http://rmarkdown.rstudio.com>.
+
+When you click the **Knit** button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
+
+```{r Data}
+library(tidyverse)
+
+fabry<-readxl::read_xlsx("./FABRY_Enrolment_nov2018_update.xlsx")
+
+glimpse(fabry[,c(13:20,24:25)])
+
+```
+
+
+```{r Metaanalysis}
+library(metafor)
+df<-read.csv("FabryPrevalence.csv")
+
+#subset and remove NA
+df2<- df%>%
+  select(Author, Year, Country, Total, Fabry) %>%
+  mutate(xi=Fabry,
+         ni=Total,
+         pi=xi/ni) %>% na.omit()
+df2 <- escalc(measure="PFT", xi=xi, ni=ni, data=df2, add=0)	
+res <- rma(yi, vi, method="REML", data=df2, slab=paste(Author))
+
+
+#forest plot
+forest(res, transf=transf.ipft.hm, targs=list(ni=df2$ni), xlim=c(-1,1),refline=0.004,cex=.8, ilab=cbind(df2$Country, df2$xi, df2$ni),
+       ilab.xpos=c(-.6,-.4,-.2),digits=3)
+
+op <- par(cex=.75, font=2)
+#vertical position 11
+text(-1.0, 11, "Author ",     pos=4)
+text(c(-.4,-.2, -.60), 11, c("Fabry", 	"Total", "Country"))
+text(.95, 11, "frequency [95% CI]", pos=2)
+par(op)
+```
+Note that the `echo = FALSE` parameter was added to the code chunk to prevent printing of the R code that generated the plot.

Fabry/FabryPrevalence.csv

@@ -0,0 +1,11 @@
+Author,Year,Country,Total,Fabry,Male,FabryMale,Female,FabryFemale,,ICHincluded,TIAincluded,RestricyedInclusion,hypertension
+,,,,,,,,,,,,,
+rolfs,2005,Germany,721,28,432,21,289,7,,,,,
+brouns,2007,Belgium,103,0,64,0,39,0,,,,,
+baptista,2010,Portugal,493,6,300,3,193,3,,,,,
+brouns,2010,Belgium,1000,10,547,2,453,3,,,,,
+wozniak,2010,USA,154,1,154,1,,,,,,,
+lanthier,2017,Canada,365,1,,1,,0,,no,yes,,
+monash,2021,Australia,205,0,192,0,13,0,,yes,no,no,
+malavera,2020,Australia,58,0,35,0,23,1,,,,no,
+bersano,2016,Italy,33,1,,,,,,yes,yes,yes,