我在着色與「猿」(慢邊着色)如何避免循環所有tip.labels(包APE)
自從我用C編程始終,我還是有些功能產生phylo樹邊發現很難以循環式的方式停止思考。
我能想到做到這一點的唯一方法是(1)循環所有tip.labels(ID),(2)找出哪些邊緣屬於它們,(3)設置所需的顏色。
這樣做是1加1,因此是非常緩慢的大樹:
tsampltime.rooted=structure(list(edge = structure(c(24L, 24L, 24L, 24L, 24L, 25L,
26L, 26L, 27L, 27L, 28L, 28L, 25L, 29L, 29L, 30L, 30L, 30L, 30L,
24L, 31L, 31L, 32L, 32L, 32L, 33L, 33L, 34L, 35L, 35L, 34L, 36L,
36L, 34L, 37L, 37L, 1L, 2L, 12L, 23L, 25L, 26L, 6L, 27L, 5L,
28L, 3L, 4L, 29L, 7L, 30L, 8L, 9L, 10L, 11L, 31L, 13L, 32L, 21L,
22L, 33L, 20L, 34L, 35L, 14L, 15L, 36L, 16L, 17L, 37L, 18L, 19L
), .Dim = c(36L, 2L)), Nnode = 14L, tip.label = c("0", "2325",
"55304", "124953", "72254", "66507", "85089", "110256", "123265",
"97350", "123721", "36770", "48692", "110612", "97224", "104337",
"124625", "128499", "120928", "88404", "73335", "75059", "17928"
), edge.length = c(0, 0.953297, 8.054944, 4.4120893, 9.173083,
1.409346, 3.752752, 0.483517, 4.620875, 0.582417, 0.510989, 12.4862723,
6.291209, 1.920329, 3.071429, 4.5027528, 5.497248, 2.777472,
5.5274749, 8.414843, 2.5467017, 3.79121, 3.824171, 3.961538,
3.804944, 2.126375, 1.75275, 1.93956, 3.3516546, 1.57418, 2.31319,
2.22528, 4.0384651, 3.898348, 2.722523, 1.87088)), .Names = c("edge",
"Nnode", "tip.label", "edge.length"), class = "phylo", order = "cladewise")
...
#distValuesPerId[,] has [LABELID,COLOR]
distValuesPerId=source('http://ubuntuone.com/5y7ZYCWfE73T5lhnUpmeXc')
...
uniqueIDs=unique(tree$tip.label)
distTrdsampledcol <-rep("black", length(tree$edge)) #init in black
for(i in uniqueIDs) { #(1)
a= c(which(tree$tip.label==i))
b= which(tree$edge[,2]== a) #(2)
distTrdsampledcol [ b ] <- distValuesPerId[i,2] #(3)
}
...
#plot(tree, edge.color=distTrdsampledcol)
誰能幫助我重新考慮這一點?這樣做更有效率嗎?
提前感謝!
j
使用'dput(樹)'一棵小樹上給我們一些樣本數據。 – nograpes 2012-07-06 19:32:53
@nograpes,在這裏它是: [dput(distValuesPerId)](http://ubuntuone.com/5y7ZYCWfE73T5lhnUpmeXc) [dput(樹)](http://ubuntuone.com/3GdyoX4ETlU5cDuRpPMfcN) [圖(樹)(http://ubuntuone.com/5ikA8OIrN9nQrpIgpzTehg) 感謝您的答覆 – lourencoj 2012-07-09 13:45:54
實際上,你可以只發布圖片你的問題裏面,你也應該複製的代碼來創建你的樹有好。我現在就爲你做。它不是張貼'data.frame' 12萬行是個好主意。你應該採取只需要指出data.frame的顏色,並張貼。 – nograpes 2012-07-09 14:34:31