Displaying 1 result from an estimated 1 matches for "mssid".
Did you mean:
msgid
2007 Jan 21
2
efficient code. how to reduce running time?
...{
F=numeric(S) #create an empty vector to store the F-values
newY=sample(y,length(y)) #permute the cancer types
newX = cbind(x, newY);
# anova for all sites
for ( i in sites )
{
a <- anova(lm(newY~factor(newX[,i])));
F[i] <- a$`F value`[1];
}
MSSid <- which (F == max(F)); # index of MSS (Most Significant Site)
maxF = cbind(maxF,max(F));
}
maxF;
}
# set the output file
sink("/tmp/R.out.3932.100")
# load the dataset
snp = read.table(file("/tmp/msoutput.3932.100"))
#print (snp);
# pi: desired proportion of varia...