首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    GEO_多组数据联合分析(去除批次效应)

    )#exp2的第三个样本有些异常,可以去掉或者用normalizeBetweenArrays标准化,把它拉回正常水平。...#把第三个样本从表达矩阵里去掉exp2 = exp2[,-3] 分别获取探针注释ids1,ids2,不同的平台,都要自行修改获得探针注释。...= trans_array(exp2,ids2)处理完之后取交集,后合并s = intersect(rownames(exp1),rownames(exp2));length(s)exp1 <- exp1...[s,]exp2 <- exp2[s,]#合并表达矩阵,cbind是按列拼接的函数exp = cbind(exp1,exp2)boxplot(exp) #合并后的表达矩阵可以看到较为明显的批次效应4 分组信息需要修改和检查分组信息...两个数据合并不用改动,如果有多个数据再往后加rep("C",ncol(exp3))...batch <- c(rep("A",ncol(exp1)),rep("B",ncol(exp2)));batch

    72320
    领券