2016-08-29 5 views
1

Я переписываю код с 2 лет назад и встречаю новую ошибку с функцией survSplit. Ошибка говорит, что мой объект не может быть найден, хотя это определенный столбец в моем фреймворке.Ошибка при использовании функции survSplit

Вот пример моего dataframe:

f12 <- data.frame(id = 1:6, 
       next.ivl= c(22.348, 1.837, 2.051,1.782,1.692, 1.730), 
       event = c(0,1,1,0,1,1), 
       enter= rep(0,6), 
       end=c(22.348, 1.837,2.051,1.782,1.629,1.730)) 

Мой предыдущий код был следующий:

cutpoints.l <- c(10/12, 1.25, 1.75, 2.25, seq(3,11)) 

f12.split <- survSplit(f12, 
        cut = cutpoints.l, 
        event = "event", 
        start = "enter", 
        end = "next.ivl", 
        episode = "ivl") 

Я пошел вперед и попытался запустить пример кода при условии here (также ниже), и я получаю еще одну ошибку («Ошибка в Surv (время, статус): временная переменная не является числовой». Эти ошибки, похоже, связаны. Что происходит? Почему мне нужно получить код ошибки?

library(survival)  
aml3<-survSplit(aml,cut=c(5,10,50),end="time",start="start", 
      event="status",episode="i") 

FYI: в настоящее время работает R версии 3.3.1 (2016-06-21)

Обновление для дальнейшего использования: правильный код должен быть

f12.split <- survSplit(Surv(next.ivl,event)~.,f12, 
       cut = cutpoints.l, 
       event = "event", 
       start = "enter", 
       end = "next.ivl", 
       episode = "ivl") 
+1

Пожалуйста, разместите свои данные, чтобы мы могли помочь –

+4

@CyrusMohammadian: Набор данных aml' является частью пакета выживания, который поставляется с каждым полным пакетом установки R. –

ответ

2

Первый аргумент должен быть формула и вы даете ему data.frame. Я понимаю, что эта ссылочная страница выглядит как недавняя (даже самая передовая) версия, смотрящая на URL, но текущая версия pkg: выживание - 2.39-5, и эта страница ссылается на версию 2.38-3. Текущая версия справочной страницы для survSplit ссылается на более раннюю версию функции, которая не имеет интерфейса формулы и имеет примеры, которые запускаются без ошибок в моем экземпляре R 3.3.1.

+0

Спасибо - похоже, что я использую старую формулировку функции. (И приношу извинения за дополнительную должность - вот мой первый вопрос!) – AliG