1

私は自分のデータに対してcoxphを実行し、次のような結果を取得しました。

> z
電話：
coxph（式= Surv（年、イベント）〜y、データ= x）


    coef exp（coef）se（coef）zp
y 0.0714 1.07 0.288 0.248 0.8


尤度比検定=0.06、1 df、p = 0.804 n = 65、イベント数= 49

保存したいだけです

y 0.0714 1.07 0.288 0.248 0.8

ファイルに。順列を実行して1000zを生成するためです。

次のようなテキストファイルに保存したいと思います。

フィン-0.37940.6840.1914 -1.983 0.0470
年齢-0.05740.9440.0220 -2.611 0.0090
レース0.31391.3690.3080 1.019 0.3100
wexp -0.1498 0.861 0.2122 -0.706 0.4800
3月-0.43370.6480.3819 -1.136 0.2600
パロ-0.08490.9190.1958 -0.434 0.6600

誰でも助けることができますか？

ありがとう！

4

1 に答える 1

2

係数は、次の方法で簡単にアクセスできます。

summary(z)[['coefficients']]

および信頼区間情報

summary(z)[['conf.int']]

summary.coxphオブジェクトの構成要素を調べるには

str(summary(z))

私のアドバイスは、順列のリストを作成することです

data_list <- list(data_1, ...., data_1000)

それから電話する

lots_models <- lapply(data_list, coxph, formula = Surv(Years, Event) ~ y)

モデルのリストを作成する

要約は次の方法で作成できます。

lots_summaries <- lapply(lots_models, summary)

係数を抽出する

all_coefficients <- lapply(lots_summaries, '[[', 'coefficients')
all_conf.int     <- lapply(lots_summaries, '[[', 'conf.int')

順列 ID 列を追加します (必要な場合)。

all_coefs_id <- lapply(seq_along(data_list), 
   function(i) cbind(all_coefficients[[i]],i))
all_ci_id <- lapply(seq_along(data_list), 
   function(i) cbind(all_conf.int[[i]],i))

次に、data.frameに結合します

all_coefs_df <- do.call(rbind, all_coefs_id)
all_ci_df <- do.call(rbind, all_ci_id)

その後、テキストファイルとして保存します

于 2012-08-09T00:34:56.377 に答える