2017-01-12 18 views
-1

私は、患者番号で索引付けされた臨床試験のデータを持っています。およびデータフレーム内の治療グループを含む。変数X1〜X7は、1〜7ヶ月目の観測値に対応します。一般的には、改善は測定パラメータの減少に相当します。時間枠を列としてデータフレームから時系列データをプロットする

Patient.No Group X1 X2 X3 X4 X5 X6 X7 
1  PX-002   1 2 24 4 82 2 2 2 
2  PX-009   5 77 6 36 33 1 2 0 
3  PX-019   1 6 7 5 4 3 2 2 
.... 
8  PX-001   54 7 7 5 5 3 2 1 
9  PX-004   2 67 3 2 3 2 1 1 
10  PX-013   3 5 24 4 2 1 1 2 
11  PX-015   2 9 4 33 63 2 1 1 

7ヶ月の間に、各患者の進化のために追跡する別々の時系列1と言うグループ-1のためにプロットするための良い方法は何でしょう。 x軸上に月がある?

enter image description here

が、私はそれを手動で行うことができますが、不格好に思える:

これは私が何を意味するかの醜いプロトタイプです。どんな良い方法?

plot(unlist(data[data$Patient.No=="PX-002",c("X1","X2","X3","X4","X5","X6","X7")]),type="b") 
lines(unlist(data[data$Patient.No=="PX-019",c("X1","X2","X3","X4","X5","X6","X7")]),type="b") 

それとも

grp1_data<-data[data$Treat..Group==1,] 

plot.ts(y = unlist(grp1_data[grp1_data$Patient.No=="PX-002",c("X1","X2","X3","X4","X5","X6","X7")]),x = c(1,2,3,4,5,6,7)) 

答えて

1

これはggplotで本当に簡単です。ところで、再現可能な例を投稿するべきです。あなたが持っていれば、私はあなたのためにこれをコード化していたでしょう。これがうまくいきたいです。

library(dplyr) 
library(tidyr) 
library(ggplot2) 

df %>% 
    gather(Month, Value, -Patient.No, -Group) %>% 
    ggplot(aes(Month, Value, color = Group)) + 
    geom_line(alpha = 0.75) + 
    ylab("Chosen Metric") + 
    ggtitle("Patient Progress by Group") 
+0

ありがとうDan!私はあなたの誇りを試しました。しかし、私はこのようなものを手に入れます。どんな考えが間違っているのか/ –

+0

Danさんに感謝します!私はあなたの提案を試みた。しかし、私はこのようなものを手に入れます。どんなアイデアが間違っているのでしょうか? http://imgur.com/a/faEH8 –

+0

あなたがあなたの質問を編集して[再現可能な例](http:// stackoverflow)でコピー/ペーストできるデータを提供すると、本当に役に立ちます。 com/questions/5963269/how-to-make-a-great-r-reproducible-example) – Zafar

関連する問題