zhoujj zhoujj2013
zhoujj2013
/ density.r
Created
July 21, 2015 07:27
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| args <- commandArgs(TRUE) | |
| prefix = args[1] | |
| #print(args[c(2:length(args))]) | |
| pdf(paste(prefix,".density.pdf",sep="")) | |
| color = rainbow(length(args) - 1) | |
| for (i in 2:length(args)){ | |
| dat = read.table(args[i],header=F) |
zhoujj2013
/ jmzeng_WGCNA.r
Created
September 17, 2017 03:43
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| setwd('WGCNA/') | |
| # 56 breast cancer cell lines were profiled to identify patterns of gene expression associated with subtype and response to therapeutic compounds. | |
| if(F){ | |
| ## https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE48213 | |
| #wget -c ftp://ftp.ncbi.nlm.nih.gov/geo/series/GSE48nnn/GSE48213/suppl/GSE48213_RAW.tar | |
| #tar -xf GSE48213_RAW.tar | |
| #gzip -d *.gz | |
| ## 首先在GSE48213_RAW目录里面生成tmp.txt文件，使用shell脚本 | |
| # awk '{print FILENAME"\t"$0}' * |grep -v EnsEMBL_Gene_ID >tmp.txt | |
| ## 然后把tmp.txt导入R语言里面用reshape2处理即可 |
OlderNewer