本人初学生信,浏览一篇帖子时,无法理解其中一段代码的意思,故向大家请教。部分代码如下:
if (!file.exists( './data/TCGA-GBM.Rdata' )) {
gzfile <- "./raw_data/TCGA-GBM.AffyU133a_log2.tsv.gz"
download.file("https://tcga.xenahubs.net/download/TCGA.GBM.sampleMap/HT_HG-U133A.gz",
destfile = gzfile)
library(R.utils)
gunzip(gzfile, remove = F)
library(data.table)
raw_data <- fread( "./raw_data/TCGA-GBM.AffyU133a_log2.tsv",
sep = ' ', header = T)
raw_data <- as.data.frame( raw_data )
raw_data[1:5, 1:6]
rownames( raw_data ) <- raw_data[, 1]
raw_data <- raw_data[, -1]
raw_data[1:5, 1:6]
raw_data <- 2^raw_data - 1
raw_data <- ceiling( raw_data )
raw_data[1:5, 1:6]
pick_row <- apply( raw_data, 1, function(x){
sum(x == 0) < 10
})
raw_data <- raw_data[pick_row, ]
dim(raw_data )
save( raw_data, file = './data/TCGA-GBM.Rdata' )
}else{
load('./data/TCGA-GBM.Rdata')
}
pick_row这行代码之中,function(x){ sum(x == 0) < 10 }做了什么操作?是何作用?希望大佬能解释一下,十分感谢