У меня есть объект экспрессии гена Seurat, который имеет слоты.
Один из слотов - @meta.data
, который является матрицей.
Я бы хотел создать столбец $ orig.ident, присвоив ему значение meta$brain.region
в качестве коэффициента. meta
- это моя оригинальная таблица метаданных.
Я делаю это для нескольких наборов данных и хотел бы сделать его обобщаемым.
Идея состоит в том, что пользователю нужно будет только ввести имя исходного объекта, и все с этого момента будет вызываться соответствующим образом.
Подсказка пользователя:
> dataset <- "path/to/gw14.RData"
> seurat.obj <- "gw14"
Затем загружается рабочее пространство, которое включает в себя объект Seurat gw14.
> load(dataset)
> seurat.obj.new <- paste0(seurat.obj, ".", 2)
Я не понимаю, почему использование get
здесь возвращает ошибку ниже:
> get(seurat.obj.new)@meta.data$orig.ident <- factor(meta$brain.region)
Error in get(seurat.obj.new)@meta.data$orig.ident = factor(meta$brain.region) :
could not find function "get<-"
В то время как его использование здесь работает, как и ожидалось:
> assign(seurat.obj.new, CreateSeuratObject(raw.data = get(seurat.obj)@raw.data,
min.cells = 0, min.genes = 0, project=age))