CNP SMNM

From Pheno Wiki
Revision as of 16:49, 24 June 2011 by Katie (Talk | contribs) (Code/Algorithms)

Jump to: navigation, search

go back to HTAC

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

Smnm bigschematic.jpg

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

Smnm all data tables.jpg

Smnm all useable tables.jpg

Smnm comp useable tables.jpg

Smnm comp useable tables nooutlier.jpg

Smnm percorr.jpg

Smnm main percorr.jpg

Smnm manip percorr.jpg

Smnm maint rt.jpg

Smnm manip rt.jpg

References

Cannon, 2005, Archives of General Psychiatry
Kim, 2004, Schizophrenia Research
Glahn, 2002, Neuroimage