## ----------------------------------------------------------------------------- 库(TSCAN)数据(lpsdata) procdata < -预处理(lpsdata ) ## ----------------------------------------------------------------------------- lpsmclust < - exprmclust (procdata ) ## ----------------------------------------------------------------------------- plotmclust (lpsmclust ) ## ----------------------------------------------------------------------------- lpsorder < - TSCANorder lpsmclust lpsorder # #----------------------------------------------------------------------------- diffval < - difftest (procdata lpsorder) #选择下的差异表达基因qvlue截止0.05头(row.names (diffval) [diffval qval < 0.05美元 ]) ## ----------------------------------------------------------------------------- STAT2expr < - log2 (lpsdata(“STAT2”)+ 1)singlegeneplot (STAT2expr,TSCANorder(lpsmclust,flip=TRUE,orderonly=FALSE)) ## ----------------------------------------------------------------------------- subpopulation <- data.frame(cell = colnames(procdata)， sub = ifelse(grepl(" un"，colnames(procdata))，0,1)， stringsAsFactors =FALSE) #比较有或没有标记基因信息的排序order1 <- TSCANorder(lpsmclust) order2 <- TSCANorder(lpsmclust，c(1、2、3)< -订单列表(order1 order2 ) ## ----------------------------------------------------------------------------- orderscore(分组人口、订单)