require(graphics)
require(utils)
require(ggplot2)
require(NeatMap)
dHN <- read.table('http://dl.dropbox.com/u/33619290/DadosHugoNeto.txt',h=T)
dendHN <-
hclust(dist(sqrt(dHN)),method='complete')
clones=c("RB92579","RB768647","Q124","RB75126","RB925211","RB835486","RB946022",
"RB912825","SP80-3280","RB855511","RB867515","SP83-5073","RB72454","SP80-1816",
"RB825336","RB855156","RB855536","SP83-2847","RB855453","TUC77-42","SP70-2233",
"RB931555","RB835054","RB855206","RB925345","RB8317","RB947501","RB93509",
"RB813804","RB855113","RB83102","SP77-5181","RB845257","RB855046","RB845197",
"RB835089","RB8495","RB825548","RB855463","RB925453","RB855595","RB945951",
"RB946903","RB855035","RB92606","RB735200","RB945961","SP70-1143","RB855036",
"IAC87-3396","CO62175","RB956911","TUC71-7","RB9350","RB966928","RB912850",
"RB845210","L60-14","SP81-3250","RB928064","RB721012","SP79-1011","SP801842",
"IAC52-326","RB935915","RB806043")
ggplot.default() +
opts(panel.background = theme_blank(),
panel.grid.major = theme_blank(),
panel.grid.minor = theme_blank(),
legend.position ='none',
axis.title.x = theme_blank(),
axis.title.y = theme_blank(),
axis.text.y = theme_text(angle=0,size=8))+
scale_x_continuous(breaks=seq(0:11),labels=seq(0,0.5,by=0.045))+
scale_y_continuous(breaks=seq(1:66),labels=clones) +
geom_vline(xintercept=2.2,colour='red',linetype=2)+
geom_text(aes(x=2.6,y=-1,label='0.06',colour='red',size=5)) +
geom_text(aes(x=3.1,y=67,label='Magic method',colour='red',size=5))+
geom_text(aes(x=5.3,y=30,label='88',colour='green',size=5))+
draw.dendrogram(dendHN, dendro.dir="right")
Allaman
(S,f,P)