分组语句targets <- data括号里的没看懂, Label = c(paste("Con", 1:4, sep=""), paste("DHT", 1:3, sep="") 这句话,R怎么能辨别 lane1 lane2 lane3 lane4 lane5 lane6 lane8 len这几列,哪个是Con哪个是DHT呢??
谢谢老师解答!
代码如下:
# 1. 载入数据 读取read count数
data <- read.delim("pnas_expression.txt", row.names=1, stringsAsFactors=FALSE);
head(data);
#输出
# lane1 lane2 lane3 lane4 lane5 lane6 lane8 len
# ENSG00000215696 0 0 0 0 0 0 0 330
# ENSG00000215700 0 0 0 0 0 0 0 2370
# ENSG00000215699 0 0 0 0 0 0 0 1842
# ENSG00000215784 0 0 0 0 0 0 0 2393
# ENSG00000212914 0 0 0 0 0 0 0 384
# ENSG00000212042 0 0 0 0 0 0 0 92
dim(data);
# [1] 37435 8
#2. 构建分组变量
#分为 Control组和DHT组 分别为4个和3个重复
targets <- data.frame(Lane = c(1:6,8), Treatment = c(rep("Control",4),rep("DHT",3)),
Label = c(paste("Con", 1:4, sep=""), paste("DHT", 1:3, sep="")));