{"id":447,"date":"2023-04-25T08:23:02","date_gmt":"2023-04-25T08:23:02","guid":{"rendered":"https:\/\/yunshangtulv.com.cn\/?p=447"},"modified":"2023-04-25T08:23:02","modified_gmt":"2023-04-25T08:23:02","slug":"r%e8%af%ad%e8%a8%80%e8%bf%9b%e9%98%b6%e7%bb%98%e5%9b%be%ef%bc%9a%e4%b8%89%e7%bb%b4pca","status":"publish","type":"post","link":"https:\/\/yunshangtulv.com.cn\/?p=447","title":{"rendered":"R\u8bed\u8a00\u8fdb\u9636\u7ed8\u56fe\uff1a\u4e09\u7ef4PCA"},"content":{"rendered":"

\u4e09\u7ef4\u4e3b\u6210\u5206\u5206\u6790\uff083D PCA\uff09\u662f\u4e00\u79cd\u7edf\u8ba1\u65b9\u6cd5\uff0c\u7528\u4e8e\u964d\u4f4e\u9ad8\u7ef4\u6570\u636e\u7684\u590d\u6742\u6027\uff0c\u540c\u65f6\u4fdd\u7559\u6570\u636e\u4e2d\u7684\u5927\u90e8\u5206\u53d8\u5f02\u3002\u901a\u8fc7\u5bf9\u539f\u59cb\u6570\u636e\u8fdb\u884c\u7ebf\u6027\u53d8\u6362\uff0cPCA\u63d0\u53d6\u51fa\u4e3b\u6210\u5206\uff0c\u8fd9\u4e9b\u4e3b\u6210\u5206\u662f\u539f\u59cb\u6570\u636e\u4e2d\u53d8\u5f02\u7684\u6b63\u4ea4\u7ebf\u6027\u7ec4\u5408\u3002\u5728\u4e09\u7ef4PCA\u4e2d\uff0c\u6211\u4eec\u63d0\u53d6\u524d\u4e09\u4e2a\u4e3b\u6210\u5206\uff0c\u5c06\u6570\u636e\u53ef\u89c6\u5316\u4e3a\u4e09\u7ef4\u56fe\u5f62\u3002<\/p>\n

\u4ee5\u4e0b\u662f\u4e00\u4e2a\u4f7f\u7528R\u8bed\u8a00\u6267\u884c\u4e09\u7ef4PCA\u7684\u793a\u4f8b\u3002\u8fd9\u91cc\u6211\u4eec\u5c06\u4f7f\u7528\u5185\u7f6e\u7684iris<\/code>\u6570\u636e\u96c6\u3002<\/p>\n

    \n
  1. \n

    \u52a0\u8f7d\u6240\u9700\u5e93\u548c\u6570\u636e\uff1a<\/h2>\n<\/li>\n<\/ol>\n
    \n
    <\/div>\n
    # \u52a0\u8f7d\u6240\u9700\u5e93<\/span>
    \nlibrary(<\/span>ggplot2)<\/span>
    \nlibrary(<\/span>plotly)<\/span><\/p>\n

    # \u52a0\u8f7d\u5185\u7f6e\u7684iris\u6570\u636e\u96c6<\/span>
    \ndata(<\/span>iris)<\/span>
    \nhead(<\/span>iris)<\/span>
    \n<\/code><\/div>\n<\/div>\n

      \n
    1. \n

      \u6570\u636e\u9884\u5904\u7406\u548cPCA\u8ba1\u7b97\uff1a<\/h2>\n<\/li>\n<\/ol>\n
      \n
      <\/div>\n
      # \u5bf9\u6570\u636e\u8fdb\u884c\u7f29\u653e<\/span>
      \nscaled_iris <-<\/span> scale(<\/span>iris[<\/span>,<\/span> -<\/span>5<\/span>]<\/span>)<\/span> # \u53bb\u9664\u7b2c\u4e94\u5217\uff08Species\uff09\uff0c\u56e0\u4e3a\u5b83\u662f\u5206\u7c7b\u53d8\u91cf<\/span><\/p>\n

      # \u8ba1\u7b97PCA<\/span>
      \npca_result <-<\/span> prcomp(<\/span>scaled_iris,<\/span> center =<\/span> TRUE<\/span>,<\/span> scale. =<\/span> TRUE<\/span>)<\/span>
      \n<\/code><\/div>\n<\/div>\n

        \n
      1. \n

        \u51c6\u5907\u7ed8\u52363D PCA\u56fe\u5f62\u7684\u6570\u636e\uff1a<\/h2>\n<\/li>\n<\/ol>\n
        \n
        <\/div>\n
        # \u63d0\u53d6PCA\u5206\u6570<\/span>
        \nscores <-<\/span> as.data.frame(<\/span>pca_result$<\/span>x[<\/span>,<\/span> 1<\/span>:<\/span>3<\/span>]<\/span>)<\/span><\/p>\n

        # \u4e3a\u6570\u636e\u6846\u6dfb\u52a0Species\u5217<\/span>
        \nscores$<\/span>Species <-<\/span> iris$<\/span>Species<\/p>\n

        # \u5c55\u793a\u524d\u51e0\u884c\u6570\u636e<\/span>
        \nhead(<\/span>scores)<\/span>
        \n<\/code><\/div>\n<\/div>\n

          \n
        1. \n

          \u4f7f\u7528plotly<\/code>\u5e93\u7ed8\u52363D PCA\u56fe\u5f62\uff1a<\/h2>\n<\/li>\n<\/ol>\n
          \n
          <\/div>\n
          # \u521b\u5efa3D\u6563\u70b9\u56fe<\/span>
          \np <-<\/span> plot_ly(<\/span>scores,<\/span> x =<\/span> ~<\/span>PC1,<\/span> y =<\/span> ~<\/span>PC2,<\/span> z =<\/span> ~<\/span>PC3,<\/span> color =<\/span> ~<\/span>Species,<\/span> type =<\/span> \"scatter3d\"<\/span>,<\/span> mode =<\/span> \"markers\"<\/span>,<\/span>
          \n marker =<\/span> list<\/span>(<\/span>size =<\/span> 5<\/span>,<\/span> opacity =<\/span> 0.8<\/span>)<\/span>)<\/span><\/p>\n

          # \u4fee\u6539\u56fe\u5f62\u5e03\u5c40<\/span>
          \np <-<\/span> p %>%<\/span> layout(<\/span>scene =<\/span> list<\/span>(<\/span>xaxis =<\/span> list<\/span>(<\/span>title =<\/span> \"PC1\"<\/span>)<\/span>,<\/span>
          \n yaxis =<\/span> list<\/span>(<\/span>title =<\/span> \"PC2\"<\/span>)<\/span>,<\/span>
          \n zaxis =<\/span> list<\/span>(<\/span>title =<\/span> \"PC3\"<\/span>)<\/span>)<\/span>)<\/span><\/p>\n

          # \u663e\u793a\u56fe\u5f62<\/span>
          \np
          \n<\/code><\/div>\n<\/div>\n

            \n
          1. \n

            \u4fdd\u5b58\u56fe\u5f62\u4e3aPDF\u6587\u4ef6\uff1a<\/h2>\n<\/li>\n<\/ol>\n
            \n
            <\/div>\n
            # \u5b89\u88c5\u548c\u52a0\u8f7dwebshot\u5e93<\/span>
            \nif<\/span> (<\/span>!<\/span>requireNamespace(<\/span>\"webshot\"<\/span>,<\/span> quietly =<\/span> TRUE<\/span>)<\/span>)<\/span> {<\/span>
            \n install.packages(<\/span>\"webshot\"<\/span>)<\/span>
            \n}<\/span>
            \nlibrary(<\/span>webshot)<\/span><\/p>\n

            # \u5c06\u56fe\u5f62\u4fdd\u5b58\u4e3aPDF\u6587\u4ef6<\/span>
            \nplotly_IMAGE(<\/span>p,<\/span> format =<\/span> \"pdf\"<\/span>,<\/span> out_file =<\/span> \"3D_PCA_plot.pdf\"<\/span>,<\/span> width =<\/span> 800<\/span>,<\/span> height =<\/span> 600<\/span>)<\/span>
            \n<\/code><\/div>\n<\/div>\n

            \u4ee5\u4e0a\u4ee3\u7801\u8be6\u7ec6\u89e3\u91ca\u5982\u4e0b\uff1a<\/p>\n

              \n
            • \u52a0\u8f7dggplot2<\/code>\u548cplotly<\/code>\u5e93\uff0c\u4ee5\u53ca\u5185\u7f6e\u7684iris<\/code>\u6570\u636e\u96c6\u3002<\/li>\n
            • \u6570\u636e\u9884\u5904\u7406\uff1a\u5bf9\u6570\u636e\u8fdb\u884c\u7f29\u653e\u4ee5\u6d88\u9664\u4e0d\u540c\u7279\u5f81\u95f4\u7684\u5c3a\u5ea6\u5dee\u5f02\uff0c\u7136\u540e\u8fdb\u884cPCA\u8ba1\u7b97\u3002<\/li>\n
            • \u63d0\u53d6PCA\u5206\u6570\u5e76\u5c06\u5b83\u4eec\u4e0e\u539f\u59cb\u6570\u636e\u96c6\u7684Species\u5217\u7ec4\u5408\u3002<\/li>\n
            • \u4f7f\u7528plotly<\/code>\u5e93\u521b\u5efa3D\u6563\u70b9\u56fe\uff0c\u5e76\u4e3a\u6bcf\u4e2a\u5750\u6807\u8f74\u6dfb\u52a0\u6807\u9898\u3002<\/li>\n
            • \u4f7f\u7528webshot<\/code>\u5e93\u5c06\u751f\u6210\u76843D PCA\u56fe\u5f62\u4fdd\u5b58\u4e3aPDF\u6587\u4ef6\u3002<\/li>\n<\/ul>\n

              \u6267\u884c\u8fd9\u4e9b\u4ee3\u7801\u540e\uff0c\u60a8\u5c06\u770b\u5230\u4e00\u4e2a\u7f8e\u89c2\u7684\u4e09\u7ef4PCA\u56fe\u5f62\uff0c\u5e76\u5c06\u5176\u4fdd\u5b58\u4e3a\u540d\u4e3a3D_PCA_plot.pdf<\/code>\u7684PDF\u6587\u4ef6<\/p>\n

               <\/p>\n

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

              \u4e09\u7ef4\u4e3b\u6210\u5206\u5206\u6790\uff083D PCA\uff09\u662f\u4e00\u79cd\u7edf\u8ba1\u65b9\u6cd5\uff0c\u7528\u4e8e\u964d\u4f4e\u9ad8\u7ef4\u6570\u636e\u7684\u590d\u6742\u6027\uff0c\u540c\u65f6\u4fdd\u7559\u6570\u636e\u4e2d\u7684\u5927\u90e8\u5206\u53d8\u5f02\u3002\u901a\u8fc7\u5bf9\u539f\u59cb […]<\/p>\n","protected":false},"author":111,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[9],"tags":[],"class_list":["post-447","post","type-post","status-publish","format-standard","hentry","category-r"],"_links":{"self":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts\/447","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/users\/111"}],"replies":[{"embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=447"}],"version-history":[{"count":0,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts\/447\/revisions"}],"wp:attachment":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=447"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=447"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=447"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}