R语言进行项目参数估计的程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请求大佬帮忙看下程序哪里出错了，小白一个论文急需啊啊啊：谢谢谢谢！！
不知道哪里错了，项目区分度难度都不应该是这么大的数字啊。数据的EXCEL放在附件了谢谢各位大佬！！！

> install.packages("ltm")
--- 在此連線階段时请选用CRAN的鏡子 ---
试开URL’https://mirrors.tongji.edu.cn/CRAN/bin/windows/contrib/4.1/ltm_1.1-1.zip'
Content type 'application/zip' length 705388 bytes (688 KB)
downloaded 688 KB

package ‘ltm’ successfully unpacked and MD5 sums checked

The downloaded binary packages are in
        C:\Users\sky\AppData\Local\Temp\Rtmp6bWWD5\downloaded_packages
> library(ltm)
载入需要的程辑包：MASS
载入需要的程辑包：msm
载入需要的程辑包：polycor
> setwd("c:/mystudy/RIRT")
> ddd<-read.fwf("data1.csv", widths=c(1,rep(1,18)))
> dim(ddd)
[1] 26 19
> ddd <- ddd[,-1]
> names(ddd) <- paste("It",1:18, sep="")
> mod.2pl <- ltm(ddd~z1, IRT.param=T)
> mod.2pl$conv
[1] 0
> summary(mod.2pl)
Error in solve.default(object$hessian) :
  Lapack例行程序dgesv: 系统正好是奇异的: U[3,3] = 0
> coef(mod.2pl)
            Dffclt           Dscrmn
It1  -1.398595e+09 4.687995e-08
It2  -2.797913e+00 1.694660e-01
It3  -1.398595e+09 4.687995e-08
It4  -5.958137e-01 5.651422e-01
It5  -1.398595e+09 4.687995e-08
It6  -2.154298e+00 1.450527e-01
It7  -1.398595e+09 4.687995e-08
It8  -2.313392e-01 3.893327e+00
It9  -1.398595e+09 4.687995e-08
It10 -2.947178e+00 3.487284e-01
It11 -1.398595e+09 4.687995e-08
It12 -3.548480e-01 3.164737e+00
It13 -1.398595e+09 4.687995e-08
It14 -4.610567e-01 7.744310e-01
It15 -1.398595e+09 4.687995e-08
It16 -1.621804e+01 3.924119e-02
It17 -1.398595e+09 4.687995e-08
It18 -1.802842e+00 2.659201e-01
> factor.scores(mod.2pl, method="EAP")
Call:
ltm(formula = ddd ~ z1, IRT.param = T)

Scoring Method: Expected A Posteriori

Factor-Scores for observed response patterns:
   It1 It2 It3 It4 It5 It6 It7 It8 It9 It10 It11 It12 It13 It14 It15 It16
1    1   0   1   0   1   0   1   1   1    0    1    0    1    0    1    1
2    1   0   1   0   1   0   1   1   1    1    1    1    1    1    1    0
3    1   0   1   0   1   1   1   0   1    1    1    0    1    0    1    0
4    1   0   1   0   1   1   1   0   1    1    1    0    1    1    1    1
5    1   0   1   1   1   0   1   1   1    1    1    1    1    0    1    1
6    1   0   1   1   1   0   1   1   1    1    1    1    1    1    1    1
7    1   0   1   1   1   1   1   0   1    0    1    0    1    0    1    1
8    1   0   1   1   1   1   1   0   1    0    1    1    1    1    1    1
9    1   0   1   1   1   1   1   1   1    1    1    0    1    1    1    1
10   1   1   1   0   1   0   1   0   1    0    1    0    1    1    1    0
11   1   1   1   0   1   0   1   0   1    0    1    1    1    0    1    0
12   1   1   1   0   1   0   1   0   1    1    1    0    1    0    1    1
13   1   1   1   0   1   0   1   0   1    1    1    1    1    1    1    1
14   1   1   1   0   1   0   1   1   1    1    1    1    1    1    1    1
15   1   1   1   0   1   1   1   1   1    1    1    1    1    0    1    0
16   1   1   1   0   1   1   1   1   1    1    1    1    1    0    1    1
17   1   1   1   1   1   0   1   0   1    1    1    0    1    0    1    1
18   1   1   1   1   1   1   1   0   1    1    1    0    1    0    1    0
19   1   1   1   1   1   1   1   0   1    1    1    0    1    1    1    1
20   1   1   1   1   1   1   1   1   1    0    1    1    1    1    1    1
21   1   1   1   1   1   1   1   1   1    0    1    1    1    1    1    1
22   1   1   1   1   1   1   1   1   1    1    1    1    1    0    1    0
23   1   1   1   1   1   1   1   1   1    1    1    1    1    1    1    0
24   1   1   1   1   1   1   1   1   1    1    1    1    1    1    1    0
25   1   1   1   1   1   1   1   1   1    1    1    1    1    1    1    1
   It17 It18 Obs   Exp     z1 se.z1
1     1    1   1 0.008 -0.351 0.438
2     1    0   1 0.038  0.533 0.586
3     1    1   1 0.078 -1.184 0.603
4     1    1   1 0.089 -0.926 0.527
5     1    0   1 0.065  0.477 0.568
6     1    1   2 0.300  0.884 0.685
7     1    0   1 0.046 -1.188 0.604
8     1    1   1 0.022 -0.208 0.433
9     1    1   1 0.043 -0.005 0.441
10    1    0   1 0.023 -1.119 0.583
11    1    0   1 0.008 -0.521 0.449
12    1    1   1 0.160 -1.161 0.596
13    1    1   1 0.053 -0.244 0.434
14    1    1   1 0.251  0.711 0.638
15    1    0   1 0.047  0.388 0.540
16    1    1   1 0.165  0.482 0.570
17    1    0   1 0.100 -1.060 0.566
18    1    1   1 0.080 -0.948 0.533
19    1    0   1 0.083 -0.804 0.495
20    1    0   1 0.114  0.750 0.649
21    1    1   1 0.229  0.868 0.681
22    1    1   1 0.166  0.669 0.626
23    1    0   1 0.218  0.888 0.687
24    1    1   1 0.453  1.020 0.719
25    1    1   1 0.891  1.040 0.724

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：项目参数 参数估计 R语言 Successfully Application