An R package for analysis of microbiome relative abundance data using zero inflated beta GAMLSS and meta-analysis across microbiome studies using random effects models
Hi,
I am using the CRAN package to analyze (meta-analysis) of an example data. I made a data frame as shown in taxtab6.rm[[5]]. All the function ran well, as per example in github. However, the function, metatab.show gives a "table with no data".
I am wondering what is the problem here. There is no error message.
Thanks a lot in advance for the reply.
I am attaching the codes below.
Hi, I am using the CRAN package to analyze (meta-analysis) of an example data. I made a data frame as shown in taxtab6.rm[[5]]. All the function ran well, as per example in github. However, the function, metatab.show gives a "table with no data". I am wondering what is the problem here. There is no error message.
Thanks a lot in advance for the reply. I am attaching the codes below.
input_1 = read.csv("Input_1.csv", header = TRUE, sep = ",") input_2 = read.csv("Input_2.csv", header = TRUE, sep = ",") input_3 = read.csv("Input_3.csv", header = TRUE, sep = ",")
taxacom.ex_1<-taxa.compare(taxtab=input_1,propmed.rel="gamlss",comvar="gender",adjustvar=c("age","smoking"),longitudinal="no") taxacom.ex_2<-taxa.compare(taxtab=input_2,propmed.rel="gamlss",comvar="gender",adjustvar=c("age","smoking"),longitudinal="no") taxacom.ex_3<-taxa.compare(taxtab=input_3,propmed.rel="gamlss",comvar="gender",adjustvar=c("age","gender"),longitudinal="no")
taxacom.ex_1$study <- "India" taxacom.ex_2$study <- "Bangladesh" taxacom.ex_3$study <- "Germany"
tabsex4<-rbind.fill(taxacom.ex_1, taxacom.ex_2, taxacom.ex_3) metab.sex<-meta.taxa(taxcomdat=tabsex4, summary.measure="RR", pool.var="id", studylab="study", backtransform=FALSE, percent.meta=0.5, p.adjust.method="fdr")
Metadat<-metatab.show(metatab=metab.sex$random,com.pooled.tab=tabsex4,tax.lev="l2",showvar="genderMale",p.cutoff.type="p", p.cutoff=1,display="data")