-1

このコードは、1 つのウィンドウに複数のプロットを表示します。各プロットは、 X1 から X13 という名前の各列に従って作成されました。as a mainSO各図に、図をプロットするために使用される列の名前を追加したいと思います。以下を参照して、どれがどれであるかを区別しここに画像の説明を入力ます。

ref= read.table("D:\\AS_asc.txt", sep="",header=TRUE)
sour1 = read.table("D:\\re.txt", sep="",header=TRUE) 
sour2= read.table("D:\\_asc.txt", sep="",header=TRUE)
columns <- paste0("X", 1:13)
par(mfrow=c(4,4))
 lapply(
 columns,
 function(column)
  {
   result1 <- (
    mean(ref[[column]]) - 
     ((sd(ref[[column]]) / sd(sour1[[column]])) * mean(sour1[[column]])) + 
     ((sd(ref[[column]]) / sd(sour1[[column]]) * sour1[[column]]))
     )   # calculate using ref and sour1
     result2 <- ((  
      mean(ref[[column]]) - 
      ((sd(ref[[column]]) / sd(sour2[[column]], na.rm=TRUE)) * mean(sour2[[column]], na.rm=TRUE)) + 
     ((sd(ref[[column]]) / sd(sour2[[column]], na.rm=TRUE) * sour2[[column]])) 
    ))  # calculate using ref and sour2
      plot(
    ref[[column]],
    result1,
    ylab = "[[column]]",
    xlab = "[[column]]",
   col  = 2
    )

   points(ref[[column]], ref[[column]], col = 'green')
  points(ref[[column]], result2, col = 'blue')
   }
   )
4

2 に答える 2

6

これをコードに適応させることができると確信していますがmain = paste(column)、プロット コマンド内に追加するのと同じくらい簡単です。

例：

columns <- paste0( "X" , 1:9 )
par(mfrow=c(3,3))
lapply( columns , function(column){ plot(1:10 , y = 1:10 , main = paste(column ) ) } )

ここに画像の説明を入力

于 2013-03-13T19:00:53.490 に答える
3

pointsおそらく呼び出しの後に、その関数内に次の行を追加するだけでよいのではないでしょうか。

title(main=column)
于 2013-03-13T18:56:33.730 に答える