Re: [R] Speeding up lots of calls to GLM

2012-03-12 Thread Davy
y doesn't click with me. Also, I have never heard of a scores test before but some googling has turned up the Lagrange multiplier test. Is this the one you mentioned. Thanks again, Davy. On 12 March 2012 20:04, Thomas Lumley wrote: > On Tue, Mar 13, 2012 at 7:11 AM, Davy wrote: > >

[R] Speeding up lots of calls to GLM

2012-03-12 Thread Davy
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 col

[R] problem with indicators for switch

2010-08-02 Thread Roy Davy
Hi, I am having some problems setting up some indicators and would appreciate any help. I have some data called 'lights' with 3 variables called x, a and b. x - is the date a - equals 1 to indicate an 'on' button is activated b - equals 1 to indicate an 'off' button is activated Essentiall