summary, tables, and correlations for data frames
dcor(data, y = NULL, x = NULL, use = "pairwise.complete.obs", ...)data("sTRACE",package="timereg")
dt<- sTRACE
dt$time2 <- dt$time^2
dt$wmi2 <- dt$wmi^2
head(dt)
#> no wmi status chf age sex diabetes time vf time2 wmi2
#> X1944 1944 1.5 9 0 84.924 1 1 1.345000 0 1.809025 2.25
#> X5783 5783 1.9 0 1 74.193 0 0 6.910000 0 47.748100 3.61
#> X784 784 0.8 9 0 78.081 0 1 0.196000 0 0.038416 0.64
#> X3763 3763 1.3 0 0 55.479 1 0 7.543000 0 56.896849 1.69
#> X2927 2927 1.6 0 1 62.997 0 0 7.126000 0 50.779876 2.56
#> X4511 4511 1.0 9 1 67.644 1 0 4.532606 0 20.544520 1.00
dcor(dt)
#> no wmi status chf age
#> no 1.0000000000 -0.009437654 0.0002042854 0.03176893 -0.08317491
#> wmi -0.0094376537 1.000000000 -0.3103788196 -0.37464791 -0.19288208
#> status 0.0002042854 -0.310378820 1.0000000000 0.34237066 0.43051203
#> chf 0.0317689348 -0.374647914 0.3423706564 1.00000000 0.34197944
#> age -0.0831749071 -0.192882080 0.4305120263 0.34197944 1.00000000
#> sex -0.1090494653 -0.023542809 -0.0671941443 -0.13205952 -0.27356308
#> diabetes 0.0103460848 -0.138602457 0.0866091232 0.08871325 0.04685823
#> time 0.0185007942 0.328582293 -0.8145783614 -0.34820213 -0.41435287
#> vf 0.0618025722 -0.089668054 0.0473290244 0.13467109 -0.03983927
#> time2 0.0091527320 0.313368613 -0.8526160479 -0.34533562 -0.41370693
#> wmi2 -0.0081611472 0.987512109 -0.2992334471 -0.36641743 -0.18629143
#> sex diabetes time vf time2
#> no -0.10904947 0.01034608 0.01850079 0.06180257 0.009152732
#> wmi -0.02354281 -0.13860246 0.32858229 -0.08966805 0.313368613
#> status -0.06719414 0.08660912 -0.81457836 0.04732902 -0.852616048
#> chf -0.13205952 0.08871325 -0.34820213 0.13467109 -0.345335622
#> age -0.27356308 0.04685823 -0.41435287 -0.03983927 -0.413706928
#> sex 1.00000000 -0.05523671 0.03630103 -0.04760863 0.037530286
#> diabetes -0.05523671 1.00000000 -0.09049143 -0.05060615 -0.087337711
#> time 0.03630103 -0.09049143 1.00000000 -0.10704592 0.976099537
#> vf -0.04760863 -0.05060615 -0.10704592 1.00000000 -0.076612883
#> time2 0.03753029 -0.08733771 0.97609954 -0.07661288 1.000000000
#> wmi2 -0.01543475 -0.12474865 0.31427694 -0.09003700 0.302021849
#> wmi2
#> no -0.008161147
#> wmi 0.987512109
#> status -0.299233447
#> chf -0.366417431
#> age -0.186291429
#> sex -0.015434752
#> diabetes -0.124748650
#> time 0.314276938
#> vf -0.090037003
#> time2 0.302021849
#> wmi2 1.000000000
dcor(dt,~time+wmi)
#> time wmi
#> time 1.0000000 0.3285823
#> wmi 0.3285823 1.0000000
dcor(dt,~time+wmi,~vf+chf)
#> vf: 0
#> chf: 0
#> time wmi
#> time 1.0000000 0.2443018
#> wmi 0.2443018 1.0000000
#> ------------------------------------------------------------
#> vf: 1
#> chf: 0
#> time wmi
#> time 1.0000000 0.8466897
#> wmi 0.8466897 1.0000000
#> ------------------------------------------------------------
#> vf: 0
#> chf: 1
#> time wmi
#> time 1.0000000 0.1903859
#> wmi 0.1903859 1.0000000
#> ------------------------------------------------------------
#> vf: 1
#> chf: 1
#> time wmi
#> time 1.000000 0.354452
#> wmi 0.354452 1.000000
dcor(dt,time+wmi~vf+chf)
#> vf: 0
#> chf: 0
#> time wmi
#> time 1.0000000 0.2443018
#> wmi 0.2443018 1.0000000
#> ------------------------------------------------------------
#> vf: 1
#> chf: 0
#> time wmi
#> time 1.0000000 0.8466897
#> wmi 0.8466897 1.0000000
#> ------------------------------------------------------------
#> vf: 0
#> chf: 1
#> time wmi
#> time 1.0000000 0.1903859
#> wmi 0.1903859 1.0000000
#> ------------------------------------------------------------
#> vf: 1
#> chf: 1
#> time wmi
#> time 1.000000 0.354452
#> wmi 0.354452 1.000000
dcor(dt,c("time*","wmi*"),~vf+chf)
#> vf: 0
#> chf: 0
#> time time2 wmi wmi2
#> time 1.0000000 0.9744619 0.2443018 0.2286891
#> time2 0.9744619 1.0000000 0.2345873 0.2233235
#> wmi 0.2443018 0.2345873 1.0000000 0.9913077
#> wmi2 0.2286891 0.2233235 0.9913077 1.0000000
#> ------------------------------------------------------------
#> vf: 1
#> chf: 0
#> time time2 wmi wmi2
#> time 1.0000000 0.9942844 0.8466897 0.8244614
#> time2 0.9942844 1.0000000 0.8793376 0.8612097
#> wmi 0.8466897 0.8793376 1.0000000 0.9968521
#> wmi2 0.8244614 0.8612097 0.9968521 1.0000000
#> ------------------------------------------------------------
#> vf: 0
#> chf: 1
#> time time2 wmi wmi2
#> time 1.0000000 0.9731838 0.1903859 0.1826388
#> time2 0.9731838 1.0000000 0.1762251 0.1697576
#> wmi 0.1903859 0.1762251 1.0000000 0.9833844
#> wmi2 0.1826388 0.1697576 0.9833844 1.0000000
#> ------------------------------------------------------------
#> vf: 1
#> chf: 1
#> time time2 wmi wmi2
#> time 1.0000000 0.9828671 0.3544520 0.3026710
#> time2 0.9828671 1.0000000 0.2827122 0.2298121
#> wmi 0.3544520 0.2827122 1.0000000 0.9842457
#> wmi2 0.3026710 0.2298121 0.9842457 1.0000000