Nothing
`acf.fnc` <-
function(dat, group="Subject", time="Trial", x = "RT", plot=TRUE, ...) {
dat = dat[order(dat[,group], dat[,time]),]
cnt=0
civec = vector(mode="numeric")
for (s in levels(dat[,group])) {
cnt = cnt+1
tmp = dat[dat[, group] == s, ]
a = stats::acf(tmp[,x], plot=FALSE)$acf[,,1]
n = nrow(tmp)
ci = -(1/n) + 2/sqrt(n)
civec[cnt]=ci
a.dfr = data.frame(Lag=0:(length(a)-1), Acf=a, Subject=rep(s, length(a)),
ci=rep(ci, length(a)))
if (cnt == 1) {
res = a.dfr
} else {
res = rbind(res, a.dfr)
}
}
dfr = res
if (plot==TRUE) {
lattice::xyplot(Acf~Lag|Subject,type="h", data=dfr, col.line="black",
panel = function(...) {
lattice::panel.abline(h=civec[lattice::panel.number()],col.line="grey")
lattice::panel.abline(h=-civec[lattice::panel.number()],col.line="grey")
lattice::panel.xyplot(...)
},
strip = lattice::strip.custom(bg="grey90"),
par.strip.text=list(cex=0.8), ...)
} else {
return(dfr)
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.