摘要: 001、 pbmc <- NormalizeData(pbmc) 002、实现过程 dat <- pbmc[["RNA"]]@counts dat <- as.data.frame(dat) for (i in 1:ncol(dat)) { dat[,i] <- log1p(dat[,i]/sum( 阅读全文
posted @ 2022-06-30 20:01 小鲨鱼2018 阅读(138) 评论(0) 推荐(0) 编辑
摘要: 001、直接输入函数 library(qqman) ## 加载包 manhattan ## 直接再终端输入函数名称回车,即可显示源代码 002、在Rstudio中,在终端中选中函数, 按f2 library(qqman) manhattan ## 选中函数名, 按f2 003、fix(函数名)或者 阅读全文
posted @ 2022-06-30 19:10 小鲨鱼2018 阅读(2521) 评论(0) 推荐(0) 编辑
摘要: 001、 a <- c(3, 4, 3, 1) b <- c(1, 2, 5, 4) c <- c(3, 1, 2, 7) dat <- data.frame(a, b, c) dat colSums(dat) ## 求和 rowSums(dat) colMeans(dat) ## 平均值 rowM 阅读全文
posted @ 2022-06-30 16:43 小鲨鱼2018 阅读(1149) 评论(0) 推荐(0) 编辑
摘要: 001、 packages <- c("data.table", "ggplot2", "qqman", "detectRUNS") ## 确保包已经安装 sapply(packages, require) sapply(packages, require, character.only = TRU 阅读全文
posted @ 2022-06-30 14:39 小鲨鱼2018 阅读(106) 评论(0) 推荐(0) 编辑
摘要: 001、应用于数据框 a <- c(2, 3, 6) b <- c(5, 2, 2) c <- c(4, 5, 1) dat <- data.frame(a, b, c) dat sapply(dat, sum) ## 对列进行循环, 求和, 返回向量 002、应用于列表 a <- c(3, 4, 阅读全文
posted @ 2022-06-30 12:31 小鲨鱼2018 阅读(639) 评论(0) 推荐(0) 编辑
摘要: 001、substr str <- "werewrt" ## 测试字符串 substr(str, 1, 3) ## 提取1-3字符 substr(str, 3) ## 必须指定起始 substr(str, 2, 2) ## 提取单个字符 002、substring str <- "werewrt" 阅读全文
posted @ 2022-06-30 12:21 小鲨鱼2018 阅读(920) 评论(0) 推荐(0) 编辑
摘要: 1、awk + >>追加实现 for i in $(seq 2 $(head -n 1 a.txt | awk '{print NF}')); do cut -f 1,$i a.txt | awk 'BEGIN{sum =0 } {if(NR != 1) sum += $2} {if($1 ~ /^ 阅读全文
posted @ 2022-06-30 10:27 小鲨鱼2018 阅读(185) 评论(0) 推荐(0) 编辑