Tag Archives: Venn

4月 05

在R中绘制韦恩图的几种方法

Posted on by

  • 使用limma包

1
2
3
4
5
6
7
8
9
10
11
12
13
14

#安装limma包
#source("http://www.bioconductor.org/biocLite.R")
#biocLite("limma")
 
#载入limma包
library(limma)
hsb2 《- read.table("http://www.ats.ucla.edu/stat/R/notes/hsb2.csv", sep=',', header=T)
attach(hsb2)
hw 《- (write《=60)
hm 《- (math 《=60)
hr 《- (read 《=60)
c3 《- cbind(hw, hm, hr)
a 《- vennCounts(c3)
vennDiagram(a, include = "both", names = c("High Writing", "High Math", "High Reading"), cex = 1, counts.col = "red")

  • 使用gplots包

1
2
3
4
5
6
7
8
9
10
11
12
13

#安装gplots包
#install.packages("gplots")
 
#载入gplots包
library(gplots)
oneName 《- function() paste(sample(LETTERS,5,replace=TRUE),collapse="")
geneNames 《- replicate(1000, oneName())
GroupA 《- sample(geneNames, 400, replace=FALSE)
GroupB 《- sample(geneNames, 750, replace=FALSE)
GroupC 《- sample(geneNames, 250, replace=FALSE)
GroupD 《- sample(geneNames, 300, replace=FALSE)
input  《-list(GroupA,GroupB,GroupC,GroupD)
venn(input)

Continue reading →

5月 30

使用Perl绘制统计图

Posted on by

注意：请将代码中的“》”（中文全角）全部替换为“>”（英文半角）。

  • Bar

代码：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

#!/use/bin/perl
 
use SVG::TT::Graph::Bar;
 
my @fields        = qw(Jan Feb Mar);
my @data_sales_02 = qw(12 45 21);
 
my $graph = SVG::TT::Graph::Bar-》new(
  {
      'height' =》 '500',
      'width'  =》 '300',
      'fields' =》 @fields,
  }
);
 
$graph-》add_data(
  {
      'data'  =》 @data_sales_02,
      'title' =》 'Sales 2002',
  }
);
 
open( my $fh, '》', "bar.svg" );
select $fh;
binmode $fh;
print $graph-》burn();
close($fh);

输出：
bar
Continue reading →