Skip to content

Instantly share code, notes, and snippets.

gset <- getGEO('GSE26253')[[1]]
pset1 <- pData(gset)
gset2 <- getGEO('GSE15460')[[1]]
gset3 <- getGEO('GSE15460')[[2]]
pset2 <- pData(gset2)
pset3 <- pData(gset3)
View(pset1)
pset1_final <- pset1[,c("pathological stage:ch1", "recurrence free survival time (month):ch1", "status (0=non-recurrence, 1=recurrence):ch1")]
View(pset1_final)
write.csv(pset1_final, file = "GSE26253PHENOFinal.csv")
gset <- getGEO('GSE26253')[[1]]
pData(gset)
gset2 <- getGEO('GSE15460')[[1]]
pData(gset2)
write.csv(pData(gset), file = "GSE26253PHENO.csv")
write.csv(pData(gset2), file = "GSE15460PHENO.csv")
##
data.matrix(gset@featureData@data$Entrez_Gene_ID)
data.matrix(gset@featureData@data$Probe_Id)
data.matrix(gset2@featureData@data$ENTREZ_GENE_ID)