注：日常工作中常遇到需要整理重复测量方差分析表，该表需结合t检验、单因素方差分析及重复测量方差分析，若分步骤处理较为麻烦，因而我整理了一个自定义函数可快速生成该表格 汇总结果表 自定义函数 library(tidyverse) #数据处理 library(broom) #整理方差分析结果 library(data.table) #数据合并 perf…
注：所用r包rcompanion、corrplot 相关性分析 library(rcompanion) cor_result <- rcorr(as.matrix(adam[,c(5,7,9,11,13)]), type = "spearman") cor_matrix <- cor_result$r #相关性矩阵 p_valu…
注：所用r包 survminer 准备工作 logrank示例数据下载 library(openxlsx) #读取数据 library(survival) #生存分析 library(survminer) #结果可视化 tt <- read.xlsx("logrank示例数据.xlsx") logrank检验和cox生存分析 log_rank…
所用R包: forestplot , forestploter 示例数据 森林图示例数据下载 单列森林图（forestplot包） 官方文档：Introduction to forest plots #forestplot包森林图代码附注释 library(forestplot) data <- read.xlsx("森林图示例数据.xlsx…
所用R包: ggplot2 1.数据准备 # 定义一个函数，用于生成每个变量的绘图所用数据 generate_test_data <- function(variable) { # 统计每列中大于 1、2、3。。。。的个数 count_values <- function(column) { c( sum(column < 0, …
所用R包: tidyverse 1.所需数据计算 t1 <- tibble(id=rep(1:25,each=4),group=1,period=rep(c(1:4),25),value=rnorm(n = 100, mean = 2, sd = 1)) t2 <- tibble(id=rep(1:25,each=4),group=2,…
注：整理了日常工作中，对于数据变量添加标签的需求 所用R包：openxlsx, Hmisc, haven R语言添加标签（stata导出可用） # 简单添加 var_labels <- () #给出标签向量 data <- imap(data, ~ { label(data[[.y]]) <- var_labels…
注：此文章整理了常用的统计图绘图代码 所用r包：ggplot2 饼图 p <- ggplot(pdata2,aes(x = factor(1),percent,fill=type))+geom_col(colour = "white")+ coord_polar(theta = "y", start = 0) + geom_text(aes(…
注：在日常使用r语言保留数据的小数位数时，常会碰到小数位数强制取整的情况，如2.00会保存为2，本文介绍了我在遇到此问题时的处理方法 所用r包：tidyverse #将所有数值型变量的小数位数统一 format_column <- function(column) { if (is.numeric(column)) { # 识别每列的原始小数…
注：本章内容介绍了excel中常用的查询填补法，该方法相较于Excel中的功能更为强大，自定义程度高，可规定按照查询到的值替换原数据或仅替换原数据的缺失值。 数据集比较 #修改数据格式使对应 bcdata <- bcdata %>% mutate(across(all_of(variables),~as.numeric(as.chara…