Rに3Dの矢印を描くにはどうすればよいですか？

Question

私はrglパッケージを使用して自分のデータの3Dプロットを作成しています。いくつかの理由（3D PCAバイプロット）のために、私はベクトルが必要です - 矢印のある線分。私は矢のように3Dの円錐を持っているので、私は立ち往生しています。

どうにかして、問題のジオメトリの周りに私の老人の心を包み込むことはできません。つまり、私はベクトルをユーザ座標系の原点から[3,3,3]へのベクトルであるベクトル

segments3d(rbind(c(0, 0, 0), c(3, 3, 3))) 

で描画します。

[3,3,3]にヒントを付けてコーンを作成したいと思います。コーンの基部は円で形成することができる。半径Rの（Y平面に垂直）XZ平面上で円を描くことは簡単である：

n <- 10 
sin.t <- sin(seq(0, 2 * pi, len= n)) 
cos.t <- cos(seq(0, 2 * pi, len= n)) 
r <- 0.1 
xv <- x + r * sin.t 
yv <- rep(y, n) 
zv <- z + r * cos.t 

しかし、どのように私は今、円が垂直ベクトルに今あるように、これらの点を変換しますか？ベクトルの方向に沿った先端から0.2の中心？私がこの変換をしたら、triangles3d関数で三角形を描画します。各三角形の先端には1つの角があり、円の点には2つの頂点があります。

これは基本的な数学であり、18歳の私には問題がないことが分かっています（28歳の私も）。どんなフック（魚とは対照的に）でも感謝します。

Answer 1

rglのデモでは、cone3d関数があります。ベースとチップを別々に取ります。いずれにせよ、あなたはこのような何か行うことができます：ここで

vec=rbind(c(0, 0, 0), c(3, 3, 3)) 
segments3d(vec) 


cone3d(base=vec[2,]-(vec[1,]+vec[2,]/6), 
    #this makes the head go 1/6th the length of the arrow 
     rad=0.5, 
     tip=vec[2,], 
     col="blue", 
     front="lines", 
     back="lines") 

をcone3d機能である：

cone3d <- function(base=c(0,0,0),tip=c(0,0,1),rad=1,n=30,draw.base=TRUE,qmesh=FALSE, 
        trans = par3d("userMatrix"), ...) { 
    ax <- tip-base 
    if (missing(trans) && !rgl.cur()) trans <- diag(4) 
    ### is there a better way? 
    if (ax[1]!=0) { 
    p1 <- c(-ax[2]/ax[1],1,0) 
    p1 <- p1/sqrt(sum(p1^2)) 
    if (p1[1]!=0) { 
     p2 <- c(-p1[2]/p1[1],1,0) 
     p2[3] <- -sum(p2*ax) 
     p2 <- p2/sqrt(sum(p2^2)) 
    } else { 
     p2 <- c(0,0,1) 
    } 
    } else if (ax[2]!=0) { 
    p1 <- c(0,-ax[3]/ax[2],1) 
    p1 <- p1/sqrt(sum(p1^2)) 
    if (p1[1]!=0) { 
     p2 <- c(0,-p1[3]/p1[2],1) 
     p2[3] <- -sum(p2*ax) 
     p2 <- p2/sqrt(sum(p2^2)) 
    } else { 
     p2 <- c(1,0,0) 
    } 
    } else { 
    p1 <- c(0,1,0); p2 <- c(1,0,0) 
    } 
    degvec <- seq(0,2*pi,length=n+1)[-1] 
    ecoord2 <- function(theta) { 
    base+rad*(cos(theta)*p1+sin(theta)*p2) 
    } 
    i <- rbind(1:n,c(2:n,1),rep(n+1,n)) 
    v <- cbind(sapply(degvec,ecoord2),tip) 
    if (qmesh) 
    ## minor kluge for quads -- draw tip twice 
    i <- rbind(i,rep(n+1,n)) 
    if (draw.base) { 
    v <- cbind(v,base) 
    i.x <- rbind(c(2:n,1),1:n,rep(n+2,n)) 
    if (qmesh) ## add base twice 
     i.x <- rbind(i.x,rep(n+2,n)) 
    i <- cbind(i,i.x) 
    } 
    if (qmesh) v <- rbind(v,rep(1,ncol(v))) ## homogeneous 
    if (!qmesh) 
    triangles3d(v[1,i],v[2,i],v[3,i],...) 
    else 
    return(rotate3d(qmesh3d(v,i,material=...), matrix=trans)) 
}