Subject: | | Error during multiple regression with covariates analysis. Part 1 |
From: | | Hector Espinos <[log in to unmask]> |
Reply-To: | | [log in to unmask][log in to unmask] <mailto:[log in to unmask]> > > <mailto:[log in to unmask] <mailto:[log in to unmask]>>> wrote: > > > > Dear Alina, > > > > > > OK, you probably read the data correctly. Could you show us > what your > > design matrix look like? To explore a bit further, you can > also use the > > "plot" button with the "contrast estimates" and "fitted responses" > > options. > > > > Best regards, > > Guillaume. > > > > > > On 03/04/17 16:20, Alina D wrote: > > > Hello, > > > > > > Thank you very much for replying. I did Extract data > betas > > this > > > voxel and the only thing that happens is that matlab returns > this : > > > > > > beta = > > > > > > 87.2718 > > > 265.9859 > > > > > > Before I tried to load the beta_0001.nii file and to get the > intensitied > > > of some voxels from there. > > > > > > Qse = 'beta_0001.nii'; > > > Wse = spm_vol(Qse); > > > Bse = spm_read_vols(Wse); > > > > > > I stored the coordinates of the voxels in a matrix of the form: > > > myindex3d = zeros(length(myindex),3) > > > > > > And then looped through the total number and ploted it : > > > scatter_data(ii) = > Bse(myindex3d(ii,1),myindex3d(ii,2),myindex3d(ii,3)) > > > > > > What would you suggest? > > > > > > Best wishes, > > > Alina > > > > > > On Mon, Apr 3, 2017 at 5:05 PM, Guillaume Flandin > <[log in to unmask] <mailto:[log in to unmask]> > <mailto:[log in to unmask] <mailto:[log in to unmask]>> > > > <mailto:[log in to unmask] <mailto:[log in to unmask]> > <mailto:[log in to unmask] <mailto:[log in to unmask]>>>> wrote: > > > > > > Dear Alina, > > > > > > How did you extract the data for a given voxel? What is > > displayed if, > > > when looking at results, you move the 'current voxel' (red > > arrow) to > > > your voxel of interest and do a right click while on the > glass > > brain and > > > select Extract data > betas > this voxel in the contextual > > menu that > > > opens? > > > > > > Best regards, > > > Guillaume. > > > > > > > > > On 03/04/17 12:51, Alina Dinu wrote: > > > > Dear all, > > > > > > > > I am interested in the beta values for certain activated > > voxels.I > > > extracted them from the beta_0001 file. I have a couple of > > questions: > > > > 1. Could someone please tell me if this is this the > correct > > file? > > > (given that I have a block design with only one condition) > > > > 2. After I plot them I notice that the values are too high > > (in the > > > order of 80-100) which is nonsensical for beta values. Could > > someone > > > please tell me if there is any kind of scalling done in > this file? > > > How are these values coded? > > > > 3. Could someone please tell me if this file also > includes the > > > confound mean predictor? > > > > > > > > Kind regards, > > > > Alina > > > > > > > > > > -- > > > Guillaume Flandin, PhD > > > Wellcome Trust Centre for Neuroimaging > > > University College London > > > 12 Queen Square > > > London WC1N 3BG > > > > > > > > > > -- > > Guillaume Flandin, PhD > > Wellcome Trust Centre for Neuroimaging > > University College London > > 12 Queen Square > > London WC1N 3BG > > > > > > -- > Guillaume Flandin, PhD > Wellcome Trust Centre for Neuroimaging > University College London > 12 Queen Square > London WC1N 3BG > >[log in to unmask] |
Date: | | Thu, 20 Apr 2017 19:35:37 +0200 |
Content-Type: | | multipart/mixed |
Parts/Attachments: |
|
|
|
|
I'm doing the second level analysis in SPM software.
I do the differences between groups. I have three groups (control,
patients with EHM and patients without EHM).
Now I have done a multiple regression with covariates for each group.
The covariates (9 covariates in total) are biochemistry parameters and
oxidate/nitrosative biomarkers.
I have done without problems the multiple regressions for the control
group and for the patients with EHM group.
Hovewer when I do the same with the other group (patients without
EHM)appear the follow error:
SPM12: spm_results_ui (v6200) 19:20:31 - 20/04/2017
=======================================================================
SPM computation : ...height
thresholdError using gammaln
Input must be nonnegative.
Error in spm_ECdensity (line 68)
EC(3,:) = a*exp(gammaln((v+k-2)/2)-b)*(k*t/v).^(1/2*(k-2))...
Error in spm_P_RF (line 63)
EC = spm_ECdensity(STAT,Z,df);
Error in spm_uc_peakFDR (line 83)
[P,p,Eu] = spm_P_RF(1,0,ui,df,STAT,R,n);
Error in wfu_spm_getSPM12 (line 866)
[up,Pp] = spm_uc_peakFDR(0.05,df,STAT,R,n,Zum,XYZum,u);
Error in spm_getSPM (line 7)
[SPM xSPM] = wfu_spm_getSPM12(varargin{:});
Error in spm_results_ui
Error while evaluating UIControl Callback
I attached the SPM generated and the script that I have used in this
case (in matlab format).
The contrast in the other groups are the same, only change the vector
values of the covariates and the scans (they are different for each
group).
Thank for help me in advance!
If anyone need more specific information please write me!
|
|
|
|