{"id":1251,"date":"2023-12-19T21:54:37","date_gmt":"2023-12-19T13:54:37","guid":{"rendered":"https:\/\/edu.secda.info\/scu12354014\/?p=1251"},"modified":"2023-12-19T21:56:55","modified_gmt":"2023-12-19T13:56:55","slug":"20231205%e7%ac%ac%e5%8d%81%e4%b8%89%e9%80%b1-%e5%a4%9a%e8%ae%8a%e9%87%8f%e4%bd%9c%e6%a5%ad","status":"publish","type":"post","link":"https:\/\/edu.secda.info\/scu12354014\/?p=1251","title":{"rendered":"20231205\u7b2c\u5341\u4e09\u9031-\u591a\u8b8a\u91cf\u4f5c\u696d"},"content":{"rendered":"<p><span style=\"color: #333399;\">####################HW###########################<\/span><br \/>\n<span style=\"color: #333399;\">#\u5f9e\u8cc7\u6599\u5eab\u4e2d\u7684 \u6d89\u5165\u5ea6\u3001\u4e92\u52d5\u3001\u54c1\u724c\u5f62\u8c61\u3001\u4fe1\u4efb \u56db\u500b\u8b8a\u9805\u4f9d\u64da \u4e3b\u6210\u5206\u5206\u6790\uff0c\u8403<\/span><br \/>\n<span style=\"color: #333399;\">#\u53d6\u7b2c\u4e00\u500b\u4e3b\u6210\u5206(\u8981\u547d\u540d)\uff0c\u7136\u5f8c\u9032\u884c\u5c0d\u8cfc\u8cb7\u610f\u9858\u7684\u5f71\u97ff\uff0c\u5206\u5225\u6aa2\u9a57<\/span><br \/>\n<span style=\"color: #333399;\">#(1) \u7dda\u6027\u8ff4\u6b78\u5206\u6790<\/span><br \/>\n<span style=\"color: #333399;\">#(2) \u4e2d\u4f4d\u6578\u8ff4\u6b78\u5206\u6790<\/span><br \/>\n<span style=\"color: #333399;\">#(3) \u662f\u5426\u6709\u4e8c\u6b21\u66f2\u7dda\u7279\u5fb5<\/span><br \/>\n<span style=\"color: #333399;\">#\u52a0\u5165\u6027\u5225\u548c\u6559\u80b2\u7a0b\u5ea6\u70ba\u63a7\u5236\u8b8a\u6578\u3002<\/span><\/p>\n<blockquote><p>Data&lt;-read.csv(&#8220;C:\/Users\/zxcvb\/Desktop\/\u7814\/\u591a\u8b8a\u91cf\/data3.csv&#8221;)<br \/>\nsummary(Data)<br \/>\nstr(Data)<br \/>\nfor(i in 2:31){<br \/>\ntemp_mean&lt;-mean(Data[,i],na.rm=T)<br \/>\nna.rows&lt;-is.na(Data[,i])<\/p>\n<p>Data[na.rows,i]&lt;-temp_mean<br \/>\n}<br \/>\nstr(Data)<br \/>\nData$\u6d89\u5165\u5ea6 &lt;-rowMeans(Data[ ,c(2:6,8)])<br \/>\nData$ \u4e92\u52d5 &lt;-rowMeans(Data[ ,c(9:13)])<br \/>\nData$ \u4fe1\u4efb &lt;-rowMeans(Data[ ,c(14:20)])<br \/>\nData$\u54c1\u724c\u5f62\u8c61 &lt;-rowMeans(Data[ ,c(21:26)])<br \/>\nData$\u8cfc\u8cb7\u610f\u9858 &lt;-rowMeans(Data[ ,c(27:31)])<\/p>\n<p>Data$High &lt;- ifelse(Data$\u8cfc\u8cb7\u610f\u9858 &gt;= median(Data$\u8cfc\u8cb7\u610f\u9858),1,0)<\/p>\n<p>Data$sex &lt;- ifelse(Data$g1==1,&#8217;Boy&#8217;,&#8217;Girl&#8217;)<br \/>\nfor (i in 1:384){<br \/>\nif( Data[i,&#8217;g2&#8242;] &lt;=3 ) {Data[i,&#8217;Edu&#8217;] &lt;- &#8216;\u9ad8\u4e2d(\u542b\u4ee5\u4e0b)&#8217;}<br \/>\nif( Data[i,&#8217;g2&#8242;] ==4 ) {Data[i,&#8217;Edu&#8217;] &lt;- &#8216;\u5927\u5b78\u3001\u5927\u5c08&#8217;}<br \/>\nif( Data[i,&#8217;g2&#8242;] ==5 ) {Data[i,&#8217;Edu&#8217;] &lt;- &#8216;\u7814\u7a76\u6240&#8217;}<br \/>\n}<\/p><\/blockquote>\n<p><span style=\"color: #333399;\">#########################PCA<\/span><br \/>\n<span style=\"color: #333399;\">##\u4e3b\u6210\u5206\u5206\u6790<\/span><br \/>\n<span style=\"color: #333399;\">###1.\u7522\u751f\u4e00\u500b\u65b0\u7684\u8b8a\u9805\uff0c\u7531\u820a\u7684\u8b8a\u9805\u7dda\u6027\u7d44\u5408(\u52a0\u6b0a)\u800c\u6210<\/span><br \/>\n<span style=\"color: #333399;\">###\u4e3b\u8981\u6e90\u81ea\u65bc \u820a\u7684\u8b8a\u9805\u6709\u9ad8\u5ea6\u76f8\u95dc<\/span><br \/>\n<span style=\"color: #333399;\">###\u6709\u4e00\u7a2e \u5f62\u6210\u6027\u8b8a\u9805 = \u4e3b\u6210\u5206<\/span><br \/>\n<span style=\"color: #333399;\">###\u76ee\u7684<\/span><br \/>\n<span style=\"color: #333399;\">###\u8907\u88fd \u820a\u7684\u8b8a\u9805\u7684\u8b8a\u7570\u6578 =&gt; \u5f62\u6210 \u4e00\u500b\u65b0\u7684\u6307\u6a19 \u7528\u5728\u5f8c\u9762\u7684\u5206\u6790<\/span><br \/>\n<span style=\"color: #333399;\">###\u628a\u53d7\u8a66\u8005\u52a0\u4ee5 \u5dee\u7570\u5316<\/span><br \/>\n<span style=\"color: #333399;\">###2.\u9664\u4e86\u5dee\u7570\u5316\u4ee5\u5916<\/span><br \/>\n<span style=\"color: #333399;\">###\u89e3\u6c7a \u9ad8\u5ea6\u5171\u7dda\u6027\u554f\u984c =&gt; \u820a\u7684\u8b8a\u9805\u9ad8\u5ea6\u76f8\u95dc<\/span><br \/>\n<span style=\"color: #333399;\">###\u505a \u4e3b\u6210\u5206\u56de\u6b78<\/span><br \/>\n<span style=\"color: #333399;\">###3.\u5728\u56e0\u7d20\u5206\u6790\u4e2d\u505a\u70ba\u8403\u53d6\u6cd5\u4e4b\u7528<\/span><\/p>\n<p><span style=\"color: #333399;\">##\u7279\u5fb5<\/span><br \/>\n<span style=\"color: #333399;\">##1.\u6709\u591a\u5c11\u500b\u820a\u7684\u8b8a\u9805 \u5c31\u6703 \u5f62\u6210\u591a\u5c11\u500b\u4e3b\u6210\u5206<\/span><br \/>\n<span style=\"color: #333399;\">##2.\u4e0a\u8ff0\u4e3b\u6210\u5206\u4e4b\u9593 \u5f7c\u6b64\u7368\u7acb<\/span><br \/>\n<span style=\"color: #333399;\">##3.\u7b2c\u4e00\u500b\u4e3b\u6210\u5206 \u7d04\u7565\u5c31\u662f \u820a\u7684\u8b8a\u9805\u7684\u5e73\u5747\u6578<\/span><br \/>\n<span style=\"color: #333399;\">##4.\u9019\u500b\u7b2c\u4e00\u4e3b\u6210\u5206\u6709\u6700\u5927\u7684\u8b8a\u7570\u6578\uff0c\u5c31\u662f\u820a\u7684\u8b8a\u9805\u7684\u6700\u5927\u7279\u5fb5\u503c<\/span><br \/>\n<span style=\"color: #333399;\">##5.\u7b2c\u4e8c\u500b\u4e3b\u6210\u5206 \u7d04\u7565\u662f \u820a\u8b8a\u9805\u4e4b\u9593\u7684\u5c0d\u6bd4(\u6b63\u8ca0)<\/span><\/p>\n<blockquote><p>PCAHW &lt;- princomp(Data[,34:37],cor = TRUE,scores= TRUE)<br \/>\nsummary(PCAHW)<\/p><\/blockquote>\n<p><span style=\"color: #333399;\">#Importance of components:<\/span><br \/>\n<span style=\"color: #333399;\">#\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 Comp.1\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 Comp.2\u00a0 \u00a0 \u00a0 \u00a0 \u00a0Comp.3\u00a0 \u00a0 \u00a0 \u00a0 Comp.4<\/span><br \/>\n<span style=\"color: #333399;\">#Standard deviation\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 1.5917706\u00a0 \u00a0 \u00a0 0.8240888\u00a0 \u00a0 \u00a00.6724596\u00a0 \u00a0 0.57874177<\/span><br \/>\n<span style=\"color: #333399;\">#Proportion of Variance\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 0.6334334\u00a0 \u00a0 \u00a00.1697806\u00a0 \u00a0 \u00a0 \u00a00.1130505\u00a0 \u00a0 0.08373551<\/span><br \/>\n<span style=\"color: #333399;\">#Cumulative Proportion\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a00.6334334\u00a0 \u00a0 0.8032140\u00a0 \u00a0 \u00a0 0.9162645\u00a0 \u00a0 1.00000000<\/span><\/p>\n<p><span style=\"color: #333399;\">#\u7b2c\u4e00\u500b\u7d44\u6210\u5206\u8ddf\u7b2c\u4e8c\u500b\uff0c\u4e00\u8d77\u89e3\u91cb\u4e8680% (0.8032140)<\/span><br \/>\n<span style=\"color: #333399;\">#\u7b2c\u4e09\u500b\u53ea\u6709\u89e3\u91cb11%\uff0c\u70ba\u751a\u9ebc\u53ef\u4ee5\u4e00\u8d77\u89e3\u91cb\uff0c\u56e0\u70ba\u4ed6\u5011\u4e92\u76f8\u7368\u7acb<\/span><\/p>\n<blockquote><p>PCAHW$loadings<\/p><\/blockquote>\n<p><span style=\"color: #333399;\">##Loadings: \u52a0\u6b0a\u4fc2\u6578<\/span><\/p>\n<p><span style=\"color: #333399;\">#\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 Comp.1\u00a0 Comp.2\u00a0 \u00a0Comp.3\u00a0 \u00a0Comp.4<\/span><br \/>\n<span style=\"color: #333399;\">#\u6d89\u5165\u5ea6\u00a0 \u00a0 \u00a0 \u00a0 \u00a00.416\u00a0 \u00a0 0.909<\/span><br \/>\n<span style=\"color: #333399;\">#\u4e92\u52d5\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a00.536\u00a0 \u00a0-0.236\u00a0 \u00a0 \u00a0-0.343\u00a0 \u00a0 \u00a0 \u00a0 0.734<\/span><br \/>\n<span style=\"color: #333399;\">#\u4fe1\u4efb\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a00.532\u00a0 \u00a0-0.223\u00a0 \u00a0 \u00a0-0.459\u00a0 \u00a0 \u00a0 \u00a0-0.675<\/span><br \/>\n<span style=\"color: #333399;\">#\u54c1\u724c\u5f62\u8c61\u00a0 \u00a0 0.506\u00a0 \u00a0-0.263\u00a0 \u00a0 \u00a0 \u00a00.819<\/span><\/p>\n<p><span style=\"color: #333399;\">#\u7b2c\u4e00\u500b\u7d44\u6210\u5206\uff0c\u76f8\u7576\u65bc\u820a\u7684\u8b8a\u9805\u7684\u5e73\u5747\u6578<\/span><br \/>\n<span style=\"color: #333399;\">#\u7b2c\u4e8c\u500b\u7d44\u6210\u5206\uff0c\u76f8\u7576\u65bc&#8217;\u5c0d\u6bd4&#8217;<\/span><\/p>\n<blockquote><p>PCA1 &lt;- PCAHW $scores[, 1]<\/p>\n<p>HW&lt;- cbind(Data,PCA1)<\/p>\n<p>summary(reg &lt;- lm(data = HW ,\u8cfc\u8cb7\u610f\u9858 ~ PCA1))<\/p><\/blockquote>\n<p><strong><span style=\"color: #333399;\">#(1) \u7dda\u6027\u8ff4\u6b78\u5206\u6790<\/span><\/strong><\/p>\n<blockquote><p>reg1 &lt;- lm(data =HW ,\u8cfc\u8cb7\u610f\u9858 ~ PCA1 + sex + Edu)<\/p><\/blockquote>\n<p><strong><span style=\"color: #333399;\">#(2) \u4e2d\u4f4d\u6578\u8ff4\u6b78\u5206\u6790<\/span><\/strong><\/p>\n<blockquote><p>library(quantreg)<\/p>\n<p>reg2 &lt;- rq(data =HW ,\u8cfc\u8cb7\u610f\u9858 ~ PCA1 + sex + Edu)<\/p><\/blockquote>\n<p><strong><span style=\"color: #333399;\">#(3) \u662f\u5426\u6709\u4e8c\u6b21\u66f2\u7dda\u7279\u5fb5<\/span><\/strong><\/p>\n<blockquote><p>HW$PCA1_squared &lt;- HW$PCA1^2<\/p>\n<p>reg3 &lt;- lm(\u8cfc\u8cb7\u610f\u9858 ~ PCA1 + PCA1_squared + sex + Edu, data = HW)<\/p>\n<p>library(stargazer)<br \/>\nstargazer(reg1,reg2,reg3,<br \/>\ntype = &#8216;html&#8217;, #html , #text<br \/>\ndigits = 2 ,<br \/>\n#header = FALSE,<br \/>\nout = &#8216;C:\/Users\/zxcvb\/Desktop\/\u7814\/\u591a\u8b8a\u91cf\/logit.html&#8217;,<br \/>\ntitle = &#8216;Regression Results&#8217;<br \/>\n)<\/p><\/blockquote>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"alignnone wp-image-1252\" src=\"https:\/\/edu.secda.info\/scu12354014\/wp-content\/uploads\/2023\/12\/\u87a2\u5e55\u64f7\u53d6\u756b\u9762-2023-12-19-213457-251x300.png\" alt=\"\" width=\"510\" height=\"610\" srcset=\"https:\/\/edu.secda.info\/scu12354014\/wp-content\/uploads\/2023\/12\/\u87a2\u5e55\u64f7\u53d6\u756b\u9762-2023-12-19-213457-251x300.png 251w, https:\/\/edu.secda.info\/scu12354014\/wp-content\/uploads\/2023\/12\/\u87a2\u5e55\u64f7\u53d6\u756b\u9762-2023-12-19-213457.png 597w\" sizes=\"(max-width: 510px) 100vw, 510px\" \/><\/p>\n","protected":false},"excerpt":{"rendered":"<p>####################HW##################<\/p>\n","protected":false},"author":16,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1,26],"tags":[],"_links":{"self":[{"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=\/wp\/v2\/posts\/1251"}],"collection":[{"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=\/wp\/v2\/users\/16"}],"replies":[{"embeddable":true,"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=1251"}],"version-history":[{"count":4,"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=\/wp\/v2\/posts\/1251\/revisions"}],"predecessor-version":[{"id":1256,"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=\/wp\/v2\/posts\/1251\/revisions\/1256"}],"wp:attachment":[{"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1251"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1251"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/edu.secda.info\/scu12354014\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1251"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}