Hi R users:
I have de following data frame (called "Sx")
Descripcion Nitratos
Cont85g 72.40
Cont85g 100.50
Cont85g 138.30
Cont80g 178.33
Cont80g 79.01
Cont80g 74.16
Cont75g 23.70
Cont75g 15.80
Cont75g 16.20
Patron80g 88.93
Patron80g 113.01
Patron80g 86.53
If I run this code without the "pdf" device
it works fine on the screen, but when I
use the "pdf" device it does not show in
the "plot1.pdf" file what I got on the screen.
What am I doing wrong?
Thank you for your help.
library(multcomp)
Sx<-data.frame(Descripcion=
factor(c("Cont85g","Cont85g","Cont85g",
"Cont80g","Cont80g","Cont80g",
"Cont75g","Cont75g","Cont75g",
"Patron80g","Patron80g","Patron80g")),
Nitratos=c(72.40,100.50, 138.30,
178.33,79.01,74.16,
23.70, 15.80,16.20,
88.93,113.01,86.53))
pdf("plot1.pdf")
m1<-aov(Nitratos~Descripcion-1,data=Sx)
vect1<-table(Sx$Descripcion)
K<-contrMat(vect1,base=4)
dnk<-glht(m1,linfct=K)
summary(dnk)
old.par<-par(no.readonly = TRUE)
par(mai=c(1,2,1.25,1),mgp=c(3,1,0))
print(plot(dnk,las=1,xlab=""))
print(abline(v=0,lty=2))
par(old.par)
dev.off()
______________________________________________
R-help@r-project.org mailing list
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.
