exp <- exprs(gset[[1]])
时间: 2023-12-24 10:04:50 浏览: 204
这行代码是用来提取一个基因表达数据集(Gene Expression Dataset)中第一个样本的基因表达值(expression values)。假设gset是一个Gene Expression Dataset对象,其中包含多个样本的基因表达值数据。exprs()函数用于提取Gene Expression Dataset对象中的基因表达值矩阵,而[[1]]用于提取矩阵的第一列,即第一个样本的基因表达值。最终,exp将包含第一个样本的基因表达值向量。
相关问题
library(here) source(here::here("code", "libraries.R")) source(here::here("code", "funcs.R")) ## function to intersect multiple vectors intersect_all <- function(a,b,...){ Reduce(intersect, list(a,b,...)) } ## read data tpm.rsem <- read.table(here::here("data", "tpm_original.txt"), sep = "\t") counts.rsem <- read.table(here::here("data", "counts_original.txt"), sep = "\t") qc <- read.table(here::here("data", "qc_original.txt"), sep = "\t", stringsAsFactors = FALSE) # quality control mappings <- readRDS(here::here("data", "mappings.RDS")) # genes info ## scater intialize min_thresh_log_tpm <- 0.1 sceset_all <- SingleCellExperiment(assays = list(counts = as.matrix(counts.rsem), tpm = as.matrix(tpm.rsem), exprs = log2(as.matrix(tpm.rsem) + 1), expressed = log2(tpm.rsem + 1) > min_thresh_log_tpm), colData = qc, rowData = mappings) ## scater quality control找出错误,使之顺利运行
这段代码是用来加载外部库和函数的。首先使用`library(here)`来加载here库。然后使用`source()`函数来加载两个外部文件,都在当前工作目录下的code文件夹中。其中一个文件名为"libraries.R",另一个文件名为"funcs.R",这些文件包含了一些自定义的函数和库。`here::here()`函数用来指定这些文件的路径,使得它们能够被正确地加载。
sceset_all <- perCellQCMetrics(sceset_all, exprs_values = "exprs", detection_limit = min_thresh_log_tpm, cell_controls = list("regular" = which(qc$pool_H12 == 0), "pool" = which(qc$pool_H12 == 1))) if (length(which(is.na(as.numeric(sceset_all$uniquely_mapped_percent)))) > 0) sceset_all <- sceset_all[, -which(is.na(as.numeric(sceset_all$uniquely_mapped_percent)))]这段代码有错误吗
抱歉,作为 AI ,我无法判断这段代码是否有错误。这是一段 R 语言代码,可以对 sceset_all 数据进行处理,但是需要更多的上下文信息才能确切地判断是否有错误。建议您检查输入的变量是否正确定义或者寻求专业人士的帮助。
阅读全文