R语言程序询问

想请各位帮忙分析一下R程序中这一句：plot(ecpe, pattern=data.ecpe$data[1,-1])  个人不是很懂。谢谢大家了，希望解答时能够 讲的详细些！

pattern可能是正则表达式，把整段代码贴出来看看吧

elvajear 发表于 2019-6-11 15:35
pattern可能是正则表达式，把整段代码贴出来看看吧

您好，谢谢您的热心帮助。我现在对于这一段代码已经理解了。不过有出现了另一个问题。在出来的数据图中，第一个图中的图例遮挡了其中的数据，所以我想将图例进行修改，但是没有成功，所以想向您寻求一下帮助，怎样修改图例，让得出的图更加的美观。谢谢您的帮助！以下是我的代码：

1. # load package
   
2. library(CDM)
   
3. 
   
4. # estimate DINA model
   
5. ecpe <- din(data.ecpe$data[,-1], data.ecpe$q.matrix)
   
6. 
   
7. # parameters and standard errors of DINA model
   
8. param <- IRT.se(ecpe, extended=TRUE)
   
9. p <- split(param, param$partype)
   
10. 
    
11. # items characteristics (cf. Table 6)
    
12. pvalues <- colMeans(data.ecpe$data[,-1], na.rm=TRUE) # item p-values
    
13. p$guess # guessing parameters
    
14. p$slip # slipping parameters
    
15. omega1 <- 1 - p$guess$est - p$slip$est # item discrimination
    
16. omega2 <- (p$guess$est + (1 - p$slip$est))/2 # item easiness
    
17. 
    
18. # skill characteristics (cf. Table 7 and 8)
    
19. p$margprobs # skill distribution Q1
    
20. p$probs # skill class distribution Q2
    
21. IRT.factor.scores(ecpe, type="MLE")[1:5,] # individual skill profile Q3
    
22. 
    
23. # plot model parameters
    
24. par(mfrow=c(2,2))
    
25. plot(ecpe, pattern=data.ecpe$data[1,-1])

复制代码

plot中的ecpe是使用din函数估计出来的结果，data.ecpe$data[1,-1]是出去第一列后的第一行

如果我没记错的话，IRT.se中计算出来的标准误是错的。建议使用我开发的包估计标准误，https://cran.r-project.org/web/packages/dcminfo/index.html

想要自己定义图的话，建议你了解一下tiff，pdf这个两个函数