data("POSTmData")
y <- as.integer(x = metadata[,"GC"] == "BV")
X <- metadata[,"mRace"]
result <- post(y = y,
X = X,
OTU = otu[,1:20],
tree = otutree,
cValues = seq(0,0.05,by=0.01))
print(x = result)
Run the code above in your browser using DataLab