plot.ssanova: A Plot Method for ssanova and ssanvoa0 Objects Tailored to...

In STAR: Spike Train Analysis with R

Description

Plot a ssanova or a ssanova0 object.

Usage

## S3 method for class 'ssanova'
plot(x, y, include, ask = FALSE, ncol = 2, nrow = 3, ...)
## S3 method for class 'ssanova0'
plot(x, y, include, ask = FALSE, ncol = 2, nrow = 3, ...)

Arguments

x	a ssanova or a ssanova0 object.
y	not used, only included for compatibility with generic method.
include	a character string with the model terms one wants to plot. If missing all terms are plotted.
ask	a logical. If TRUE terms are plotted (on a common y scale) one after the other and the user is invited to hit the enter key to generate the next plot. If FALSE (default) all terms are drawn on a suitable number of X11 devices. The number of terms on each device is controlled by arguments ncol and nrow.
ncol	the number of columns of the display matrix used on each device when ask is set to FALSE.
nrow	the number of rows of the display matrix used on each device when ask is set to FALSE.
...	not used only there for method definition compatibility.

Value

Nothing returned. The method is used for its side effect, plots are generated.

Note

The designed is inspired by the plot method for gam objects in package mgcv.

Author(s)

Christophe Pouzat [e-mail]

See Also

quickPredict, plot.quickPredict

Examples

## Not run: 
data(e060824spont)
DFA <- subset(mkGLMdf(e060824spont,0.004,0,59),neuron==1)
DFA <- within(DFA,i1 <- isi(DFA,lag=1))
DFA <- DFA[complete.cases(DFA),]
m2u1 <- mkM2U(DFA,"lN.1",0,29)
m2ui <- mkM2U(DFA,"i1",0,29,maxiter=200)
DFA <- within(DFA,e1t <- m2u1(lN.1))
DFA <- within(DFA,i1t <- m2ui(i1))
with(DFA,plot(ecdf(e1t[time>29]),pch="."))
abline(a=0,b=1,col=2,lty=2)
with(DFA,plot(ecdf(i1t[time>29]),pch="."))
abline(a=0,b=1,col=2,lty=2)
m1.fit <- gssanova(event~e1t*i1t, data=subset(DFA,time>29), family="binomial", seed=20061001)
inter.pred <- m1.fit %qp% "e1t:i1t"
contour(inter.pred,what="mean",nlevels=10,col=2,lwd=2)
contour(inter.pred,what="sd",nlevels=5,col=1,lwd=1,lty=2,add=TRUE)
inter.predN <- changeScale(inter.pred,attr(m2u1,"qFct"),attr(m2ui,"qFct"))
contour(inter.predN,what="mean",nlevels=5,col=2,lwd=1)
contour(inter.predN,what="sd",nlevels=3,col=1,lwd=1,lty=2,add=TRUE)
plot(m1.fit,nr=3,nc=1)

## End(Not run)

STAR documentation built on May 2, 2019, 11:44 a.m.