Displaying 20 results from an estimated 10000 matches similar to: "how to calculate left kronecker product?"
2012 Feb 09
1
Row-wise kronecker product with Matrix package
I'm trying to calculate the row-wise kronecker product A \Box B of two
sparse matrices A and B, and am struggling to find a quick way to do this
that takes advantage of sparseness. I thought a good idea would be to use
"rep" to construct 2 matrices of the same dimension of the end product, and
multiply these two together:
library(Matrix)
A<-Matrix(c(1,0,0,0,0,1,2,0), 2, 4)
2005 Jul 13
2
Kronecker matrix product
Hi
I want to write a little function that takes a matrix X of size
m-by-n, and a list L of length "m", whose elements are matrices all
of which have
the same number of columns but possibly a different number of rows.
I then want to get a sort of dumbed-down kronecker product in which
X[i,j] is replaced by X[i,j]*L[[j]]
where L[[j]] is the j-th of the "m" matrices. For
2008 Oct 22
4
A matrix automation problem
[I am really sorry if it is double posted, I doubt me previous post could not
reach forum due to some problem with net]
Suppose I have a matrix :
a = matrix(1:9, 3)
>From this matrix, I construct 9 additional matrices :
i = 1:9
bi = a * i
Now combining all those 9 new matrices, I construct a final metrix as :
c =
b1 b4 b7
b2 b5 b8
b3 b6 b8
I want to automate this procedure for any
2010 Nov 25
1
Request: kronecker to get a sep= argument
kronecker, with make.dimnames=TRUE uses a hardwired sep=":" in the line
tmp <- outer(dnx[[i]], dny[[i]], FUN = "paste", sep = ":")
For an application in which dimnames arise from an n-way array, where
different dimensions have
different roles, and I would like to be able to use kronecker in the form
kronecker(A, B, make.dimnames=TRUE,
2005 Dec 08
1
kronecker(... , make.dimnames=TRUE)
Hi
I'm using kronecker() with a matrix and a vector. I'm interested in
the column names that kronecker() returns:
> a <- matrix(1:9,3,3)
> rownames(a) <- letters[1:3]
> colnames(a) <- LETTERS[1:3]
> b <- c(x=1,y=2)
> kronecker(a,b,make.dimnames=TRUE)
A: B: C:
a:x 1 4 7
a:y 2 8 14
b:x 2 5 8
b:y 4 10 16
c:x 3 6 9
c:y 6 12 18
>
The
2009 Feb 16
2
solve.QP with box and equality constraints
Dear list,
I am trying to follow an example that estimates a 2x2 markov transition
matrix across several periods from aggregate data using restricted least
squares.
I seem to be making headway using solve.QP(quadprog) as the unrestricted
solution matches the example I am following, and I can specify simple
equality and inequality constraints. However, I cannot correctly specify a
constraint
2012 Jun 09
2
Matrix package loading problem "Error : object ‘kronecker’ is not exported by 'namespace:methods'"
Hi R users all ,
I have a clean install of R-2.15.0 in a win32 machine and a problem loading
Matrix 1.0-6 that looks like this:
> library(Matrix)
Loading required package: lattice
Error : object ?kronecker? is not exported by 'namespace:methods'
Error: package/namespace load failed for ?Matrix?
>
I understand that kronecker() now has an S4 generic in package methods. Is
that a
2012 Mar 12
3
how to calculate a variance and covariance matrix for a vector
Hello,
I have a vector {a, b1, b2, b3, b4}. How can I calculate the following
matrix:
var(a) cov(a, b1) cov(a, b2) cov(a, b3) cov(a, b4)
cov(a, b1) var(b1) cov(a, b2) cov(a, b3) cov(a, b4)
...
...
cov(a, b1) cov(a, b2) cov(a, b3) cov(a, b4) var(b4)
I would very appreciate your inputs. Thank you very much.
Sincerely,
Jialin Huang
[[alternative HTML version deleted]]
2011 Jun 05
2
kronecker sum
Dear All,
Could someone please suggest how to find the Kronecker sum of two 2x2
matrices,
i.e. given two matrices:
-A A
a -a
and
-B B
b -b
I need:
-A-B A B 0
a -a-B 0 B
b 0 -A-b A
0 b a -a-b
Many thanks,
Lara
[[alternative HTML version deleted]]
2008 Nov 12
2
Outer, kronecker, etc.
`outer` (and related functions like kronecker) require that their
functional argument operate elementwise on arrays. This means for
example that
outer( 1:2, 3:4, list)
or
outer(1:2,3:4,function(a,b){1})
gives an error.
Is there a version of `outer`/`kronecker`/etc. that takes arbitrary
functions and does its own elementwise mapping? In the first example
above, I'd expect the
2003 Nov 18
1
aov with Error and lme
Hi
I searched in the list and only found questions
without answers e.g.
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/19955.html
: Is there a way to get the same results with lme as
with aov with Error()?
Can anybody reproduce the following results with lme:
id<-c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5)
2010 Jan 11
2
sparseM and kronecker product_R latest version
Dear all,
I just installed the new version of R, 2.10.1, and I am currently
using the package sparseM. (I also use a 64 bit windows version)
I got a problem that I never had: when I try to multiply with a
kronecker product (%x%) two sparse matrixes I get the following
message:
Error in dim(x) <- length(x) : invalid first argument
I never had this problem with previous versions of R.
May
2003 May 21
1
callNextMethod
Hi,
I don't understand why this code doesn't work (f(b2)):
/////////////////
setClass("B0", representation(b0 = "numeric"))
setClass("B1", representation("B0", b1 = "character"))
setClass("B2", representation("B1", b2 = "logical"))
f <- function(x) class(x)
setMethod("f", "B0",
2006 Sep 30
1
Gradient problem in nlm
Hello everyone!
I am having some trouble supplying the gradient function to nlm in R for
windows version 2.2.1.
What follows are the R-code I use:
fredcs39<-function(a1,b1,b2,x){return(a1+exp(b1+b2*x))}
loglikcs39<-function(theta,len){
value<-sum(mcs39[1:len]*fredcs39(theta[1],theta[2],theta[3],c(8:(7+len))) -
pcs39[1:len] * log(fredcs39(theta[1],theta[2],theta[3],c(8:(7+len)))))
2009 Nov 29
1
Plotting observed vs. fitted values
Dear Wiza[R]ds,
I am very grateful to Duncan Murdoch for his assistance with this problem.
His help was invaluable. However, the problem has become a little more
complicated for me. Now, in each plot, I need to plot the observed and
fitted values of a supine and upright posture experiment. Here is what I
have and how far I got.
# tritiated (3H)-Norepinephrine(NE) disappearance from plasma
#
2012 May 13
2
Discrete choice model maximum likelihood estimation
Hello,
I am new to R and I am trying to estimate a discrete model with three
choices. I am stuck at a point and cannot find a solution.
I have probability functions for occurrence of these choices, and then I
build the likelihood functions associated to these choices and finally I
build the general log-likelihood function.
There are four parameters in the model, three of them are associated to
2007 Feb 25
1
Repeated measures logistic regression
Dear all,
I'm struggling to find the best (set of?) function(s) to do repeated
measures logistic regression on some data from a psychology experiment.
An artificial version of the data I've got is as follows. Firstly,
each participant filled in a questionnaire, the result of which is a
score.
> questionnaire
ID Score
1 1 6
2 2 5
3 3 6
4 4 2
...
2009 Nov 29
3
Plotting observed vs. Predicted values, change of symbols
Dear Wiz[R]ds,
I am deeply grateful for the help from Duncan Murdoch, Gray Calhoun, and
others. We are almost there. For whatever reason, I can't change the symbol
from a circle to a triangle in the upright posture plots. Any ideas? I have
included the problem in full.
# tritiated (3H)-Norepinephrine(NE) disappearance from plasma
# concentrations supine and upright
# supine
datasu <-
2013 Nov 01
2
computation of hessian matrix
below is a code to compute hessian matrix , which i need to generate 29 number of different matrices for example first element in x1 and x2 is use to generate let say matrix (M1) and second element in x1 and x2 give matrix (M2) upto matrix (M29) corresponding to the total number of observations and b1 and b2 are constant.
can some one guide me or help to implement this please. I did not
2008 Jun 05
1
nls() newbie convergence problem
I'm sure this must be a nls() newbie question, but I'm stumped.
I'm trying to do the example from Draper
and Yang (1997). They give this snippet of S-Plus code:
Specify the weight function:
weight < - function(y,x1,x2,b0,b1,b2)
{
pred <- b0+b1*x1 + b2*x2
parms <- abs(b1*b2)^(1/3)
(y-pred)/parms
}
Fit the model
gmfit < -nls(~weight(y,x1,x2,b0,b1,b2),