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