Extact and calculate p-values of bootstrap GLM fit
# S3 method for cytoglm 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)) 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