摘要: @ARGV既然以@开头,标明这是一个list,你当做数组看就行了。含义是包含了程序从命令行得到的所有参数。比如运行a.pl 10 20 30 a b c则,$ARGV[0]是10$ARGV[3]是a 阅读全文
posted @ 2013-01-15 10:04 蒋蒋 阅读(604) 评论(0) 推荐(0) 编辑
摘要: 1. windows安装(在命令提示窗口输入以下自动下载)ppm install Config::Std 阅读全文
posted @ 2013-01-10 11:12 蒋蒋 阅读(226) 评论(0) 推荐(0) 编辑
摘要: Q-Q图是一种散点图,对应于正态分布的Q-Q图,就是由标准正态分布的分位数为横坐标,样本值为纵 坐标的散点图. 要利用QQ图鉴别样本数据是否近似于正态分布,只需看QQ图上的点是否近似地在一条直线附近,而且该直线的斜率为标准差,截距为均值. 用QQ图还可获得样本偏度和峰度的粗略信息. Q-Q图可以用于检验数据的分布,所不同的是,Q-Q图是用变量数据分布的分位数与所指定分布的分位数之间的关系曲线来进行检验的。P-P图和Q-Q图的用途完全相同,只是检验方法存在差异 由于P-P图和Q-Q图的用途完全相同,只是检验方法存在差异。要利用QQ图鉴别样本数据是否近似于正态分布,只需看QQ图上的点是否近似... 阅读全文
posted @ 2013-01-08 17:05 蒋蒋 阅读(1229) 评论(0) 推荐(0) 编辑
摘要: 摘自:吴喜之:《非参数统计》(第二版),中国统计出版社,2006年10月:P164-1651、ks.test() 例如零假设为N(15,0.2),则ks.test(x,"pnorm",15,0.2)。如果不是正态分布,还可以选"pexp", "pgamma"等。2、shapiro.test() 可以进行关于正态分布的Shapiro-Wilk检验。3、nortest包 lillie.test()可以实行更精确的Kolmogorov-Smirnov检验。 ad.test()进行Anderson-Darling正态性检验。 cvm.test 阅读全文
posted @ 2013-01-08 13:07 蒋蒋 阅读(22124) 评论(0) 推荐(0) 编辑
摘要: Linux多数发行版自带的sort程序,非常强大,在此只说多字段排序 sort 有个参数-k,可以指定字段,有比较复杂的语法,不在文本范围内。 一下为一段数据(从基因中得到,仅仅作为demo),文件名为 data chr13 3008566 3008677 chr9 3024384 3024515 chr19 3157071 3157172 chr5 3236386 3236476 chr13 3041044 3041191 chr12 3045343 3045532 chr6 3087308 3087625 chr5 3109870 3110091 chr9 3115454 3115531 阅读全文
posted @ 2013-01-07 13:36 蒋蒋 阅读(1055) 评论(0) 推荐(0) 编辑
摘要: 例:d <- data.frame(name=c("zhao","qian","sun","li"),weight=c(62,58,79,60),height=c(178,169,180,173))x <- d$namey1 <- d$weighty2 <- d$height用包plotrix#png("2y_1.png")library(plotrix)xpos <- 1:4twoord.plot(xpos,y1,xpos,y2,xlim=c(0,5),lylim 阅读全文
posted @ 2013-01-07 10:51 蒋蒋 阅读(896) 评论(0) 推荐(0) 编辑
摘要: slices <- c(10, 12, 4, 16, 8)lbls <- c("US", "UK", "Australia", "Germany", "France")pct <- round(slices/sum(slices)*100)lbls <- paste(lbls, pct) # add percents to labelslbls <- paste(lbls,"%",sep="") # ad % to label 阅读全文
posted @ 2013-01-07 10:51 蒋蒋 阅读(1780) 评论(0) 推荐(0) 编辑
摘要: for i in `find . -name "countlist.txt.num"`;do echo $i&&cat $i&& echo;done 阅读全文
posted @ 2013-01-04 10:05 蒋蒋 阅读(461) 评论(0) 推荐(0) 编辑
摘要: 如果是有进程名为abcd,abcde等有包含"abc"字符的,那么判断将会不准备,如果精确的匹配到abc呢,当时不知道,只能很诚实的回答,不清楚其实答案很简单,用grep –w "abc" 或者是grep "\<abc\>"都可以实现-w, --word-regexp 强制 PATTERN 仅完全匹配字词 阅读全文
posted @ 2012-12-31 10:13 蒋蒋 阅读(87658) 评论(0) 推荐(2) 编辑
摘要: $0代表程序名,包含完整路径 阅读全文
posted @ 2012-12-28 09:37 蒋蒋 阅读(199) 评论(0) 推荐(0) 编辑