# load library
# library(variancePartition)
# load simulated data:
data(varPartData)
form <- ~ Age + (1|Individual) + (1|Tissue)
res <- fitVarPartModel( geneExpr[1:10,], form, info )
# evaluate the collinearity score on the first model fit
# this reports the correlation matrix between coefficients estimates
# for fixed effects
# the collinearity score is the maximum absolute correlation value
# If the collinearity score > .999 then the variance partition
# estimates may be problematic
# In that case, a least one variable should be omitted
colinearityScore(res[[1]])
Run the code above in your browser using DataLab