Displaying 1 result from an estimated 1 matches for "isnp".
Did you mean:
isn
2005 Mar 04
0
Is aggregate() what I need here?
...snps<-t(snps)
texp<-t(exp)
sink(paste("output.", sep="", Sys.getenv("LSB_JOBINDEX")))
#loops below are hardwired for 5 gene-expression levels (some genes
have two
#probes, and those are treated as separate genes for now) and 100 SNPs.
for (iexp in 1:5){
for (isnp in 1:100){
genotype<-factor(tsnps[,isnp])
# make sure there is more than one genotypic class before doing
ANOVA
if (length(unique(genotype))>1) {
expression<-texp[,iexp]
stuff<-anova(lm(expression ~ genotype))
qq=c(iexp,isnp)
print (qq)
p...