DB design log (9)

May 15

One of the major part of today's work is finishing the shopping-cart and confirm-order components. The shopping-cart component need communication between jsp pages. And every item in the shopping-cart need a special name. It's a key point to find a available method to send and receive the variable. I realized this by using a enumeration class, which contains the name of the process-waiting variables.

After some adjusting effort, this method works. Although a teammate found a traditional way, using hidden input class, I just got some exposure to that method without any revising, because my method also works.

The most challenging part is the confirm-order part. In this process the error occured clarified me that there could only exist ONE ResultSet. This problem made me entangled for quite a long time, however, finally, I managed to rule out the problem and solved the error.

您可以使用R语言中的以下包来生成RNA-seq项目log2(ratios)折线图并对不同基因集群进行聚类: 1. edgeR 2. DESeq2 3. clusterProfiler 以下是一个简单的R代码示例,向您展示如何使用这些包来完成任务: #加载必要的包 library(edgeR) library(DESeq2) library(clusterProfiler) #读取RNA-seq数据 counts <- read.table("path/to/counts.csv", header=TRUE, row.names=1, sep=',') #使用edgeR包来标准化数据 y <- DGEList(counts) y <- calcNormFactors(y) design <- model.matrix(~group) y <- estimateDisp(y, design) fit <- glmQLFit(y, design) qlf <- glmQLFTest(fit, coef=2) log2_ratios <- qlf$table$logFC #绘制折线图 plot(log2_ratios, type="l") #使用DESeq2进行差异表达分析和基因聚类 dds <- DESeqDataSetFromMatrix(countData=countdata,colData=coldata,design=~batch+condition) dds <- DESeq(dds) res <- results(dds, contrast=c("condition", "treated", "control")) qvalue <- res[,"padj"] log2FC <- res[,"log2FoldChange"] heatmap.2(expr, dendrogram="both", scale="row", trace="none", labRow=rownames(expr), Colv=FALSE, hclustfun=hclust, distfun=function(x) as.dist(1-cor(t(x),method="spearman")), ColSideColors=color_bar) #使用clusterProfiler包进行GO富集分析 res_ordered <- res[order(res$pvalue),] res20 <- res_ordered[1:20,] background <- rownames(countdata) ids <- rownames(res20) enrich <- enrichGO(ids=ids,OrgDb="org.Hs.eg.db",ont="BP",pvalueCutoff=0.05,readable=T, gene2Symbol= T,background=background) barplot(enrich,showCategory=15,font.size=7,cex.names=0.7)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值