Run_JTree(treelet)
Run_JTree()所属R语言包：treelet
Treelet basis/hierarchical tree construction
Treelet基础/层次结构树的构建
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Returns information on the simultaneous construction of the treelet orthonormal basis and hierarchical tree, including which nodes were merged at each step and the basis at each specified step of the construction.
返回有关标准正交的基础和分层树treelet，包括被合并的节点在每个指定的步骤，每个步骤和基础建设的同步建设。
用法----------Usage----------
Run_JTree(X, maxlev, whichsave)
参数----------Arguments----------
参数：X
the covariance matrix of the data. For example, if using this function on genetics data to improve estimates of heritability, as in the Crossett et al arXiv paper, this argument will be the estimated additive genetic relationship matrix \hat{A}.
数据的协方差矩阵。例如，如果使用此功能的遗传学资料，以提高估计的遗传，克罗塞特等arXiv纸，这种说法是估计加性遗传关系矩阵\hat{A}。
参数：maxlev
the maximum height of the tree. This must be an integer between 1 and nrow(X)-1.
树的最大高度。这必须是一个整数1和nrow(X)-1之间。
参数：whichsave
a vector containing the levels of the tree, specified as integers between 1 and maxlev, for which you want to save the basis functions and the covariance matrix.
一个向量，包含水平的树，指定为整数1和maxlev，你要保存的基础功能和协方差矩阵。
Details
详细信息----------Details----------
This function serves as a wrapper for the functions Build_JTree and JTree_Basis, which build the hierarchical tree and calculate the basis and covariance matrix at each level, respectively.
此功能作为一个包装的功能Build_JTree和JTree_Basis，建立层次结构树和计算的基础和协方差矩阵，分别在每个级别。
值----------Value----------
a list with components
与组件的列表
参数：basis
This is a list with maxlev elements. Only those elements that are specified in the whichsave argument will be non-null entries in the list. For the non-null entries, this is the orthonormal treelet basis calculated at that level of the tree.
这是maxlev元素的列表。只有那些元素中指定的whichsave参数非空列表中的条目。对于非空条目，这是的标准正交treelet基础的树在这一水平计算。
参数：Zpos
A matrix of dimension maxlev x 2. Each row records which two nodes/clusters of the tree were combined at each step in its construction.
的矩阵尺寸maxlev×2。的每一行的记录，其中两个节点/合并树簇在其结构中的每个步骤。
参数：T
This is a list with maxlev elements, where each element is a 2x2 Jacobi rotation matrix for each step of the treelet algorithm.
这是一个maxlev元素，其中每个元素是一个2×2的Jacobi旋转矩阵每个步骤的treelet算法的列表。
参数：PCidx
A matrix of dimension maxlev x 2, where each row is a permutation of (1,2) indicating which of the two nodes/clusters merged at that step is the sum variable (value of 1) and which is the difference (value of 2).
A矩阵的维度maxlev×2，其中的每一行是一个置换(1,2)表示合并在该步骤中的两个节点/聚类的总和变量（值1）和它的区别是（值为2）。
参数：all_nodes
A matrix of dimension maxlev x nrow(X) giving node/cluster labels at each step of the treelet algorithm. A label of zero indicates a node/cluster that was merged with another node/cluster and was the difference variable.
维的矩阵maxlevXnrow(X)给节点/聚类标签在每一步的treelet算法。零标签表示节点/聚类与另一个节点/聚类和合并的差异变量。
参数：TreeCovs
This is a list with maxlev elements. Only those elements that are specified in the whichsave argument will be non-null entries in the list. For the non-null entries, this is the covariance matrix calculated at that level of the tree. The covariances in this matrix are those between the weights (orthogonal projections onto local basis vectors) in the basis expansion of the data vector.
这是maxlev元素的列表。只有那些元素中指定的whichsave参数非空列表中的条目。对于非空条目，这是在该级别的树的协方差矩阵计算。该矩阵中的协方差的权重之间（当地基矢量的正交投影）的数据矢量的基础上扩大。
（作者）----------Author(s)----------
Trent Gaugler <a href="mailto:tgaugler@andrew.cmu.edu">tgaugler@andrew.cmu.edu</a>
参考文献----------References----------
参见----------See Also----------
Build_JTree, JTree_Basis, TCS
Build_JTree，JTree_Basis，TCS
实例----------Examples----------
data(Ahat)
out=Run_JTree(Ahat,49,49)
#The information in out$Zpos[1,] and out$all_nodes[1,][的信息了$ Zpos [1]和出all_nodes上[1]]
#both show which two individuals were the first merged [两个人都显示的第一个合并]
#in the tree. The remaining rows give information[在树中。其余行提供信息]
#on subsequent merges in the tree.[在随后的合并树中。]
basis=out$basis[[49]]
cov=out$TreeCovs[[49]]
temp=basis
#This is how you can use the basis and cov output[这是如何使用的基础和覆盖输出]
#to reconstruct the estimated relationship matrix.[重建的估计关系矩阵。]
#See how close temp and the original Ahat are:[查看密切温度和原Ahat：]
Ahat1=round(Ahat,14)
temp1=round(temp,14)
sum(Ahat1!=temp1)
#In this example, we do start seeing discrepancies in the 15th digit and beyond.[在这个例子中，我们开始看到在第15位和超越的差异。]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-1 11:50:20
