R语言-工作目录

Python015

R语言-工作目录,第1张

新建方法

1.进入Rsudio界面就选择打开file-New Project-New Directory-New Project-设置名称,保存文件

2..建立一个文件夹,复制之前的.Rproj,修改自己想要的名字,进入Rsudio界面就选择打开file-open Project,选择刚刚新建的.Rproj

最最重点是后面所有文件读取,需要读取的文件放到工作目录,再用需要的函数可以直接读取,默认为工作目录下操作,我要用ReadAffy读取CEL文件,就把CEL文件放入我的工作目录,然后直接读取就可以了,见下图

4.脚本、数据、图片的默认保存位置,也是文件读取的默认位置—工作目录

最后谢谢曾老师和小洁老师耐心指导,菜鸟上线,很容易一学就忘,但是不要紧,办法总比困难多!

杀杀

前两天一个同学问了一个批量读数据(文件)的问题

感觉现在分析测序数据的话,在R中一次性读取多个文件的功能还是挺需要用到的

决定总结一下下

总不能写n行read.csv…balabala…去读取,那太麻烦了

于是决定使用批量读数据

首先我们要获取这个文件夹里所有的文件名字以便读取

接下来我们需要构建读文件的路径

ps: 在做这种批量操作时还是有必要先花时间测试一个样本是否能够成功,然后再去做批量的循环,这样也能方便自己找bug

因此我们先试着读取一个文件

这里因为例子中的文件格式比较特殊,所以使用的读取函数是read.delim,如果是.csv或.txt等格式,也可以替换成read.csv, read.table等函数。

好了测试通过,接下来开始批量读文件

需要注意的是,合并文件的时候,需要注意你是要取所有文件中行名的交集还是并集,或者最后你是否有一个统一的行名来规整所有文件。(当然如果你确定每个文件的行名都是相同的,就可以直接cbind),我遇到过需要取交集的情况,可以写一个循环对每个文件的行名做交集,用最终的交集文件去match出所有你要的行然后合并。

但是这个例子中,有的文件的行名是1-4000,有的是5000-8000,存在非常大的差异,而需求是最后整合成一个行名为0-8000的矩阵,因此我先建立一个0-8000的列作为行名,把每个文件和这列进行一个match操作,然后未match到的填充为0。

其实没有很难的部分,主要是需要读取文件夹中的所有文件名,然后循环读取就行了。