We present an algorithm for multilinear decomposition that allows for arbitrary shifts along one modality. The method is applied to neural activity arranged in the three modalities space, time, and trial. Thus, the algorithm models neural activity as a linear superposition of components with a fixed time course that may vary across either trials or space in its overall intensity and latency. Its utility is demonstrated on simulated data as well as actual EEG, and fMRI data. We show how shift-invariant multilinear decompositions of multiway data can successfully cope with variable latencies in data derived from neural activity - a problem that has caused degenerate solutions especially in modeling neuroimaging data with instantaneous multilinear decompositions. Our algorithm is available for download at www.erpwavelab.org.