th <- c(.14, .17, .17, .21, .27, .27, .67, .67, .67, .67, .69, .74)
dat <- c(.14,.14, .19, .35, .14, .25, .39, .46, .58, .69, .64, .69)
dat8 <- c(.25, .32, .21, .43, .14, .21, .35, .50, .38, .63, .44, .57)
dat15 <- c(.20,.19, .18, .29, .14, .25, .37, .70, .58, .79, .66, .53)
datpl <- c(.14, .11, .18, .36, .11, .26, .40, .47, .57, .71, .64, .73)

grip <- function(da){
	test <- abs(da - th)
	print(mean(test))
}

gripper <- function(){
	grip(dat)
	grip(dat8)
	grip(dat15)
	grip(datpl)
}

gripper()

th <- dat

gripper()

#phi = 0.4
th <- c(.17,.2, .2, .25, .31, .31, .63, .20, .43, .32, .17, .69)

gripper()

#phi = 0.32

th <- c(0.13, 0.16, 0.16, 0.66, 0.26, 0.26, 0.76, 0.68, 0.69, 0.69, 0.70, 0.75)

gripper()


