-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathkrus96train.R
84 lines (67 loc) · 3.5 KB
/
krus96train.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
krus96train <-
function(blocks = 15, subjs = 56, ctxt = TRUE, seed = 1) {
set.seed(seed)
## note "x7"=bias cue, which is always on
sr <- rbind(
##I1, PC1, PR1, I2, PC2, PR2, X, C1, R1, C2, R2
c( 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0), # I1.PC1
c( 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0), # I1.PC1
c( 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0), # I1.PC1
c( 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0), # I1.PR1
c( 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0), # I2.PC2
c( 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0), # I2.PC2
c( 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0), # I2.PC2
c( 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1) # I2.PR2
)
trainingitems <- data.frame(cbind("", sr))
colnames(trainingitems) <- c("stim", "x1", "x2", "x3", "x4", "x5",
"x6", "x7", "t1", "t2", "t3", "t4")
trainingitems[,"stim"] <- c(rep("I1.PC1", 3), "I1.PR1",
rep("I2.PC2", 3), "I2.PR2")
testitems <- rbind(
##I1, PC1, PR1, I2, PC2, PR2, X, C1, R1, C2, R2
c( 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0), # I1
c( 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0), # PC1
c( 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0), # PR1
c( 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0), # I2
c( 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0), # PC2
c( 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0), # PR2
c( 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0), # PC1.PR1
c( 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0), # PC2.PR2
c( 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0), # I1.PC1.PR1
c( 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0), # I2.PC2.PR2
c( 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0), # I1.PC2
c( 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0), # I1.PR2
c( 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0), # I2.PC1
c( 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0), # I2.PR1
c( 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0), # PC1.PR2
c( 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0), # PC2.PR1
c( 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0), # I1.PC1.PR2
c( 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0) # I2.PC2.PR1
)
teststim <- c("I1", "PC1", "PR1", "I2", "PC2", "PR2", "PC1.PR1",
"PC2.PR2", "I1.PC1.PR1", "I2.PC2.PR2", "I1.PC2",
"I1.PR2", "I2.PC1", "I2.PR1", "PC1.PR2", "PC2.PR1",
"I1.PC1.PR2", "I2.PC2.PR1")
bigtr <- NULL
for(subj in 1:subjs) {
tr <- data.frame(matrix(0, ncol=14, nrow = (blocks * nrow(sr))))
for (i in 1:blocks){
samp <- sample(1:nrow(sr), nrow(sr))
tr[(1:nrow(sr)) + nrow(sr) * (i-1), 2] <- i
tr[(1:nrow(sr)) + nrow(sr) * (i-1), 3] <- trainingitems[samp,
"stim"]
tr[(1:nrow(sr)) + nrow(sr) * (i-1), 4:14] <- sr[samp,]
}
tr[1,1]<-1
colnames(tr) <- c("ctrl", "block", colnames(trainingitems))
testrials <- data.frame(2, blocks + 1, teststim, testitems)
colnames(testrials) <- colnames(tr)
traintrials<-nrow(tr)
tr <- rbind(tr, testrials)
tr[(traintrials+1):(nrow(tr)),"ctrl"]<-2
bigtr <- rbind(bigtr, tr)
}
if(!ctxt) bigtr <- bigtr[,c(1:9, 11:14)]
return(bigtr)
}