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