[解決済み] 相対生存率解析のための率表の作成
2022-02-15 13:06:21
質問内容
relsurvパッケージは相対生存率を計算することができます。
入力は、データ(私の場合は大腸癌と診断された患者)と、次の年齢区分までの生存率の変化(一般的には性別と診断年別)を含む率表です。
私の場合、年齢、性別、診断年、社会経済的スコアに応じた相対生存率を得たいのです。1つのデータ・フレームに率表に入るべきすべてのデータを持っています。
> head(lifetable)
sex seifa age 2007 2008 2009 2010 2011 2012 2013
1 female SEIFA1 0 0.9947100 0.9953400 0.9956200 0.9959200 0.9964100 0.9966000 0.9966800
2 female SEIFA1 1 0.9996625 0.9996613 0.9997038 0.9996489 0.9996425 0.9996667 0.9997556
3 female SEIFA1 2 0.9997750 0.9997968 0.9998105 0.9998119 0.9998350 0.9998167 0.9998656
4 female SEIFA1 3 0.9998375 0.9998374 0.9998460 0.9998621 0.9998625 0.9998667 0.9999022
5 female SEIFA1 4 0.9998750 0.9998645 0.9998697 0.9998871 0.9998900 0.9999000 0.9999267
6 female SEIFA1 5 0.9998827 0.9999077 0.9998887 0.9998956 0.9999037 0.9999160 0.9999280
これを「ratetable」に変換する手っ取り早い方法を探しています。いろいろ検索してみたのですが、何も見つかりませんでした...。どなたか教えてください。
どのように解決するのですか?
このトピックを閉じるために、私は解決策を見つけました。
- 標準レート表(年齢、性別、年別)を5つ作成し、各世帯のカテゴリーごとに1つずつ作成する。
for (i in seq_along(levels(lifetable$seifa))){ # Create standard ratetables (by age,sex,year) for each seifa category lt_ss_f <- subset(lifetable, seifa == paste("SEIFA", i, sep="") & sex =="female"); lt_ss_f <- lt_ss_f[-c(1:3)]; lt_ss_f <- as.matrix(lt_ss_f) lt_ss_m <- subset(lifetable, seifa == paste("SEIFA", i, sep="") & sex =="male") ; lt_ss_m <- lt_ss_m[-c(1:3)]; lt_ss_m <- as.matrix(lt_ss_m) assign(paste("SEIFA", i, sep=""), transrate(lt_ss_m,lt_ss_f, yearlim=c(2003,2013), int.length=1))}
- joinrate()を使用して5つの料金テーブルを結合し、4次元(年齢、性別、年、seifa)を含む完全な料金テーブルを作成します。
FULLrateTable <- joinrate(list(SEIFA1=SEIFA1, SEIFA2=SEIFA2, SEIFA3=SEIFA3, SEIFA4=SEIFA4, SEIFA5=SEIFA5),dim.name="SEIFA")
関連
-
[解決済み】"置換する項目数が置換長の倍数でない "と表示されるのはなぜですか?
-
[解決済み】R: 複数行の ggplot2 コマンドで「単項演算子エラー」が発生する。
-
[解決済み] 変数の型(リスト)が無効です
-
[解決済み] テスト
-
[解決済み】Rで立方根と対数変換をする
-
[解決済み】R Markdown - html出力でフォントサイズとフォントタイプを変更する
-
[解決済み] lmer エラー: グループ化係数は観測数未満でなければならない
-
[解決済み】bstTreeの予測値の混乱行列、エラー:'The data must contain some levels that overlap reference'.
-
[解決済み] 因子を日付形式に変換するにはどうすればいいですか?
-
[解決済み】行列式で「数値/複雑な行列/ベクトルの引数を必要とする」というエラーが発生する?
最新
-
nginxです。[emerg] 0.0.0.0:80 への bind() に失敗しました (98: アドレスは既に使用中です)
-
htmlページでギリシャ文字を使うには
-
ピュアhtml+cssでの要素読み込み効果
-
純粋なhtml + cssで五輪を実現するサンプルコード
-
ナビゲーションバー・ドロップダウンメニューのHTML+CSSサンプルコード
-
タイピング効果を実現するピュアhtml+css
-
htmlの選択ボックスのプレースホルダー作成に関する質問
-
html css3 伸縮しない 画像表示効果
-
トップナビゲーションバーメニュー作成用HTML+CSS
-
html+css 実装 サイバーパンク風ボタン
おすすめ
-
[解決済み】 'height' はベクトルか行列でなければならない barplot のエラー
-
[解決済み】scale_color_manual()が動作しない件
-
[解決済み】「次のオブジェクトは 'package:xxx' からマスクされています」とはどういう意味ですか?
-
[解決済み】reshape2 meltの警告メッセージ
-
[解決済み】Rで立方根と対数変換をする
-
[解決済み] lmer エラー: グループ化係数は観測数未満でなければならない
-
[解決済み】R: predict() の数値 'envir' 引数が長さ1でない。)
-
[解決済み】R 置換する項目数が置換長の倍数でない/しかし結果は正しい
-
[解決済み】Rで結果の行数がベクトル長(arg 2)の倍数でない件
-
[解決済み】R4DSのエラー比較(1)は、アトミック型とリスト型でのみ可能です