데이터 전처리

R에서 limma로 하는 차등 발현 분석

John Blischak

Instructor

독소루비신 유발 심독성 기전

  • 독소루비신 작용 기전 연구를 위한 2x2 설계:
    • 유전자형 2종: 야생형 (wt), Top2b 결손 (top2b)
    • 처리 2종: PBS (pbs), 독소루비신 (dox)
    • Zhang et al. 2012
dim(eset)
Features  Samples 
   29532       12
table(pData(eset)[, c("genotype", "treatment")])
        treatment
genotype dox pbs
   top2b   3   3
   wt      3   3
R에서 limma로 하는 차등 발현 분석

특성 검사

plotDensities(eset, 
    group = pData(eset)[, "genotype"], 
    legend = "topright")

R에서 limma로 하는 차등 발현 분석

전처리 단계

  • 로그 변환

  • 분위수 정규화

  • 필터링

R에서 limma로 하는 차등 발현 분석

검증: Top2b 박스플롯

boxplot(<y-axis> ~ <x-axis>, main = "<title>")
boxplot(<gene expression> ~ <phenotype>, main = "<feature>")
boxplot(<Top2b expression> ~ <genotype>, main = "<Top2b info>")
R에서 limma로 하는 차등 발현 분석

변동 원인 확인

  • 주성분 분석

  • limma 함수 plotMDS

  • 예상되는 클러스터 수는?

R에서 limma로 하는 차등 발현 분석

연습해 봅시다!

R에서 limma로 하는 차등 발현 분석

Preparing Video For Download...