Difference between revisions of "CNP SMNM"
(→Cleaning Rules) |
(→Code/Algorithms) |
||
Line 28: | Line 28: | ||
===Code/Algorithms=== | ===Code/Algorithms=== | ||
+ | //ALL SUBJECTS <br/> | ||
+ | //get breakdown of numbers in each group:<br/> | ||
+ | bysort status: summarize ptid<br/> | ||
+ | bysort la5ccontrol: summarize ptid<br/> | ||
+ | bysort gender: summarize ptid<br/> | ||
+ | |||
+ | //see how many complete cases there are.<br/> | ||
+ | count if smnm_main_mn !=.<br/> | ||
+ | |||
+ | //exclusions<br/> | ||
+ | |||
+ | generate float smnm_percent_corr=(smnm_manip_mn + smnm_main_mn)/2<br/> | ||
+ | count if smnm_percent_corr<=.5<br/> | ||
+ | count if smnm_trialcount != 40 & smnm_main_tt !=.<br/> | ||
+ | count if (smnm_main_mnrt >6000 | smnm_manip_mnrt <6000) & smnm_main_tt !=.<br/> | ||
+ | |||
+ | //create variable summing exclusion criteria<br/> | ||
+ | generate float useable_smnm= 1 if smnm_trialcount==40 & smnm_percent_corr>.5 & smnm_main_mnrt< 6000 & smnm_manip_mnrt<6000<br/> | ||
+ | count if useable_smnm==1<br/> | ||
+ | |||
+ | //SUBJECTS WITH USEABLE DATA<br/> | ||
+ | bysort status: summarize ptid if useable_smnm==1<br/> | ||
+ | bysort la5ccontrol: summarize ptid if useable_smnm==1<br/> | ||
+ | bysort gender: summarize ptid if useable_smnm==1<br/> | ||
+ | summarize age smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1<br/> | ||
+ | bysort status: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1<br/> | ||
+ | bysort la5ccontrol: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1<br/> | ||
+ | bysort gender: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1<br/> | ||
+ | |||
+ | //COMPLETERS WITH USEABLE DATA<br/> | ||
+ | generate float comp_useable_smnm= 1 if smnm_trialcount==40 & smnm_percent_corr>.5 & smnm_main_mnrt< 6000 & smnm_manip_mnrt<6000 &status==2<br/> | ||
+ | count if comp_useable_smnm==1<br/> | ||
+ | bysort status: summarize ptid if comp_useable_smnm==1<br/> | ||
+ | bysort la5ccontrol: summarize ptid if comp_useable_smnm==1<br/> | ||
+ | bysort gender: summarize ptid if comp_useable_smnm==1<br/> | ||
+ | summarize age smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1<br/> | ||
+ | bysort status: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1<br/> | ||
+ | bysort la5ccontrol: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1<br/> | ||
+ | bysort gender: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1<br/> | ||
+ | |||
+ | |||
+ | //OUTLIERS<br/> | ||
+ | sktest smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt<br/> | ||
+ | |||
+ | grubbs smnm_manip_mn if comp_useable_smnm==1, gen (grubb_smnm_manip_mn)<br/> | ||
+ | grubbs smnm_main_mn if comp_useable_smnm==1, gen(grubb_smnm_main_mn)<br/> | ||
+ | grubbs smnm_manip_mnrt if comp_useable_smnm==1, gen (grubb_smnm_manip_mnrt)<br/> | ||
+ | grubbs smnm_main_mnrt if comp_useable_smnm==1, gen(grubb_smnm_main_mnrt)<br/> | ||
+ | grubbs smnm_percent_corr if comp_useable_smnm==1, gen(grubb_smnm_percorr)<br/> | ||
+ | |||
+ | generate float comp_useable_nooutlier_smnm=1 if comp_useable_smnm==1 & grubb_smnm_manip_mn !=1 & grubb_smnm_main_mn !=1 & grubb_smnm_manip_mnrt <br/>!=1 & grubb_smnm_main_mnrt !=1 & grubb_smnm_percorr !=1<br/> | ||
+ | count if comp_useable_nooutlier_smnm==1<br/> | ||
+ | bysort status: summarize ptid if comp_useable_nooutlier_smnm==1<br/> | ||
+ | bysort la5ccontrol: summarize ptid if comp_useable_nooutlier_smnm==1<br/> | ||
+ | bysort gender: summarize ptid if comp_useable_nooutlier_smnm==1<br/> | ||
+ | summarize age smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1<br/> | ||
+ | bysort status: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1<br/> | ||
+ | bysort la5ccontrol: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1<br/> | ||
+ | bysort gender: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1<br/> | ||
+ | |||
+ | histogram smnm_percent_corr if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(SMNM Overall Percent Correct)<br/> | ||
+ | histogram smnm_main_mn if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Maintenance Percent Correct)<br/> | ||
+ | histogram smnm_manip_mn if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Manipulation Percent Correct)<br/> | ||
+ | histogram smnm_main_mnrt if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Maintenance Response Time (Outliers >1650))<br/> | ||
+ | histogram smnm_manip_mnrt if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Manipulation Response Time (Outliers >1800))<br/> | ||
===Data Distributions=== | ===Data Distributions=== |
Revision as of 16:49, 24 June 2011
go back to HTAC
Contents
Basic Task Description
This is a working memory task with conditions challenging both maintenance and manipulation of items being held in working memory.
Task Procedure
Task Structure Detail
Task Schematic
Task Parameters Table
Stimuli
Dependent Variables
Cleaning Rules
1. If any of the calculated variables are missing, that subject should be listed for exclusion. 2. If Trial Count is not =40, then that subject should be listed for exclusion. 3. If overall performance (average percent correct) is below 50% (chance, in this task), that subject should be listed for exclusion. 4. If average RT is greater than 6000ms that subject should be listed for exclusion.
It is worth considering dropping people with performance >50% on either maintenance or manipulation (rather than just on the average).
If you want to know if your group looks approximately like it should, in every sample run so far, accuracy is lower in manipulation than maintenance, and reaction time is greater in manipulation than maintenance. This should be the first group level analysis, to quickly verify everything is correct. It would be highly unexpected to have a deviation from this pattern, even though the position along the y axis (overall accuracy or RT) may vary, especially in patient groups.
Code/Algorithms
//ALL SUBJECTS
//get breakdown of numbers in each group:
bysort status: summarize ptid
bysort la5ccontrol: summarize ptid
bysort gender: summarize ptid
//see how many complete cases there are.
count if smnm_main_mn !=.
//exclusions
generate float smnm_percent_corr=(smnm_manip_mn + smnm_main_mn)/2
count if smnm_percent_corr<=.5
count if smnm_trialcount != 40 & smnm_main_tt !=.
count if (smnm_main_mnrt >6000 | smnm_manip_mnrt <6000) & smnm_main_tt !=.
//create variable summing exclusion criteria
generate float useable_smnm= 1 if smnm_trialcount==40 & smnm_percent_corr>.5 & smnm_main_mnrt< 6000 & smnm_manip_mnrt<6000
count if useable_smnm==1
//SUBJECTS WITH USEABLE DATA
bysort status: summarize ptid if useable_smnm==1
bysort la5ccontrol: summarize ptid if useable_smnm==1
bysort gender: summarize ptid if useable_smnm==1
summarize age smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1
bysort status: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1
bysort la5ccontrol: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1
bysort gender: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if useable_smnm==1
//COMPLETERS WITH USEABLE DATA
generate float comp_useable_smnm= 1 if smnm_trialcount==40 & smnm_percent_corr>.5 & smnm_main_mnrt< 6000 & smnm_manip_mnrt<6000 &status==2
count if comp_useable_smnm==1
bysort status: summarize ptid if comp_useable_smnm==1
bysort la5ccontrol: summarize ptid if comp_useable_smnm==1
bysort gender: summarize ptid if comp_useable_smnm==1
summarize age smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1
bysort status: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1
bysort la5ccontrol: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1
bysort gender: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_smnm==1
//OUTLIERS
sktest smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt
grubbs smnm_manip_mn if comp_useable_smnm==1, gen (grubb_smnm_manip_mn)
grubbs smnm_main_mn if comp_useable_smnm==1, gen(grubb_smnm_main_mn)
grubbs smnm_manip_mnrt if comp_useable_smnm==1, gen (grubb_smnm_manip_mnrt)
grubbs smnm_main_mnrt if comp_useable_smnm==1, gen(grubb_smnm_main_mnrt)
grubbs smnm_percent_corr if comp_useable_smnm==1, gen(grubb_smnm_percorr)
generate float comp_useable_nooutlier_smnm=1 if comp_useable_smnm==1 & grubb_smnm_manip_mn !=1 & grubb_smnm_main_mn !=1 & grubb_smnm_manip_mnrt
!=1 & grubb_smnm_main_mnrt !=1 & grubb_smnm_percorr !=1
count if comp_useable_nooutlier_smnm==1
bysort status: summarize ptid if comp_useable_nooutlier_smnm==1
bysort la5ccontrol: summarize ptid if comp_useable_nooutlier_smnm==1
bysort gender: summarize ptid if comp_useable_nooutlier_smnm==1
summarize age smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1
bysort status: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1
bysort la5ccontrol: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1
bysort gender: summarize smnm_percent_corr smnm_main_mn smnm_manip_mn smnm_main_mnrt smnm_manip_mnrt if comp_useable_nooutlier_smnm==1
histogram smnm_percent_corr if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(SMNM Overall Percent Correct)
histogram smnm_main_mn if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Maintenance Percent Correct)
histogram smnm_manip_mn if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Manipulation Percent Correct)
histogram smnm_main_mnrt if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Maintenance Response Time (Outliers >1650))
histogram smnm_manip_mnrt if comp_useable_smnm==1, frequency bfcolor(gs11) blcolor(black) normal normopts( clcolor(black) ) xtitle(Manipulation Response Time (Outliers >1800))
Data Distributions
References
Cannon, 2005, Archives of General Psychiatry
Kim, 2004, Schizophrenia Research
Glahn, 2002, Neuroimage