[BioC] Is it possible to fit a hierarchical linear model in LIMMA?

[pingzhao Hu]

Hi,
Say two groups A and B. Each group has 10 individuals.
Assume that we have two variables: Y and X.
Therefore, the hierarchical linear model will like:

Yij=B0i+B1Xij+ERRORij, (i=A,B; j=1,2,...,10)

We are interested in testing the significance of B1.
Although some R packages can do this job, we are interested in
using the shrinkage property in LIMMA or other R packages.

Thank you for your help in advance!

Pingzhao Hu
The Center for Applied Genomics
The Hospital for Sick Children
Toronto, Canada