ParkerICI / vite

R package for analysis of single-cell data using graphs
GNU General Public License v3.0
10 stars 4 forks source link

Error in run_scaffold_analysis #13

Open karoljacek opened 3 years ago

karoljacek commented 3 years ago

Hello, Thank you a lot for creating and maintaining this package. I have a problem when I try to run scaffold_analysis:

library(vite) files_list <- list.files(path = ".", pattern = "*.fcs", full = TRUE, ignore.case = TRUE) input.files <- c("group1.clustered.txt","group2.clustered.txt") col.names <- as.character(fs@frames$683_myelo.fcs@parameters@data$desc) landmarks.data <- load_landmarks(files.list = files_list, asinh.cofactor = 5, transform.data = T) run_scaffold_analysis(input.files, ref.file = input.files[1], landmarks.data = landmarks.data, col.names = col.names) ###############
Error in dimnames(x) <- dn : length of 'dimnames' [2] not equal to array extent

I have found that the problem is inside function get_distances_from_landmarks, in line: image

colnames(dd) <- tab$cellType

When I looked inside tab and dd files i found this: tab: image dd: image

So program is trying to sign one element vector tab$cellType="myelo" into colnames(dd) which requires 50 of them.

I would like to ask you for help with this issue.

I'm using R 4.0.5, Rstudio server 1.4.1006 and ubuntu 20.04 You can download files from: https://www.dropbox.com/s/dcia2xlhy11prjw/scaffold_issue.tar.gz?dl=0

Best, Karol Jacek