[R] Help vectorizing data generation for IRT graded response model

Hello, I have code that will generate data for a 5 category IRT graded response model. However, the code could be improved through vectorizing. Here is the code below: ###Inputs N <- 100 #Number of people taking test n <- 10 #Number of items nCat <- 5 #number of categories ###Generate Item pa

[R] Kernel Density Estimation in R

Hello, I am wanting to run a simulation study in R comparing several different bandwidth selection methods for data simulated from several different distribution types (normal, lognormal, bimodal, etc.) and wanted to know how to calculate the mean integrated square errors for the optimal smoothin