1

を使用してRstudioいます。ここggsurvのチュートリアルから、データのカプラン・マイヤー曲線を描画するためにパッケージの関数を使用GGallyしています（生存分析用）。凡例を単独で処理するため、代わりに使用しています。plotggsurv

リンクに示されているように、複数の曲線は色で区別されます。線種で区別したい。チュートリアルにはそのためのオプションがないようです。以下は私のコマンドです：

surv1 <- survfit(Surv(DaysOfTreatment,Survived)~AgeOnFirstContactGroup)
print(ggsurv(surv1, lty.est = 3)+ ylim(0, 1))

lty.est=3（または2）すべての行に同じ破線が表示されます。行ごとに異なる破線が必要です。を使用lty=typeするとエラーが発生します: object 'type' not found。lty=typeで動作しますggplotが、プロットggplotを直接処理しません。survfit

ggsurvどちらかまたは単純な線種で曲線を区別する方法を教えてください（ただし、凡例を処理するため、plot好みます）ggsurv

4

1 に答える 1

9

ggsurv のドキュメントから

lty.est: 生存曲線の線種。ベクトルの長さは、1 または階層の数に等しい必要があります。

したがって、層ごとに異なる線種を取得するには、プロットする線の数と同じ長さのベクトルに等しい lty.est を設定します。各値は異なる線種に対応します。

たとえばsurvival、パッケージの肺データを使用して

library(GGally)
library(survival)
data(lung)
surv1 <- survfit(Surv(time,status) ~ sex, data = lung)
ggsurv(surv1, lty.est=c(1,2), surv.col = 1)

次のプロットを与える

ここに画像の説明を入力

ggplot テーマやその他の ggplot 要素をプロットに追加することもできます。たとえば、次のように Cowplot テーマを使用して外観を改善できます。

library(ggplot2)
library(cowplot)
ggsurv(surv1, lty.est=c(1,2), surv.col = 1) + theme_cowplot()

ここに画像の説明を入力

線種で区別した後に凡例ラベルを変更する必要がある場合は、この方法で行うことができます

ggsurv(surv1, lty.est=c(1,2), surv.col = 1) +
  guides(colour = FALSE) +
  scale_linetype_discrete(name   = 'Sex', breaks = c(1,2), labels = c('Male', 'Female'))

ここに画像の説明を入力

于 2016-05-28T22:47:05.673 に答える