Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/69.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
子集在R中不起作用_R_Subset - Fatal编程技术网

子集在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

我刚开始使用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,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
基因列表