In qqmath, how would one go about having 'distribution' change with
panel.number? I've tried
set.seed(1)
mydata <- data.frame(ind = factor(rep(2:4, each = 100)))
mydata$val <- rt(300, df=rep(2:4, each = 100))
plot<-qqmath(~ val | ind,
layout=c(3,1),
data = mydata,
prepanel = function(x, distribution, ...) {
mydist<-function(p) qt(p, df = panel.number() + 1)
prepanel.qqmathline(x, distribution=mydist,...)
},
panel = function(x, distribution, ...) {
mydist<-function(p) qt(p, df = panel.number() + 1)
panel.qqmathline(x, distribution=mydist,...)
panel.qqmath(x, distribution=mydist,...)
})
But no luck.
Ben
______________________________________________
R-help@stat.math.ethz.ch 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.
