Kronecker products.

Article View: aus.stats.s

Article #16079

#16079

From: Rolf Turner
Date: Thu, 05 Mar 1992 03:04

16 lines
844 bytes

Doug Nychka asks ``Is there an easy way to do Kronecker products in S?''

Yes:

kron<-function(a,b){
	matrix(aperm(a %o% b, c(3, 1, 4, 2)),
	nrow = nrow(a) * nrow(b), ncol = ncol(a) * ncol(b))
}


				cheers,
	
					Rolf Turner
					rolf@math.unb.ca
D

Technical Details

Message-ID: <9203050219.AA19215@halmos.math.UNB.CA>

Path: arc1-ca.pugleaf.net!sparky!uunet!munnari.oz.au!uniwa!dmsperth!dmssyd.syd.dms.CSIRO.AU!news

Show Full Headers

Complete Article Headers

X-Google-Language: ENGLISH,ASCII-7-bit
X-Google-Thread: fa931,98cab5f953b9f573,start
X-Google-Attributes: gidfa931,public
X-Google-ArrivalTime: 1992-03-05 00:04:31 PST
Path: arc1-ca.pugleaf.net!sparky!uunet!munnari.oz.au!uniwa!dmsperth!dmssyd.syd.dms.CSIRO.AU!news
From: rolf@math.unb.ca (Rolf Turner)
Newsgroups: aus.stats.s
Subject: Kronecker products.
Message-ID: <9203050219.AA19215@halmos.math.UNB.CA>
Date: 5 Mar 92 03:04:46 GMT
Sender: daemon@syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Lines: 15
X-Gateway: dmssyd.syd.dms.CSIRO.AU