- 阅读权限
- 255
- 威望
- 1 级
- 论坛币
- 132048 个
- 通用积分
- 9219.7014
- 学术水平
- 925 点
- 热心指数
- 1073 点
- 信用等级
- 703 点
- 经验
- 130943 点
- 帖子
- 3353
- 精华
- 1
- 在线时间
- 3498 小时
- 注册时间
- 2009-7-18
- 最后登录
- 2025-5-18
|
大概是这个样子
- x <- data.frame(v1 = rnorm(10), v2 = rnorm(10), v3 = rnorm(10), v4 = rnorm(10), v5 = rnorm(10))
- n <- ncol(x)
- y <- matrix(NA, n, n, byrow = TRUE)
- y[upper.tri(y)] <- 0
- y[lower.tri(y, diag = TRUE)] <- 1:(n*(n+1)/2)
- y <- t(y)
- pdf("out.pdf")
- par(mar=c(2, 2, 2, 2) + 0.1)
- nf <- layout(y)
- for (i in 1:n) {
- for (j in i:n) {
- if (i == j) {
- plot(1, xlim = c(0, 1), ylim = c(0, 1), type = "n", xaxt = 'n', yaxt = 'n', ann = FALSE)
- text(0.5, 0.5, colnames(x)[i], cex = 2)
- } else {
- plot(x[, i], x[, j], xaxt = 'n', yaxt = 'n', ann = FALSE)
- if (i == 1) axis(side = 3, at = seq(-2, 2, 0.5))
- if (j == n) axis(side = 4, at = seq(-2, 2, 0.5))
- }
- }
- }
- dev.off()
复制代码
|
-
总评分: 学术水平 + 1
热心指数 + 1
查看全部评分
|