Closed tdhock closed 8 years ago
hey @anujkhare I ran your code on the ovarian data set
data(ovarian) X <- with(ovarian, cbind(age=age, residual.disease=resid.ds-1, treatment=rx-1, ecog.ps=ecog.ps-1)) y_l <- ovarian$futime y_r <- ovarian$futime y_r[ovarian$fustat == 0] <- NA y_surv <- log(cbind(y_l, y_r)) fit <- iregnet(X, y_surv, family="gaussian")
currently we have
> fit$beta[, 1:7] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 6.772081 6.772118 6.772103 6.772112 6.772108 6.77211 6.772109 [2,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.00000 0.000000 [3,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.00000 0.000000 [4,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.00000 0.000000 [5,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.00000 0.000000 >
Can you please add row names to the beta matrix? (Intercept) and column names of input X matrix, same as glmnet.
(Intercept)
hey @anujkhare I ran your code on the ovarian data set
currently we have
Can you please add row names to the beta matrix?
(Intercept)
and column names of input X matrix, same as glmnet.