Extact and calculate p-values of bootstrap GLM fit

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

Arguments

object

A cytoglm 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))
glm_fit <- CytoGLMM::cytoglm(df,
                             protein_names = protein_names,
                             condition = "condition",
                             group = "donor",
                             num_boot = 10) # in practice >=1000
summary(glm_fit)
#> # A tibble: 10 × 3
#>    protein_name pvalues_unadj pvalues_adj
#>    <chr>                <dbl>       <dbl>
#>  1 m01                    0.2         0.5
#>  2 m02                    0.2         0.5
#>  3 m03                    0.2         0.5
#>  4 m04                    0.2         0.5
#>  5 m05                    0.4         0.5
#>  6 m06                    0.4         0.5
#>  7 m08                    0.4         0.5
#>  8 m09                    0.4         0.5
#>  9 m07                    0.6         0.6
#> 10 m10                    0.6         0.6