03 2017 档案

R 画structure图
摘要:id percent k1_B04_WL-1.fs_1 0.021 k31_B04_WL-1.fs_1 0.624 k21_B04_WL-1.fs_1 0.355 k1 K=3 数据输入如上面所示: library(ggplot2) ggplot(b,aes(id,percent,fill=k))+ 阅读全文

posted @ 2017-03-29 16:27 yangyzh 阅读(879) 评论(0) 推荐(0)

perl ExtUtils::Manifest
摘要:Can't locate ExtUtils/Manifest.pm in @INC (you may need to install the ExtUtils::Manifest module) yum install perl-devel 阅读全文

posted @ 2017-03-27 12:58 yangyzh 阅读(499) 评论(0) 推荐(0)

解决libstdc++.so.5问题
摘要:./bowtie2-buildbowtie2-build: error while loading shared libraries: libtinfo.so.5: cannot open shared object file: No such file or directory 系统为 fedor 阅读全文

posted @ 2017-03-23 12:05 yangyzh 阅读(2288) 评论(0) 推荐(0)

R ggplot2 线性回归
摘要:摘自 http://f.dataguru.cn/thread-278300-1-1.html library(ggplot2) x=1:10y=rnorm(10)a=data.frame(x= x, y= y) # plot(y~x,data=a) lm_eqn = function(df){m = 阅读全文

posted @ 2017-03-09 19:19 yangyzh 阅读(1693) 评论(0) 推荐(0)

perl 计算方差中值平均数 Statistics::Descriptive;
摘要:http://search.cpan.org/~shlomif/Statistics-Descriptive-3.0612/lib/Statistics/Descriptive.pmuse Statistics::Descriptive;my @a=(100,88,77,18,79,96,57,88 阅读全文

posted @ 2017-03-09 16:14 yangyzh 阅读(2085) 评论(0) 推荐(0)

R ggplot2 翻转坐标
摘要:p <- ggplot(mpg, aes(class, hwy)) p + geom_boxplot() p + geom_boxplot() + coord_flip() 阅读全文

posted @ 2017-03-09 15:51 yangyzh 阅读(967) 评论(0) 推荐(0)

R ggplot2 改变颜色
摘要:p<-ggplot(iris,aes(Petal.Length,Petal.Width,color=Species))+geom_point()cols=c("red","green","blue") scale_colour_manual和scale_fill_manual用法相同更改颜色和leg 阅读全文

posted @ 2017-03-08 21:56 yangyzh 阅读(933) 评论(0) 推荐(0)

SNP密度分布模式
摘要:1. window=100k,step=2k 统计每个window的snp密度,然后用mixtools的normalmixEM(两个组分的混合模型)统计snp的分布模式。 R command: library(mixtools) SNPdensity=read.table("snp.density. 阅读全文

posted @ 2017-03-08 21:47 yangyzh 阅读(1096) 评论(0) 推荐(0)