Rosa Oliveira wrote:
> Dear Jim,
>
> when I run your code (even the one you send me, not in my data), I get:
>
> Don't know how to automatically pick scale for object of type function.
> Defaulting to continuous
> Error in data.frame(x = c(0.1, 0.2, 0.1, 0.2, 0.1, 0.2, 0.1, 0.2, 0.1, :
> arguments imply differing number of rows: 24, 0
Well, let's agree on the data first. Using your original dataset:
my.data<-read.table(text="region sample factora factorb factorc
0.1 10 0.895 0.903 0.378
0.2 10 0.811 0.865 0.688
0.1 20 0.735 0.966 0.611
0.2 20 0.777 0.732 0.653
0.1 30 0.600 0.778 0.694
0.2 30 0.466 174.592 0.461
0.1 40 0.446 0.432 0.693
0.2 40 0.392 0.294 0.686",header=TRUE)
library(plotrix)
par(tcl=-0.1)
gap.plot(x=rep(seq(10,45,by=5),3),
y=unlist(my.data[,c("factora","factorb","factorc")]),
main="A plot of factorial mystery",
gap=c(1.1,174),ylim=c(0,175),ylab="factor score",xlab="Group",
xticlab=c(" \n0.1\n10"," \n0.2\n10"," \n0.1\n20"," \n0.2\n20",
" \n0.1\n30"," \n0.2\n30"," \n0.1\n40"," \n0.2\n40"),
ytics=c(0,0.5,1,174.59),pch=rep(1:3,each=8),col=rep(c(4,2,3),each=8))
mtext(c("Region","Sample"),side=1,at=6,line=c(0,1))
lines(seq(10,45,by=5),my.data$factora,col=4)
lines(seq(10,45,by=5),my.data$factorb[c(1:5,NA,7,8)],col=2)
lines(seq(10,45,by=5),my.data$factorc,col=3)
legend(18,1.8,c("factora","factorb","factorc"),pch=1:3,col=c(4,2,3))
This produces a plot, and I realize that it is not the one you
describe. As before, if you can let us know what is wrong with it,
maybe we can fix it.
Jim
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
