From minc-users@bic.mni.mcgill.ca Fri Apr 1 05:07:04 2005 From: minc-users@bic.mni.mcgill.ca (=?ISO-8859-1?Q?Christina_Ro=DFmanith?=) Date: Fri Apr 1 05:07:04 2005 Subject: [MINC-users] New dcm2mnc test code for Linux In-Reply-To: References: Message-ID: <424D1DA0.5000500@neuro.ma.uni-heidelberg.de> Hi, I've tested your software for our data. Do I have to worry about the following WARNING messages? Christina Rossmanith ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Checking file types... assuming files are Syngo DICOM (CD/Export). Parsing 568 files |<--------------------------------------------------->| Sorting 568 files... Done sorting files. Processing files, one series at a time... -Parsing series info |<--------------------------------------------------->| -Creating minc file |<--------------------------------------------------->| -Parsing series info |<--------------------------------------------------->| -Creating minc file |<--------------------------------------------------->| -Parsing series info |<--------------------------------------------------->| WARNING: calculated slice width (-1.0000000000) disagrees with file's slice width (1.0000000000) -Creating minc file |<--------------------------------------------------->| -Parsing series info |<--------------------------------------------------->| WARNING: calculated slice width (-1.2500000000) disagrees with file's slice width (1.2500000000) -Creating minc file |<--------------------------------------------------->| -Parsing series info |<--------------------------------------------------->| -Creating minc file |<--------------------------------------------------->| -Parsing series info |<--------------------------------------------------->| -Creating minc file |<--------------------------------------------------->| Done processing files. Robert VINCENT wrote: >Hi all, > >I've got a new binary release of my dcm2mnc converter available for >Linux at: > >http://www.bic.mni.mcgill.ca/~bert/dcm2mnc-2.0.002-i386-linux.tar.gz > >If you have any "problem data", please give this new code a try. I have >added some support for Philips and GE proprietary fields to this version. >It should handle functional data from either manufacturer, and may handle >PET as well as MRI data in some cases. > >Please let me know of any interesting successes or problems! > >I'm still working towards a source release of this code sometime >reasonably soon. Ultimately it will be folded into the MINC >distributions. > >Thanks, > > -bert > > >_______________________________________________ >MINC-users@bic.mni.mcgill.ca >http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > > From minc-users@bic.mni.mcgill.ca Fri Apr 1 21:02:05 2005 From: minc-users@bic.mni.mcgill.ca (Roger Luechinger) Date: Fri Apr 1 21:02:05 2005 Subject: [MINC-users] Dicom to Minc converter for Phillips ? In-Reply-To: References: <1110281812181263.19332@webmail> <42411FDC.1010605@biomed.ee.ethz.ch> Message-ID: <42490BA6.6030607@biomed.ee.ethz.ch> Hi Kuan, Kuan H. Kho wrote: > I think you're right about the slice info: that's very important to > most imagers. As we make whole 3D scans instead of slices (both > structural and functional images), the simple rawtominc conversion > suffices. But I indeed did not realize that there may be a bunch of > people that indeed do need that information! Your right for most people the minc-File without the geometry or at least with identical orientation and origin for the fmri and Anatomical Volume is fine. > About a rectominc tool: I can work with the rawtominc, but I heard > that on the newest Phillips scanner software (release ..? I can't keep > count anymore), the PAR/REC format changed a bit. I could be wrong, > but one may have to check that. There is a minor change for the par/rec-File between Scanner Software-Releaase 9.x and 10.x (or in other word if you scanner host computer switch from OpenVMS (<=9) to Windows (>=10)) Rel. <=9.x can export Par Version 3 Files Rel >=10 export Par Version 4 Files But only the par-File with the additional information changes, not the raw-Data in the rec-File. Therefore the rawtominc will work fine. (Up to Rel 1.2/Rel11 which is current). In the par-Files some information like recon-Size moved from the upper (scan-related part) to the lower part with slice related infos. Regards Roger -- ^ ^ o o ------------------------OOO----(_)----OOO--------------------- Dr. Roger Luechinger E-Mail: luechinger@biomed.ee.ethz.ch Universitaetsspital MR-Zentrum VMR30 www/PGP: http://www.mr.ethz.ch/~rluchin/ CH-8091 Zuerich Tel: ++41 1 255 30 64 Fax: ++41 1 255 45 06 -------------------------oooO-------Oooo---------------------- From minc-users@bic.mni.mcgill.ca Wed Apr 6 16:31:04 2005 From: minc-users@bic.mni.mcgill.ca (Roy Yves) Date: Wed Apr 6 15:31:04 2005 Subject: [MINC-users] register, regions of interests and Tailarach coordinates Message-ID: Hello: Is there a tool, or functionality in register or elsewhere that will allow me to: 1- Is there a tool that, given a set of regions of interests identified using their names or Broadmann number, and a threshold for the t value, will search a t-stat map and return the Tailarach/MNI coordinates of the peaks that are present in those regions, together with the value of the peak? 2- Is there a tool that takes as input a region of interest identified using their names or Broadmann number and return the boundaries of that regions in Tailarach/MNI coordinates? 3- Is there a tool that takes as input Tailarach/MNI coordinates and return their names and/or Broadmann number? 4- Is there a functionality in register that will facilitate that task of performing 1, 2 and 3. If there is no such tool, anybody would be interested in collaborating in developing one? Many thanks. Yves Roy From minc-users@bic.mni.mcgill.ca Thu Apr 7 00:54:04 2005 From: minc-users@bic.mni.mcgill.ca (Jeni Chen) Date: Wed Apr 6 23:54:04 2005 Subject: [MINC-users] register, regions of interests and Tailarach coordinates In-Reply-To: Message-ID: Hi Yves, On Wed, 6 Apr 2005, Roy Yves wrote: > Hello: > > Is there a tool, or functionality in register or elsewhere that will allow me to: > > 1- Is there a tool that, given a set of regions of interests identified using their names or Broadmann number, and a threshold for the t value, will search a t-stat map and return the Tailarach/MNI coordinates of the peaks that are present in those regions, together with the value of the peak? > > 2- Is there a tool that takes as input a region of interest identified using their names or Broadmann number and return the boundaries of that regions in Tailarach/MNI coordinates? Wouldn't Keith's "stat_summary" in the fmristat do this? There's also 'find_peaks' that will give you all the peaks above a given threshold in a .tag file. You can then load the tag file in Register to see these peaks. Usage: find_peaks [options] .mnc .tag find_peaks -h > 3- Is there a tool that takes as input Tailarach/MNI coordinates and return their names and/or Broadmann number? There is a program called Talairach Daemon (http://ric.uthscsa.edu/projects/talairachdaemon.html) It will read tab or space delimited records from text files containing lists of Talairach coordinates arranged in x-y-z order, and gives you the name or Broadman area number in a text file. There is also an option of doing a single search by entering the coordinates in a dialog window. > 4- Is there a functionality in register that will facilitate that task of performing 1, 2 and 3. I guess you can use 'find_peaks' to get the tag file (a text file) and import it into Talairach Daemon to get the names that correspond to those coordinates. cheers, Jeni ***************************** Jen-I Chen Agente de recherche Université de Montréal Tel: 514-343-6111 ext. 4359 Fax: 514-343-2111 ***************************** From minc-users@bic.mni.mcgill.ca Thu Apr 7 18:55:05 2005 From: minc-users@bic.mni.mcgill.ca (roberto) Date: Thu Apr 7 17:55:05 2005 Subject: [MINC-users] classify_clean Message-ID: <42541E6E.6020401@gmail.com> Hi there! (this is my first e-mail to the list... ) 1. how does classify_clean work with multiple input volumes? (a t1, t2 and pd) it would be nice if there was an example of it in the wiki... 2. also, I'm having problems running nu_correct (v1.10) with some volumes (which appear to have long headers). Any idea? thanks! roberto From k.h.kho@azu.nl Thu Apr 7 18:56:02 2005 From: k.h.kho@azu.nl (Kuan H. Kho) Date: Thu Apr 7 17:56:02 2005 Subject: [MINC-users] register, regions of interests and Tailarach coordinates In-Reply-To: References: Message-ID: Hi, I'm not sure about such a function in register, but if I see your questions, you'd probably want to check out the "Talairach Daemon". See http://ric.uthscsa.edu/projects/talairachdaemon.html But beware, there are some pitfalls concerning MNI/Talairach - (see e.g. Brett, M., Johnsrude, I. S., & Owen, A. M. (2002). The problem of functional localization in the human brain. Nat Rev Neurosci, 3(3), 243-249. They also have some info on their CBU website) Good luck, Kuan ---- Kuan H. Kho, MD Depts. of Psychiatry & Neurosurgery, G03.124 University Medical Center Utrecht Heidelberglaan 100 P.O.Box 85500 3508 GA Utrecht The Netherlands P +31-30-2507969/2507953 F +31-30-2542100 E k.h.kho_'at'_azu.nl On Apr 6, 2005 9:30 PM, Roy Yves wrote: > Hello: > > Is there a tool, or functionality in register or elsewhere that will allow me to: > > 1- Is there a tool that, given a set of regions of interests identified using their names or Broadmann number, and a threshold for the t value, will search a t-stat map and return the Tailarach/MNI coordinates of the peaks that are present in those regions, together with the value of the peak? > > 2- Is there a tool that takes as input a region of interest identified using their names or Broadmann number and return the boundaries of that regions in Tailarach/MNI coordinates? > > 3- Is there a tool that takes as input Tailarach/MNI coordinates and return their names and/or Broadmann number? > > 4- Is there a functionality in register that will facilitate that task of performing 1, 2 and 3. > > If there is no such tool, anybody would be interested in collaborating in developing one? > > Many thanks. > > Yves Roy > > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > From minc-users@bic.mni.mcgill.ca Fri Apr 8 15:58:04 2005 From: minc-users@bic.mni.mcgill.ca (Jason Lerch) Date: Fri Apr 8 14:58:04 2005 Subject: [MINC-users] register, regions of interests and Tailarach coordinates In-Reply-To: References: Message-ID: <4f4cc754634cd1256188d3f1ef9e16dc@bic.mni.mcgill.ca> On Apr 6, 2005, at 3:30 PM, Roy Yves wrote: > Hello: > > Is there a tool, or functionality in register or elsewhere that will > allow me to: > > 1- Is there a tool that, given a set of regions of interests > identified using their names or Broadmann number, and a threshold for > the t value, will search a t-stat map and return the Tailarach/MNI > coordinates of the peaks that are present in those regions, together > with the value of the peak? 'find_peaks', I believe, will do the job that you want. Except for the Brodmann area part, which is next to impossible since we don't have real maps of most Brodmann Areas. > > 2- Is there a tool that takes as input a region of interest identified > using their names or Broadmann number and return the boundaries of > that regions in Tailarach/MNI coordinates? That would be a rather long list of numbers for a boundary, would it not? I'm not quite sure I understand what you are after. > > 3- Is there a tool that takes as input Tailarach/MNI coordinates and > return their names and/or Broadmann number? You could get the tag file produced by find_peaks and look at it in register together with a probability map generated for a certain area, or, alternately, the map of structures used by ANIMAL to do its segmentation to get some idea about the area each peak is in. > > 4- Is there a functionality in register that will facilitate that task > of performing 1, 2 and 3. The combination of tag files of peaks together with anatomical data and structural probability maps does the job pretty well for me anyways. Cheers, Jason From minc-users@bic.mni.mcgill.ca Fri Apr 8 16:02:04 2005 From: minc-users@bic.mni.mcgill.ca (Jason Lerch) Date: Fri Apr 8 15:02:04 2005 Subject: [MINC-users] classify_clean In-Reply-To: <42541E6E.6020401@gmail.com> References: <42541E6E.6020401@gmail.com> Message-ID: <6da69a4d359ba5ed3dd4531388c5a0a2@bic.mni.mcgill.ca> On Apr 6, 2005, at 1:37 PM, roberto wrote: > Hi there! (this is my first e-mail to the list... ) Welcome! > > 1. how does classify_clean work with multiple input volumes? (a t1, t2 > and pd) > it would be nice if there was an example of it in the wiki... Theoretically? Your best bet is taking a look at Zijdenbos, A. P., Forghani, R. and Evans, A. C. (2002) Automatic "pipeline" analysis of 3-D MRI data for clinical trials: application to multiple sclerosis. IEEE Trans Med Imaging 21, 1280-1291. If you mean practically, then simply append them to the command line, i.e. classify_clean -clean_tags t1.mnc t2.mnc pd.mnc output_classified.mnc > > 2. also, I'm having problems running nu_correct (v1.10) with some > volumes (which appear to have long headers). Any idea? N3's header handling has always been rather dumb - so unless someone fixed that recently, you probably have to remove the dicom part of the header. Create a script like the following #!/usr/local/bin/perl5 -w open IN, "$ARGV[0]"; open OUT, ">/tmp/mincedittmp_$$.txt"; while () { if ($_ =~ /dicom_0x.+/) { # do nothing } else { print OUT $_; } } close IN; close OUT; system("mv /tmp/mincedittmp_$$.txt $ARGV[0]"); and run 'mincedit mincfile.mnc above_script.pl' before running N3. You might of course have to adapt that script to your needs. Jason > > thanks! > roberto > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users From minc-users@bic.mni.mcgill.ca Sat Apr 9 19:52:03 2005 From: minc-users@bic.mni.mcgill.ca (Jamila Ahdidan) Date: Sat Apr 9 18:52:03 2005 Subject: [MINC-users] glim_image Message-ID: <20050409225130.62397.qmail@web60807.mail.yahoo.com> Hi all, I'm trying to understand what is exactly being done with glim_imge. I'm running the command: glim_image -t_stat "CBF_results/results_""$n"".mnc" t_stat voxel column 1 "Matrices/matrix_norm_""$n"".glim" with a glim matrix containing only the filenames and the group indication. I had previously run the same command but with column 2 with a matrix with 2 columns one for the intercept and another one for the group assignation. This last way of doing gives me bad results, with a lot of noise and a very small range of value, where nothing is significant. The other way gives me cleaner results that look like what I expected, and higher significance. My dilema is whether I should just keep my good results and forget about the intercept in the matrix, or I should stick to the intercept and conclude that nothing is interpretable from my results. I'm really desperate to have a clear explanation of what is happening and what I should choose and why. Many thanks to you, Jamila __________________________________ Do you Yahoo!? Yahoo! Small Business - Try our new resources site! http://smallbusiness.yahoo.com/resources/ From minc-users@bic.mni.mcgill.ca Sat Apr 9 21:25:04 2005 From: minc-users@bic.mni.mcgill.ca (Jason Lerch) Date: Sat Apr 9 20:25:04 2005 Subject: [MINC-users] glim_image In-Reply-To: <20050409225130.62397.qmail@web60807.mail.yahoo.com> References: <20050409225130.62397.qmail@web60807.mail.yahoo.com> Message-ID: <59eb99d481ef3ae1fa94c175f11fbe23@bic.mni.mcgill.ca> On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > My dilema is whether I should just keep > my good results and forget about the intercept in the > matrix, or I should stick to the intercept and > conclude that nothing is interpretable from my > results. Is there any reason to force the slope to be 0 at x=0? If you have standard VBM density data, then that is an invalid assumption, since there is every reason to allow the y to take on an arbitrary value at x=0, so you would include that column of ones for your intercept term. If you have different data then this assumption might be valid - something that is the case, for example, when looking at asymmetry VBM. But by and large you will want an intercept. Hope this helps, Jason From minc-users@bic.mni.mcgill.ca Sun Apr 10 04:42:05 2005 From: minc-users@bic.mni.mcgill.ca (Jamila Ahdidan) Date: Sun Apr 10 03:42:05 2005 Subject: [MINC-users] glim_image In-Reply-To: 6667 Message-ID: <20050410074117.98583.qmail@web60809.mail.yahoo.com> --0-460192004-1113118877=:98142 Content-Type: text/plain; charset=us-ascii Hi Jason, Well I have to say that I don't use glim_image to perform a VBM study, but to perform a t test at each voxel to assess the difference between my group of patients and my group of controls. So, I don't really know whether I want to force to be 0 at x=0. (I don't really know what that means!). Do you think I'm using the wrong minc command, and if yes do you have another idea? Thanks a lot, Jamila Jason Lerch wrote: On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > My dilema is whether I should just keep > my good results and forget about the intercept in the > matrix, or I should stick to the intercept and > conclude that nothing is interpretable from my > results. Is there any reason to force the slope to be 0 at x=0? If you have standard VBM density data, then that is an invalid assumption, since there is every reason to allow the y to take on an arbitrary value at x=0, so you would include that column of ones for your intercept term. If you have different data then this assumption might be valid - something that is the case, for example, when looking at asymmetry VBM. But by and large you will want an intercept. Hope this helps, Jason _______________________________________________ MINC-users@bic.mni.mcgill.ca http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users --------------------------------- Do you Yahoo!? Yahoo! Mail - Find what you need with new enhanced search. Learn more. --0-460192004-1113118877=:98142 Content-Type: text/html; charset=us-ascii
Hi Jason,
Well I have to say that I don't use glim_image to perform a VBM study, but to perform a t test at each voxel to assess the difference between my group of patients and my group of controls. So, I don't really know whether I want to force to be 0 at x=0. (I don't really know what that means!).
Do you think I'm using the wrong minc command, and if yes do you have another idea?
 
Thanks a lot,
Jamila

Jason Lerch <jason@bic.mni.mcgill.ca> wrote:

On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote:

> My dilema is whether I should just keep
> my good results and forget about the intercept in the
> matrix, or I should stick to the intercept and
> conclude that nothing is interpretable from my
> results.

Is there any reason to force the slope to be 0 at x=0? If you have
standard VBM density data, then that is an invalid assumption, since
there is every reason to allow the y to take on an arbitrary value at
x=0, so you would include that column of ones for your intercept term.
If you have different data then this assumption might be valid -
something that is the case, for example, when looking at asymmetry VBM.
But by and large you will want an intercept.

Hope this helps,

Jason

_______________________________________________
MINC-users@bic.mni.mcgill.ca
http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users

Do you Yahoo!?
Yahoo! Mail - Find what you need with new enhanced search. Learn more. --0-460192004-1113118877=:98142-- From minc-users@bic.mni.mcgill.ca Mon Apr 11 10:32:04 2005 From: minc-users@bic.mni.mcgill.ca (Jason Lerch) Date: Mon Apr 11 09:32:04 2005 Subject: [MINC-users] glim_image In-Reply-To: <20050410074117.98583.qmail@web60809.mail.yahoo.com> References: <20050410074117.98583.qmail@web60809.mail.yahoo.com> Message-ID: <8c0e01bd01076e17a3be24a16110b2be@bic.mni.mcgill.ca> --Boundary_(ID_+WcXmHmYo1zpkvYSkQNFBw) Content-type: multipart/alternative; boundary="Boundary_(ID_r58iT5LGFw7366+ZOeQ0zA)" --Boundary_(ID_r58iT5LGFw7366+ZOeQ0zA) Content-type: text/plain; charset=US-ASCII; format=flowed Content-transfer-encoding: 7BIT Greetings again, I've attached a graph which shows the difference - the data was generated using the following function: y = 2 + 0.5x + Error where x is a simple sequence between 0 and 10. There are two regression lines through the graph. The one in blue fits a model including the intercept term, the one in red fits a model without an intercept. You can see that the red line takes on a value of 0 when x=0 - and that therefore the fit is not as accurate as the blue line. None of this stuff is specific to glim_image - this is all standard linear model statistics. One of the best online references that I know about is here: http://www.itl.nist.gov/div898/handbook/ though there surely are others as well. Good luck, Jason --Boundary_(ID_r58iT5LGFw7366+ZOeQ0zA) Content-type: text/enriched; charset=US-ASCII Content-transfer-encoding: 7BIT Greetings again, I've attached a graph which shows the difference - the data was generated using the following function: y = 2 + 0.5x + Error where x is a simple sequence between 0 and 10. There are two regression lines through the graph. The one in blue fits a model including the intercept term, the one in red fits a model without an intercept. You can see that the red line takes on a value of 0 when x=0 - and that therefore the fit is not as accurate as the blue line. None of this stuff is specific to glim_image - this is all standard linear model statistics. One of the best online references that I know about is here: http://www.itl.nist.gov/div898/handbook/ though there surely are others as well. Good luck, Jason --Boundary_(ID_r58iT5LGFw7366+ZOeQ0zA)-- --Boundary_(ID_+WcXmHmYo1zpkvYSkQNFBw) Content-type: application/pdf; x-unix-mode=0644; name=sample-graph.pdf Content-transfer-encoding: quoted-printable Content-disposition: inline; filename=sample-graph.pdf %PDF-1.1=0A%=81=E2=81=E3=81=CF=81=D3 1=200=20obj=0A<<=0A/CreationDate=20(D:20050411092246)=0A/ModDate=20= (D:20050411092246)=0A/Title=20(R=20Graphics=20Output)=0A/Producer=20(R=20= 1.8.0)=0A/Creator=20(R)=0A>>=0Aendobj=0A2=200=20obj=0A<<=0A/Type=20= /Catalog=0A/Pages=203=200=20R=0A>>=0Aendobj=0A4=200=20obj=0A<<=0A= /ProcSet=20[/PDF=20/Text]=0A/Font=20<<=20/F1=206=200=20R=20/F2=207=200=20= R=20/F3=208=200=20R=20/F4=209=200=20R=20/F5=2010=200=20R=20/F6=2011=200=20= R=20>>=0A>>=0Aendobj=0A5=200=20obj=0A<<=0A/Type=20/Encoding=0A= /BaseEncoding=20/PDFDocEncoding=0A/Differences=20[=2045/minus=20= 96/quoteleft=0A144/dotlessi=20/grave=20/acute=20/circumflex=20/tilde=20= /macron=20/breve=20/dotaccent=0A/dieresis=20/.notdef=20/ring=20/cedilla=20= /.notdef=20/hungarumlaut=20/ogonek=20/caron=20/space]=0A>>=0Aendobj=0A6=20= 0=20obj=0A<<=0A/Type=20/Font=0A/Subtype=20/Type1=0A/Name=20/F1=0A= /BaseFont=20/Helvetica=0A/Encoding=205=200=20R=0A>>=0Aendobj=0A7=200=20= obj=0A<<=0A/Type=20/Font=0A/Subtype=20/Type1=0A/Name=20/F2=0A/BaseFont=20= /Helvetica-Bold=0A/Encoding=205=200=20R=0A>>=0Aendobj=0A8=200=20obj=0A<<=0A= /Type=20/Font=0A/Subtype=20/Type1=0A/Name=20/F3=0A/BaseFont=20= /Helvetica-Oblique=0A/Encoding=205=200=20R=0A>>=0Aendobj=0A9=200=20obj=0A= <<=0A/Type=20/Font=0A/Subtype=20/Type1=0A/Name=20/F4=0A/BaseFont=20= /Helvetica-BoldOblique=0A/Encoding=205=200=20R=0A>>=0Aendobj=0A10=200=20= obj=0A<<=0A/Type=20/Font=0A/Subtype=20/Type1=0A/Name=20/F5=0A/BaseFont=20= /Symbol=0A>>=0Aendobj=0A11=200=20obj=0A<<=0A/Type=20/Font=0A/Subtype=20= /Type1=0A/Name=20/F6=0A/BaseFont=20/ZapfDingbats=0A>>=0Aendobj=0A12=200=20= obj=0A<<=0A/Type=20/Page=0A/Parent=203=200=20R=0A/Contents=2013=200=20R=0A= /Resources=204=200=20R=0A>>=0Aendobj=0A13=200=20obj=0A<<=0A/Length=2014=20= 0=20R=0A>>=0Astream 1=20J=201=20j=2010=20M=20q=0AQ=20q=2059.04=2073.44=20342.72=20299.52=20= re=20W=20n=0A0.000=200.000=200.000=20RG=0A0.75=20w=0A[]=200=20d=0ABT=0A= /F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2068.77=20188.61=20Tm=20(l)=20= Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2070.36=20= 151.65=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=2071.94=20219.65=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=2073.53=20161.40=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=2075.12=20152.65=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2076.70=20136.05=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2078.29=20= 156.30=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=2079.88=20181.38=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=2081.46=20142.68=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=2083.05=20147.71=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2084.64=20174.46=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2086.22=20= 176.28=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=2087.81=20147.43=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=2089.40=20211.56=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=2090.98=20170.88=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2092.57=20189.66=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=2094.16=20= 158.76=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=2095.74=20131.88=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=2097.33=20171.83=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=2098.92=20147.40=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20100.50=20142.13=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20102.09=20= 155.20=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20103.68=20135.64=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20105.26=20187.33=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20106.85=20186.49=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20108.44=20148.07=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20110.02=20= 174.73=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20111.61=20194.26=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20113.20=20208.46=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20114.78=20175.02=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20116.37=20200.11=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20117.96=20= 155.81=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20119.54=20194.04=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20121.13=20178.17=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20122.72=20191.88=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20124.30=20148.65=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20125.89=20= 169.98=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20127.48=20223.09=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20129.06=20160.63=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20130.65=20192.78=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20132.24=20170.60=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20133.82=20= 188.05=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20135.41=20174.16=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20137.00=20188.15=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20138.58=20154.46=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20140.17=20148.70=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20141.76=20= 177.18=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20143.34=20219.02=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20144.93=20193.97=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20146.52=20184.64=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20148.10=20199.72=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20149.69=20= 197.99=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20151.28=20164.75=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20152.86=20202.67=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20154.45=20225.96=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20156.04=20240.06=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20157.62=20= 168.44=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20159.21=20183.85=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20160.80=20226.05=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20162.38=20221.75=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20163.97=20187.12=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20165.56=20= 200.03=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20167.14=20211.80=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20168.73=20197.20=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20170.32=20239.27=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20171.90=20184.06=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20173.49=20= 248.66=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20175.08=20192.48=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20176.66=20172.39=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20178.25=20186.51=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20179.84=20196.31=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20181.42=20= 240.00=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20183.01=20224.24=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20184.60=20176.02=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20186.18=20241.60=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20187.77=20224.31=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20189.36=20= 235.82=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20190.94=20143.91=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20192.53=20218.00=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20194.12=20184.92=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20195.70=20228.22=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20197.29=20= 219.12=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20198.88=20207.88=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20200.46=20227.31=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20202.05=20161.94=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20203.64=20181.04=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20205.22=20= 290.24=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20206.81=20174.36=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20208.40=20232.69=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20209.98=20185.89=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20211.57=20251.06=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20213.16=20= 229.19=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20214.74=20238.69=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20216.33=20200.55=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20217.92=20233.32=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20219.50=20173.82=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20221.09=20= 235.35=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20222.68=20189.08=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20224.26=20189.28=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20225.85=20222.06=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20227.44=20236.34=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20229.02=20= 235.45=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20230.61=20230.51=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20232.20=20223.49=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20233.78=20216.62=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20235.37=20215.90=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20236.96=20= 195.91=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20238.54=20245.16=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20240.13=20233.09=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20241.72=20245.87=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20243.30=20190.45=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20244.89=20= 216.34=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20246.48=20223.09=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20248.06=20170.06=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20249.65=20219.92=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20251.24=20281.03=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20252.82=20= 235.62=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20254.41=20238.66=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20256.00=20241.41=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20257.58=20228.76=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20259.17=20225.62=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20260.76=20= 228.58=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20262.34=20213.33=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20263.93=20217.44=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20265.52=20207.56=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20267.10=20235.69=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20268.69=20= 231.52=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20270.28=20236.44=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20271.86=20234.20=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20273.45=20232.63=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20275.04=20252.53=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20276.62=20= 248.15=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20278.21=20248.39=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20279.80=20190.69=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20281.38=20240.24=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20282.97=20268.41=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20284.56=20= 224.82=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20286.14=20237.65=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20287.73=20257.25=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20289.32=20249.21=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20290.90=20273.89=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20292.49=20= 251.27=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20294.08=20271.37=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20295.66=20210.01=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20297.25=20236.89=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20298.84=20210.78=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20300.42=20= 242.84=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20302.01=20241.85=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20303.60=20244.74=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20305.18=20277.56=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20306.77=20296.69=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20308.36=20= 249.09=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20309.94=20237.33=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20311.53=20263.62=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20313.12=20254.62=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20314.70=20249.23=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20316.29=20= 294.05=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20317.88=20281.20=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20319.46=20251.98=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20321.05=20175.48=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20322.64=20275.23=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20324.22=20= 270.96=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20325.81=20255.76=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20327.40=20241.04=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20328.98=20284.75=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20330.57=20251.33=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20332.16=20= 266.85=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20333.74=20310.67=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20335.33=20231.28=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20336.92=20274.38=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20338.50=20254.67=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20340.09=20= 233.18=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20341.68=20317.28=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20343.26=20298.15=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20344.85=20312.18=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20346.44=20247.68=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20348.02=20= 293.10=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20349.61=20280.74=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20351.20=20273.46=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20352.78=20249.17=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20354.37=20275.89=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20355.96=20= 291.75=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20357.54=20292.13=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20359.13=20239.51=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20360.72=20295.09=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20362.30=20282.07=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20363.89=20= 248.71=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20365.48=20243.81=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20367.06=20265.56=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20368.65=20257.62=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20370.24=20277.04=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20371.82=20= 282.68=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20373.41=20274.24=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20375.00=20295.98=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20376.58=20285.51=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20378.17=20273.14=20Tm=20= (l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20379.76=20= 267.43=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=20= 7.48=20381.34=20311.39=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20Tf=201=20Tr=20= 7.48=200=200=207.48=20382.93=20303.41=20Tm=20(l)=20Tj=200=20Tr=0A/F6=201=20= Tf=201=20Tr=207.48=200=200=207.48=20384.52=20285.02=20Tm=20(l)=20Tj=200=20= Tr=0A/F6=201=20Tf=201=20Tr=207.48=200=200=207.48=20386.10=20332.88=20Tm=20= (l)=20Tj=200=20Tr=0AET=0AQ=20q=0A0.000=200.000=200.000=20RG=0A0.75=20w=0A= []=200=20d=0A71.73=2073.44=20m=20389.07=2073.44=20l=20S=0A71.73=2073.44=20= m=2071.73=2066.24=20l=20S=0A135.20=2073.44=20m=20135.20=2066.24=20l=20S=0A= 198.67=2073.44=20m=20198.67=2066.24=20l=20S=0A262.13=2073.44=20m=20= 262.13=2066.24=20l=20S=0A325.60=2073.44=20m=20325.60=2066.24=20l=20S=0A= 389.07=2073.44=20m=20389.07=2066.24=20l=20S=0ABT=0A0.000=200.000=200.000=20= rg=0A/F1=201=20Tf=2012.00=200.00=20-0.00=2012.00=2068.40=2047.52=20Tm=20= (0)=20Tj=0A/F1=201=20Tf=2012.00=200.00=20-0.00=2012.00=20131.86=2047.52=20= Tm=20(2)=20Tj=0A/F1=201=20Tf=2012.00=200.00=20-0.00=2012.00=20195.33=20= 47.52=20Tm=20(4)=20Tj=0A/F1=201=20Tf=2012.00=200.00=20-0.00=2012.00=20= 258.80=2047.52=20Tm=20(6)=20Tj=0A/F1=201=20Tf=2012.00=200.00=20-0.00=20= 12.00=20322.26=2047.52=20Tm=20(8)=20Tj=0A/F1=201=20Tf=2012.00=200.00=20= -0.00=2012.00=20382.39=2047.52=20Tm=20(10)=20Tj=0AET=0A59.04=20109.75=20= m=2059.04=20361.87=20l=20S=0A59.04=20109.75=20m=2051.84=20109.75=20l=20S=0A= 59.04=20160.17=20m=2051.84=20160.17=20l=20S=0A59.04=20210.59=20m=2051.84=20= 210.59=20l=20S=0A59.04=20261.02=20m=2051.84=20261.02=20l=20S=0A59.04=20= 311.44=20m=2051.84=20311.44=20l=20S=0A59.04=20361.87=20m=2051.84=20= 361.87=20l=20S=0ABT=0A/F1=201=20Tf=200.00=2012.00=20-12.00=200.00=20= 41.76=20106.41=20Tm=20(0)=20Tj=0A/F1=201=20Tf=200.00=2012.00=20-12.00=20= 0.00=2041.76=20156.83=20Tm=20(2)=20Tj=0A/F1=201=20Tf=200.00=2012.00=20= -12.00=200.00=2041.76=20207.26=20Tm=20(4)=20Tj=0A/F1=201=20Tf=200.00=20= 12.00=20-12.00=200.00=2041.76=20257.68=20Tm=20(6)=20Tj=0A/F1=201=20Tf=20= 0.00=2012.00=20-12.00=200.00=2041.76=20308.11=20Tm=20(8)=20Tj=0A/F1=201=20= Tf=200.00=2012.00=20-12.00=200.00=2041.76=20355.19=20Tm=20(10)=20Tj=0AET=0A= 59.04=2073.44=20m=0A401.76=2073.44=20l=0A401.76=20372.96=20l=0A59.04=20= 372.96=20l=0A59.04=2073.44=20l=0AS=0AQ=20q=0ABT=0A0.000=200.000=200.000=20= rg=0A/F1=201=20Tf=2012.00=200.00=20-0.00=2012.00=20208.89=2018.72=20Tm=20= (X-scale)=20Tj=0A/F1=201=20Tf=200.00=2012.00=20-12.00=200.00=2012.96=20= 201.69=20Tm=20(Y-scale)=20Tj=0AET=0AQ=20q=2059.04=2073.44=20342.72=20= 299.52=20re=20W=20n=0A0.000=200.000=201.000=20RG=0A0.75=20w=0A[]=200=20d=0A= 59.04=20155.98=20m=20401.76=20292.67=20l=20S=0A1.000=200.000=200.000=20= RG=0A59.04=20101.61=20m=20401.76=20321.20=20l=20S=0AQ=0Aendstream=0A= endobj=0A14=200=20obj=0A13098=0Aendobj=0A3=200=20obj=0A<<=0A/Type=20= /Pages=0A/Kids=20[=0A12=200=20R=0A]=0A/Count=201=0A/MediaBox=20[0=200=20= 432=20432]=0A>>=0Aendobj=0Axref=0A0=2015=0A0000000000=2065535=20f=20=0A= 0000000020=2000000=20n=20=0A0000000162=2000000=20n=20=0A0000014412=20= 00000=20n=20=0A0000000211=2000000=20n=20=0A0000000328=2000000=20n=20=0A= 0000000584=2000000=20n=20=0A0000000680=2000000=20n=20=0A0000000781=20= 00000=20n=20=0A0000000885=2000000=20n=20=0A0000000993=2000000=20n=20=0A= 0000001071=2000000=20n=20=0A0000001155=2000000=20n=20=0A0000001237=20= 00000=20n=20=0A0000014390=2000000=20n=20=0Atrailer=0A<<=0A/Size=2015=0A= /Info=201=200=20R=0A/Root=202=200=20R=0A>>=0Astartxref=0A14496=0A%%EOF=0A= --Boundary_(ID_+WcXmHmYo1zpkvYSkQNFBw) Content-type: multipart/alternative; boundary="Boundary_(ID_FqB2HAv6dQB6GGTiiimTqw)" --Boundary_(ID_FqB2HAv6dQB6GGTiiimTqw) Content-type: text/plain; charset=ISO-8859-1; format=flowed Content-transfer-encoding: quoted-printable On Apr 10, 2005, at 3:41 AM, Jamila Ahdidan wrote: > Hi Jason, > Well I have to say that I don't use glim_image to perform a VBM study,=20= > but to perform a t test at each voxel to assess the difference between=20= > my group of patients and my group of controls. So, I don't really know=20= > whether I want to force to be 0 at x=3D0. (I don't really know what = that=20 > means!). > Do you think I'm using the wrong minc command, and if yes do you have=20= > another idea? > =A0 > Thanks a lot, > Jamila > > Jason Lerch wrote: > > On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > > > My dilema is whether I should just keep > > my good results and forget about the intercept in the > > matrix, or I should stick to the intercept and > > conclude that nothing is interpretable from my > > results. > > Is there any reason to force the slope to be 0 at x=3D0? If you have > standard VBM density data, then that is an invalid assumption, since > there is every reason to allow the y to take on an arbitrary value at > x=3D0, so you would include that column of ones for your intercept = term. > If you have different data then this assumption might be valid - > something that is the case, for example, when looking at asymmetry = VBM. > But by and large you will want an intercept. > > Hope this helps, > > Jason > > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > > Do you Yahoo!? > Yahoo! Mail - Find what you need with new enhanced search. Learn = more.= --Boundary_(ID_FqB2HAv6dQB6GGTiiimTqw) Content-type: text/enriched; charset=ISO-8859-1 Content-transfer-encoding: quoted-printable On Apr 10, 2005, at 3:41 AM, Jamila Ahdidan wrote: Hi Jason, Well I have to say that I don't use glim_image to perform a VBM study, but to perform a t test at each voxel to assess the difference between my group of patients and my group of controls. So, I don't really know whether I want to force to be 0 at x=3D0. (I don't really know what that means!). Do you think I'm using the wrong minc command, and if yes do you have another idea? =A0 Thanks a lot, Jamila Jason Lerch < wrote: On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > My dilema is whether I should just keep > my good results and forget about the intercept in the > matrix, or I should stick to the intercept and > conclude that nothing is interpretable from my > results. Is there any reason to force the slope to be 0 at x=3D0? If you have standard VBM density data, then that is an invalid assumption, since there is every reason to allow the y to take on an arbitrary value at=20= x=3D0, so you would include that column of ones for your intercept term. If you have different data then this assumption might be valid -=20 something that is the case, for example, when looking at asymmetry VBM. But by and large you will want an intercept. Hope this helps, Jason _______________________________________________ MINC-users@bic.mni.mcgill.ca http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users Do you Yahoo!? Yahoo! Mail - Find what you need with new enhanced search. 0000,0000,EEEELearn more.= --Boundary_(ID_FqB2HAv6dQB6GGTiiimTqw)-- --Boundary_(ID_+WcXmHmYo1zpkvYSkQNFBw)-- From minc-users@bic.mni.mcgill.ca Mon Apr 11 12:46:04 2005 From: minc-users@bic.mni.mcgill.ca (Nameeta Lobo) Date: Mon Apr 11 11:46:04 2005 Subject: [MINC-users] Error running AAL software Message-ID: <1113234347.425a9bab0b0d4@netmail.bsd.uchicago.edu> Hello all, I had sent this e-mail sometime ago. I am interested in using the AAL software. I use AFNI for my data analysis though. From what I have read AAL is compatible with SPM alone and so downloaded SPM2 and then followed the instructions on the AAL software README file to download AAL http://www.cyceron.fr/freeware/. Whenever I open Matlab, type in aal and then hit any of the labeling procedure buttons, I get the following error ??? Undefined function or variable 'spm'. Error in ==> /home/people/nlobo/SPM2/spm2/toolbox/aal/spm_getSPM.m On line 145 ==> SPMid = spm('SFnBanner',mfilename,SCCSid); Error in ==> /home/people/nlobo/SPM2/spm2/toolbox/aal/gin_dlabels.m On line 16 ==> [SPM,xSPM] = spm_getSPM; ??? Error while evaluating uicontrol Callback. I dont seem to know what the error is. Also, is it possible to use these ROIs with AFNI directly by just resampling them or something. These ROIs are in the ANALYZE format and so conversion to other formats could be quite simple unless I am totally wrong about this. I have carried out all my analysis in AFNI and dont have any knowledge of SPM. Any help will be really appreciated. Thanks a lot for your time . Nameeta From minc-users@bic.mni.mcgill.ca Mon Apr 11 16:59:03 2005 From: minc-users@bic.mni.mcgill.ca (Jamila Ahdidan) Date: Mon Apr 11 15:59:03 2005 Subject: [MINC-users] glim_image In-Reply-To: 6667 Message-ID: <20050411195849.59630.qmail@web60805.mail.yahoo.com> --0-1277789368-1113249529=:59368 Content-Type: text/plain; charset=us-ascii Dear Jason, Thanks for your explanation. It help a lot! I have one (or 2) more question. The results from a glim_image with a matrix an intercept and the one without intercept are the same, except from the range of t values that is bigger without intercept, which allows me to define highly significant areas (even after the very conservative bonferroni correction). So, I wonder whether the results from the glim without intercept can be used to give an estimate of what could be found if we had a bigger sample size. Do you think I can use the results from the glim without intercept in some way? and if yes how? I hope you'll find the time to answer! Many regards, Jamila Jason Lerch wrote: Greetings again, I've attached a graph which shows the difference - the data was generated using the following function: y = 2 + 0.5x + Error where x is a simple sequence between 0 and 10. There are two regression lines through the graph. The one in blue fits a model including the intercept term, the one in red fits a model without an intercept. You can see that the red line takes on a value of 0 when x=0 - and that therefore the fit is not as accurate as the blue line. None of this stuff is specific to glim_image - this is all standard linear model statistics. One of the best online references that I know about is here: http://www.itl.nist.gov/div898/handbook/ though there surely are others as well. Good luck, Jason On Apr 10, 2005, at 3:41 AM, Jamila Ahdidan wrote: > Hi Jason, > Well I have to say that I don't use glim_image to perform a VBM study, > but to perform a t test at each voxel to assess the difference between > my group of patients and my group of controls. So, I don't really know > whether I want to force to be 0 at x=0. (I don't really know what that > means!). > Do you think I'm using the wrong minc command, and if yes do you have > another idea? > > Thanks a lot, > Jamila > > Jason Lerch wrote: > > On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > > > My dilema is whether I should just keep > > my good results and forget about the intercept in the > > matrix, or I should stick to the intercept and > > conclude that nothing is interpretable from my > > results. > > Is there any reason to force the slope to be 0 at x=0? If you have > standard VBM density data, then that is an invalid assumption, since > there is every reason to allow the y to take on an arbitrary value at > x=0, so you would include that column of ones for your intercept term. > If you have different data then this assumption might be valid - > something that is the case, for example, when looking at asymmetry VBM. > But by and large you will want an intercept. > > Hope this helps, > > Jason > > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > > Do you Yahoo!? > Yahoo! Mail - Find what you need with new enhanced search. Learn more. --------------------------------- Do you Yahoo!? Yahoo! Small Business - Try our new resources site! --0-1277789368-1113249529=:59368 Content-Type: text/html; charset=us-ascii
Dear Jason,
Thanks for your explanation. It help a lot!
I have one (or 2) more question. The results from a glim_image with a matrix an intercept and the one without intercept are the same, except from the range of t values that is bigger without intercept, which allows me to define highly significant areas (even after the very conservative bonferroni correction). So, I wonder whether the results from the glim without intercept can be used to give an estimate of what could be found if we had a bigger sample size. Do you think I can use the results from the glim without intercept in some way? and if yes how?
 
I hope you'll find the time to answer!
 
Many regards,
Jamila


Jason Lerch <jason@bic.mni.mcgill.ca> wrote:
Greetings again,

I've attached a graph which shows the difference - the data was
generated using the following function:

y = 2 + 0.5x + Error

where x is a simple sequence between 0 and 10.

There are two regression lines through the graph. The one in blue fits
a model including the intercept term, the one in red fits a model
without an intercept. You can see that the red line takes on a value of
0 when x=0 - and that therefore the fit is not as accurate as the blue
line.

None of this stuff is specific to glim_image - this is all standard
linear model statistics. One of the best online references that I know
about is here:

http://www.itl.nist.gov/div898/handbook/

though there surely are others as well.

Good luck,

Jason



On Apr 10, 2005, at 3:41 AM, Jamila Ahdidan wrote:

> Hi Jason,
> Well I have to say that I don't use glim_image to perform a VBM study,
> but to perform a t test at each voxel to assess the difference between
> my group of patients and my group of controls. So, I don't really know
> whether I want to force to be 0 at x=0. (I don't really know what that
> means!).
> Do you think I'm using the wrong minc command, and if yes do you have
> another idea?
>  
> Thanks a lot,
> Jamila
>
> Jason Lerch wrote:
>
> On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote:
>
> > My dilema is whether I should just keep
> > my good results and forget about the intercept in the
> > matrix, or I should stick to the intercept and
> > conclude that nothing is interpretable from my
> > results.
>
> Is there any reason to force the slope to be 0 at x=0? If you have
> standard VBM density data, then that is an invalid assumption, since
> there is every reason to allow the y to take on an arbitrary value at
> x=0, so you would include that column of ones for your intercept term.
> If you have different data then this assumption might be valid -
> something that is the case, for example, when looking at asymmetry VBM.
> But by and large you will want an intercept.
>
> Hope this helps,
>
> Jason
>
> _______________________________________________
> MINC-users@bic.mni.mcgill.ca
> http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users
>
> Do you Yahoo!?
> Yahoo! Mail - Find what you need with new enhanced search. Learn more.

Do you Yahoo!?
Yahoo! Small Business - Try our new resources site! --0-1277789368-1113249529=:59368-- From minc-users@bic.mni.mcgill.ca Mon Apr 11 17:11:04 2005 From: minc-users@bic.mni.mcgill.ca (Jason Lerch) Date: Mon Apr 11 16:11:04 2005 Subject: [MINC-users] glim_image In-Reply-To: <20050411195849.59630.qmail@web60805.mail.yahoo.com> References: <20050411195849.59630.qmail@web60805.mail.yahoo.com> Message-ID: <047734a7c69a8f60cfc7dd3a8e37613a@bic.mni.mcgill.ca> --Apple-Mail-16--348434962 Content-Transfer-Encoding: quoted-printable Content-Type: text/plain; charset=ISO-8859-1; format=flowed Greetings again, no, you cannot use the results of the regression without an intercept=20 unless you have a good a priori reason of eliminating the intercept=20 from your model - which it doesn't seem like you have. Sorry! Jason On Apr 11, 2005, at 3:58 PM, Jamila Ahdidan wrote: > Dear Jason, > Thanks for your explanation. It help a lot! > I have one (or 2) more question. The results from a glim_image with a=20= > matrix an intercept and the one without intercept are the same, except=20= > from the range of t values that is bigger without intercept, which=20 > allows me to define highly significant areas (even after the very=20 > conservative bonferroni correction). So, I wonder whether the results=20= > from the glim without intercept can be used to give an estimate of=20 > what could be found if we had a bigger sample size. Do you think I can=20= > use the results from the glim without intercept in some way? and if=20 > yes how? > =A0 > I hope you'll find the time to answer! > =A0 > Many regards, > Jamila > > > Jason Lerch wrote: > Greetings again, > > I've attached a graph which shows the difference - the data was > generated using the following function: > > y =3D 2 + 0.5x + Error > > where x is a simple sequence between 0 and 10. > > There are two regression lines through the graph. The one in blue fits > a model including the intercept term, the one in red fits a model > without an intercept. You can see that the red line takes on a value=20= > of > 0 when x=3D0 - and that therefore the fit is not as accurate as the = blue > line. > > None of this stuff is specific to glim_image - this is all standard > linear model statistics. One of the best online references that I = know > about is here: > > http://www.itl.nist.gov/div898/handbook/ > > though there surely are others as well. > > Good luck, > > Jason > > > > On Apr 10, 2005, at 3:41 AM, Jamila Ahdidan wrote: > > > Hi Jason, > > Well I have to say that I don't use glim_image to perform a VBM=20 > study, > > but to perform a t test at each voxel to assess the difference=20 > between > > my group of patients and my group of controls. So, I don't really=20= > know > > whether I want to force to be 0 at x=3D0. (I don't really know what=20= > that > > means!). > > Do you think I'm using the wrong minc command, and if yes do you = have > > another idea? > > =A0 > > Thanks a lot, > > Jamila > > > > Jason Lerch wrote: > > > > On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > > > > > My dilema is whether I should just keep > > > my good results and forget about the intercept in the > > > matrix, or I should stick to the intercept and > > > conclude that nothing is interpretable from my > > > results. > > > > Is there any reason to force the slope to be 0 at x=3D0? If you have > > standard VBM density data, then that is an invalid assumption, since > > there is every reason to allow the y to take on an arbitrary value = at > > x=3D0, so you would include that column of ones for your intercept=20= > term. > > If you have different data then this assumption might be valid - > > something that is the case, for example, when looking at asymmetry=20= > VBM. > > But by and large you will want an intercept. > > > > Hope this helps, > > > > Jason > > > > _______________________________________________ > > MINC-users@bic.mni.mcgill.ca > > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > > > > Do you Yahoo!? > > Yahoo! Mail - Find what you need with new enhanced search. Learn=20 > more. > > Do you Yahoo!? > Yahoo! Small Business - Try our new resources site!= --Apple-Mail-16--348434962 Content-Transfer-Encoding: quoted-printable Content-Type: text/enriched; charset=ISO-8859-1 Greetings again, no, you cannot use the results of the regression without an intercept unless you have a good a priori reason of eliminating the intercept from your model - which it doesn't seem like you have. Sorry! Jason On Apr 11, 2005, at 3:58 PM, Jamila Ahdidan wrote: Dear Jason, Thanks for your explanation. It help a lot! I have one (or 2) more question. The results from a glim_image with a matrix an intercept and the one without intercept are the same, except from the range of t values that is bigger without intercept, which allows me to define highly significant areas (even after the very conservative bonferroni correction). So, I wonder whether the results from the glim without intercept can be used to give an estimate of what could be found if we had a bigger sample size. Do you think I can use the results from the glim without intercept in some way? and if yes how? =A0 I hope you'll find the time to answer! =A0 Many regards, Jamila Jason Lerch < wrote: Greetings again, I've attached a graph which shows the difference - the data was generated using the following function: y =3D 2 + 0.5x + Error where x is a simple sequence between 0 and 10. There are two regression lines through the graph. The one in blue fits a model including the intercept term, the one in red fits a model without an intercept. You can see that the red line takes on a value of=20 0 when x=3D0 - and that therefore the fit is not as accurate as the blue line. None of this stuff is specific to glim_image - this is all standard linear model statistics. One of the best online references that I know about is here: http://www.itl.nist.gov/div898/handbook/ though there surely are others as well. Good luck, Jason On Apr 10, 2005, at 3:41 AM, Jamila Ahdidan wrote: > Hi Jason, > Well I have to say that I don't use glim_image to perform a VBM study, > but to perform a t test at each voxel to assess the difference between > my group of patients and my group of controls. So, I don't really know > whether I want to force to be 0 at x=3D0. (I don't really know what that > means!). > Do you think I'm using the wrong minc command, and if yes do you have > another idea? > =A0 > Thanks a lot, > Jamila > > Jason Lerch wrote: > > On Apr 9, 2005, at 6:51 PM, Jamila Ahdidan wrote: > > > My dilema is whether I should just keep > > my good results and forget about the intercept in the > > matrix, or I should stick to the intercept and > > conclude that nothing is interpretable from my > > results. > > Is there any reason to force the slope to be 0 at x=3D0? If you have > standard VBM density data, then that is an invalid assumption, since > there is every reason to allow the y to take on an arbitrary value at > x=3D0, so you would include that column of ones for your intercept term. > If you have different data then this assumption might be valid - > something that is the case, for example, when looking at asymmetry VBM. > But by and large you will want an intercept. > > Hope this helps, > > Jason > > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > > Do you Yahoo!? > Yahoo! Mail - Find what you need with new enhanced search. Learn more. Do you Yahoo!? Yahoo! Small Business - 0000,0000,EEEETry our new resources site!= --Apple-Mail-16--348434962-- From Andrew Janke Mon Apr 11 22:07:03 2005 From: Andrew Janke (Andrew Janke) Date: Mon Apr 11 21:07:03 2005 Subject: [MINC-users] Error running AAL software In-Reply-To: <1113234347.425a9bab0b0d4@netmail.bsd.uchicago.edu> References: <1113234347.425a9bab0b0d4@netmail.bsd.uchicago.edu> Message-ID: On Apr 12, 2005 1:45 AM, Nameeta Lobo wrote: > Hello all, > > I had sent this e-mail sometime ago. I am interested in using the AAL software. > I use AFNI for my data analysis though. From what I have read AAL is compatible > with SPM alone and so downloaded SPM2 and then followed the instructions on the > AAL software README file to download AAL http://www.cyceron.fr/freeware/. > Whenever I open Matlab, type in aal and then hit any of the labeling procedure > buttons, I get the following error > > ??? Undefined function or variable 'spm'. > > Error in ==> /home/people/nlobo/SPM2/spm2/toolbox/aal/spm_getSPM.m > On line 145 ==> SPMid = spm('SFnBanner',mfilename,SCCSid); > > Error in ==> /home/people/nlobo/SPM2/spm2/toolbox/aal/gin_dlabels.m > On line 16 ==> [SPM,xSPM] = spm_getSPM; > > ??? Error while evaluating uicontrol Callback. > > I dont seem to know what the error is. > > Also, is it possible to use these ROIs with AFNI directly by just resampling > them or something. These ROIs are in the ANALYZE format and so conversion to > other formats could be quite simple unless I am totally wrong about this. I > have carried out all my analysis in AFNI and dont have any > knowledge of SPM. Nameeta I suspect you meant to send this to the SPM mailing list not the MINC mailing list... :) -- Andrew Janke (a.janke@gmail.com || www.cmr.uq.edu.au/~rotor) Australia->Brisbane H: +61 7 3390 6332 || M: +61 4 2138 8581 From Andrew Janke Sat Apr 16 09:54:04 2005 From: Andrew Janke (Andrew Janke) Date: Sat Apr 16 08:54:04 2005 Subject: [MINC-users] Re: Display + Cygwin In-Reply-To: References: <42381FF0.9060007@drcmr.dk> <423999EF.9070402@drcmr.dk> <4252986C.1040505@drcmr.dk> <425E4CAD.2040807@drcmr.dk> Message-ID: Just in case anyone else out there needs a fix to the Display backup file naming problem under Cygwin, here is a binary version that should tide you over untill the next release of MINC (1.4) http://packages.bic.mni.mcgill.ca/cygwin/Display-cygwin-binary-2005-04-15.exe.gz a From minc-users@bic.mni.mcgill.ca Mon Apr 18 16:54:04 2005 From: minc-users@bic.mni.mcgill.ca (Roy Yves) Date: Mon Apr 18 15:54:04 2005 Subject: [MINC-users] mincinfo -attvalue df:df blah_t.mnc problem Message-ID: Hello: mincinfo -attvalue df:df blah_t.mnc gives me the following error: ncvarid: ncid 3: Variable nout found. Should't the df be a variable present in a fmristat statistical output file ? Merci Yves From minc-users@bic.mni.mcgill.ca Mon Apr 18 17:16:04 2005 From: minc-users@bic.mni.mcgill.ca (Sylvain MILOT) Date: Mon Apr 18 16:16:04 2005 Subject: [MINC-users] mincinfo -attvalue df:df blah_t.mnc problem In-Reply-To: Message-ID: Not sure if it is set or not but you could try (assuming it has df in its name) mincheader blah_t.mnc | grep -i df to see what turns up. HTH, Sylvain On Mon, 18 Apr 2005, Roy Yves wrote: > Hello: > > mincinfo -attvalue df:df blah_t.mnc > > gives me the following error: > > ncvarid: ncid 3: Variable nout found. > > Should't the df be a variable present in a fmristat statistical output file ? > > Merci > Yves > > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > --- Sylvain Milot (sylvain@bic.mni.mcgill.ca) (trinity@bic.mni.mcgill.ca) Brain Imaging Centre Montreal Neurological Institute Webster 2B, Room 208 Montreal, Qc., Canada, H3A 2B4 Phone : (514) 398-4965, Fax: 398-8948 Mobile : (514) 712-1768 Office : 527 Av Des Pins O., Room 204 Montreal, Qc., H2W 1S4 From minc-users@bic.mni.mcgill.ca Wed Apr 20 20:24:03 2005 From: minc-users@bic.mni.mcgill.ca (minc-users@bic.mni.mcgill.ca) Date: Wed Apr 20 19:24:03 2005 Subject: [MINC-users] error running mincmath Message-ID: <1114039381.4266e455c959d@web.mail.umich.edu> Command: (executed during nu_evaluate) mincmath -clobber -verbose -copy_header -zero -div in.mnc.gz in_nuc_field.mnc out.mnc Error message: Files in_nuc_field.mnc and in.mnc.gz have different coordinate spacings (zspace) nu_evaluate: crashed while running mincmath (termination status=256) When I run mincinfo on the two volumes, I get the same output: image dimensions: zspace yspace xspace dimension name length step start -------------- ------ ---- ----- zspace 112 1.5 -86.396 yspace 256 -1.48438 178.728 xspace 256 -1.48438 190 any ideas on what causing? Thank you. From minc-users@bic.mni.mcgill.ca Thu Apr 21 12:18:09 2005 From: minc-users@bic.mni.mcgill.ca (D. Louis Collins) Date: Thu Apr 21 11:18:09 2005 Subject: [MINC-users] error running mincmath In-Reply-To: <1114039381.4266e455c959d@web.mail.umich.edu> References: <1114039381.4266e455c959d@web.mail.umich.edu> Message-ID: <94d149057afdb0f1c504973f2c726b9d@mcgill.ca> Haytham, If the mincinfo output for the two volumes is the same, (ie, same starts, steps and number of elements _and_ same dimension ordering), then I suspect that the two volumes have different direction cosines. Check with 'mincheader file.mnc | grep direction_cosine'. you could rerun with '-nocheck_dimensions' - this will force the code to do voxel-to-voxel manipulations, without checking against the dimension information. -Louis On Apr 20, 2005, at 7:23 PM, haytham@umich.edu wrote: > Command: (executed during nu_evaluate) > > mincmath -clobber -verbose -copy_header -zero -div in.mnc.gz > in_nuc_field.mnc > out.mnc > > > Error message: > Files in_nuc_field.mnc and in.mnc.gz have different coordinate > spacings (zspace) > nu_evaluate: crashed while running mincmath (termination status=256) > > > When I run mincinfo on the two volumes, I get the same output: > > image dimensions: zspace yspace xspace > dimension name length step start > -------------- ------ ---- ----- > zspace 112 1.5 -86.396 > yspace 256 -1.48438 178.728 > xspace 256 -1.48438 190 > > > any ideas on what causing? > > Thank you. > _______________________________________________ > MINC-users@bic.mni.mcgill.ca > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > From minc-users@bic.mni.mcgill.ca Thu Apr 21 19:56:04 2005 From: minc-users@bic.mni.mcgill.ca (Robert VINCENT) Date: Thu Apr 21 18:56:04 2005 Subject: [MINC-users] dcm2mnc 2.0.03 available In-Reply-To: <41C545A2921542ykuma-tky@umin.ac.jp> Message-ID: Hello again, I've released an updated binary version of the dcm2mnc conversion tool for Intel Linux users interested in testing it: http://www.bic.mni.mcgill.ca/~bert/dcm2mnc-2.0.03-i386-linux.tar.gz A new option, "-minmax", has been added. This option controls whether or not the converter uses the DICOM fields 0028,0106 and 0028,0107 to set the valid minimum and maximum voxel values. This option may be needed to obtain proper scaling on some PET systems. A number of new warning messages have been added, and some of the more spurious messages present in the past releases have been suppressed. In general, this converter should fix most of the bugs that have been reported to me, or at least those I was able to reproduce. I hope to have a source release available before the end of the month. -bert From cjb@pet.auh.dk Mon Apr 25 08:53:04 2005 From: cjb@pet.auh.dk (Christopher Bailey) Date: Mon Apr 25 07:53:04 2005 Subject: [MINC-users] DEPICT read minc? Message-ID: <1114429974.8500.44.camel@kafka.pet.auh.dk> Dear list, I would like to try out Dr. Gunn et al.'s DEPICT parametric imaging software on some dynamic MINC data. I am only missing some m-files for reading/writing MINC, which I'd rather be grateful for than author of ;)! I have EMMA (0.9.7) installed. Sincerely, Chris -- Christopher Bailey PET Centre and Center for Functionally Integrative Neuroscience Aarhus University Hospital, Denmark http://www.cfin.au.dk/ From cjb@pet.auh.dk Wed Apr 27 06:06:05 2005 From: cjb@pet.auh.dk (Christopher Bailey) Date: Wed Apr 27 05:06:05 2005 Subject: [MINC-users] DEPICT read minc? In-Reply-To: References: Message-ID: <1114592697.11279.7.camel@kafka.pet.auh.dk> Hi Bert, Dr. Gunn was kind enough to send me the file I needed: DEPICT_read_minc.m Basically, it uses EMMA's miinquire and mireadimages to get the image data into a struct that DEPICT can then process. I'm having some problems, though, related to the rest of the code being rather hard-coded to Analyze. In fact the MINC reading seems not to be functional as it is (perhaps written for an earlier version of EMMA), so I'm looking into what I need to change. I'm also still struggling with the exact format of the input file, so no parametric images yet! -Chris On Tue, 2005-04-26 at 19:03 -0400, Robert VINCENT wrote: > Hi Christopher, > > Has anyone replied to this yet? I'm not familiar with the DEPICT code or > I would've replied sooner. I'm not certain what m-files you would need, > EMMA should include all of the m-files to perform basic MINC I/O in > Matlab. > > -bert > > On Mon, 25 Apr 2005, Christopher Bailey wrote: > > > > > Dear list, > > > > I would like to try out Dr. Gunn et al.'s DEPICT parametric imaging > > software on some dynamic MINC data. I am only missing some m-files for > > reading/writing MINC, which I'd rather be grateful for than author > > of ;)! I have EMMA (0.9.7) installed. > > > > Sincerely, > > > > Chris > > -- > > Christopher Bailey > > PET Centre and > > Center for Functionally Integrative Neuroscience > > Aarhus University Hospital, Denmark > > http://www.cfin.au.dk/ > > > > > > > > _______________________________________________ > > MINC-users@bic.mni.mcgill.ca > > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users > > >