library(utils)
rforge <- "http://r-forge.r-project.org"
install.packages("estimate", repos=rforge, dependencies=TRUE)
library(estimate)
OvarianCancerExpr <- system.file("extdata", "sample_input.txt",
package="estimate")
read.table(OvarianCancerExpr)[1:4,1:4]
filterCommonGenes(input.f=OvarianCancerExpr,
output.f="OV_10412genes.gct",
id="GeneSymbol")
estimateScore(input.ds = "OV_10412genes.gct",
output.ds="OV_estimate_score.gct",
platform="affymetrix")
plotPurity(scores="OV_estimate_score.gct", samples="s516",
platform="affymetrix")
scores=read.table("OV_estimate_score.gct",skip = 2,header = T)
rownames(scores)=scores[,1]
scores=t(scores[,3:ncol(scores)])
scores
estimate 计算肿瘤纯度
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 泛癌研究之使用MutSigCV看肿瘤突变异质性 发表于June 2013:Mutational heterogen...
- Hackl, H., Charoentong, P., Finotello, F., & Trajanoski, ...
- 全文1232字,阅读时长5分钟左右;买足保额,保险才真的有用。 买保险其实非常重要的一点就是: 保额一定要足够!!...