爲什麼加載時的函數沒有給出任何輸出？

Question

這裏我已經定義了功能CpG.T.Test，在裏面我運行了t -test。但最後，t和p都不存在。

CpG.T.Test <- function(betam, pheno){ 
    R = 50000 
    t = numeric(R) 
    p = numeric(R) 

    for (i in 1:R){ 
    beta.v <- betam[i, ] 
    pheno.v <- pheno 

    test <- t.test (beta.v ~ pheno.v) 
    t[i] = test$statistic 
    p[i] = test$p.value 
    } 
} 

Answer 1

因爲你需要返回2個變量，你可以將它們組合成一個單一的數據幀，並按如下退回它

CpG.T.Test <- function(betam, pheno){ 
    R = 50000 
    t = numeric(R) 
    p = numeric(R) 

    for (i in 1:R){ 
    beta.v <- betam[i, ] 
    pheno.v <- pheno 

    test <- t.test (beta.v ~ pheno.v) 
    t[i] = test$statistic 
    p[i] = test$p.value 
    } 

    cbind(t,p) # the last variable used in the scope is usually returned. 
    ## OR 
    ## return cbind(t,p) # or you can explicitly return the variable 
}