zmia 大佬有话说 :
请问有没有会R语言的大佬 一个脚本小修改
本帖最后由 zmia 于 2020-7-3 18:06 编辑
#### table1 —-
setwd(tabPath)
# https://cran.r-project.org/web/packages/table1/vignettes/table1-examples.html
# 变量名设置
labels = list(
# variables = list(age_at_initial_pathologic_diagnosis = "Age (year)",
# gender = "Gender",
# relative_family_cancer_history = "Family history of cancer",
# pathologic_stage = "TNM stage",…),
variables = labs_list, # 等同于上示设置
groups = list("", paste0(symbol, " Expression"), "")) # 表头大组
# 设置函数用于添加p值列
rndr = function(x, name, …) {
# 将data_for_table1替换为真实的数据
if (length(x) == 0) {
y <- data_for_table1[]
s <- rep("", length(render.default(x=y, name=name, …)))
if (is.numeric(y)) {
p <- t.test(y ~ data_for_table1[,symbol])$p.value# symbol为用于分组的基因名
} else {
p <- chisq.test(table(y, droplevels(data_for_table1[,symbol])))$p.value
}
s <- sub("<", "<", format.pval(p, digits=3, eps=0.001))
s
} else {
render.default(x=x, name=name, …)
}
}
rndr.strat <- function(label, n, …) {
ifelse(n==0, label, render.strat.default(label, n, …))
}
# 设置分组
levels(data_for_table1[,symbol]) = c("High","Low","P-value")
strata = c(list(Total = data_for_table1), split(data_for_table1, data_for_table1[,symbol]))
# 显示表格(可直接拷贝到word里修改, 或在Viewer窗口Export -> Save as Web Page…)
tab1 = table1(strata, labels, groupspan = c(1, 2, 1), droplevels = F, render = rndr, render.strat = rndr.strat)
需要从 https://cran.r-project.org/web/packages/table1/vignettes/table1-examples.html 中拷贝与原脚本相同的 rndr 和 rndr.strat 函数替换掉上面代码中的函数,并修改函数体内相应的变量名称。
loc大佬多,特来求助,本人对代码一窍不通,10元小红包不成敬意
heihai 大佬有话说 :
我也对代码一窍不通https://cdn.jsdelivr.net/gh/hishis/forum-master/public/images/patch.gif
minsc 大佬有话说 :
Rust?
营养快线 大佬有话说 :
别说钱的事儿,会的看到了肯定会琢磨, loc的技术大佬从来不含糊。你加上十块钱这题就无解了