å©ç¨Hmiscå
ä¸çrcorr.censå½æ° å±éï¼ - åªè½å¤çä¸ä¸ªé¢æµåé - 对è¶
è¿2åç±»çåç±»åéå¤çç²ç³
# å è½½å
åçææ°æ®æ¡,è¿éçææ°æ®æ¡ä¸»è¦æ¯ä¸ºäºæ¹ä¾¿å¤§å®¶ç解ï¼å 为大家é常é½æ¯å°Excelçæ°æ®è¯»è¿Rï¼åå¨ä¸ºæ°æ®æ¡æ ¼å¼
library(survival)
library(Hmisc)
age <- rnorm(200, 50, 10)
bp <- rnorm(200,120, 15)
d.time <- rexp(200)
cens <- runif(200,.5,2)
death <- d.time <= cens
os <- pmin(d.time, cens)
sample.data <- data.frame(age = age,bp = bp,os = os,death = death)
#让æ们çä¸ä¸çæçä¾åæ°æ®çå6è¡
head(sample.data)
## age bp os death
## 1 33.18822 114.6965 1.106501 FALSE
## 2 41.86970 123.2265 1.365944 FALSE
## 3 50.41484 124.9522 0.867119 FALSE
## 4 45.66936 127.3237 1.155765 TRUE
## 5 39.79024 134.8846 1.257501 TRUE
## 6 31.89088 140.9382 1.125504 FALSE
rcorr.censç代ç åç»æï¼ç¬¬ä¸ä¸ªå¼å°±æ¯Cææ°ï¼åæ¶ä¹æDxyçå¼
rcorr.cens(sample.data$age, Surv(sample.data$os, sample.data$death))
## C Index Dxy S.D. n missing
## 4.528492e-01 -9.430156e-02 5.565299e-02 2.000000e+02 0.000000e+00
## uncensored Relevant Pairs Concordant Uncertain
## 1.290000e+02 3.172800e+04 1.436800e+04 8.072000e+03
rcorrcensç代ç åç»æ,注ærcorrcensçåæ³æ¯åæformulaï¼å
¬å¼ï¼çå½¢å¼ï¼è¾ä¸ºæ¹ä¾¿ï¼èrcorr.censç åæ³æ¯åªè½å¨åé¢åä¸ä¸ä¸ª
èªåé,并ä¸ä¸æ¯ædata = ...çåæ³ï¼æç¹ç¹çãè¾ä¸ºéæ¾çæ¯è¿ä¸¤ç§æ¹æ³å¾å°çCææ°ç
æ å误éè¦éè¿S.D./2é´æ¥å¾å°ã
r <- rcorrcens(Surv(os, death) ~ age + bp,data = sample.data)
r
## Somers' Rank Correlation for Censored Data Response variable:Surv(os, death)
##
## C Dxy aDxy SD Z P n
## age 0.453 -0.094 0.094 0.056 1.69 0.0902 200
## bp 0.498 -0.003 0.003 0.054 0.06 0.9517 200