Displaying 1 result from an estimated 1 matches for "sys_pheno".
Did you mean:
sp_pheno
2009 Oct 12
1
Speed up and limit memory usage of lm()
...and memory usage becomes a huge issue. I have therefore been trying to speed the process and limit the memory usage.
Here follows part of the code do give better understanding of what I am doing:
modela <- lm(RSSYS10 ~ RS_AGE + RS_AGESQ + SEX + RS_BMI)
frssys <- function(SNP_A1,data=sys_pheno, model=modela){
modelb <- lm(RSSYS10 ~ RS_AGE + RS_AGESQ + SEX + RS_BMI + SNP_A1,data=data)
modelc <- lm(RSSYS10 ~ RS_AGESQ + SEX + RS_BMI + SNP_A1 * RS_AGE,data=data)
p1 <- anova(model,modelb)$P[2]
p2 <- anova(model,modelc)$P[2]
p3 <- anova(modelb,modelc)$P[2]...