蒋蒋 - 博客园

2012年12月14日

摘要： mean（trim）去掉最大值和最小值后的平均数（Trim mean）> s<-c(.25, 0.4, 1.00,2.00,3.00,4.00,5.00,8.00,12.00,50.00) > s [1] 0.25 0.40 1.00 2.00 3.00 4.00 5.00 8.00 12.00 50.00> mean(s,trim=1/10) [1] 4.425#####也就是去除s中的最大值最小值等价于> s<-c(.25,1.00,2.00,3.00,4.00,5.00,8.00,12.00) > mean(s) [1] 4.40625##### 阅读全文

posted @ 2012-12-14 12:26 蒋蒋阅读(477) 评论(0) 推荐(0) 编辑

excel快速下拉函数到指定的位置

摘要： excel中怎样快速下拉函数到指定的位置?例如：设置B1=A1*2，然后把B列往下拉时，就分别得到B2+A2*2,B3=A3*2,且自动填充数值，可是，如果行数非常多，达到几万个，则要消耗很长时间，怎样能快速的拉到最后一行，或者拉到指定的位置呢？在名称框（编辑栏最左边的框）里输入要填充公式的区域 B1:B20000回车输入公式=A1*2按CTRL+ENTER 阅读全文

posted @ 2012-12-14 12:23 蒋蒋阅读(1689) 评论(0) 推荐(0) 编辑

2012年12月13日

order与sort函数

摘要： sort函数是对数列本身进行排序，（升）order函数是对数列元素排序后在原数列中的位置进行确定。> j=scan()1: 22: 43: 54: 65: 16: 3> sort(j)[1] 1 2 3 4 5 6> order(j)[1] 5 1 6 2 3 4关系：sort(j)=j [order(j)] ##右边意义：对j中的位置进行检索，返回真值阅读全文

posted @ 2012-12-13 19:58 蒋蒋阅读(1984) 评论(0) 推荐(0) 编辑

Adjust P-values

摘要： Adjust P-values for Multiple ComparisonsDescriptionGiven a set of p-values, returns p-values adjusted using one of several methods.Usagep.adjust(p, method = p.adjust.methods, n = length(p))p.adjust.methods# c("holm", "hochberg", "hommel", "bonferroni", "B 阅读全文

posted @ 2012-12-13 19:53 蒋蒋阅读(1339) 评论(0) 推荐(0) 编辑

R语言常用函数参考

摘要：基本一、数据管理vector：向量 numeric：数值型向量 logical：逻辑型向量 character；字符型向量 list：列表 data.frame：数据框 c：连接为向量或列表 length：求长度 subset：求子集 seq，from:to，sequence：等差序列 rep：重复 NA：缺失值 NULL：空对象 sort，order，unique，rev：排序 unlist：展平列表 attr，attributes：对象属性 mode，typeof：对象存储模式与类型 names：对象的名字属性二、字符串处理character：字符型向量 nchar：字符数 substr：阅读全文

posted @ 2012-12-13 19:46 蒋蒋阅读(55927) 评论(1) 推荐(0) 编辑

cummeRbund包

摘要： cummeRbundhttp://www.bioconductor.org/packages//2.10/bioc/html/cummeRbund.htmlAnalysis, exploration, manipulation, and visualization of Cufflinks high-throughput sequencing data.Bioconductor version: Release (2.10)Allows for persistent storage, access, exploration, and manipulation of Cufflinks high 阅读全文

posted @ 2012-12-13 19:40 蒋蒋阅读(1099) 评论(0) 推荐(0) 编辑

输入复制内容&交&并

摘要： y<-read.table("clipboard",header=T)union(x, y) # 并集 iintersect(x, y) # 交集 setdiff(x, y) # 差集 setequal(x, y) # 集合是否相等阅读全文

posted @ 2012-12-13 19:38 蒋蒋阅读(165) 评论(0) 推荐(0) 编辑

R输出到文件

摘要： > sink("intersect.txt") ####建立文件名 > length(intersect(dgeseq1[,1],cuff1[,1])) ####要输出的内容 > sink() ####输出要是输出的格式不是窗口显示的形式时：write.table(inter,"inter.txt",col.names=F,row.names=F,quote=F,sep="\t")、#inter:输出的内容#inter.txt:输出的文件名阅读全文

posted @ 2012-12-13 19:36 蒋蒋阅读(3801) 评论(0) 推荐(0) 编辑

perl读取跳过第一行

摘要： open READ,"<","G:/project-xiaomai/yourfile.txt"; readline READ; # skip the first line while (<READ>) { print $_; } close READ; 阅读全文

posted @ 2012-12-13 15:36 蒋蒋阅读(2206) 评论(0) 推荐(0) 编辑

grep

摘要：几个Perl grep函数的示例1.统计匹配表达式的列表元素个数$num_apple=grep/^apple$/i,@fruits;在标量上下文里，grep返回匹配中的元素个数；在列表上下文里，grep返回匹配中的元素的一个列表。所以，上述code返回apple单词在@fruits数组中存在的个数。因为$num_apple是个标量，它强迫grep结果位于标量上下文里。2.从列表里抽取唯一元素@unique=grep{++$count{$_}<2} qw(abacddefgfhh); print"@unique\n"; 上述code运行后会返回：abcdefgh 即qw 阅读全文

posted @ 2012-12-13 15:34 蒋蒋阅读(1260) 评论(0) 推荐(0) 编辑

蒋蒋

公告