子集在R中不起作用
我刚开始使用R进行一些基本分析。我试图使用子集在R中不起作用,r,subset,R,Subset,我刚开始使用R进行一些基本分析。我试图使用子集从我的一个txt文件中提取一些数据。 我的数据文件的头部如下所示 head(gps2) logFC Pval Gene 1 -0,155887003 0,872819236 Cidec 2 -2,853534804 0,158917555 Hsd3b5 3 1,031753975 0,120636931 Cyp2a22 4 0,323919709 0,292010361 Osbpl3 5 -2,842
子集
从我的一个txt文件中提取一些数据。
我的数据文件的头部如下所示
head(gps2)
logFC Pval Gene
1 -0,155887003 0,872819236 Cidec
2 -2,853534804 0,158917555 Hsd3b5
3 1,031753975 0,120636931 Cyp2a22
4 0,323919709 0,292010361 Osbpl3
5 -2,842982917 0,039520301 Elovl3
6 -2,223250518 6,1759E-06 Sult5a1
我试图用Pval<0.05
和logFC<-0.5
得到子集。
我编写了如下命令:
genelist=read.table('data.txt',header=T,stringsAsFactors = F)
genelist_down=subset(genelist,logFC< -0.5&Pval<0.05)
head(genelist_down)
logFC Pval Gene
1 -0,155887003 0,872819236 Cidec
8 -0,033531014 0,933301176 Stap1
20 -0,034793814 0,965019851 Snora61
24 -0,291441914 0,63454618 Ar
27 -0,118891524 0,810291984 Ddit4
28 -0,090525284 0,8405535 1810011O10Rik
很明显,它不起作用。。。
有人能帮我解决这个问题吗…因为您的文件格式是
,
表示一个小数点,您的列都被读取为字符,因此子集给出的结果不正确。将dec=“,”
添加到您的read.table()
调用中,然后再次尝试子集()
genelist <- read.table(
"data.txt",
header = TRUE,
stringsAsFactors = FALSE,
dec = ","
)
subset(genelist, logFC < -0.5 & Pval < 0.05)
# logFC Pval Gene
# 5 -2.842983 3.95203e-02 Elovl3
# 6 -2.223251 6.17590e-06 Sult5a1
genelist由于您的文件格式为,
表示小数点,因此您的列都被读取为字符,因此子集给出了错误的结果。将dec=“,”
添加到您的read.table()
调用中,然后再次尝试子集()
genelist <- read.table(
"data.txt",
header = TRUE,
stringsAsFactors = FALSE,
dec = ","
)
subset(genelist, logFC < -0.5 & Pval < 0.05)
# logFC Pval Gene
# 5 -2.842983 3.95203e-02 Elovl3
# 6 -2.223251 6.17590e-06 Sult5a1
基因列表