llrdata <-
function(n = 100, q=5)
{
# Generates data for loglinear regression.
#
y <- 0 * 1:n
beta <- 0 * 1:q
beta[1:3] <- 1
alpha <- -2.5
x <- matrix(rnorm(n * q), nrow = n,
ncol = q)
x <- 0.5*x + 1
SP <- alpha + x%*%beta
y <- rpois(n,lambda=exp(SP))
list(x=x,y=y)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.