笔记—R语言做相关气泡图

zz/2024/5/21 21:58:59

library(corrplot)
data <- read.table(file.choose(), header = T,sep = '\t')
new_data <- data[,-1]
?cor
??par

pr <- cor(new_data, method = "pearson")
pr1 <- cor(x= new_data[1:10],y = new_data[11:18], method = "pearson")
pr2 <- cor(x = data[1:3],y = data[4:9])
View(pr1)
res <- cor.mtest(data, method = "pearson",conf.level = 0.95)
View(res)
col1 <- colorRampPalette(c("#7F0000", "red", "#FF7F00", "yellow", "white",
                           "cyan", "#007FFF", "blue","#00007F"))
col2 <- colorRampPalette(c( "#053061","#2166AC","#4393C3","#92C5DE","#D1E5F0",
                           "#FFFFFF","#FDDBC7","#F4A582","#D6604D","#B2182B","#67001F" ))
col3 <- colorRampPalette(c("red", "white", "blue"))
col4 <- colorRampPalette(c("#7F0000", "red", "#FF7F00", "yellow", "#7FFF7F",
                           "cyan", "#007FFF", "blue", "#00007F"))

#添加显著性
#首先构建cor.mtest函数
cor.mtest <- function(mat, ...) {
  mat <- as.matrix(mat)
  n <- ncol(mat)
  p.mat<- matrix(NA, n, n)
  diag(p.mat) <- 0
  for (i in 1:(n - 1)) {
    for (j in (i + 1):n) {
      tmp <- cor.test(mat[, i], mat[, j], ...)
      p.mat[i, j] <- p.mat[j, i] <- tmp$p.value
    }
  }
  colnames(p.mat) <- rownames(p.mat) <- colnames(mat)
  p.mat
}
pv <- cor.mtest(new_data)
pv1<- pv[c(1:10),c(11:18)]

corrplot(pr1, tl.cex = 0.75,method = "ellipse", tl.col = "black",
         cl.ratio    = 0.1,mar = c(1.1,1.1,2.1,2.1), p.mat=pv1, sig.level = c(0.001,0.01,0.05), insig = "label_sig",
         pch.cex = 1.5, tl.srt = 90,col = col2(100),)

 


http://www.ngui.cc/zz/2700893.html

相关文章

如何查看 安卓证书 的签名

如何查看 安卓证书 的签名 自有安卓证书的签名查看方法 1&#xff09;通过命令查看 电脑上要装有Java 找见Java目录下的keytool.exe 打开运行&#xff0c;输入cmd&#xff0c;打开命令提示符&#xff0c;进入Java所在的盘 通过 cd 命令进入keytool.exe所在的文件夹 输入keyto…

vue 使用 swiper 实现轮播的那些事

首先运行 npm下载 npm install swiper --save-dev在需要用到的页面中 <template><div class"banner"><div class"swiper-container"><div class"swiper-wrapper"><div class"swiper-slide"><img s…

scroll-behavior属性:滚动框指定滚动行为及锚点

1、scroll-behavior属性:auto(滚动条立即滚动) / smooth(窗口平稳滚动) / inherit / initial / unset 2、锚点 <!DOCTYPE html> <html lang="en"> <head><meta charset="UTF-8"><title>锚点平滑跳转</title>&l…

如何搭建vue3.0脚手架

一 创建vue3.0脚手架 1、如果已经安装过vue2.0,要先卸载vue2.0,没有安装过的话可以忽略这一步(也可以不卸载vue2.0,这样就可以创建不同版本的项目啦!) npm uninstall vue-cli -g2、安装vue3.0,注意:Node 版本升级至 8.9 及以上 npm install -g @vue/cli3、vue3.0创建…

ThinkPHP下搭建后台(H-ui.admin)

一、下载 下载这里就不多说了&#xff0c;自行去官网下载即可 二、安装 ThinkPHP直接解压到&#xff0c;网页根目录&#xff08;我的是www目录下&#xff09; H-ui.admin先解压在一个文件夹中&#xff0c;稍后慢慢接入 三、搭建 首先要保证ThinkPHP搭建成功&#xff0c;如果是在…

MapReduce中词频统计简单实现

输入数据 test1 test2 test2 test3 test3 test3 test5 test5 test4 test5 test5 test4 test5 test4 test4 Maven必须配置 注意:Windos本地运行需要确定本地有Hadoop依赖并确保和Pom配置文件中版本一致,WordCountDriver中第6点输入输出需要自行修改 <?xml version=&qu…

MapReduce 订单和商品pid替换(Map Join)基本实现

思路简介 首先在Map阶段初始化时读取替换表的内容存入内存中,然后在Map阶段map方法时将其进行替换。所以只需要Map阶段即可,不需要Reduce阶段。 输入数据 order.txt 订单表数据(间隔:\t) 订单id 商品id 数量 1001 01 1 1002 02 2 1003 03 3 1004 01 4 1005 02 5 1006 0…

c++移动语义、万能引用、引用折叠、完美转发

想了解标题的内容&#xff0c;我们可以先简单了解下c中的左值和右值。 左值和右值 左值&#xff1a;能对表达式取地址&#xff0c;就是左值右值&#xff1a;不能对表达式取地址&#xff0c;就是右值将亡值&#xff1a;和右值引用相关的表达式&#xff0c;这样的表达式通常是将…

Ubuntu14.04安装CMake3.5.1

1、下载cmake-3.5.1.tar.gz&#xff1a;https://cmake.org/download/ 2、把 cmake-3.5.1.tar.gz放到任意临时目录(因为Cmake的安装路径默认在&#xff1a;/usr/local/中)&#xff0c;接着查看 cmake-3.5.1.tar.gz的权限&#xff08;如果权限不足&#xff0c;就不能顺利操作&…

ASP.NET 实现邮箱发送随机验证码功能(详解)

一&#xff1a;前端页面添加一个邮件发送按钮&#xff08;button控件&#xff09; <div> <asp:Button ID"Button1" runat"server" Text"发送验证码" OnClick"Button1_Click" /> </div> 二&#xff1a;后台对butt…