######################################### # Itens item1<--1.9 item2<--0.6 item3<--0.25 item4<-0.30 item5<-0.45 veros<-function(theta=0) { px1<-1/(1+exp(-1*(theta-item1))) px2<-1/(1+exp(-1*(theta-item2))) px3<-1/(1+exp(-1*(theta-item3))) px4<-1/(1+exp(-1*(theta-item4))) px5<-1/(1+exp(-1*(theta-item5))) px1;px2;px3;px4;px5; L1<-1*log(px1)+(1-1)*log(1-px1);L1 L2<-1*log(px2)+(1-1)*log(1-px2);L2 L3<-0*log(px3)+(1-0)*log(1-px3);L3 L4<-0*log(px4)+(1-0)*log(1-px4);L4 L5<-0*log(px5)+(1-0)*log(1-px5);L5 Lv<-L1+L2+L3+L4+L5; Lv } L_1<-veros(theta=-3) L_2<-veros(theta=-2) L_3<-veros(theta=-1) L_4<-veros(theta=0) L_5<-veros(theta=1) L_6<-veros(theta=2) L_7<-veros(theta=3) log_veros<-c(L_1,L_2,L_3,L_4,L_5,L_6,L_7) plot(seq(-4,4,len=7),log_veros,xlab='theta',type='l') #######################