Модель Fine и Grey в R с вмененными наборами данных

Question

У меня есть длинный (вертикально сложенный) набор данных, содержащий 10 импутаций (переменная «импутация» определяет номер импутации). Вменение было сделано в SAS, но я хотел бы вычислить некоторую c-статистику, используя R.

Я знаю, как рассчитать c-stats, используя функцию cindex и функцию FGR для одного вмененного набора данных. Я не уверен, как я повторю это в вертикально сложенном наборе данных. Я пытался использовать функцию «с», но безуспешно.

Вот мои коды:

fgr.model <- FGR(Hist(time, outcome) ~ x1 + x2 + x3, data=mydata1, cause=1)

cscore <- cindex(list(fgr.model), forumula=Hist(time, outcome)~1,
          cens.model="marginal", data=mydata1, eval.time=c(1826), cause=1)

Как рассчитать c-stats с использованием функции cindex и функции FGR в наборе данных с вертикальным стеком?

Answer 1

Основываясь на деталях в комментариях: я понял, что вам нужно рассчитать статистику по разделу исходного кадра данных mydata1 - т.е. вам нужно выбрать только те строки, которые соответствуют одному «вменению» (этоимя, которое вы используете, ссылаетесь ли вы с неправильным словом на «входной сеанс»? просто любопытно)

Сначала вы должны создать новый фрейм данных, содержащий только данные для одного «вменения», в следующемВ примерах мы рассмотрим операцию для вменения № 4. Существуют различные способы выполнения работы.

Первый способ, который работает, если имена столбцов верны:

mydata1portion = mydata1[mydata1$ImputationNo==4,]

Второй, который работает, еслиПоложение / порядок столбцов правильный:

mydata1portion = mydata1[mydata1$V1==4,]

Третий способ в случае упорядочения вменений / строк данных.

mydata1portion = mydata1[(100*(4-1)+1):(100*4),]

Первые два вы используете значение для столбца ImputationNo чтобы отфильтровать фрейм данных, в последнем вы вырезаете фрейм данных на основе положения строк.

Наконец, вы можете вычислить статистикуна полученных mydata1portion и НЕ на полных mydata1.