摘要：
library(ggplot2)library(grid)grid.newpage()pushViewport(viewport(layout = grid.layout(1,2)))vplayout = function(x, y) viewport(layout.pos.row = x, lay 阅读全文
2017年3月6日
摘要：
###manhattan plots in qqplot2library(ggplot2)setwd("~/ncbi/zm/XPCLR/")read.table("LW.gene.xpclr.position.txt",header=F,sep=',')->ILhead(IL)#V1 V2 V3 V 阅读全文
2017年1月13日
摘要：
import sysinfile = sys.argv[1]outfile = sys.argv[2]gff_list = []fh = open(infile)for x in fh: if x[0] == '#': continue y = x.strip().split('\t') gff_l 阅读全文
2016年12月4日
2016年12月3日
摘要：
#using R #Prepare simulation data: a <- 1:5 b <- a^2 #predicting model model <- lm(b~poly(a,2,raw=T)) predicted <- predict(model,data.frame(x=a),inter 阅读全文