search for: snp_a1

Displaying 2 results from an estimated 2 matches for "snp_a1".

Did you mean: snp_a
2007 Nov 16
1
constraint matrices in vglm (VGAM package)
Hello R users, I am performing a multinomial logit regression and would like to constrain a few model coefficients to be equal. Here is my model: multi <- vglm(case123con ~ SNP_A1+SNP_A2+age, multinomial, work.analy) where case123con is a four level categorical variable (case 1, case 2, case 3, control) and SNP_A1 and SNP_A2 are indicator functions (yes/no). The output of this model is: Coefficients: Value Std. Error t value (Intercept):1...
2009 Oct 12
1
Speed up and limit memory usage of lm()
...cution time 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,m...