Kysymys:
Voinko mallintaa teknisiä kopioita DESeq2: ssa?
Konrad Rudolph
2017-05-24 14:32:27 UTC
view on stackexchange narkive permalink

Käytän tavallisesti collapseReplicates (tai romahtaa ylävirtaan) teknisten kopioiden käsittelemiseksi.

Kuitenkin nykyisessä RNA-seq-kokeellisessa suunnittelussa näytteet sekvensoitiin kahdesti käyttämällä erilaisia ​​kirjastojen valmistusprotokollia, mikä johtaa merkittäviin eroihin tuloksena olevissa laskenta-arvioissa (itse asiassa eri protokollan valinta selittää suurimman osan variaatiosta PCA: n mukaisessa massadatassa). Siksi haluaisin mallintaa nämä erot lisäämällä ne kovariaattina DESeq2-malliin.

Voisin vain ajatella, että voisin lisätä suunnittelukauteen toisen tekijän, joka vastaa DESeq2-vinjetissä oleva pasilla-esimerkki, joka lisää tekijän type kirjastotyypille (yksi pää vs. pariliitos). pasilla-vinjettissä todetaan kuitenkin, että nämä erilaiset kirjastot ovat itse asiassa itsenäisiä biologisia kopioita, eivät teknisiä kopioita, kuten olin aiemmin olettanut.

Olen nyt huolestunut siitä, että teknisten jäljennösten suunnittelu voi lisätä keinotekoisesti tehoa. Voiko malli olla pelastettavissa?

On syytä huomata, että meillä myös on biologisia jäljitelmiä ja täydellinen sijoitus; ts. jokainen yhdistelmä hoitoa ja kirjaston valmistelua biologisten kopioiden kanssa.

Kaksi vastused:
Devon Ryan
2017-05-24 15:20:26 UTC
view on stackexchange narkive permalink

Käytännössä teknisiä kopioita ei ole mahdollista sisällyttää kyseiseen suunnitteluun (ainakin DESeq2: een). Huolesi tehon paisuttamisesta on täsmälleen oikea ja ainoa tapa torjua sitä olisi lisätä paritustekijä, kuten tapaus- tai kasvain-normaalitutkimuksissa. Eli jotain:

  ryhmäkirjastoPrep näyte1 WT A 12 WT B 13 MUT A 24 MUT B 25 WT A 36 WT B 37 MUT A 48 MUT B 4  

Tässä näyte yhdistää tekniset kopiot yhteen. Tämä on kuitenkin joko puutteellista tai lopulta informatiivisempaa.

A_Skelton73
2017-05-24 16:09:19 UTC
view on stackexchange narkive permalink

Jos ne ovat todella teknisiä kopioita, niin ei ole mitään tapaa mallintaa niitä DESeq2: n avulla *, kuten olet viitannut toimintoon collapseReplicates . DESeq2 / Mike Loven yleinen suositus collapseReplicates -toiminnolla on vain lisätä lukemat yhteen teknisiä kopioita varten.

Jos haluat mallintaa ne romahtamisen sijaan, voit voom muuntaa tietosi ja seurata Limma-käyttöoppaan.

(Koodiesimerkki kohdasta 11.3):

Suunnittelu

  | Tiedostonimi | Cy3 | Cy5 || --------- |: -----: | ----: || Tiedosto1 | wt1 | mu1 || Tiedosto2 | wt1 | mu1 || Tiedosto3 | wt2 | mu2 || Tiedosto4 | wt2 | mu2 |  

Esimerkki

  > biolrep <- c (1, 1, 2, 2) > corfit <- duplicateCorrelation (MA, ndups = 1, lohko = biolrep) > fit <-lmFit (MA, block = biolrep, cor = corfit $ consensus) > fit <- eBayes (fit) > topTable (fit,  

* Muokkaa: Kuten @Devon Ryan mainitsee, voit suunnitella pariliitetyn mallin ryhmässä DESeq2 , mutta varo tekemästä tätä täyttä sijoitusta.

Tämä vastaus saa minut ajattelemaan, että tarvitsemme taulukon asettelutukea tällä sivustolla: ominaisuuksien / laskelmien /… taulukot ovat täällä erittäin yleisiä.
@KonradRudolph - Hyväksyi, Oli sääli, kun huomasin, että markdown-taulukoita ei tuettu!


Tämä Q & A käännettiin automaattisesti englanniksi.Alkuperäinen sisältö on saatavilla stackexchange-palvelussa, jota kiitämme cc by-sa 3.0-lisenssistä, jolla sitä jaetaan.
Loading...