初学R，练个手

> data <- outer(1:6, 1:6, "+")
> data
[,1] [,2] [,3] [,4] [,5] [,6]
[1,]    2    3    4    5    6    7
[2,]    3    4    5    6    7    8
[3,]    4    5    6    7    8    9
[4,]    5    6    7    8    9   10
[5,]    6    7    8    9   10   11
[6,]    7    8    9   10   11   12
> gailv <- function(x) length(x)/length(data)
> fc <- factor(data)
> tapply(data, fc, gailv)
2          3          4          5          6          7          8
0.02777778 0.05555556 0.08333333 0.11111111 0.13888889 0.16666667 0.13888889
9         10         11         12
0.11111111 0.08333333 0.05555556 0.02777778

> fr <- table(data)
> fr
data
2  3  4  5  6  7  8  9 10 11 12
1  2  3  4  5  6  5  4  3  2  1
> plot(as.numeric(names(fr)), fr, type="h", xlab="Determinant", ylab="Frequency")