今天主要是关于向量、数据框、向量及数据框的提取介绍。
Part1
一 1 生成向量
x1 <- c(1,4,"apple")
x2 <- seq(from = 3,to = 15,by = 3)
x3 <- rnorm(10)
x4 <- rep("gene",times=3)
x1;x2;x3;x4
[1] "1" "4" "apple"
[1] 3 6 9 12 15
[1] 0.3145981 0.5347617 0.3309096 -0.1297031 -0.9895892 0.0303523 0.8205070 -0.5117744 1.4034537 -0.6964412
[1] "gene" "gene" "gene"
一 2 提取向量
# (1)根据位置
x1[3]
# [1] "apple"
#(2) 根据逻辑值判断(逻辑值个数必须等于变量长度)
x2[x2>5]
# [1] 6 9 12 15
Part2
二 1 生成数据框
X <- data.frame(
gene=paste0("gene",1:5),
val=rnorm(5,mean = 1)
)
二 2 查看数据框属性
str(X)
# 'data.frame': 5 obs. of 2 variables:
# $ gene: Factor w/ 5 levels "gene1","gene2",..: 1 2 3 4 5
# $ val : num 1.18 0.474 0.267 1.597 1.741
dim(X)
# [1] 5 2
##查看行名和列名
colnames(X)
rownames(X)
二 3 数据框导出
write.table(X,file = "huahua.txt",sep = "\t",quote = F)
二 4 变量保存
save(X,file = "huahua.txt")
二 5提取元素
X[1,] #第一行
X[,1] #第一列
X[2,2] #第二行二列
X$gene #取gene这一列
Part3
三 选修部分
a <-data.frame(case=paste0("S",1:50),values=runif(50))
plot(a$case,a$values)
`有强迫症的花花说她不喜欢反复出现a`
## 解决办法 法一
attach(a)
plot(case,values)
detach(a)
## 法二
with(a,{
plot(case,values)
x<<-summary(values) #求和并赋值给x,<<的意思是作为全局变量,也就是出了大括号仍有效。
})
x #运行完后打印x
问题 save(X,file="test.RData")这句代码如果报错X not found,是为什么,应该怎么解决?
没有X这个环境变量,或者X没有被赋值。
网友评论