\begin{Schunk}
\begin{Sinput}
 N <- nrow(dat.kj)
 x <- cbind(dat.kj$v1,
             dat.kj$v2,
             dat.kj$s1,
             dat.kj$s2)
 y <- cbind(dat.kj$c1,
             dat.kj$c2)
 Kx <- ncol(x)
 Ky <- ncol(y)
 datstan <- list(N, Kx, Ky, y, x)
\end{Sinput}
\end{Schunk}
