Diff for "FAQ/singcase" - CBU statistics Wiki
location: Diff for "FAQ/singcase"
Differences between revisions 28 and 29
Revision 28 as of 2007-03-19 14:41:09
Size: 1691
Editor: PeterWatson
Comment:
Revision 29 as of 2007-03-19 14:41:27
Size: 1691
Editor: PeterWatson
Comment:
Deletions are marked like this. Additions are marked like this.
Line 34: Line 34:
 * [:FAQ/singcase/abnorm Confidence Intervals for Abnormality]  * [:FAQ/singcase/abnorm: Confidence Intervals for Abnormality]
Line 36: Line 36:
 * [http://www.abdn.ac.uk/~psy086/dept/singslope.htm: Downloadable Program for comparing individual with Controls]  * [http://www.abdn.ac.uk/~psy086/dept/singslope.htm Downloadable Program for comparing individual with Controls]

Comparing a within subjects group difference to that of a single case

Pairwise case

Suppose a group of controls and a single case do a pair of conditions (c1, c2, say) and we are interested in seeing if the relationship between the conditions in the controls is the same as that in the single case.

Since we have a single case we can treat its difference as a constant. So we end up with a one sample t-test on the control group difference minus the case difference. Ie we test if the mean difference in controls - difference in case equals 0.

There is a one-sample t-test option under analyze:means in SPSS. The two-sample (independent samples) may also be performed on the difference in SPSS. Alternatively the below example will work using R. Here c1 and c2 are the two conditions and sub denotes if the subject is a control or a case.

Note that these tests assume that the slopes in the control and patient groups have the same variance. [attachment:crawford.pdf A detailed discussion of the methods is available.]

c1 <- c(2,1,2,3,1,2,2,1,1,3,3,1)
c2 <- c(3,4,1,1,6,1,3,2,2,5,2,2)
sub <- c(1,1,1,1,1,1,1,1,1,1,1,2)
diff <- c1–c2
diffn <- diff[sub == 1]
diffp <- diff[sub == 2]
diffy <- diffn - diffp 
const <- gl(1,length(diffn))
id <- gl(length(diffn),1)
t.test(diffy, mu=0)
df <- length(diffn)-1
toutc <- (mean(diffn)-diffp) / sqrt(var(diffn)*((nc+1)/nc))
pv <- 2*pt(-abs(toutc), df)

None: FAQ/singcase (last edited 2017-04-04 15:57:39 by PeterWatson)