R_Tips0 のバックアップの現在との差分(No.11)

追加された行はこの色です。
削除された行はこの色です。
[[Rの備忘録]]

R で2種類のグラフィックスを作成するとする．片方にはカテゴリA, C, D があり，片方にはA, B, C  があるとする．ここで両方に共通するカテゴリを同じ色と指定したい．もっとも簡単な方法は?もちろん，手作業でアドホックに処理するのはなしで．


#ref(colorPlot.png)


思い出さないなぁ．
もっと簡単な方法があったはずだが，とりあえず，こんな方法でも実現できる．

まず二つの組に登場するカテゴリを一緒くたにしてベクトル z として扱い，それぞれの組のカテゴリ x, y が，この共通ベクトル z の何番目にあるかを調べるという方策である．
 x <- c("A","C", "D")
 y <- c("A","B","C")
 z <- unique(c(x,y))

#y の各要素がzの何番目にあるか知りたい

 y %in% z #意味無し

# 添字番号を取るためにループする?こんな厄介な処理が必要なのかなぁ?
 xcol <- NULL
 ycol <-  NULL
 
 for( i in 1:length(x)){
  for(j in 1:length(z)){
    if(x[i] == z[j]){
      cat(j, "\n")
      xcol <- c(xcol, j)
    }
  }
 }
 for( i in 1:length(y)){
  for(j in 1:length(z)){
    if(y[i] == z[j]){
      cat(j, "\n")      
      ycol <- c(ycol, j)
    }
  }
 }

 xx <- sample(x, 100, rep = T)
 xx.n <- sample(1:5, 100, rep = T)
 xx.l <- sample(x,100, rep = T)
 
 yy <- sample(y,100, rep = T)
 yy.n <- sample(1:5, 100, rep = T)
 yy.l <- sample(y,100, rep = T)

 par(mfrow = c(1,2))
 matplot(prop.table(table(xx.n, xx.l)), type = "l",
          col = xcol, lwd = 3)
 legend(1,.1, legend = x, col =  xcol,  lwd = 3)
 
 matplot(prop.table(table(yy.n, yy.l)), type = "l",
           col = ycol, lwd = 3)
 legend(1,.1, legend = y, col =  ycol,  lwd = 3)
# いや，絶対もっと簡単な方法があったはず

#article
**QGNOaMIppzo [#x84d87a0]
>[[askazupw]] (2010-03-16 (火) 17:50:04)~
~
KXmJJs  <a href="http://cipxoliwsiob.com/">cipxoliwsiob</a>, [url=http://mzzpghriorje.com/]mzzpghriorje[/url], [link=http://xjvavnocrjrc.com/]xjvavnocrjrc[/link], http://vatkcfdgkavj.com/~

//

#comment

**refNOuarGJEvh [#k3238b54]
>[[xjedfpm]] (2010-03-16 (火) 07:17:00)~
~
bfUgpi  <a href="http://jeyijfblxmlp.com/">jeyijfblxmlp</a>, [url=http://xbfjnrjanjuy.com/]xbfjnrjanjuy[/url], [link=http://vixcegaomvxr.com/]vixcegaomvxr[/link], http://akrttxrmozsy.com/~

//

#comment

**obUjrEGKnqjKtl [#o3921d98]
>[[nprqdaz]] (2010-03-15 (月) 21:58:14)~
~
dGVOcC  <a href="http://cinfcnjzwvbv.com/">cinfcnjzwvbv</a>, [url=http://mqteiivuonpp.com/]mqteiivuonpp[/url], [link=http://zzfwyiwelrct.com/]zzfwyiwelrct[/link], http://fzeezzpptxuo.com/~

//

#comment

**gPihbjxQMpGJ [#w6ffeaed]
>[[akbwlnojt]] (2010-03-15 (月) 14:44:47)~
~
KVRr1f  <a href="http://talfxsaxivmb.com/">talfxsaxivmb</a>, [url=http://bshbwijrlgvt.com/]bshbwijrlgvt[/url], [link=http://gjivjituejzc.com/]gjivjituejzc[/link], http://vmrtvbnvknmp.com/~

//

#comment

**xLuhfRRBpTHdiw [#r4b32939]
>[[tgyeshlvbj]] (2010-03-15 (月) 08:12:59)~
~
ZzaD7F  <a href="http://tydklqezvvcm.com/">tydklqezvvcm</a>, [url=http://uagjffbjplrh.com/]uagjffbjplrh[/url], [link=http://uuxryokkqyye.com/]uuxryokkqyye[/link], http://gkzpokodupzu.com/~

//

#comment

**xZeGRydPkntgIVjA [#cec0790a]
>[[sztstyyl]] (2010-03-15 (月) 06:03:45)~
~
HUeAc4  <a href="http://ahlxgvmezevr.com/">ahlxgvmezevr</a>, [url=http://vfrrsxreqzql.com/]vfrrsxreqzql[/url], [link=http://dutcvyiijblb.com/]dutcvyiijblb[/link], http://eexrfqzfylec.com/~

//

#comment

**iXKuVMATuTy [#q54348eb]
>[[ixpdgh]] (2010-03-11 (木) 05:35:42)~
~
KPVjm7  <a href="http://wrjotfpavtjk.com/">wrjotfpavtjk</a>, [url=http://upjofzuucuoh.com/]upjofzuucuoh[/url], [link=http://szduypfeokkj.com/]szduypfeokkj[/link], http://chjtcswcciaq.com/~

//

#comment

**KffyJSMCOWLk [#kb90eece]
>[[cvxzjkz]] (2010-03-09 (火) 09:31:31)~
~
ea0uDo  <a href="http://wvbfvsockewc.com/">wvbfvsockewc</a>, [url=http://szyhpyncwmsl.com/]szyhpyncwmsl[/url], [link=http://tfgrernjxars.com/]tfgrernjxars[/link], http://becrpjojakjr.com/~

//

#comment

**PLYFlHrEJKrJ [#wf68aac3]
>[[rzaleavbcj]] (2010-02-26 (金) 01:30:00)~
~
HqlNeZ  <a href="http://yomwpsstsuxl.com/">yomwpsstsuxl</a>, [url=http://bxgokflgsnob.com/]bxgokflgsnob[/url], [link=http://gtypaijmnvyz.com/]gtypaijmnvyz[/link], http://owbhxrykshts.com/~

//

#comment
アールメカブ

R_Tips0 のバックアップの現在との差分(No.11)
