Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error: in normalize.quantiles(dataset0): vector types do not match in copyVector #51

Open
rln0005 opened this issue May 3, 2023 · 5 comments

Comments

@rln0005
Copy link

rln0005 commented May 3, 2023

I am having an issue using Scissor to annotate bulk RNA sequencing (of paired tumor & normal samples) with single-cell tumor sequencing.

Here is my code:

## import bulk RNA seq data
bulk <- as.matrix(read.csv("gene_count_matrix.csv"))
# convert the bulk dataset to matrix
bulk_dataset <- as.matrix(bulk_dataset)

## import bulk phenotype data
pheno <- read.table("Scissor_pheno.txt")

## import single-cell multiome seq data
sc <- Read10X_h5('C:/path/to/filtered_feature_bc_matrix.h5')
sc <- sc$`Gene Expression`
# perform seurat pre-processing 
sc_dataset <- Seurat_preprocessing(sc)

## run Scissor function 
phenotype <- pheno
tag <- c('Normal', 'Tumor')
info <- Scissor(bulk_dataset, sc_dataset, phenotype, tag=tag, alpha=0.5, family="binomial")

This results in the following error:

Error in normalize.quantiles(dataset0) : 
  vector types do not match in copyVector

I have tried several different things (including manual installation of preprocessCore, using BiocManager::install("preprocessCore", configure.args="--disable-threading", force = TRUE)
...but nothing has worked. I also had issues with the normalize.quantiles() before I coerced the bulk dataset into a matrix.

Can someone advise on if I've made an error, and/or how to mitigate this issue?

@Able-Andy
Copy link

Have you figured this problem? I got the same problem as you.

@rln0005
Copy link
Author

rln0005 commented May 8, 2023

No, I haven't been able to figure this out.

Have you figured this problem? I got the same problem as you.

@Able-Andy
Copy link

Hi,I‘ve solved this problem, maybe the forms of input bulk_dataset and phenotype are wrong, the form of phenotype is "named num".
Also, my problem is my bulk_dataset ,maybe there are some errors when I got my bulk.
You can see this 'https://www.zhblog.net/qa/object-class-r.html' ,hoping it will help.

No, I haven't been able to figure this out.

Have you figured this problem? I got the same problem as you.

@NIKO6666666
Copy link

Hi,I‘ve solved this problem, maybe the forms of input bulk_dataset and phenotype are wrong, the form of phenotype is "named num". Also, my problem is my bulk_dataset ,maybe there are some errors when I got my bulk. You can see this 'https://www.zhblog.net/qa/object-class-r.html' ,hoping it will help.

No, I haven't been able to figure this out.

Have you figured this problem? I got the same problem as you.

I can't open this link, it seems to indicate a risk of privacy breach. Can you tell me the correct phenotype format and bulk matrix format? Thank you very much

@DingtaoHu
Copy link

I am having an issue using Scissor to annotate bulk RNA sequencing (of paired tumor & normal samples) with single-cell tumor sequencing.

Here is my code:

## import bulk RNA seq data
bulk <- as.matrix(read.csv("gene_count_matrix.csv"))
# convert the bulk dataset to matrix
bulk_dataset <- as.matrix(bulk_dataset)

## import bulk phenotype data
pheno <- read.table("Scissor_pheno.txt")

## import single-cell multiome seq data
sc <- Read10X_h5('C:/path/to/filtered_feature_bc_matrix.h5')
sc <- sc$`Gene Expression`
# perform seurat pre-processing 
sc_dataset <- Seurat_preprocessing(sc)

## run Scissor function 
phenotype <- pheno
tag <- c('Normal', 'Tumor')
info <- Scissor(bulk_dataset, sc_dataset, phenotype, tag=tag, alpha=0.5, family="binomial")

This results in the following error:

Error in normalize.quantiles(dataset0) : 
  vector types do not match in copyVector

I have tried several different things (including manual installation of preprocessCore, using BiocManager::install("preprocessCore", configure.args="--disable-threading", force = TRUE) ...but nothing has worked. I also had issues with the normalize.quantiles() before I coerced the bulk dataset into a matrix.

Can someone advise on if I've made an error, and/or how to mitigate this issue?

I have solved the problem with the following code: bulk_dataset[] <- lapply(bulk_dataset, function(x) as.numeric(as.character(x)))

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants