set.seed(123)
x <- sample(1:4, 250, replace=TRUE)
pr1 <- pram(x)
length(which(pr1$x == x))
x2 <- sample(1:4, 250, replace=TRUE)
length(which(pram(x2)$x == x2))
data(free1)
marstatPramed <- pram(free1[,"MARSTAT"])
Run the code above in your browser using DataLab