今天主要是關(guān)于向量、數(shù)據(jù)框亿絮、向量及數(shù)據(jù)框的提取介紹耕渴。
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)根據(jù)位置
x1[3]
# [1] "apple"
#(2) 根據(jù)邏輯值判斷(邏輯值個數(shù)必須等于變量長度)
x2[x2>5]
# [1] 6 9 12 15
Part2
二 1 生成數(shù)據(jù)框
X <- data.frame(
gene=paste0("gene",1:5),
val=rnorm(5,mean = 1)
)
二 2 查看數(shù)據(jù)框?qū)傩?/p>
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 數(shù)據(jù)框?qū)С?/p>
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)
`有強迫癥的花花說她不喜歡反復(fù)出現(xiàn)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副硅,是為什么姥宝,應(yīng)該怎么解決?
沒有X這個環(huán)境變量恐疲,或者X沒有被賦值腊满。