Winnie09
commented
4 months ago Dear Marta,
Thank you for reaching out and for your interest in using the Lamian package for your analysis. Based on the error message you've encountered (Error in tcrossprod(x, y) : requires numeric/complex matrix/vector arguments), it seems the issue arises from the input format or type of the arguments you're passing to the lamian_test() function, particularly those involved in the tcrossprod() operation within the function. This error typically occurs when tcrossprod() is called with non-numeric or incorrectly structured arguments. Here are a few suggestions to troubleshoot and resolve this issue:
Check the Expression Matrix Format: Ensure that the expression matrix (expr) you're passing is a numeric matrix. The error might stem from using data@assays[["RNA"]]@Counts directly without ensuring it's in the correct format. You've correctly converted it to a matrix with as.matrix(), but it's worth double-checking that this matrix contains only numeric values.
Verify the Design Matrix: Your design matrix construction seems correct, but ensure that it is indeed numeric and properly formatted. Sometimes, data conversion or subsetting operations might unintentionally alter the data type. It could be helpful to explicitly convert the intercept and group to numeric types if they are not already.
Pseudotime Vector: Ensure that the pseudotime vector is purely numeric and correctly named with cell names. The naming of cells in the pseudotime vector should precisely match those in the expression matrix and cell annotation dataframe.
Cell Annotation Data Frame: Double-check the cellanno dataframe to ensure that it's structured correctly with cell names and sample/group identifiers. The cell names (Cell) should match exactly those used in the expression matrix and pseudotime vector.
Debugging Tips:
Try running the lamian_test() function with a smaller subset of your data to see if the error persists. This could help isolate the issue. Print out or inspect the structures of expr, cellanno, pseudotime, and design just before calling lamian_test() to confirm they're all in the expected format. Expression Matrix in expr Parameter: There seems to be a discrepancy in your code regarding which expression matrix is used. Initially, you create exp from data@assays[["SCT"]]@DaTa, but for the lamian_test() function, you switch to data@assays[["RNA"]]@CountS. Ensure you're using the intended matrix for your analysis.
If after these checks the issue remains unresolved, it could be helpful to directly inspect the lamian_test() function's handling of the input arguments, particularly where tcrossprod() is called, to ensure compatibility with the inputs. Sometimes, directly debugging with a smaller set of data can reveal where the mismatch occurs. For example, you might also check if the values in the pseudotime are unique.
I hope these suggestions help you move forward with your analysis. Please let me know how it works. Thanks!
MartaSanchezCarbonell
Dear Lamian team,
Thank you very much for creating this tool, it looks really promising and useful for my analysis. However, I am facing an error. I am studying the gene expression over time from two conditions (WT and KO) and I would like to apply your XDE test from module 3. I performed the pseudotime analysis with Slingshot and I am trying to use the function lamian_test() with test.type = 'variable' to know the genes that are differentially expressed between conditions along pseudotime.
My working object is a Seurat object where I have the pseudotime calculated with Slingshot in a slot in metadata. I have tried to supply all the arguments following your tutorial https://winnie09.github.io/Wenpin_Hou/pages/Lamian.html and trying to give them in the correct format (either data frames, matrices or vectors). But, I am obtaining the following error: Error in tcrossprod(x, y) : requires numeric/complex matrix/vector arguments
I am pasting my code here:
Checking if all cells have pseudotime value
summary(data$slingshot_pseudotime_separated)
Creating exp matrix
exp<-as.matrix(data@assays[["SCT"]]@data)
Creating design
intercept<-c("1","1") group<-c("1","0") Genotype<-c("WT","KO") design<-data.frame(intercept,group) row.names(design) <- Genotype design<-as.matrix(design)
Creating cell annotation
Geno<-as.data.frame(data$Genotype) Cell<-row.names(Geno) Geno_1<-Geno$
data$Genotypecellanno<-data.frame(Cell,Geno_1)Creating pseudotime numeric vector with cell names
pseudotime<-data$slingshot_pseudotime_separated names(pseudotime) <-Cell
XDE test
Res <- lamian_test( expr = as.matrix(data@assays[["RNA"]]@counts), cellanno = cellanno, pseudotime = pseudotime, design = design, test.type = 'variable', testvar = 2, permuiter = 100,
This is for permutation test only.
We suggest that users use default permuiter = 100.
Alternatively, we can use test.method = 'chisq' to swich to the chi-square test.
ncores = 1 )
Could you please help me to find where is the error?
Thank you very much in advance and I wish you a nice day!
Marta