cur*_*wed 2 plot r bar-chart lattice reshape
我在表(从 csv 导入)中有一堆数据,格式如下:
date classes score
9/1/11 french 34
9/1/11 english 34
9/1/11 french 34
9/1/11 spanish 34
9/2/11 french 34
9/2/11 english 34
9/3/11 spanish 34
9/3/11 spanish 34
9/5/11 spanish 34
9/5/11 english 34
9/5/11 french 34
9/5/11 english 34
Run Code Online (Sandbox Code Playgroud)
忽略分数列,这不重要。
我需要根据日期统计参加英语或西班牙语或法语课程的学生总数,即。我需要首先按日期对其进行分组,然后根据语言将每一天划分为更多的块,并将其绘制为堆叠条形图,如下所示。每个条形代表一个日期,条形的每个横截面代表一种语言。
一旦我以矩阵形式获取数据,其中每一行代表一个日期,每列代表一个属性(或语言),我就想出了如何做到这一点。所以我假设数据是 csv 格式的:
ie french english spanish
9/1/11 2 1 1
9/2/11 1 1 0
9/3/11 0 0 2
9/5/11 1 2 1
Run Code Online (Sandbox Code Playgroud)
那么我可以这样做:
directory<-"C:\\test\\language.csv"
ourdata6<-read.csv(directory)
language<-as.matrix(ourdata6)
barchart(prop.table(language), horizontal=FALSE, auto.key = list(space='right',cex=.5,border=T,points=F, lines=F,lwd=5,text=c('french','spanish','enligsh'),cex=.6), main = list(label="Distribution of classes 10",cex=2.5), ylab = list(", cex=1.7),xlab.top=list("testing",cex=1.2))
Run Code Online (Sandbox Code Playgroud)
挑战是将数据从原始格式转换为我需要的格式。
我试过
a<-count(language, c("date", "classes"))
Run Code Online (Sandbox Code Playgroud)
它给了我按两者排序的计数,但它是垂直形式的
ie
9/1/11 french 2
9/1/11 english 1
9/1/11 spanish 1
etc...
Run Code Online (Sandbox Code Playgroud)
我需要旋转它,以便它成为每个日期的一行。此外,如果其中一些可能为零,那么我需要它们的占位符,即。第一列必须对应于法语,第二列必须对应于我当前的设置工作的英语。
关于如何做到这一点的任何想法,或者我使用 matrix + prop.table 的方法是否正确?有没有更简单的方法来做到这一点?
假设您的数据位于名为 的数据框中df,您可以借助dplyr和tidyr包来实现:
library(dplyr)
library(tidyr)
wide <- df %>% select(date,classes) %>%
group_by(date,classes) %>%
summarise(n=n()) %>% # as @akrun said, you can also use tally()
spread(classes, n, fill=0)
Run Code Online (Sandbox Code Playgroud)
使用您提供的示例数据,这会产生以下数据框:
date english french spanish
9/1/11 1 2 1
9/2/11 1 1 0
9/3/11 0 0 2
9/5/11 2 1 1
Run Code Online (Sandbox Code Playgroud)
现在您可以使用以下内容制作lattice情节:
barchart(date ~ english + french + spanish, data=wide, stack = TRUE,
main = list(label="Distribution of language classes",cex=1.6),
xlab = list("Number of classes", cex=1.1),
ylab = list("Date", cex=1.1),
auto.key = list(space='right',cex=1.2,text=c('Enligsh','French','Spanish')))
Run Code Online (Sandbox Code Playgroud)
这给出了以下情节:

编辑:您也可以使用ggplot2,而不是使用点阵图,这(至少在我看来)更容易理解。一个例子:
# convert the wide dataframe to a long one
long <- wide %>% gather(class, n, -date)
# load ggplot2
library(ggplot2)
# create the plot
ggplot(long, aes(date, n, fill=class)) +
geom_bar(stat="identity", position="stack") +
coord_flip() +
theme_bw() +
theme(axis.title=element_blank(), axis.text=element_text(size=12))
Run Code Online (Sandbox Code Playgroud)
这使:

| 归档时间: |
|
| 查看次数: |
3957 次 |
| 最近记录: |