calc.iS.X(SpatioTemporal)
calc.iS.X()所属R语言包：SpatioTemporal
Matrix Multiplication with Block Matrices
分块矩阵的矩阵乘法
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Computes either the product between a block diagonal, square matrix iS and a block matrix X; the quadratic form of a block diagonal, square matrix, t(X)*iS*X; or a block matrix multiplied by a vector, X*alpha.
计算块的产品之间的对角线，正方形矩阵iS和分块矩阵X;二次形式的块对角方矩阵，t(X)*iS*X;或块矩阵乘以一个向量，X*alpha的“。
用法----------Usage----------
calc.iS.X(X, iS, dim)
calc.X.iS.X(X, iS.X, dim)
calc.mu.B(dim, X, alpha)
参数----------Arguments----------
参数：X
A list of matrices with which to form the block matrix. Usually <br> mesa.data.model$X, where mesa.data.model is obtained from <br> create.data.model. The list should contain dim$m matrices, with the i:th matrix being of size dim$n - by - dim$p[i].
与形成的块矩阵的矩阵的列表。通常<BR> mesa.data.model$X，其中mesa.data.model获得<BR> create.data.model。列表应包含dim$m矩阵，与i：dim$n - - dim$p[i]个矩阵的大小。
参数：iS
A block diagonal, square matrix, with dim$m blocks each of size <br> dim$n - by - dim$n.
块对角方阵，与dim$m块，每块的大小<BR> dim$n - BY - dim$n。
参数：iS.X
Matrix containing the product of iS and X. Output from calc.iS.X.
矩阵产品的iS和X。输出从calc.iS.X。
参数：dim
Information regarding parameter dimension, output from loglike.dim.
有关参数尺寸，输出从loglike.dim。
参数：alpha
A list of dim$m vectors, with the i:th vector being of length dim$p[i].
列表，与idim$m矢量：次矢量的长度dim$p[i]。
值----------Value----------
Returns iS*X, X'*iS*X, or X*alpha.
返回的是* X，X*是* X，或X *α。
（作者）----------Author(s)----------
Johan Lindstr枚m
参见----------See Also----------
Block matrices can be created by make.sigma.B, make.sigma.B.full and <br> make.sigma.nu.
分块矩阵可以创建make.sigma.B，make.sigma.B.full和<br> make.sigma.nu。
Other block matrix functions makeCholBlock, block.mult, calc.tF.times.mat, <br> dot.prod, sumLogDiag.
其他分块矩阵的功能makeCholBlock，block.mult，calc.tF.times.mat，参考dot.prod，sumLogDiag。
This function is called by loglike.
此功能被称为loglike。
实例----------Examples----------
## Create a block diagonal matrix, ...[＃创建一个块对角矩阵，...]
iS <- rbind(c(2,1,0,0),c(1,3,0,0),
c(0,0,3,2),c(0,0,2,4))
## ... a block matrix ...[＃...分块矩阵...]
X <- list(matrix(c(1,2)),matrix(c(2,2,3,4),2,2))
## ... with alternative form, ...[＃...另一种形式，...]
Xt <- rbind(cbind(X[[1]],matrix(0,2,2)),
cbind(matrix(0,2,1),X[[2]]))
## ... and a vector alpha.[＃...和矢量α-。]
alpha <- list(c(1),c(-2,1))
## Create matching dimension structure.[＃创建匹配的维度结构。]
dim <- list(p=c(1,2), m=2, n=2)
## Compute iS * X[＃计算* X]
iS.X <- calc.iS.X(X,iS,dim)
## or[＃]
iS %*% Xt
## Compute X'* iS * X[＃计算X*是* X]
calc.X.iS.X(X, iS.X, dim)
## or[＃]
t(Xt) %*% iS %*% Xt
## Compute X* alpha[＃计算X *阿尔法。]
calc.mu.B(dim, X, alpha)
## or[＃]
cbind(X[[1]] %*% alpha[[1]], X[[2]] %*% alpha[[2]])
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-30 13:01:51
