差分

このページの2つのバージョン間の差分を表示します。

この比較画面へのリンク

次のリビジョン		前のリビジョン
r:survival:plots [2018/12/11 15:58] – created watalur:survival:plots [2018/12/11 16:13] (現在) – watalu
行 1: 行 1:
-==== 信頼性データ図 ====+==== 信頼性データ図 Reliability Data Diagram ====
  
 横軸に経過時間を取り、縦軸に個人、製品、IDなどを取って、寿命データの中の個々の対象の寿命の長さを横線で表した図を、信頼性データ図というらしい。これを、データの並び順と、小さい順の２種類の順序、描けるようにした関数が次のplot.reliability.data.diagramである。 横軸に経過時間を取り、縦軸に個人、製品、IDなどを取って、寿命データの中の個々の対象の寿命の長さを横線で表した図を、信頼性データ図というらしい。これを、データの並び順と、小さい順の２種類の順序、描けるようにした関数が次のplot.reliability.data.diagramである。
行 6: 行 6:
  
 <code> <code>
-plot.reliability.data.diagram <- function(x,sort=TRUE,...) {+plot.reliability.data.diagram <- function(x,sort=TRUE,xlab=NULL,ylab=NULL,cex.plot=1,...) {
   if(is.numeric(x)==TRUE) {   if(is.numeric(x)==TRUE) {
     n <- length(x)     n <- length(x)
行 41: 行 41:
     }     }
   }   }
 +  if(is.null(xlab)) {
 +    xlab = "Time"
 +  }
 +  if(is.null(ylab)) {
 +    ylab = "#"
 +  }
 +  par(cex=cex.plot)
   plot(x.time,c(1:n),   plot(x.time,c(1:n),
        xlim=c(0,max(x.time)),        xlim=c(0,max(x.time)),
        type="n",        type="n",
-       ylab="#", +       ylab=ylab, 
-       xlab="Lifetime data", +       xlab=xlab, 
-       yaxt="n") +       yaxt="n", 
-  axis(2,1:n,labels=c(n:1)) +       cex=cex.plot) 
-  points(x.time,c(1:n),pch=c(1,4)[x.status+1],...)+  axis(2,1:n,labels=c(n:1), 
 +       cex=cex.plot) 
 +  points(x.time,c(1:n),pch=c(1,4)[x.status+1], 
 +       cex=cex.plot,...)
   for( i in c(1:n) ) {   for( i in c(1:n) ) {
-    lines(c(0,x.time[i]),c(i,i),...)+    lines(c(0,x.time[i]),c(i,i), 
 +       cex=cex.plot,...)
   }   }
 } }
 </code> </code>
  
 +使い方は
 +<code>
 +X = data.frame(time=c(1.5,5.3,2.5,7.8,6.2,4.1,8.2,10.3,1.5,13.4),
 +                status=c(1,1,1,1,1,1,1,1,1,1))
 +plot.reliability.data.diagram(X,sort=TRUE)
 +</code>
 +
 +これで次のグラフが描かれる。
 +
 +{{:r:survival:reliability-data-diagram-1.png|}}
 +
 +もしくは
 +<code>
 +X = data.frame(time=c(1.5,5.3,2.5,7.8,6.2,4.1,8.2,10.3,1.5,13.4),
 +                status=c(1,1,1,1,1,1,1,1,1,1))
 +plot.reliability.data.diagram(X,sort=FALSE)
 +</code>
 +
 +データの並びが変わる。
 +
 +{{:r:survival:reliability-data-diagram-2.png|}}
 +
 +軸ラベルを指定したり、グラフの文字や打点の大きさを変えることもできる。
 +<code>
 +X = data.frame(time=c(1.5,5.3,2.5,7.8,6.2,4.1,8.2,10.3,1.5,13.4),
 +                status=c(1,1,1,1,1,1,1,1,1,1))
 +plot.reliability.data.diagram(X,sort=TRUE,cex.plot=0.8,xlab="Time (days)",ylab="No.")
 +</code>
 +
 +{{:r:survival:reliability-data-diagram-3.png|}}
 +
 +==== レキシスダイヤグラム Lexis Diagram ====
 +
 +いつか加筆予定。