5.1 Linear Regression
limma_a_b is used to perform linear regression.
lm_res <- limma_a_b(eset = m, model.str = "~ AGE", coef.str = "AGE")
head(arrange(lm_res, adj.P.Val)) # top 6 rows sorted by adjusted p-value## logFC AveExpr t P.Value adj.P.Val
## NP_001077077.1 -0.10468905 -4.163336e-17 -4.186740 0.0001865011 0.5982300
## NP_001138679.1 0.01185948 7.806256e-18 3.821284 0.0002953129 0.5982300
## NP_003286.1 0.01335415 -6.938894e-18 3.870673 0.0002292102 0.5982300
## NP_055336.1 -0.01819954 1.233581e-18 -3.986807 0.0002306581 0.5982300
## NP_000026.2 -0.01689006 -5.703200e-18 -2.847600 0.0056747708 0.8201268
## NP_000176.2 -0.02059760 -3.802134e-18 -3.055901 0.0031033317 0.8201268
## B
## NP_001077077.1 0.001863506
## NP_001138679.1 -1.049384099
## NP_003286.1 -0.880645496
## NP_055336.1 -0.598496581
## NP_000026.2 -3.895795433
## NP_000176.2 -3.338485703Below is a graphical representation of the results for a specific feature.
