search for: yongchuan

...ble. How are 1.1, 1.2, 1.3 duplicate row names? Thx. > table <- read.table('latestWithNumber.txt', header=T) Error in "row.names<-.data.frame"(`*tmp*`, value = c("1.1", "1.2", "1.3", : duplicate 'row.names' are not allowed Yongchuan

...model$frail returns the group level frailty terms. Does this mean this is the average of the frailty values for the respective groups? Also, if I'm fitting it to a gamma frailty, how do I extract the rate and scale parameters for the different gamma distributions fit to each group? Thanks. Yongchuan

...coxph object looks more correct. > survfit(resulting) Call: survfit.coxph(object = resulting) n events median 0.95LCL 0.95UCL 115 15 Inf Inf Inf Is there a limit to the size of the data set that I read in? Or am I just doing something silly above? Thanks much. Yongchuan (this is the coxph regression: resultag <- coxph(Surv(Start,Stop,PrepayDate)~modBalance + closingCoupon+lienPosition +originalFICO,table)

Question: When survfit() function is used upon a coxph object, the 'n' returned is vastly smaller (n=6) than the number of distinct loans in the dataset used. I am trying to estimate a Cox proportional hazards model for a set of loans (over 6000) using using time varying covariates. For this 6000+ loans, I have some 62,000 different vectors representing the loans at different periods of