时间戳

首页 > 解决方案 > 在 RStudio 中，这行代码出现两个错误。有人可以帮我找出问题所在吗？我给了代码表格开始

r - 在 RStudio 中，这行代码出现两个错误。有人可以帮我找出问题所在吗？我给了代码表格开始

问题描述

##--> Compute Z-score (Mean expression = 0 and SD = 1 for each gene)

Z_score_transcriptomic_data <- data_rma_volcano
Z_score_transcriptomic_data  <- Z_score_transcriptomic_data [match(row.names(sig_limma_cont), row.names(Z_score_transcriptomic_data )),]
Z_score_transcriptomic_data  <- unique(Z_score_transcriptomic_data)

center_rowmeans <- function(x) {
    mean_x = rowMeans(x)
    sd_x = sd(x)
    (x - mean_x) / sd_x 
}

data_transcriptomic_HM <- cbind(center_rowmeans(Z_score_transcriptomic_data[,c(1:6)]))
    center_rowmeans(Z_score_transcriptomic_data[,c(7:12)]),
    center_rowmeans(Z_score_transcriptomic_data[,c(13:18)]),
    center_rowmeans(Z_score_transcriptomic_data[,c(19:24)]),

colnames(data_transcriptomic_HM) <- c(c("N_1", "N_2", "N_3", "N_4", "N_5", "N_6","CM_1", "CM_2", "CM_3", "CM_4", "CM_5", "CM_6", "EM_1", "EM_2", "EM_3", "EM_4", "EM_5", "EM_6", "EMRA_1", "EMRA_2", "EMRA_3", "EMRA_4", "EMRA_5", "EMRA_6"))

##--> Hierarchical clustering
clustering_hc <- hclust(as.dist(1-cor(t(data_transcriptomic_HM[,c(1:24)]))), method = "ward.D2")
clusters <- cutree(tree = clustering_hc , k = 10)

    ##--> Hierarchical clustering
    clustering_hc <- hclust(as.dist(1-cor(t(data_transcriptomic_HM[,c(1:24)]))), method = "ward.D2")
    clusters <- cutree(tree = clustering_hc , k = 10)

在 RConsole 中运行此代码后，我收到以下错误：

   Error in base::rowMeans(x, na.rm = na.rm, dims = dims, ...) : 
  'x' must be numeric
Called from: base::rowMeans(x, na.rm = na.rm, dims = dims, ...)
Browse[1]>     center_rowmeans(Z_score_transcriptomic_data[,c(7:12)]),
Error during wrapup: unexpected ',' in "    center_rowmeans(Z_score_transcriptomic_data[,c(7:12)]),"
Error: no more error handlers available (recursive errors?); invoking 'abort' restart
Browse[1]>     center_rowmeans(Z_score_transcriptomic_data[,c(13:18)]),
Error during wrapup: unexpected ',' in "    center_rowmeans(Z_score_transcriptomic_data[,c(13:18)]),"
Error: no more error handlers available (recursive errors?); invoking 'abort' restart
Browse[1]>     center_rowmeans(Z_score_transcriptomic_data[,c(19:24)]),
Error during wrapup: unexpected ',' in "    center_rowmeans(Z_score_transcriptomic_data[,c(19:24)]),"
Error: no more error handlers available (recursive errors?); invoking 'abort' restart
Browse[1]> 
> colnames(data_transcriptomic_HM) <- c(c("N_1", "N_2", "N_3", "N_4", "N_5", "N_6","CM_1", "CM_2", "CM_3", "CM_4", "CM_5", "CM_6", "EM_1", "EM_2", "EM_3", "EM_4", "EM_5", "EM_6", "EMRA_1", "EMRA_2", "EMRA_3", "EMRA_4", "EMRA_5", "EMRA_6"))
Error in colnames(data_transcriptomic_HM) <- c(c("N_1", "N_2", "N_3",  : 
  object 'data_transcriptomic_HM' not found
> 
> ##--> Hierarchical clustering
> clustering_hc <- hclust(as.dist(1-cor(t(data_transcriptomic_HM[,c(1:24)]))), method = "ward.D2")
Error in h(simpleError(msg, call)) : 
  error in evaluating the argument 'x' in selecting a method for function 'cor': error in evaluating the argument 'x' in selecting a method for function 't': object 'data_transcriptomic_HM' not found
> clusters <- cutree(tree = clustering_hc , k = 10)
Error in nrow(tree$merge) : object 'clustering_hc' not found

如果我遗漏了您需要回答的任何细节，我很乐意回答，非常适合这个新网站，所以我不太了解如何很好地提问。

标签： rbioinformaticsbioconductor

解决方案

推荐阅读