Question: degPatterns: Error in lower.to.upper.tri.inds(n) : 'n' must be >= 2
0
gravatar for salamandra
2.2 years ago by
salamandra330
salamandra330 wrote:

When running this:

clusters <- degPatterns(cluster_rlog, metadata = TableBJ, time = "condition", col=NULL)

I get this error:

Error in lower.to.upper.tri.inds(n) : 'n' must be >= 2

and also get a warning:

Large number of genes given. Please,make sure is not an error. NormallyOnly DE genes are useful for this function.

How to solve the error?

My data looks like this:

head(cluster_rlog)
                    BJ_1      BJ_2      BJ_3 d15_CD49f+_1 d15_CD49f+_2 d25_CD49f+_1 d25_CD49f+_2 d25_CD49f+_3
ENSG00000000003  8.416135  8.458166  8.366318     7.870605     7.977721     7.790931     7.857108     8.043796
ENSG00000000457  7.374256  7.478103  7.333263     7.574723     7.376089     7.797221     7.623953     7.572987
ENSG00000001084  8.051259  8.078769  8.149138     8.368629     8.486787     8.241352     8.401051     8.431956

                d25_CD49f+CD34+_1 d25_CD49f+CD34+_2 d25_CD49f+CD34+_3
ENSG00000000003          8.005091          8.215726          8.231341
ENSG00000000457          7.805268          7.804865          7.776632
ENSG00000001084          8.335795          8.529618          8.769456

TableBJ
          sampleName         fileName time        celltype           condition
1               BJ_1 1A_ATCACG_counts   d0              BJ               d0_BJ
2               BJ_2 2A_CAGATC_counts   d0              BJ               d0_BJ
3               BJ_3 3A_ATCACG_counts   d0              BJ               d0_BJ
4       d15_CD49f+_1 2B_ACTTGA_counts  d15 CD49fposCD34neg d15_CD49fposCD34neg
5       d15_CD49f+_2 3B_CGATGT_counts  d15 CD49fposCD34neg d15_CD49fposCD34neg
6       d25_CD49f+_1 1C_TTAGGC_counts  d25 CD49fposCD34neg d25_CD49fposCD34neg
7       d25_CD49f+_2 2C_GATCAG_counts  d25 CD49fposCD34neg d25_CD49fposCD34neg
8       d25_CD49f+_3 3C_TTAGGC_counts  d25 CD49fposCD34neg d25_CD49fposCD34neg
9  d25_CD49f+CD34+_1 1E_ACAGTG_counts  d25 CD49fposCD34pos d25_CD49fposCD34pos
10 d25_CD49f+CD34+_2 2E_GGCTAC_counts  d25 CD49fposCD34pos d25_CD49fposCD34pos
11 d25_CD49f+CD34+_3 3F_GCCAAT_counts  d25 CD49fposCD34pos d25_CD49fposCD34pos

Complete code follows:

library("DESeq2")
sampleNames <- c("BJ_1", "BJ_2", "BJ_3", "d15_CD49f+_1", "d15_CD49f+_2", "d25_CD49f+_1", "d25_CD49f+_2", "d25_CD49f+_3", "d25_CD49f+CD34+_1", "d25_CD49f+CD34+_2", "d25_CD49f+CD34+_3")

sampleFiles <- c("1A_ATCACG_counts", "2A_CAGATC_counts", "3A_ATCACG_counts", "2B_ACTTGA_counts", "3B_CGATGT_counts", "1C_TTAGGC_counts", "2C_GATCAG_counts", "3C_TTAGGC_counts", "1E_ACAGTG_counts", "2E_GGCTAC_counts", "3F_GCCAAT_counts")

time <- factor(c(rep('d0',3), rep('d15',2),rep('d25',6)))

celltype <- factor(c(rep('BJ',3),rep('CD49fposCD34neg',5),rep('CD49fposCD34pos',3)))

condition <- factor(paste0(time,'_',celltype))

TableBJ <- data.frame(sampleName = sampleNames, fileName = sampleFiles, time = time, celltype = celltype, condition=condition)

ddsBJ <- DESeqDataSetFromHTSeqCount(sampleTable = TableBJ, design= ~ condition)
ddsHTSeqBJ <- ddsBJ[rowSums(counts(ddsBJ)) > 1, ]
rld <- rlog(ddsHTSeqBJ, blind=TRUE)
dds_lrt <- DESeq(ddsHTSeqBJ, test="LRT", reduced = ~ 1)
dds_res <- results(dds_lrt, alpha = 0.05)
ddsdatres <- as.data.frame(dds_res)
ddsdatres <- ddsdatre[!is.na(ddsdatre$padj),] 
expgenes <- ddsdatres
res.sig <- ddsdatres[ddsdatres$padj < 0.05,] # select just genes < 0.05
res.sig <- res.sig[order(res.sig$padj),] # and then sort this table by p-value (smaller p-values on top):
rld_mat <- assay(rld)
cluster_rlog <-subset(rld_mat, row.names(rld_mat)%in%row.names(res.sig))

library('DEGreport')
rownames(TableBJ) <- TableBJ[,1]
clusters <- degPatterns(cluster_rlog, metadata = TableBJ, time = "condition", col=NULL)
clustering rna-seq • 969 views
ADD COMMENTlink modified 2.1 years ago • written 2.2 years ago by salamandra330
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 859 users visited in the last hour