1

我希望制作一个图表，显示对许多协议声明的响应。该图表应允许对不同组的受访者和陈述项目进行比较。

我基本上借鉴了 Kastellec & Leoni 提供的 R 代码（图 5；http ://tables2graphs.com/doku.php?id=03_descriptive_statistics#figure_5 ）。

与它们相比，我希望 x 轴的范围为 -5 到 5，表格的格式应为 2x2。

此代码应生成类似于我使用的示例数据：

mydata<-expand.grid(
col1=c('item1', 'item2', 'item3', 'item4'),
col2=c('0', '1', '3', '4'),
col3=c('T1', 'T2', 'C1', 'C2'))
mydata$col4=sapply(rnorm(64,0,1), function(x) {round(x,2)})

注： - col1：陈述项（K&L数据集中的“变量”） - col2：答案类别（K&L中的“句号”） - col3：受访者组（K&L中的“立法机关”） - col4：比例（“比例” " 在 K&L)

这是代码：

library(lattice)
    library(car)

    ltheme <- canonical.theme(color = FALSE) 
    ltheme$strip.background$col <- "lightgrey"  
    lattice.options(default.theme = ltheme)     

    mydata$col2<-factor(mydata$col2,
    levels=c(0,1,3,4),
    labels=c("strongly disagree", "disagree", "agree", "strongly agree"), ordered=TRUE)

    mydata$col3<-factor(mydata$col3,
    levels=c("T1", "C1", "T2", "C2"), ordered=TRUE)
    levels(mydata$col3)<-c("treatment group 1", "control group 1", "treatment group 2",  "control group 2")

    mydata$col1<-factor(mydata$col1,
    levels=c("item1", "item2", "item3", "item4"),
    labels=c("item 4", "item 3", "item 2", "item 1"), ordered=TRUE)

    prop.vec<-c(-5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5)

    plot<-dotplot(mydata$col1~mydata$col4|mydata$col3, xlab="levels of agreement",
    data=mydata,
    groups=mydata$col2,
    layout=c(2,2),
    scales=list(cex=0.65,
    x=list(at=prop.vec),
    alternating=3),
    par.strip.text=list(lines=2.5,cex=0.65),
    panel=function(...){
       panel.abline(v=prop.vec, col="lightgrey")
       panel.abline(h=1:11, col="lightgrey", lty=2)
       panel.xyplot(...)},
    as.table=TRUE,
    par.settings=simpleTheme(pch=c(19,1,2,17), cex=0.7), 
    auto.key=list(space="bottom", column=4, cex=0.65)
    )

    trellis.device(file="figure.pdf", device="pdf", color=FALSE, width=6, height=8)

    print(plot, newpage=FALSE)              

    dev.off()

在此处输入图像描述

4

1 回答 1

0

问题已解决 - 我更正了代码，以便它生成一个漂亮的图，显示对声明协议的响应，如上所示。如果您愿意，请随时进一步改进它。非常感谢您的帮助！

于 2012-08-20T09:26:33.960 回答