2017-12-30 20 views
-2

Rの単純な折れ線グラフで、モンテカルロシミュレーションの結果をプロットしようとしています。考え方はx軸S0とStとYの値を各ランで取得することです。モンテカルロプロットの結果

Run #  S0 St 
Call Run 1 100 141,952 
Call Run 2 100 106,711 
Call Run 3 100 97,0053 
Call Run 4 100 107,015 
Call Run 5 100 64,1486 

私が直面している問題は、私がプロットするとき、StはY軸、S0はX軸とみなされることです。誰かが私がこれをどのように動かすことができるかについてのアドバイスを提供することができますか? PS

:私のMCのsimuationsをcsv形式で(1000回のシミュレーションを中心に)C++で行われ、私はデータを入れてはR.

にロードするためにファイルをエクセルこの、tidyverseアプローチを使用して、事前

+0

あなたが間違ったコードを共有しているかどうか、あなたが期待しているプロットの例を –

+0

これは、http://www.goddardconsulting.ca/image-files/MonteCarloPathsManyを探しているグラフィックのようなものです.jpg Imは非常に新鮮な初心者のRプロットで、私は基本的なプロット関数を使用しました。 –

+0

あなたの出力には出発点と終了点しかないようです。中間の反復に対応する値はありません。 –

答えて

0

にありがとうコードです。あなたがgeom_lineggplotを実行するために、あなたのデータフレームを再構築する必要があります。

mydf <- data.frame(
    run=c(paste0("Run",1:5)), 
    S0=c(rep(x=100,5)), 
    St=c(141.952,106.711,97.0053,107.015,64.1486) 
) 

mydf_long <- mydf %>% 
    gather(key=Time_sim,value=value,c(S0,St)) %>% 
    mutate(Time_sim=str_replace_all(string=Time_sim,pattern="S",replacement=""), 
      Time_sim=str_replace_all(string=Time_sim,pattern="t",replacement="1")) %>% 
    arrange(run,Time_sim) 

mydf_long$Time_sim <- as.numeric(mydf_long$Time_sim) 

ggplot(data=mydf_long,aes(y = value, x = Time_sim, colour = run)) + 
    geom_line() 

これはプロットである:

enter image description here

正直に言うと、私はから始まる、あなたがRであなたの分析を複製することをお勧めMCのシミュレーション。

+0

ありがとうございました。これは私が探していたものです。私は時間1,2と3を言うことができたら、私はこのコードを編集してより多くのステップを追加すべきですか? –

+0

私のコードを実行するときは、長形式のデータフレームを探索します( 'view(mydf_long)'を使用して)。データの構造を理解する必要があると思います。 –