Displaying 1 result from an estimated 1 matches for "getinteractions2".
2012 Mar 12
1
Speeding up lots of calls to GLM
...to but I'd say it won't
be more than 8. In addition I have been advised that the best way to
account for multiple testing here is by permutation test, which in this
context becomes almost unfeasible, since I would have to run each
interaction ~10,000 times.
Thanks in advance, Cheers
Davy.
getInteractions2 = function(data, fSNPcol, ccCol)
{
#fSNPcol is the number of the column that contains the first SNP
#ccCol is the number of the column that contains the outcome variable
require(lmtest)
a = data.frame()
snps = names(data)[-1:-(fSNPcol-1)]
names(data)[ccCol] = "PHENOTYPE"
t...