library(hgu133plus2.db)
library(GEOquery)
a=getGEO(‘GSE103091’,destdir = ‘.’)
metadata=pData(a[[1]])[,c(2,11,12,13,14,15)]
datTraits = data.frame(gsm=metadata[,1],
age =trimws(sapply(as.character(metadata
c
h
a
r
a
c
t
e
r
i
s
t
i
c
s
c
h
1.1
)
,
f
u
n
c
t
i
o
n
(
x
)
s
t
r
s
p
l
i
t
(
x
,
"
:
"
)
[
[
1
]
]
[
2
]
)
)
,
m
e
t
a
=
t
r
i
m
w
s
(
s
a
p
p
l
y
(
a
s
.
c
h
a
r
a
c
t
e
r
(
m
e
t
a
d
a
t
a
characteristics_ch1.1),function(x) strsplit(x,":")[[1]][2])), meta=trimws(sapply(as.character(metadata
characteristicsch1.1),function(x)strsplit(x,":")[[1]][2])),meta=trimws(sapply(as.character(metadatacharacteristics_ch1.2),function(x) strsplit(x,"😊[[1]][2])),
mfs = trimws(sapply(as.character(metadata
c
h
a
r
a
c
t
e
r
i
s
t
i
c
s
c
h
1.3
)
,
f
u
n
c
t
i
o
n
(
x
)
s
t
r
s
p
l
i
t
(
x
,
"
:
"
)
[
[
1
]
]
[
2
]
)
)
,
d
e
a
t
h
=
t
r
i
m
w
s
(
s
a
p
p
l
y
(
a
s
.
c
h
a
r
a
c
t
e
r
(
m
e
t
a
d
a
t
a
characteristics_ch1.3),function(x) strsplit(x,":")[[1]][2])), death=trimws(sapply(as.character(metadata
characteristicsch1.3),function(x)strsplit(x,":")[[1]][2])),death=trimws(sapply(as.character(metadatacharacteristics_ch1.4),function(x) strsplit(x,"😊[[1]][2])),
os=trimws(sapply(as.character(metadataKaTeX parse error: Undefined control sequence: \t at position 194: …trix.txt",sep="\̲t̲",header=T,chec…probe_id <- rownames(exprSet)
exprSet_symbol <- merge(probe2symbol, exprSet, by = “probe_id”)
dim(exprSet) #查看合并前的探针数目
dim(exprSet_symbol) #查看合并后的探针数目。可以见到除去了部分没有注释的探针。
exprSet_symbolKaTeX parse error: Expected 'EOF', got '#' at position 18: …obe_id <- NULL #̲删掉不要的probe_id e…symbol), FUN = median)
rownames(exprSet_symbol) <- exprSet_symbol[,1]
exprSet_symbol <- exprSet_symbol[,-1]
mydata <- exprSet_symbol
save(mydata,datTraits,file = ‘GSE103091-wgcna-input.RData’)
芯片注释
最新推荐文章于 2023-09-27 21:13:52 发布