240 發(fā)簡信
IP屬地:廣西壯族自治區(qū)
  • 120
    小的PPI網(wǎng)絡(luò)圖,R也可以搞定它

    一般來說網(wǎng)絡(luò)圖都是R語言導(dǎo)出數(shù)據(jù)放入cytoscape里面繪制捣染,今天想試試在R語言內(nèi)部去畫牵咙,確實有相應(yīng)的R包(igraph)可以做葛家,不是專門針對PPI網(wǎng)絡(luò)來的蟀伸,自己組織一下數(shù)...

  • 120
    多分組芯片數(shù)據(jù)的差異分析-升級版

    1.背景知識 今天的目的是為了從一個5分組數(shù)據(jù)中篩選出兩個分組读跷,去做后續(xù)分析。擴展了一下禾唁,5個分組的數(shù)據(jù)其實也可以直接一步到位作差異分析的效览。 數(shù)據(jù)編號:GSE54238實驗設(shè)...

  • R讀取excel文件的最佳方式

    背景 眾所周知,excel生成的文件后綴是xls或者xlsx荡短,早期讀取這樣的文件丐枉,要么就另存為csv格式,要么就用readxls或者xlsx這樣的R包掘托,可能會遇到Rjava報...

  • 文章寫得不錯瘦锹,但是有個問題

    res = results(dds2, contrast=c("condition", "control", "treat")) 這個地方應(yīng)該是存在問題的

    下面是help文件中result中contrast參數(shù)的描述

    this argument specifies what comparison to extract from the object to build a results table. one of either:

    a character vector with exactly three elements: the name of a factor in the design formula, the name of the numerator level for the fold change, and the name of the denominator level for the fold change (simplest case)

    a list of 2 character vectors: the names of the fold changes for the numerator, and the names of the fold changes for the denominator. these names should be elements of resultsNames(object). if the list is length 1, a second element is added which is the empty character vector, character(). (more general case, can be to combine interaction terms and main effects)

    a numeric contrast vector with one element for each element in resultsNames(object) (most general case)

    If specified, the name argument is ignored.

    可見,當contrast參數(shù)有三個元素時依次代表 數(shù)據(jù)中作為treat/control標簽的列的名稱, fold change的分子闪盔, 以及fold change的分母弯院。
    fold change 一般算的時treat/control,作者也在他的描述中說他算的是treat versus control泪掀,所以這里正常應(yīng)該是
    res = results(dds, contrast=c("condition", "treat", "control"))
    在作者后面分析部分
    > diff_gene_deseq2 <-subset(res, padj < 0.05 & abs(log2FoldChange) > 1)
    > head(diff_gene_deseq2)
    log2 fold change (MLE): condition control vs treat
    也展示了他這里算的log2 fold change是control vs treat

  • 120
    單細胞轉(zhuǎn)錄組分析中的各種數(shù)據(jù)結(jié)構(gòu)

    單細胞分析世界里數(shù)據(jù)結(jié)構(gòu)多種多樣量九,主流的四種數(shù)據(jù)結(jié)構(gòu)分別是Bioconductor主導(dǎo)的SingleCellExperiment适掰,Seurat中的SeuratObject格式...

  • 120
    使用R語言展示我們生信技能樹全國巡講的征程

    話說有那么一天。全能打雜·小潔需要寫一個全國巡講推廣軟文稚伍,打開原來的推文看到這樣一個路徑圖 (之前做的弯予,但是忘記該怎么做了) 兩個月過去了,又去了幾個城市个曙,圖該更新了锈嫩,這可咋...

亚洲A日韩AV无卡,小受高潮白浆痉挛av免费观看,成人AV无码久久久久不卡网站,国产AV日韩精品