差分

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

この比較画面へのリンク

両方とも前のリビジョン前のリビジョン
次のリビジョン		前のリビジョン
r:maintenance:hardtime [2018/12/10 09:16] – [ブロック取り替え] watalur:maintenance:hardtime [2018/12/10 09:18] (現在) – [時間計画保全] watalu
行 5: 行 5:
  
 <code> <code>
-select.distribution <- function(dist="weibull") {+select.distribution = function(dist="weibull") {
   if(dist %in% c("weibull","gamma","lnorm")) {   if(dist %in% c("weibull","gamma","lnorm")) {
     distrname <- dist     distrname <- dist
行 87: 行 87:
 <code> <code>
 # 確率密度関数 # 確率密度関数
-f <- function(x) {+f = function(x) {
   return(model$d(x,model.parameters[1],model.parameters[2]))   return(model$d(x,model.parameters[1],model.parameters[2]))
 } }
行 95: 行 95:
 <code> <code>
 # 生存関数 # 生存関数
-F.bar <- function(x) {+F.bar = function(x) {
   return(1-model$p(x,model.parameters[1],model.parameters[2]))   return(1-model$p(x,model.parameters[1],model.parameters[2]))
 } }
行 106: 行 106:
 <code> <code>
 # コストレートの式 (integrateは1変数関数の数値積分をしてくれるRの関数) # コストレートの式 (integrateは1変数関数の数値積分をしてくれるRの関数)
-g <- function(x) {+g = function(x) {
   return((Cc*integrate(f,0,x)$value+Cp*(1-integrate(f,0,x)$value))/(integrate(F.bar,0,x)$value))   return((Cc*integrate(f,0,x)$value+Cp*(1-integrate(f,0,x)$value))/(integrate(F.bar,0,x)$value))
 } }
行 134: 行 134:
 <code> <code>
 # コストレートのグラフを描く # コストレートのグラフを描く
-g.list <- NULL+g.list = NULL
 for( i in c(1:500) ) { for( i in c(1:500) ) {
-  g.list <- append(g.list,g(i/10))+  g.list = append(g.list,g(i/10))
 } }
 plot(c(1:500)/10,g.list,type="b") plot(c(1:500)/10,g.list,type="b")