2012-10-26 45 views
-2

我的代碼顯示如下:如何繪製與非分類x軸(格)的柱狀圖

read.table("GO.txt",header=T,sep=',')->go 
library(lattice) 
barchart(go[,1]~go[,2],horiz=F,ylim=c(30,29666), 
layout=c(1,1),stack=F, 
auto.key=list(space='right'), 
ylab="Yes", 
scales=list(x=list(rot=45),y=list(log = T))) 

不過,我不希望繪製與排序X的圖表,又該我做? 格子排序我的數據,我不需要。

數據:

number_of_A-Unigene Class Ontology 
    45 biological adhesion biological_process 
    5328 biological regulation biological_process 
    8 carbon utilization biological_process 
    6 cell killing biological_process 
    119 cell proliferation biological_process 
    3248 cellular component organization or biogenesis biological_process 
    16594 cellular process biological_process 
    244 death biological_process 
    4354 developmental process biological_process 
    3547 establishment of localization biological_process 
    917 growth biological_process 
    429 immune system process biological_process 
    3716 localization biological_process 
    30 locomotion biological_process 
    15726 metabolic process biological_process 
    4182 multicellular organismal process biological_process 
    1673 multi-organism process biological_process 
    877 negative regulation of biological process biological_process 
    10 nitrogen utilization biological_process 
    11 pigmentation biological_process 
    640 positive regulation of biological process biological_process 
    4808 regulation of biological process biological_process 
    2462 reproduction biological_process 
    2437 reproductive process biological_process 
    7812 response to stimulus biological_process 
    190 rhythmic process biological_process 
    2001 signaling biological_process 
    2 sulfur utilization biological_process 
    44 viral reproduction biological_process 
    19852 cell cellular_component 
    1763 cell junction cellular_component 
    19852 cell part cellular_component 
    31 extracellular matrix cellular_component 
    9 extracellular matrix part cellular_component 
    686 extracellular region cellular_component 
    30 extracellular region part cellular_component 
    3698 macromolecular complex cellular_component 
    9829 membrane cellular_component 
    3432 membrane part cellular_component 
    1439 membrane-enclosed lumen cellular_component 
    15252 organelle cellular_component 
    6024 organelle part cellular_component 
    1753 symplast cellular_component 
    216 antioxidant activity molecular_function 
    15917 binding molecular_function 
    15103 catalytic activity molecular_function 
    2 channel regulator activity molecular_function 
    433 electron carrier activity molecular_function 
    319 enzyme regulator activity molecular_function 
    8 metallochaperone activity molecular_function 
    522 molecular transducer activity molecular_function 
    708 nucleic acid binding transcription factor activity molecular_function 
    26 nutrient reservoir activity molecular_function 
    102 protein binding transcription factor activity molecular_function 
    6 protein tag molecular_function 
    613 receptor activity molecular_function 
    1039 structural molecule activity molecular_function 

5翻譯調節活性molecular_function 2478轉運活性molecular_function

+0

請張貼'dput(去)',而不是剪切和粘貼您的打印數據。 – BenBarnes

回答

0

「格子」 把你的因子變量和默認的(字母)順序排列它們。您可以更改訂單。也許:

go[,2] <- factor(go[,2] unique(as.character(go[,2]))) 

這應該設置在同一順序的水平,他們將第一次出現,如果你是通過翻閱。另一個等值,雖然更巴洛克,轉讓將是:

go[,2] <- factor(go[,2] as.character(go[,2][!duplicated(go[,2])])) 

然後做你的情節。 (我沒有試圖輸入你的數據,他們看起來太雜亂了,實際上,我確實給了它一個選項卡作爲分隔符,如果他們曾經有過標籤,那麼他們會丟失它們,並且元素的數量是可變的如果白色空間是分隔符)

在重複的問題,使用數據標籤進行排序:

go[[2]] <- factor(go[[2]], levels=unique(as.character(go[[2]]))) 

+0

3Q,但你的建議是沒用的......它仍然被排序。 – Sandy

+0

恐怕什麼是'無用'就是你的榜樣。我解釋了爲什麼它不能用read.table(...,header = TRUE)加載。如果您提供輸入代碼,則可能有進展的可能性。 –

+0

好吧,非常感謝,我的數據代碼在這裏發佈。 http://goo.gl/2laUX – Sandy