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...