{"id":620,"date":"2017-07-06T13:36:48","date_gmt":"2017-07-06T05:36:48","guid":{"rendered":"https:\/\/www.beneplot.com\/?p=620"},"modified":"2020-07-18T17:28:15","modified_gmt":"2020-07-18T09:28:15","slug":"%e9%a2%98%e9%87%8f%e4%b8%8e%e6%a8%a1%e5%9e%8b%e8%ae%a1%e7%ae%97%e7%b2%be%e5%ba%a6%e7%9a%84%e5%85%b3%e7%b3%bb-r-markdown-%e8%87%aa%e5%8a%a8%e5%8c%96%e7%94%9f%e6%88%90pdf%e5%88%86%e6%9e%90%e6%8a%a5","status":"publish","type":"post","link":"https:\/\/www.beneplot.com\/?p=620","title":{"rendered":"\u9898\u91cf\u4e0e\u6a21\u578b\u8ba1\u7b97\u7cbe\u5ea6\u7684\u5173\u7cfb (R Markdown \u81ea\u52a8\u5316\u751f\u6210PDF\u5206\u6790\u62a5\u544a)"},"content":{"rendered":"

\u5728\u81ea\u52a8\u5316\u62a5\u544a\u7b97\u6cd5 (R\u4e0eLatex)<\/a><\/strong><\/span>\u4e00\u6587\u4e2d\uff0c\u8bb2\u8ff0\u4e86\u5982\u4f55\u7528R\u4e0eLatex\u751f\u6210\u81ea\u52a8\u5316\u62a5\u544a\uff0cMarkdown\u4e5f\u662f\u4e00\u79cd\u81ea\u52a8\u5316\u751f\u6210\u62a5\u544a\u7684\u8bed\u8a00\uff0c\u76f8\u6bd4Latex\uff0c\u8bed\u53e5\u66f4\u7b80\u5355\u3002<\/span><\/h3>\n

\u5728Rstudio\u4e2d\uff0c\u53ef\u4ee5\u76f4\u63a5\u8fdb\u884cR Markdown\uff0c\u5176\u6587\u4ef6\u540d\u540e\u7f00\u4e3a.rmd\u3002<\/span><\/h3>\n

\u672c\u6587\u901a\u8fc7\u5206\u6790\u9898\u91cf\u4e0eIRT\u6a21\u578b\u8ba1\u7b97\u7cbe\u5ea6\u7684\u5173\u7cfb\uff0c\u5229\u7528R\u8bed\u53e5\u4e0eMarkdown\u7ed3\u5408\u751f\u6210\u81ea\u52a8\u5316\u7684\u5206\u6790\u62a5\u544a\u3002<\/span><\/h3>\n

*\u6ce8\u610f\uff01\u5982\u679c\u4e0d\u901a\u8fc7\u4e00\u4e9b\u8bbe\u5b9a\uff0c\u65e0\u8bba\u662flatex\u8fd8\u662fmarkdown\u90fd\u4e0d\u8bc6\u522b\u4e2d\u6587\uff0c\u7ed3\u679c\u4f1a\u62a5\u9519\u3002\u5728\u672c\u4f8b\u4e2d\uff0c\u9700\u786e\u4fdd\u5df2\u5b89\u88c5rticles\u5305\uff0c\u65b0\u5efarmarkdown\u6587\u4ef6\u65f6\u4eceFrom Template\u4e2d\u9009\u62e9Ctex Document\u3002<\/span><\/h3>\n

\u00a0<\/h3>\n
\u4ee5\u4e0b\u662f\u4ee3\u7801\uff1a<\/span><\/span><\/span><\/code><\/div>\n
\u00a0<\/div>\n
---<\/span><\/code><\/div>\n
title: \"\u9898\u91cf\u4e0e\u6a21\u578b\u8ba1\u7b97\u7cbe\u5ea6\u7684\u5173\u7cfb\"<\/span><\/code><\/div>\n
author: \"Bene\"<\/span><\/code><\/div>\n
documentclass: ctexart<\/span><\/code><\/div>\n
output:<\/span><\/code><\/div>\n
\u00a0 rticles::ctex:<\/span><\/code><\/div>\n
\u00a0 pdf_document: default<\/span><\/code><\/div>\n
\u00a0 html_document: default<\/span><\/code><\/div>\n
---<\/span><\/code><\/div>\n
\u00a0<\/div>\n
```{r setup, results='hide',include=FALSE}<\/span><\/code><\/div>\n
knitr::opts_chunk$set(echo = TRUE)<\/span><\/code><\/div>\n
```<\/span><\/code><\/div>\n
\u00a0<\/div>\n
```{r, echo=FALSE, results='hide', message = FALSE, warning = FALSE}<\/span><\/code><\/div>\n
rm(list = ls())\u00a0 #\u6e05\u7a7a\u6240\u6709\u5185\u5b58#<\/span><\/code><\/div>\n
options(warn=-1) ##\u8be5\u8bed\u53e5\u7528\u4e8e\u7981\u6b62\u62a5\u544a\u8b66\u544a\uff0c\u4e0e\u4e0a\u65b9warning = FALSE\u610f\u4e49\u76f8\u540c##<\/span><\/code><\/div>\n
\u00a0<\/div>\n
##\u8f7d\u5165\u5206\u6790\u5305###<\/span><\/code><\/div>\n
library(TAM)<\/span><\/code><\/div>\n
library(sirt)<\/span><\/code><\/div>\n
library('ggplot2')<\/span><\/code><\/div>\n
library(knitr)<\/span><\/code><\/div>\n
\u00a0<\/div>\n
options(digits=3)##\u8bbe\u7f6e\u5c0f\u6570\u70b9\u4fdd\u75593\u4f4d##<\/span><\/code><\/div>\n
set.seed(123)##\u8bbe\u7f6e\u968f\u673a\u7ed3\u679c\u56fa\u5b9a\u5316\uff0c\u76ee\u7684\u5728\u4e8e\u53ef\u4ee5\u91cd\u73b0\u5206\u6790\u7ed3\u679c##<\/span><\/code><\/div>\n
\u00a0<\/div>\n
##\u6784\u5efa\u6a21\u62df\u7684\u5206\u6790\u6570\u636e###<\/span><\/code><\/div>\n
N = 100 ##\u6837\u672c\u6570\u91cf##<\/span><\/code><\/div>\n
i = 50\u00a0 ##\u9898\u91cf##<\/span><\/code><\/div>\n
b = stats::rnorm( i ) ##\u968f\u673a\u751f\u621050\u9053\u9898\u7684\u9898\u76ee\u53c2\u6570##<\/span><\/code><\/div>\n
d= c(3:i) ##\u7528\u4e8e\u6784\u5efa\u542b\u6709\u4e0d\u540c\u9898\u91cf\u7684\u6a21\u578b##<\/span><\/code><\/div>\n
y = stats::rnorm( N ) ##\u968f\u673a\u751f\u6210100\u4e2a\u6837\u672c\u53c2\u6570##<\/span><\/code><\/div>\n
\u00a0<\/div>\n
#####\u6784\u5efaIRT\u80fd\u529b\u5206\u6790function####<\/span><\/code><\/div>\n
PersonAbility = function(d)<\/span><\/code><\/div>\n
{<\/span><\/code><\/div>\n
\u00a0 return( tam.wle(TAM::tam.mml(resp = sirt::sim.raschtype( y, sample(b,d)), constraint = \"items\"))$theta )<\/span><\/code><\/div>\n
}<\/span><\/code><\/div>\n
\u00a0<\/div>\n
###\u5728\u4e0d\u540c\u9898\u91cf\u4e0b\uff0c\u6c42\u6d4b\u8bd5\u8005\u7684IRT\u80fd\u529b\u503c##<\/span><\/code><\/div>\n
din = as.data.frame(sapply(d, PersonAbility))<\/span><\/code><\/div>\n
\u00a0<\/div>\n
##\u4ee5\u5747\u65b9\u6839\u8bef\u5dee(RMSE)\u4f5c\u4e3a\u8bef\u5dee\u8861\u91cf\u6307\u6807\uff0c\u6784\u5efaRMSE\u7684\u5206\u6790function##<\/span><\/code><\/div>\n
rmse = function(y, h)<\/span><\/code><\/div>\n
{<\/span><\/code><\/div>\n
\u00a0 return(sqrt(mean((y - h) ^ 2)))<\/span><\/code><\/div>\n
}<\/span><\/code><\/div>\n
\u00a0<\/div>\n
```<\/span><\/code><\/div>\n
\u00a0<\/div>\n
## \u5206\u6790\u7ed3\u679c<\/span><\/code><\/div>\n
### \u5728\u672c\u4f8b\u4e2d, RMSE = sqrt(mean((\u539f\u59cb\u89c2\u6d4b\u503c - \u9884\u6d4b\u503c) ^ 2)), \u5176\u4e2d\u9884\u6d4b\u503c\u4e3a\u7ebf\u6027\u6a21\u578b\u4e2d\u7684\u9884\u6d4b\u503c\u3002\u4e0b\u8868\u4e3a\u4e0d\u540c\u9898\u91cf\u4e0b\u7684\u8bef\u5dee\u3002<\/span><\/code><\/div>\n
\u00a0<\/div>\n
```{r, echo=FALSE}<\/span><\/code><\/div>\n
\u00a0<\/div>\n
###\u6784\u5efa\u9898\u91cf\u4e0e\u8bef\u5dee\u5bf9\u5e94\u7684\u5206\u6790\u6a21\u578b###<\/span><\/code><\/div>\n
\u00a0<\/div>\n
performance = data.frame()<\/span><\/code><\/div>\n
\u00a0<\/div>\n
for (d in 1:length(din))<\/span><\/code><\/div>\n
{<\/span><\/code><\/div>\n
\u00a0 fit = lm(y ~ din[,d])<\/span><\/code><\/div>\n
\u00a0\u00a0<\/span><\/code><\/div>\n
\u00a0 performance = rbind(performance,<\/span><\/code><\/div>\n
\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0data.frame(Degree = d,<\/span><\/code><\/div>\n
\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 RMSE = rmse(din[,d], predict(fit))))<\/span><\/code><\/div>\n
\u00a0\u00a0<\/span><\/code><\/div>\n
}<\/span><\/code><\/div>\n
performance$Degree = performance$Degree + 2<\/span><\/code><\/div>\n
\u00a0<\/div>\n
###\u4ee5\u4e0b\u51e0\u53e5\u4ee3\u7801\u7528\u4e8e\u5c06\u957f\u8868\u683c\u8f6c\u5316\u4e3a\u77ed\u8868\u683c\uff0c\u4ee5\u4fbf\u663e\u793a\u7f8e\u89c2###<\/span><\/code><\/div>\n
#quantile(performance$Degree)<\/span><\/code><\/div>\n
#\u00a0 0%\u00a0 25%\u00a0 50%\u00a0 75% 100%\u00a0<\/span><\/code><\/div>\n
#\u00a0 3.0 14.8 26.5 38.2 50.0\u00a0<\/span><\/code><\/div>\n
\u00a0<\/div>\n
sub1= subset(performance, performance$Degree >= 3 & performance$Degree <= 14)<\/span><\/code><\/div>\n
sub2= subset(performance, performance$Degree >= 15 & performance$Degree <= 26)<\/span><\/code><\/div>\n
sub3= subset(performance, performance$Degree >= 27 & performance$Degree <= 38)<\/span><\/code><\/div>\n
sub4= subset(performance, performance$Degree >= 39 & performance$Degree <= 50)<\/span><\/code><\/div>\n
xxx = cbind(sub1,sub2,sub3,sub4)<\/span><\/code><\/div>\n
\u00a0<\/div>\n
kable(xxx,col.names=c(\"\u9898\u91cf\",\"\u5206\u6790\u8bef\u5dee\",\"\u9898\u91cf\",\"\u5206\u6790\u8bef\u5dee\",\"\u9898\u91cf\",\"\u5206\u6790\u8bef\u5dee\",\"\u9898\u91cf\",\"\u5206\u6790\u8bef\u5dee\"),align='cccccccc')<\/span><\/code><\/div>\n
\u00a0<\/div>\n
\u00a0<\/div>\n
```<\/span><\/code><\/div>\n
\u00a0<\/div>\n
## \u56fe\u5f62<\/span><\/code><\/div>\n
### \u5bf9\u7ed3\u679c\u4f5c\u56fe\uff0c\u663e\u793a\u6d4b\u8bd5\u9898\u91cf\u4e0e\u5206\u6790\u8bef\u5dee\u7684\u5173\u7cfb\u3002<\/span><\/code><\/div>\n
\u00a0<\/div>\n
```{r, echo=FALSE}<\/span><\/code><\/div>\n
ggplot(performance, aes(x = Degree, y = RMSE)) +<\/span><\/code><\/div>\n
\u00a0 labs(x=\"Item number\",y=\"Analytical error\")+<\/span><\/code><\/div>\n
\u00a0 geom_point() +<\/span><\/code><\/div>\n
\u00a0 geom_line()<\/span><\/code><\/div>\n
```<\/span><\/code><\/div>\n
\u00a0<\/div>\n

\u6700\u7ec8\u751f\u6210\u7684PDF\u5982\u4e0b\uff1a<\/span><\/h3>\n

\u81ea\u52a8\u5316\u751f\u6210\u7684PDF\u62a5\u544a<\/a><\/strong><\/span><\/h3>\n

 <\/p>\n

 <\/p>\n","protected":false},"excerpt":{"rendered":"

\u5728\u81ea\u52a8\u5316\u62a5\u544a\u7b97\u6cd5 (R\u4e0eLatex)\u4e00\u6587\u4e2d\uff0c\u8bb2\u8ff0\u4e86\u5982\u4f55\u7528R\u4e0eLatex\u751f\u6210\u81ea\u52a8\u5316\u62a5\u544a\uff0cMarkdown\u4e5f\u662f\u4e00\u79cd […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":[],"categories":[10],"tags":[],"yoast_head":"\n\n\n\n\n\n\n\n\n\n\n\n\n\n