simRasch <- function(n,v,l,R=666,standard=F){
set.seed(R)
if (standard==F){
sigma <- seq(length=n,-3,3); theta <- seq(length=v,-3,3)
} else {
sigma <- rnorm(n,0,1); theta <- rnorm(v,0,1)
}
E <- outer(sigma, theta, function(r,c) (r-c) )
eta <- sapply(1:ncol(E), function(x) plogis(E[,x]))
Rasch <- eta
for (i in 1:nrow(Rasch)) {
for (j in 1:ncol(Rasch)) {
Rasch[i,j] <- rbinom(1,(l-1),eta[i,j])
}
}
colnames(Rasch) <- sapply(1:ncol(Rasch), function(x) paste("Item_",x,sep=""))
rownames(Rasch) <- sapply(1:nrow(Rasch), function(x) paste("Ind_",x,sep=""))
Result <- list("data"=Rasch,"abil"=sigma,"diff"=theta)
return(Result)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.