機械がまだ動いている確率をグラフにする

信頼性の低い機械がある工場の監督者だとしたら、時間の経過とともに機械が動き続ける確率を把握したいはずです。この演習では、最初の30日間にわたって機械がまだ動いている確率をプロットします。

この演習はコースの一部です

Rで学ぶ確率の基礎

1日目から30日目まで、各日に機械が「まだ動いている」確率のベクトルを計算し、still_working として保存します。これは、最初の引数に数値のベクトルを渡して pgeom() を1回呼び出すだけで行えます。機械は毎日10%の確率で故障します。
コマンド qplot(still_working) を実行して、最初の30日それぞれで機械がまだ動いている確率をグラフ化します。x軸に日数、y軸に確率が表示されます。

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Calculate the probability of machine working on day 1-30
still_working <-

# Plot the probability for days 1 to 30
qplot(1:30, still_working)