关于survival包里的survdiff函数的源程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助大神，我在查看survival包里的survdiff源代码时，发现其中有一个survdiff.fit函数，但我在R里面调用它时却提示没有这个对象，在网上搜索也搜不到相关信息，请问有哪位大神知道这个函数是什么吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Survival Diff SUR DIF 源程序 survival 源程序

R是调用S语言的函数，survdiff.fit是S语言的函数。代码在survival包里找到的，如下：

1. survdiff.fit <- function(y, x, strat, rho=0) {
   
2.     #
   
3.     # This routine is almost always called from survdiff
   
4.     #  If called directly, remember that it does no error checking
   
5.     #
   
6.     n <- length(x)
   
7.     if (ncol(y) !=2) stop ("Invalid y matrix")
   
8.     if (nrow(y) !=n | length(x) !=n) stop("Data length mismatch")
   
9. 
   
10.     ngroup <- length(unique(x))
    
11.     if (ngroup <2) stop ("There is only 1 group")
    
12.     if (is.category(x)) x <- as.numeric(x)
    
13.     else x <- match(x, unique(x))
    
14. 
    
15.     if (missing(strat)) strat <- rep(1,n)
    
16.     else strat <- as.numeric(as.factor(strat))
    
17.     nstrat <- length(unique(strat))
    
18.     if (length(strat) !=n) stop("Data length mismatch")
    
19. 
    
20.     ord <- order(strat, y[,1], -y[,2])
    
21.     strat2 <- c(1*(diff(strat[ord])!=0), 1)
    
22. 
    
23.     xx <- .C(Csurvdiff2, as.integer(n),
    
24.                    as.integer(ngroup),
    
25.                    as.integer(nstrat),
    
26.                    as.double(rho),
    
27.                    as.double(y[ord,1]),
    
28.                    as.integer(y[ord,2]),
    
29.                    as.integer(x[ord]),
    
30.                    as.integer(strat2),
    
31.                    observed = double(ngroup*nstrat),
    
32.                    expected = double(ngroup*nstrat),
    
33.                    var.e    = double(ngroup * ngroup),
    
34.                    double(ngroup), double(n))
    
35. 
    
36.     if (nstrat==1)  list(expected = xx$expected,
    
37.                          observed = xx$observed,
    
38.                          var      = matrix(xx$var.e, ngroup, ngroup))
    
39.     else            list(expected = matrix(xx$expected, ngroup),
    
40.                          observed = matrix(xx$observed, ngroup),
    
41.                          var      = matrix(xx$var.e, ngroup, ngroup))
    
42.     }

复制代码