Skip to content

Commit 8541ca2

Browse files
updated
1 parent dcbdc61 commit 8541ca2

File tree

3 files changed

+15
-23
lines changed

3 files changed

+15
-23
lines changed

MCA/analysis_library.R

Lines changed: 13 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -210,14 +210,14 @@ calculate_cluster <- function(cluster_index,
210210
}
211211

212212
# plot the bun values
213-
plot_bun_value <- function(bun_value,annotation_info){
213+
plot_dun_value <- function(dun_value,annotation_info){
214214

215215
# Parameter in the function
216-
# bun_value: the bun values
216+
# dun_value: the bun values
217217
# annotation_info: the cell annotation information
218218

219219
# calculate the log2 values
220-
index_tmp = -log2(bun_value)
220+
index_tmp = -log2(dun_value)
221221

222222
# prepare the data for the boxplot
223223
longData = melt(as.matrix(index_tmp))
@@ -303,28 +303,28 @@ pdf_dun_tsne <- function(tissue_name){
303303

304304
cluster_n = unique(knn5)
305305

306-
bun_value = c()
306+
dun_value = c()
307307

308308
for(i in c(1:dim(annotation_info)[1])){
309309

310310
if(i %in% cluster_n){
311311

312312
tmp_val = calculate_cluster(i,as.numeric(knn5),data_tsne)
313313

314-
bun_value = cbind(bun_value,tmp_val[[1]])
314+
dun_value = cbind(dun_value,tmp_val[[1]])
315315

316316
}else{
317317

318-
bun_value = cbind(bun_value,c(1,1,1,1,1,1))
318+
dun_value = cbind(dun_value,c(1,1,1,1,1,1))
319319

320320
}
321321

322322
}
323323

324-
pp = plot_bun_value(bun_value,annotation_info)
324+
pp = plot_dun_value(dun_value,annotation_info)
325325

326326
ggsave(
327-
filename = paste0(tissue_name, "_tsne_bun_index.pdf"),
327+
filename = paste0(tissue_name, "_tsne_dun_index.pdf"),
328328
plot = pp,
329329
width = 8,
330330
height = 6
@@ -389,28 +389,28 @@ pdf_dun_tsne1 = function(tissue_name){
389389

390390
cluster_n = unique(knn5)
391391

392-
bun_value = c()
392+
dun_value = c()
393393

394394
for(i in c(1:dim(annotation_info)[1])){
395395

396396
if(i %in% cluster_n){
397397

398398
tmp_val = calculate_cluster(i,as.numeric(knn5),data_tsne)
399399

400-
bun_value = cbind(bun_value,tmp_val[[1]])
400+
dun_value = cbind(dun_value,tmp_val[[1]])
401401

402402
}else{
403403

404-
bun_value = cbind(bun_value,c(1,1,1,1,1,1))
404+
dun_value = cbind(dun_value,c(1,1,1,1,1,1))
405405

406406
}
407407

408408
}
409409

410-
pp = plot_bun_value(bun_value,annotation_info)
410+
pp = plot_dun_value(dun_value,annotation_info)
411411

412412
ggsave(
413-
filename = paste0(tissue_name, "_tsne_bun_index.pdf"),
413+
filename = paste0(tissue_name, "_tsne_dun_index.pdf"),
414414
plot = pp,
415415
width = 8,
416416
height = 6

MCA/analysis_main.R

Lines changed: 0 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -834,13 +834,3 @@ for(i in c(1:6)){
834834
}
835835

836836
# ---------------------------------------------------------------------------------
837-
838-
839-
840-
841-
842-
843-
844-
845-
846-

README.md

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,5 @@
1+
[![DOI](https://zenodo.org/badge/173669168.svg)](https://zenodo.org/badge/latestdoi/173669168)
2+
13
# SCRABBLE_PAPER
24

35
Data and analysis for the paper "Single-cell RNA-Seq imputation constrained by bulk RNA-Seq data".

0 commit comments

Comments
 (0)