要显示的数据来自实验(T1-T8代表大脑的不同部分),如下:
[[Block1]]
sum
[T1,] 6
[T2,] 6
[T3,] 4
[T4,] 5
[T5,] 8
[T6,] 9
[T7,] 8
[T8,] 6
[[Block2]]
sum
[T1,] 3
[T2,] 3
[T3,] 4
[T4,] 5
[T5,] 4
[T6,] 2
[T7,] 1
[T8,] 5
[[Block3]]
sum
[T1,] 3
[T2,] 3
[T3,] 4
[T4,] 2
[T5,] 4
[T6,] 8
[T7,] 3
[T8,] 1
[[Block4]]
sum
[T1,] 6
[T2,] 5
[T3,] 4
[T4,] 3
[T5,] 9
[T6,] 8
[T7,] 2
[T8,] 6
[[Block5]]
sum
[T1,] 8
[T2,] 3
[T3,] 4
[T4,] 5
[T5,] 7
[T6,] 6
[T7,] 2
[T8,] 2
[[Block6]]
sum
[T1,] 10
[T2,] 9
[T3,] 6
[T4,] 8
[T5,] 9
[T6,] 4
[T7,] 6
[T8,] 7
依旧..超过100个街区..
我想以下列方式可视化数据以查看整体价值 在每个地区都非常阻挡..
对于一个区块,我得到如下所示的线图:
但是对于100个块可视化相同是很乏味的。使用R..I尝试使用热图来查看它作为单个图的最佳方法是什么。我宁愿将它们视为图形。
最后它应该是这样的(我有一个粗略的图形)..我不知道如何在R中为单个图中的几个块或其他更好的方式来实现它:
答案 0 :(得分:7)
这是使用lattice xyplot
的替代方案。数据示例是现实的矩阵(100x8)。我试图删除条带以优化绘图区域。我认为结果只对获得数据的全局概念或主要趋势有用。
dat <- matrix(sample(1:10,100*8,rep=TRUE),nrow=8,
dimnames=list(paste0('T',1:8),paste0('Block',1:100)))
library(reshape2)
dat.m <- melt(dat)
xyplot(value~Var1|Var2,
data=dat.m,type=c('l','p'),
strip =FALSE,layout = c(10,10))
答案 1 :(得分:2)
在我看来,这基本上是ggplot2
的用途。这是一个重新创建您的数据,以及一个非常基本的情节。
# Recreate your data.
data<-c(6,6,4,5,8,9,8,6,3,3,4,5,4,2,1,5,3,3,4,2,4,8,3,1,6,5,4,3,9,8,2,6,8,3,4,5,7,6,2,2,10,9,6,8,9,4,6,7)
list<-split(data,rep(1:6,each=8))
names(list)<-paste0('Block',1:6)
library(ggplot2)
library(reshape2)
dat<-melt(list)[2:1]
names(dat)<-c('Block','Value')
dat$brain.section<-rep(1:8,6)
ggplot(dat,aes(x=brain.section,y=Value,group=Block)) + geom_line() + facet_grid(Block~.)
您可以对颜色和布局非常满意,但如果您不了解ggplot2
,可以使用它作为开始使用的东西。
以下是相同数据的热图如下所示
ggplot(dat,aes(x=brain.section,fill=Value,y=Block)) + geom_tile()
答案 2 :(得分:1)
这是一种或多或少匹配所需结果的替代方案。我想,由于要显示大量的块,因此比例并不重要。
## Recreate the data
my.data <- c(6,6,4,5,8,9,8,6,3,3,4,5,4,2,1,5,3,3,4,2,4,8,3,1,6,5,4,3,9,8,2,6,8,3,4,5,7,6,2,2,10,9,6,8,9,4,6,7)
n.block <- 6
n.sect <- 8
my.list <- split(my.data, rep(1:n.block, each = n.sect))
names(my.list) <- paste0("Block", 1:n.block)
sect.name <- paste0("T", 1:n.sect)
## Plot
scale.fact <- max(my.data)
plot(my.list[[1]], type = "n", axes = FALSE, ylim = c(1, n.block + 1), xlab = "", ylab = "")
for (i in seq(along = my.list)){
lines(i + my.list[[i]]/scale.fact)
}
axis(1, at = 1:n.sect, labels = sect.name, tick = TRUE)
axis(2, at = 1:n.block + sapply(my.list, function(x) x[[1]][1])/scale.fact,
labels = names(my.list), tick = TRUE, las = 1)