rm(list = ls())
setwd("D:/IH")
getwd()
library( "GEOquery" )
g989 =read.csv("34989.csv")
m136 =read.csv("69136.csv")
g=grep("[linc]",g989$GeneName)
str(g)
gglinc989 = g989[g,]
write.csv(mmm136,"mmm136.csv")
write.csv(gglinc989,"gglinc989.csv")
index1<-duplicated(mm136$miRNA_ID)
levels(index1)
mm136<-mm136[!index1,]
index2<-duplicated(gg989$GeneName)
levels(index2)
gg989<-gg989[!index2,]
mm136=read.csv("mm136.csv")
gg989=read.csv("gg989.csv")
mmm136 <- mm136[mm136$P.Val < 0.05, ]
ggg989 <- gg989[gg989$P.Val < 0.05, ]
library(dplyr)
mmm136 <- subset(mm136, mm136$adj.P.Val < 0.05)