cor.unbalance(CORREP)
cor.unbalance()所属R语言包:CORREP
Multivariate Correlation Estimator (Unequal Number of Replicates)
多元相关估计(复制不等数)
译者:生物统计家园网 机器人LoveR
描述----------Description----------
cor.unbalance estimates correlation from replicated data of unequal number of replicates. different from cor.balance, cor.unbalance takes a pair of variables at a time because of unequal number of replicates. the variance of each row of the data MUST equal to 1 (see example below)
cor.unbalance估计从数据复制复制不平等的相关性。不同的cor.balance,cor.unbalance因为不平等的数量,时间发生在一对变量复制。每一行数据的方差必须等于1(见下面的例子)
用法----------Usage----------
cor.unbalance(x, m1, m2)
参数----------Arguments----------
参数:x
data matrix, column represents samples (conditions), and row represents variables (genes), see example below for format information
数据矩阵,列代表样本(条件),行代表变量(基因),见下面的例子的格式信息
参数:m1
number of replicates for one variable (gene)
数复制一个变量(基因)
参数:m2
number of replicates for another variable (gene)
数复制另一个变量(基因)
Details
详情----------Details----------
The multivariate correlation estimator assumes replicated omics data are iid samples from the multivariate normal distribution. It is derived by maximizing the likelihood function. Note that the off-diagonal elements in the returned correlation matrix (G by G) is the average of off-diagonals of MLE of correlation matrix of a pair of variables (m1+m2 by m1+m2).
多元相关估计,假设复制组学数据是独立同分布多元正态分布的样本。它是派生的似然函数最大化。请注意,在返回的相关矩阵(由G G)的非对角线元素是一对变量(M1 + M2由M1 + M2)的相关矩阵的MLE的对角线平均。
值----------Value----------
A correlation matrix containing only one distinct correlation coefficient for the pair of variables (genes)
对变量(基因)包含只有一个明显的相关系数相关矩阵
作者(S)----------Author(s)----------
Dongxiao Zhu and Youjuan Li
参考文献----------References----------
参见----------See Also----------
cor.unbalance, cor
cor.unbalance,cor
举例----------Examples----------
library("CORREP")
d0 <- NULL
for(l in 1:10)
d0 <- rbind(d0, rnorm(8))
## The simulated data corresponds to the real-world data of 2 genes and 10 conditions, gene expression[#模拟数据对应2基因在现实世界中的数据和10条件下,基因表达]
## profiles were replicated 3 and 5 times. [#配置文件被复制的3倍和5倍。]
## Note this function can only take calculate correlation matrix between two genes at a time.[#注意:此功能只能采取两个基因之间的相关性矩阵计算一次。]
d0<- t(d0)
## This step is to make the standard deviation of each replicate equal to 1[#这一步是让每个复制的标准差等于1]
## so that we can model the covariance matrix as correlation matrix.[#所以,我们可以模拟相关矩阵的协方差矩阵。]
d0.std <- apply(d0, 1, function(x) x/sd(x))
M <- cor.unbalance(t(d0.std), m1=3, m2=5)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|