library(foreign) t1 <- read.spss("weller.sav") d1 <- as.data.frame(t1) d1[,16] <- (d1[,16]=="Risky Choice")-0 Months <- d1[,2] Adult <- is.na(d1$AGE) Age <- (Months>90 | Adult)+Adult Cage <- Months>90 Sex <- d1[,1] ns <- nrow(d1) R <- as.matrix(d1[,16:69]) Rsq <- rep(NA,ns) Cups <- rep(c(2,3,5,2,3,5),rep(9,6)) # p is 1/Cups Ps <- 1/Cups # doesn't work as well as Cups Qs <- rep(c(2,2,2,3,3,3,5,5,5),6) # quarters, max outcome GL <- rep(c(1,-1),c(27,27)) # gain/loss EV <- (GL*Qs/Cups)-GL # EV(risk) compared to no-risk sEV <- sign(EV) # works better than EV, but signif only for adults Regg <- matrix(NA,ns,3) Regl <- matrix(NA,ns,3) colnames(Regl) <- colnames(Regg) <- c("Int","Qs","Cups") for (i in 1:ns) {Regg[i,] <- coef(lm(R[i,1:27]~Qs[1:27]+Cups[1:27])) Regl[i,] <- coef(lm(R[i,28:54]~Qs[1:27]+Cups[1:27])) } RelWtG <- (Regg[,2]+Regg[,3]) RelWtL <- (Regl[,2]+Regl[,3]) var.test(RelWtL ~ Adult) var.test(RelWtG ~ Adult) by(RelWtL,Adult,var,na.rm=T) by(RelWtL,Adult,var,na.rm=T) # Dom means i dominates j, so this is violations of dominance Dom <- matrix(NA,27,27) Errg <- rep(0,ns) Errl <- rep(0,ns) for (i in 1:27) for (j in 1:27) {Dom[i,j] <- sign(Cups[j]-Cups[i])+sign(Qs[i]-Qs[j]) if (Dom[i,j]>0) {Errg <- Errg+(R[,i]R[,j])} if (Dom[i,j]>0) {Errl <- Errl+(R[,i+27]>R[,j+27])} if (Dom[i,j]<0) {Errl <- Errl+(R[,i+27]