Extact and calculate p-values of GLMM fit
# S3 method for cytoglmm summary(object, method = "BH", ...)
object | A |
---|---|
method | Multiple comparison adjustment method |
... | Other parameters |
tibble
data frame
set.seed(23) df = generate_data() protein_names = names(df)[3:12] df = dplyr::mutate_at(df, protein_names, function(x) asinh(x/5)) glmm_fit = CytoGLMM::cytoglmm(df, protein_names = protein_names, condition = "condition", group = "donor") #> Warning: moment-based covariance matrix estimate is not positive semi-definite; using projection summary(glmm_fit) #> # A tibble: 10 × 3 #> protein_name pvalues_unadj pvalues_adj #> <chr> <dbl> <dbl> #> 1 m03 0.0000682 0.000682 #> 2 m02 0.000163 0.000813 #> 3 m01 0.000681 0.00227 #> 4 m06 0.155 0.388 #> 5 m05 0.371 0.646 #> 6 m08 0.388 0.646 #> 7 m04 0.512 0.731 #> 8 m09 0.859 0.954 #> 9 m10 0.874 0.954 #> 10 m07 0.954 0.954