私の目的は、lapplyまたはそれに類似したmapplyを使って関数geom_point2への引数としてのパスリストです。同様の状況で、私はのようにgeom_cladelabelするリスト(またはリスト)を渡す成功を収めて:mapplyを使ってgeom_point2に引数を渡す
mapply(function (x,y,z,w,v,u,t,s) geom_cladelabel(node=x, label=y,
align=F, etc. # Where x y z etc are lists.
問題がgeom_point2内部aes
の使用に関連しています。 (geom_cladelabelにはない):
geom_point2の場合、ノード情報はaes
であり、できません。通常、エラーメッセージは表示されませんが、動作しません。
目的はこの例を動作させることですが、geom_point2
を2回書き込む代わりにmapplyを使用します。ここで
# source("https://bioconductor.org/biocLite.R")
# biocLite("ggtree")
library(ggtree)
library(ape)
#standard code
newstree<-rtree(10)
node1<-getMRCA(newstree,c(1,2))
node2<-getMRCA(newstree,c(3,4))
ggtree(newstree)+
geom_point2(aes(subset=(node ==node1)), fill="black", size=3, shape=23)+
geom_point2(aes(subset=(node ==node2)), fill="green", size=3, shape=23)
#desire to substitute the geom_point2 layers with mapply or lapply:
#lapply(c(node1,node2), function (x) geom_point2(aes(subset=(node=x)))))
サンプル入力データを使って[再現可能な例](https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)を提供すると、コードを実行してテストすることができます。 – MrFlick
'getMRCA'関数はどこから来ますか? – MrFlick
パッケージエイプ – Ferroao