Выдержанный тест с одной выборкой - как я могу использовать смещение в формуле?

Question

Я хотел бы использовать функцию Survdiff для теста с одним образцом.

Я не понял, куда мне включать в фурмулу смещение.

surdiff должно выглядеть следующим образом:

survdiff(formula, data, subset, na.action, rho=0)

Например, в тесте k-samples формула будет иметь вид: Surv(futime, fustat), а предикторы: rx,И мы будем использовать Surdiff следующим образом:

survdiff(Surv(futime, fustat) ~ rx,data=ovarian) 

(Этот пример взят из примера R)

Что я не понимаю, так это то, как я должен включать смещениетермин в предикторах в случае, если я хочу выполнить тест с одной выборкой.В объяснении R они говорят, что я должен использовать смещение (sp), где sp - вектор вероятностей.

Кто-нибудь может сказать мне, каков правильный синтаксис для этого, и они имеют в виду «включить»?

Заранее спасибо !!Alley

Answer 1

Вы включаете offset(var) как только вещь в RHS формулы модели, и var должно быть вероятность выживания , последняя может быть вычислена с использованием survexp function.

Если бы вы прочитали немного дальше в ?survdiff, вплоть до примеров, вы бы видели это:

## Expected survival for heart transplant patients based on
## US mortality tables
expect <- survexp(futime ~ ratetable(age=(accept.dt - birth.dt),
                  sex = 1, year = accept.dt, race = "white"), jasa,
                  cohort = FALSE, ratetable = survexp.usr)

## actual survival is much worse (no surprise)
survdiff(Surv(jasa$futime, jasa$fustat) ~ offset(expect))

, который показывает вам точно как включить один offset термин в формулу survdiff.