[R] Survival model error
Meghna Govil
megov31 at yahoo.com
Mon Nov 6 14:45:35 CET 2017
Thanks David. Could you show me how to do that in my example ?
Thanks,
Meghna
> On Nov 6, 2017, at 12:58 AM, David Winsemius <dwinsemius at comcast.net> wrote:
>
> You should stop trying to use matrices on the RHS and using separate vectors to Surv. Instead use a data argument and have the names in your formula refer to column names.
>
> —
> David
>
> Sent from my iPhone
>
>> On Nov 5, 2017, at 7:21 PM, Meghna Govil via R-help <r-help at r-project.org> wrote:
>>
>>
>> Hi - Below is my code and then the error when I run the last line.
>>
>> time_np <- train1_na$tte
>> event_np <- train1_na$censored
>>
>>
>> X_np <- cbind(
>> train1_na$AMT,
>> train1_na$DISCOUNT_AMT,
>> train1_na$high_price_pcnt,
>> train1_na$EM_RECEIVED,
>> train1_na$DM_RECEIVED,
>> train1_na$TXN_WITH_RINGCODE,
>> train1_na$WEB,
>> train1_na$clearance_pcnt,
>> train1_na$bts_pcnt,
>> train1_na$sales_pcnt,
>> train1_na$holiday_pcnt,
>> train1_na$TXN,
>> train1_na$REDEEMED_REWARDS
>>
>> )
>>
>> # Kaplan-Meier non-parametric analysis
>>
>> kmsurvival_np <- survfit(Surv(time_np,event_np) ~ X_np)
>>
>> Error in `[.default`(y, who, 1) : (subscript) logical subscript too long
>>
>> Any ideas?
>> I have tried several things and still get this error.
>>
>>
>> Thanks,
>> Meghna
>>
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
More information about the R-help
mailing list