Displaying 1 result from an estimated 1 matches for "upregulated_list".
2012 Jun 12
1
Analyzing large files faster
...;2.172342" "Trp53inp1"
"1422557_s_at" "1.58e-09" "2.9569254" "Mt1"
etc.
using the following code:
muscle = read.table(file="/Users/bob/Desktop/Muscle/musclesmall.txt", header
= TRUE, colClasses = "character", fill = TRUE)
upregulated_list = c()
downregulated_list = c()
nochange = c()
p_thresh = 6.51e-06
x=1
while (x <= nrow(muscle)) {
this_pval = muscle[x,"adj.P.Val"]
this_M = muscle[x, "logFC"]
if (muscle[x, "Gene.symbol"] == "") {
x= x +1
}
else {if ((this_M >= 1.0) & (this_p...