Aufgrund der Zweifel an der Vergleichbarkeit der Kalibrierungen in unterschiedlichen Sehentfernungen versuchen wir, die Eichparameter statistisch aufzubereiten.
Die benötigten Daten sind zunächst aus Einzeldateien für jeden Versuch zusammenzuführen. Dazu wird zunächst eine Liste mit Dateinamen erstellt, z.B. mit Hilfe der Shell. Da außerdem Dateien wie OVRS_D36_ps1_eich.dat und OVRS_D36_ps1_SD.dat, sollen nur die ersten Namensteile wie OVRS_D36_ps1 stehen bleiben.
Die Liste sei in der Datei files.txt enthalten. Aus den eich.dat-Dateien entnehmen wir 24 Zahlen aus den Zeilen 2-7:
liste <- readLines("files.txt")
n = length(liste)
arr = data.frame(dim=c(n,28))
for( i in 1:n ) arr[i,1:24] <- scan(paste(liste[[i]],"_eich.dat",sep=""),skip=1,nlines=6)
Die SD-Daten werden zugefügt:
for( i in 1:n ) arr[i,24+1:4] <- scan(paste(liste[[i]],"_SD.dat",sep=""))
Für die eingelesenen Daten wählen wir geeignete Variablennamen
colnames(arr)[24+1:4] <- c("SD.hl","SD.hr","SD.vl","SD.vr")
colnames(arr)[1:12] <- c(outer( c("hl","vl","hr","vr"),c("I","g.h","g.v"),FUN=paste,sep="."))
Die Versuchsbedingungen und Probanden werden aus dem Dateinamen extrahiert:
arr$VP <- substring(rownames(arr),10,12) arr$dist <- factor(as.integer(substring(rownames(arr),7,7))*10) arr$n <- as.integer(substring(rownames(arr),8,8))