Extact and calculate p-values of GLMM fit

# S3 method for cytoglmm
summary(object, method = "BH", ...)

Arguments

object

A cytoglmm class

method

Multiple comparison adjustment method

...

Other parameters

Value

tibble data frame

Examples

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