kassambara / survminer

Survival Analysis and Visualization
https://rpkgs.datanovia.com/survminer/
504 stars 162 forks source link

Cannot compute pvalue when using interval censoring #311

Open sbihorel opened 6 years ago

sbihorel commented 6 years ago

Expected behavior

p-value is calculated when the data contain interval censoring.

Actual behavior

Errors are returned no matter which log.rank.weights value is set

### for instance, default log.rank.weights value
Error in survival::survdiff(eval(fit$call$formula), data = data) : 
  Right censored data only

Steps to reproduce the problem

require(survival)
require(survminer)
fit <- survfit(Surv(time=start, time2=stop, event = event, type = 'interval')~transplant, data = heart)
ggsurvplot(fit = fit, data = heart, pval = TRUE)

session_info()

R version 3.4.4 (2018-03-15)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Linux Mint 18.1

Matrix products: default
BLAS: /usr/lib/libblas/libblas.so.3.6.0
LAPACK: /usr/lib/lapack/liblapack.so.3.6.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8    LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] survival_2.42-3 survminer_0.4.2 ggpubr_0.1.6    magrittr_1.5    ggplot2_2.2.1  

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.11        pillar_1.2.1        compiler_3.4.4      plyr_1.8.4          bindr_0.1           tools_3.4.4        
 [7] tibble_1.4.2        nlme_3.1-137        gtable_0.2.0        lattice_0.20-35     pkgconfig_2.0.1     rlang_0.2.0        
[13] Matrix_1.2-14       psych_1.7.8         cmprsk_2.2-7        yaml_2.1.14         parallel_3.4.4      bindrcpp_0.2       
[19] gridExtra_2.2.1     knitr_1.17          dplyr_0.7.4         stringr_1.2.0       survMisc_0.5.4      grid_3.4.4         
[25] data.table_1.10.4-3 glue_1.2.0          R6_2.2.2            KMsurv_0.1-5        km.ci_0.5-2         foreign_0.8-69     
[31] tidyr_0.8.0         purrr_0.2.4         reshape2_1.4.2      scales_0.4.1        splines_3.4.4       assertthat_0.2.0   
[37] mnormt_1.5-5        xtable_1.8-2        colorspace_1.3-2    labeling_0.3        stringi_1.1.5       lazyeval_0.2.0     
[43] munsell_0.4.3       broom_0.4.3         zoo_1.8-1 
dchiu911 commented 5 years ago

any updates with this issue?

sbkanaan commented 2 years ago

Still having the same issue with the pval when using interval censoring. any updates?

richardyzx commented 1 year ago

any updates with this issue?

antoine4ucsd commented 7 months ago

any update? thank you