R语言日常笔记（4）修改基因最大表达值

rm(list=ls()) 

setwd('D:\\work\\F1\\mut')

load('mRNA_exprSet.Rdata')


library(dplyr)
library(tidyr)

library(tcltk)



u <- 1:2000
plot.new()
pb <- tkProgressBar("进度","已完成 %",  0, 100)


for (i in 1:nrow(mRNA_exprSet) ) { 
  mRNA_exprSet[i,which.max(mRNA_exprSet[i,])] <- NA
  info <- sprintf("已完成 %d%%", round(i*100/nrow(mRNA_exprSet)))  
  #进度条就是根据循环里面的i来看看循环到哪一步了
  setTkProgressBar(pb,
  i*100/nrow(mRNA_exprSet),
  sprintf("进度 (%s)", info), info)
}



u <- 1:2000
plot.new()
pb <- tkProgressBar("进度","已完成 %",  0, 100)


for (i in 1:nrow(mRNA_exprSet) ) { 
  mRNA_exprSet[i,which.max(mRNA_exprSet[i,])] <- NA
  info <- sprintf("已完成 %d%%", round(i*100/nrow(mRNA_exprSet)))  
  #进度条就是根据循环里面的i来看看循环到哪一步了
  setTkProgressBar(pb,
                   i*100/nrow(mRNA_exprSet),
                   sprintf("进度 (%s)", info), info)
}




u <- 1:2000
plot.new()
pb <- tkProgressBar("进度","已完成 %",  0, 100)


for (i in 1:nrow(mRNA_exprSet) ) { 
  mRNA_exprSet[i,which.max(mRNA_exprSet[i,])] <- NA
  info <- sprintf("已完成 %d%%", round(i*100/nrow(mRNA_exprSet)))  
  #进度条就是根据循环里面的i来看看循环到哪一步了
  setTkProgressBar(pb,
                   i*100/nrow(mRNA_exprSet),
                   sprintf("进度 (%s)", info), info)
}




library(tcltk)
u <- 1:2000
plot.new()
pb <- tkProgressBar("进度","已完成 %",  0, 100)

for(i in 1:nrow(mRNA_exprSet) ){
  mRNA_exprSet[ i, 
                is.na(mRNA_exprSet[i,])] <- rowMeans(mRNA_exprSet,
                                                     na.rm = T)[i]
  info <- sprintf("已完成 %d%%", round(i*100/nrow(mRNA_exprSet)))  
  #进度条就是根据循环里面的i来看看循环到哪一步了
  setTkProgressBar(pb,
                   i*100/nrow(mRNA_exprSet),
                   sprintf("进度 (%s)", info), info)
}