[GenABEL-dev] fixing memory leak in MixABEL::FastMixedModel -- remaining issues

[Yurii Aulchenko]

William,

I have just committed the code which you have contributed in order to
fix memory leak in MixABEL::FastMixedModel (which you probably have
seen on the list already)

The package does compile, but there is one problem -- the test
statistics is uniformly zero

see example:

[code]
library(MixABEL)
example(FastMixedModel)
# betas are quite the same from GWFGLS and FMM
cor(mms[,"effB"],fmm$coef)^2
plot(mms[,"effB"],fmm$coef)
# but chi2's are clearly not
cor(mms[,"chi2.1df"],fmm$chi.sq)^2
# cause FMM produce 0's
summary(fmm$chi.sq)
[/code]

William, can you please have a look at this? I will send instructions
on how to get dev-version working in next e-mail.

best regards,
Yurii