setwd("D:\\")
library("CTT")
library("TAM")
d1<-read.csv("h2.csv",header = T)
mykey<-c(4,2,4,1,1,1,2,3,3,2,
         3,4,4,4,4,1,3,3,2,3,
         3,3,2,1,1,2,2,3,2,4,
         4,2,1,1,3,4,3,2,2,2)

         
         

myScore=score(d1,mykey,output.scored = T)
score<-data.frame(myScore$scored)
my1PL<-tam(myScore$scored)
myCTT<-itemAnalysis(myScore$scored)
write.csv(myCTT$itemReport,"h2_result.csv")



par(mai=c(0.1,0.1,0.1,0.1),mfrow=c(8,5))
for (i in 1:40) {
  plot(my1PL,item=i)  
}