computeStat.RdComputes the regression coefficients.
computeStat(ps, main_factor, time_var, subjectID_var, b, compStatParallel = FALSE)
| ps | phyloseq object. |
|---|---|
| main_factor | Character string. The name of the covariate variable. |
| time_var | Character string. The name of the time variable. |
| subjectID_var | Character string. The name of the subject ID variable. |
| b | A numeric. The block size to account for the dependence within-subject. |
| compStatParallel | A logical. True is used to use parallel in |
dataframe with the last column corresponds to taxa names and other columns are for estimated coefficients for each covariate using GEE, including Intercept.
Computes the library normalization factors using the median-ratio method, computes the observation-level weight using the mean-variance relationship (modification to voom), computes the residuals using glm glm.nb, computes the working correlation based on the given block size, computes the regression coefficient using the generalized estimating equation geem.