This definitely saves me quite some digging!
Thank you Donald!!
Bruno
On 28/01/2014 17:42, MCLAREN, Donald wrote:
> Not sure if anyone has the code pre-compiled, but all you need is:
>
> xX.xKXs = spm_sp('Set',spm_filter(xX.K,W*xX.X)); % KWX
>
> xX.xKXs.X = full(xX.xKXs.X);
>
> xX.pKX = spm_sp('x-',xX.xKXs);
>
> KWY = spm_filter(xX.K,W*Y);
>
> beta = xX.pKX*KWY;
>
>
> Hope this helps.
>
>
> Best Regards, Donald McLaren
> =================
> D.G. McLaren, Ph.D.
> Research Fellow, Department of Neurology, Massachusetts General Hospital and
> Harvard Medical School
> Postdoctoral Research Fellow, GRECC, Bedford VA
> Website: http://www.martinos.org/~mclaren
> Office: (773) 406-2464
> =====================
> This e-mail contains CONFIDENTIAL INFORMATION which may contain PROTECTED
> HEALTHCARE INFORMATION and may also be LEGALLY PRIVILEGED and which is
> intended only for the use of the individual or entity named above. If the
> reader of the e-mail is not the intended recipient or the employee or agent
> responsible for delivering it to the intended recipient, you are hereby
> notified that you are in possession of confidential and privileged
> information. Any unauthorized use, disclosure, copying or the taking of any
> action in reliance on the contents of this information is strictly
> prohibited and may be unlawful. If you have received this e-mail
> unintentionally, please immediately notify the sender via telephone at (773)
> 406-2464 or email.
>
>
> On Tue, Jan 28, 2014 at 12:35 PM, Bruno L. Giordano <[log in to unmask]
> <mailto:[log in to unmask]>> wrote:
>
> Hello,
>
> I am planning to run some permutation analysis by permuting
> condition labels in the first-level model.
>
> I wonder whether someone managed to hack the first-level model
> routines so that they take in input and give in output variables
> stored in the workspace rather than stored on the disk (I would
> actually need only the parameter estimates, and not the stats, which
> should make things much easier). I suspect that in the case of a
> permutation analysis disk access is going to be a considerable
> limiting factor for the speed of the computations, and I wouldn't
> like to overload the parallel grid with unnecessary overheads.
>
> Thank you for any tip.
>
> Bruno
>
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~__~~
> Bruno L. Giordano, PhD
> Institute of Neuroscience and Psychology
> 58 Hillhead Street, University of Glasgow
> Glasgow, G12 8QB, Scotland
> T +44 (0) 141 330 5484 <tel:%2B44%20%280%29%20141%20330%205484>
> Www: http://www.brunolgiordano.net
>
>
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Bruno L. Giordano, PhD
Institute of Neuroscience and Psychology
58 Hillhead Street, University of Glasgow
Glasgow, G12 8QB, Scotland
T +44 (0) 141 330 5484
Www: http://www.brunolgiordano.net
|