我有一个矩阵Expr
,其中行代表变量和列样本。我有一个分类向量groups
(包含“A”、“B”或“C”)我想测试哪些变量 'Expr' 可以通过样本属于 a 的事实来解释group
。
我的策略是使用广义加法模型(具有负二项分布)对问题进行建模。然后我想以变量明智的方式使用似然比测试来获得每个变量的 ap 值。我愿意:
require(VGAM)
m <- vgam(Expr ~ group, family=negbinomial)
m_alternative <- vgam(Expr ~ 1, family=negbinomial)
接着:
lr <- lrtest(m, m_alternative)
最后一步是错误的,因为它测试的是两个模型的整体似然比,而不是变量。我想获得每个变量的 p 值向量,而不是单个 p 值。
我该怎么做?(我对 R 很陌生,所以请原谅我的愚蠢)