1,120
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search→Implementing Additional Statistical Tests
## MISSING VALUES : IGNORED
## MISSING VALUES : IGNORED
single.q.lm <- function() {
single.q.lm <- function() {
cname <- c("BETA","SEBETA","TSTAT")
cname <- c("BETA","SEBETA","TSTAT") # column names for additional variables in the EPACTS output
m <- nrow(genos)
m <- nrow(genos)
p <- rep(NA,m)
p <- rep(NA,m)
if ( m > 0 ) {
if ( m > 0 ) {
for(i in 1:m) {
for(i in 1:m) {
r <- summary(lm(pheno~genos[i,]+cov-1))$coefficients[1,]
r <- summary(lm(pheno~genos[i,]+cov-1))$coefficients[1,] # run simple linear regression
p[i] <- r[4]
p[i] <- r[4] # store p-value to p[i]
add[i,] <- r[1:3]
add[i,] <- r[1:3] # store additional variables to add[i,]
}
}
}
}
