#CODIGO CRM
rSamejimaIRT=function(N,nit,alpha,beta,delta,theta){
  y=matrix(,N,nit)
  for (i in 1:N)
    for (j in 1:nit){{
	  y[i,j]=plogis(alpha[i]*theta[i]+beta[j]+rnorm(1,0,sqrt(delta[j])))
    }}
  return(t)
}


#CODIGO CRM INFLAZIONADO

rSamejimaIRTinf=function(N,nit,alpha,beta,delta,tau,b0,b1){
  t=matrix(,N,nit)
  for (p in 1:N)
    for (i in 1:nit){{
      C=sample(0:2,1,prob=c(
	  plogis(b1[i]-alpha[i]*tau[p])-plogis(b0[i]-alpha[i]*tau[p]),
	  plogis(b0[i]-alpha[i]*tau[p]),
	  1-plogis(b1[i]-alpha[i]*tau[p])))
	  
      if(C==0) t[p,i]=plogis(alpha[i]*tau[p]+beta[i]+rnorm(1,0,sqrt(delta[i])))
      else t[p,i]=C-1
    }}
  return(t)
}