行列plot用変換関数

こんなん探せばpackageにあるだろたぶん…

plot.converter<-                                # 変換しますよ〜
function(A){                                    # 変数は行列
  x<-  t(t(c(rep(1:nrow(A),ncol(A)))))          # xを用意
  y<- c()                                       # yからの…
  for(i in 1:(T+1)){
    y<- c(y,rep(i,N+1))
  }
  y<- t(t(y))                                   # yを用意
  convert.matrix<- cbind(x,y,c(A)) 
  return(convert.matrix)                        # 3列の行列完成。1列目からx,y,zに対応
}

通りすがりさんが使える関数を伝授してくださった。

expand.grid(1:3,1:4)
   Var1 Var2
1     1    1
2     2    1
3     3    1
4     1    2
5     2    2
6     3    2
7     1    3
8     2    3
9     3    3
10    1    4
11    2    4
12    3    4

というわけで、これを使うと

matrix.converter<- 
function(A){
  return(cbind(expand.grid(1:nrow(A),1:ncol(A)),c(A)))
}

先生のを使うと

ryamada.converter<- 
function(A){
  return(cbind(which(A<=max(A),arr.ind=TRUE),c(A)))
}