作者Neisseria (Neisseria)
看板R_Language
标题[问题] 有一行 R 程式码看不懂...
时间Fri Feb 24 17:01:01 2017
[问题类型]:
程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
使用者(已经有用R 做过不少作品)
[问题叙述]:
以下程式码在做 quantile normalization
但 temp <- temp$ix 这行看不懂
而且我觉得怪怪的，怎麽没有除 median 的动作，是不是我误会了什麽 @@
[程式范例]:
# Perform quantile normalization
# expression.data is a matrix
running.sort <- NULL
for (i in 1:dim(expression.data)[2]) {
temp <- sort(expression.data[,i])
running.sort <- cbind(running.sort,temp)
lines(sort(expression.data[,i]),col=col.2.use[i])
}
running.median <- apply(running.sort,1,median)
running.normalized <- NULL
for (i in 1:dim(expression.data)[2]) {
temp <- expression.data[,i]
temp <- sort(temp,index.return = TRUE)
temp <- temp$ix
names(temp) <- as.numeric(unlist(running.median))
temp <- sort(temp)
temp <- as.numeric(unlist(names(temp)))
running.normalized <- cbind(running.normalized,temp)
rm(temp)
}
这程式码节录自一段分析数据的脚本，我试过，本身应该没错
如果有错，大概是我贴错了
[环境叙述]:
N/A
[关键字]:
quantile normalization
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.171.141.237
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1487926864.A.DE6.html
※ 编辑: Neisseria (1.171.141.237), 02/24/2017 17:37:34
1F：→ x88776544pc: ix是排序完的值在排序前的index位置 02/24 22:14
2F：→ x88776544pc: Code利用ix把median or mean 代回原本对应的位置 02/24 22:16