Package 'OPI'

Description

Global environment for OPI to hold machine specific constants, etc.

Usage

.opi_env

Format

An object of class environment of length 10.

Description

Given a value in cd/$\mbox{m}^2$, return the equivalent dB value. Default is to use HFA units, so maximum stimulus is 10000 apostilbs.

Usage

cdTodb(cd, maxStim = 10000/pi)

Arguments

Value

A dB value for cd cd/$\mbox{m}^2$.

Examples

# candela to decibels
dB <- cdTodb(10000/pi)  # 0 dB
dB <- cdTodb(1000/pi)   # 10 dB
dB <- cdTodb(100/pi)    # 20 dB
dB <- cdTodb(10/pi)     # 30 dB
dB <- cdTodb(1/pi)      # 40 dB
dB <- cdTodb(0.1/pi)    # 50 dB

Description

It should be called before any other OPI functions.

Usage

chooseOPI(machine = NULL)

chooseOpi(machine = NULL)

Arguments

Value

NULL on success or list of machines otherwise.

Description

Given a value in dB, return the cd/$\mbox{m}^2$ equivalent. Default is to use HFA units, so maximum stimulus is 10000 apostilbs.

Usage

dbTocd(db, maxStim = 10000/pi)

Arguments

Value

cd/$\mbox{m}^2$ value for db dB.

Examples

# decibels to candela
cd <- dbTocd(0)   # 10000/pi
cd <- dbTocd(10)  # 1000/pi
cd <- dbTocd(20)  # 100/pi
cd <- dbTocd(30)  # 10/pi
cd <- dbTocd(40)  # 1/pi

Description

Convert degrees to pixels for machine 'machine'

Usage

degTopix(xy, machine = "compass")

Arguments

Value

xy converted to pixels (top-left is (0,0)) for the machine or NA if machine is unknown

Examples

degTopix(c(0, 0), machine="compass")  # c(960, 960) pixels
degTopix(c(-15, 2)) # c(495, 898) pixels

Description

fourTwo is a 4-2 dB staircase beginning at level est terminating after two reversals. The final estimate is the average of the last two presentations. It also terminates if the minStimulus is not seen twice, or the maxStimulus is seen twice.

Usage

fourTwo.start(est = 25, instRange = c(0, 40), verbose = FALSE, makeStim, ...)

fourTwo.step(state, nextStim = NULL)

fourTwo.stop(state)

fourTwo.final(state)

Arguments

Details

This is an implementation of a 4-2 1-up 1-down staircase. The initial staircase starts at est and proceeds in steps of 4 dB until the first reversal, and 2dB until the next reversal. The mean of the last two presentations is taken as the threshold value. Note this function will repeatedly call opiPresent for a stimulus until opiPresent returns NULL (ie no error occurred). If more than one fourTwo is to be interleaved (for example, testing multiple locations), then the fourTwo.start, fourTwo.step, fourTwo.stop and fourTwo.final calls can maintain the state of the fourTwo after each presentation, and should be used. See examples below.

Value

Multiple locations

fourTwo.start returns a list that can be passed to fourTwo.step, fourTwo.stop, and fourTwo.final. It represents the state of a fourTwo at a single location at a point in time and contains the following.

• name fourTwo

• A copy of all of the parameters supplied to fourTwo.start: startingEstimate=est, minStimulus=instRange[1], maxStimulus=instRange[2], makeStim, and opiParams=list(...)

• currentLevel The next stimulus to present.

• lastSeen The last seen stimulus.

• lastResponse The last response given.

• stairResult The final result if finished (initially NA).

• finished "Not" if staircase has not finished, or one of "Rev" (finished due to 2 reversals), "Max" (finished due to 2 maxStimulus seen), "Min" (finished due to 2 minStimulus not seen)

• numberOfReversals Number of reversals so far.

• currSeenLimit Number of times maxStimulus has been seen.

• currNotSeenLimit Number of times minStimulus not seen.

• numPresentations Number of presentations so far.

• stimuli Vector of stimuli shown at each call to fourTwo.step.

• responses Vector of responses received (1 seen, 0 not) received at each call to fourTwo.step.

• responseTimes Vector of response times received at each call to fourTwo.step.

fourTwo.step returns a list containing

• state The new state after presenting a stimuli and getting a response.

• resp The return from the opiPresent call that was made.

fourTwo.stop returns TRUE if the staircase is finished (2 reversals, or maxStimulus is seen twice or minStimulus is not seen twice).

fourTwo.final returns the final estimate of threshold (mean of last two reversals). This issues a warning if called before the staircase has finished, but still returns a value.

See Also

dbTocd, opiPresent, FT

Examples

# Stimulus is Size III white-on-white as in the HFA
makeStim <- function(db, n) {
  s <- list(x=9, y=9, level=dbTocd(db), size=0.43, color="white",
            duration=200, responseWindow=1500)
  class(s) <- "opiStaticStimulus"
  return(s)
}
chooseOpi("SimHenson")
if (!is.null(opiInitialize(type="C", cap=6)$err))
  stop("opiInitialize failed")

##############################################
# This section is for multiple fourTwos
##############################################
makeStimHelper <- function(db,n, x, y) {  # returns a function of (db,n)
  ff <- function(db, n) db+n
  body(ff) <- substitute({
    s <- list(x=x, y=y, level=dbTocd(db), size=0.43, color="white",
              duration=200, responseWindow=1500)
    class(s) <- "opiStaticStimulus"
    return(s)}, list(x=x,y=y))
  return(ff)
}
# List of (x, y, true threshold) triples
locations <- list(c(9,9,30), c(-9,-9,32), c(9,-9,31), c(-9,9,33))

# Setup starting states for each location
states <- lapply(locations, function(loc) {
  fourTwo.start(makeStim=makeStimHelper(db,n,loc[1],loc[2]),
                tt=loc[3], fpr=0.03, fnr=0.01)})

# Loop through until all states are "stop"
while(!all(st <- unlist(lapply(states, fourTwo.stop)))) {
  i <- which(!st)                         # choose a random,
  i <- i[runif(1, min=1, max=length(i))]  # unstopped state
  r <- fourTwo.step(states[[i]])               # step it
  states[[i]] <- r$state                  # update the states
}

finals <- lapply(states, fourTwo.final)    # get final estimates of threshold
for(i in 1:length(locations)) {
  cat(sprintf("Location (%+2d,%+2d) ",locations[[i]][1], locations[[i]][2]))
      cat(sprintf("has threshold %4.2f\n", finals[[i]]))
}

if (!is.null(opiClose()$err))
  warning("opiClose() failed")

Description

FT begins with a 4-2dB staircase beginning at level est. If the final estimate (last seen) is more than 4dB away from est, a second 4-2 staircase is completed beginning at the estimate returned from the first

Usage

FT(est = 25, instRange = c(0, 40), verbose = FALSE, makeStim, ...)

FT.start(est = 25, instRange = c(0, 40), makeStim, ...)

FT.step(state, nextStim = NULL)

FT.stop(state)

FT.final(state)

Arguments

Details

This is an implementation of a 4-2 1-up 1-down staircase as implemented in the first Humphrey Field Analyzer. The initial staircase starts at est and proceeds in steps of 4 dB until the first reversal, and 2dB until the next reversal. The last seen stimulus is taken as the threshold value. If, after the first staircase, the threshold is more than 4 dB away from the starting point, then a second staircase is initiated with a starting point equal to the threshold found with the first staircase.

Note this function will repeatedly call opiPresent for a stimulus until opiPresent returns NULL (ie no error occurred)

If more than one FT is to be interleaved (for example, testing multiple locations), then the FT.start, FT.step, FT.stop and FT.final calls can maintain the state of the FT after each presentation, and should be used. If only a single FT is required, then the simpler FT can be used. See examples below

Value

Single location

Returns a list containing

• npres Total number of presentations

• respSeq Response sequence stored as a list of (seen,dB) pairs

• first First staircase estimate in dB

• final Final threshold estimate in dB

Multiple locations

FT.start returns a list that can be passed to FT.step, FT.stop, and FT.final. It represents the state of a FT at a single location at a point in time and contains the following.

• name FT

• A copy of all of the parameters supplied to FT.start: startingEstimate=est, minStimulus=instRange[1], maxStimulus=instRange[2], makeStim, and opiParams=list(...).

• currentLevel The next stimulus to present.

• lastSeen The last seen stimulus.

• lastResponse The last response given.

• firstStairResult The result of the first staircase (initially NA).

• secondStairResult The result of the first staircase (initially NA, and could remain NA).

• finished TRUE if staircase has finished (2 reversals, or max/min seen/not-seen twice).

• numberOfReversals Number of reversals so far.

• currSeenLimit Number of times maxStimulus has been seen.

• currNotSeenLimit Number of times minStimulus not seen.

• numPresentations Number of presentations so far.

• stimuli Vector of stimuli shown at each call to FT.step.

• responses Vector of responses received (1 seen, 0 not) received at each call to FT.step.

• responseTimes Vector of response times received at each call toFT.step.

FT.step returns a list containing

• state The new state after presenting a stimuli and getting a response.

• resp The return from the opiPresent call that was made.

FT.stop returns TRUE if the first staircase has had 2 reversals, or maxStimulus is seen twice or minStimulus is not seen twice and the final estimate is within 4 dB of the starting stimulus. Returns TRUE if the second staircase has had 2 reversals, or maxStimulus is seen twice or minStimulus is not seen twice

FT.final returns the final estimate of threshold based on state, which is the last seen in the second staircase, if it ran, or the first staircase otherwise

FT.final.details returns a list containing

• final The final threshold.

• first The threshold determined by the first staircase (might be different from final).

• stopReason Either Reversals, Max, or Min which are the three ways in which FT can terminate.

• np Number of presentation for the whole procedure (including both staircases if run).

References

A. Turpin, P.H. Artes and A.M. McKendrick. "The Open Perimetry Interface: An enabling tool for clinical visual psychophysics", Journal of Vision 12(11) 2012.

H. Bebie, F. Fankhauser and J. Spahr. "Static perimetry: strategies", Acta Ophthalmology 54 1976.

C.A. Johnson, B.C. Chauhan, and L.R. Shapiro. "Properties of staircase procedures for estimating thresholds in automated perimetry", Investigative Ophthalmology and Vision Science 33 1993.

See Also

dbTocd, opiPresent, fourTwo.start

Examples

# Stimulus is Size III white-on-white as in the HFA
makeStim <- function(db, n) {
  s <- list(x=9, y=9, level=dbTocd(db), size=0.43, color="white",
            duration=200, responseWindow=1500)
  class(s) <- "opiStaticStimulus"
  return(s)
}
chooseOpi("SimHenson")
if (!is.null(opiInitialize(type="C", cap=6)$err))
  stop("opiInitialize failed")

result <- FT(makeStim=makeStim, tt=30, fpr=0.15, fnr=0.01)
if (!is.null(opiClose()$err))
  warning("opiClose() failed")

##############################################
# This section is for multiple FTs
##############################################
makeStimHelper <- function(db,n, x, y) {  # returns a function of (db,n)
  ff <- function(db, n) db+n
  body(ff) <- substitute({
    s <- list(x=x, y=y, level=dbTocd(db), size=0.43, color="white",
              duration=200, responseWindow=1500)
    class(s) <- "opiStaticStimulus"
    return(s)
  }, list(x=x,y=y))
  return(ff)
}

# List of (x, y, true threshold) triples
locations <- list(c(9,9,30), c(-9,-9,32), c(9,-9,31), c(-9,9,33))
# Setup starting states for each location
states <- lapply(locations, function(loc) {
  FT.start(makeStim=makeStimHelper(db,n,loc[1],loc[2]),
           tt=loc[3], fpr=0.03, fnr=0.01)})

# Loop through until all states are "stop"
while(!all(st <- unlist(lapply(states, FT.stop)))) {
  i <- which(!st)                         # choose a random,
  i <- i[runif(1, min=1, max=length(i))]  # unstopped state
  r <- FT.step(states[[i]])               # step it
  states[[i]] <- r$state                  # update the states
}

finals <- lapply(states, FT.final)    # get final estimates of threshold
for(i in 1:length(locations)) {
  cat(sprintf("Location (%+2d,%+2d) ",locations[[i]][1], locations[[i]][2]))
  cat(sprintf("has threshold %4.2f\n", finals[[i]]))
}

if(!is.null(opiClose()$err))
  warning("opiClose() failed")

Description

Implementation of opiPresent for the Kowa AP7000 machine. Version for opiKineticStimulus.

This is for internal use only. Use opiPresent() with stim as class opiStaticStimulus and you will get the Value back.

Arguments

Details

stim is a list containing at least the following 3 elements:

• path, a list containing x 2 x-coordinate in degrees (floating point) (range $[-80,80]$) and y, list of 2 y-coordinate in degrees (floating point) (range $[-70,65]$).

• sizes list of 1 size; one of .opi_env$KowaAP7000$SIZES_DEGREES.

• colors list of 1 color; one of .opi_env$KowaAP7000$COLOR_WHITE, .opi_env$KowaAP7000$COLOR_GREEN, .opi_env$KowaAP7000$COLOR_BLUE, or .opi_env$KowaAP7000$COLOR_RED

• levels list of 1 level; luminance in cd/$\mbox{m}^2$, and is rounded to the nearest whole dB for display (range 0 to 50). 0dB is 10000aps.

• speeds list of 1 speed; degrees per second range $[3, 5]$.

• duration of stimulus on in milliseconds (range $[100, 1200]$).

• responseWindow from start of stimulus presentation in milliseconds (max 5000).

Value

A list containing

• err String message or NULL for no error.

• seen 1 if seen, 0 otherwise.

• time Reaction time (if seen).

• x Location of button press in degrees.

• y Location of button press in degrees.

Description

Implementation of opiPresent for the Kowa AP7000 machine. Version for opiStaticStimulus.

This is for internal use only. Use opiPresent() with stim as class opiStaticStimulus and you will get the Value back.

Arguments

Details

stim is a list containing at least the following 3 elements:

• x, x-coordinate in degrees (floating point) (range $[-80,80]$).

• y, y-coordinate in degrees (floating point) (range $[-70,65]$).

• level is luminance in cd/$\mbox{m}^2$, and is rounded to the nearest whole dB for display (range 0 to 50). 0dB is 10000aps.

• duration of stimulus on in milliseconds (range $[100, 1200]$).

• responseWindow from start of stimulus presentation in milliseconds (max 5000).

• size one of .opi_env$KowaAP7000$SIZES_DEGREES.

• color one of .opi_env$KowaAP7000$COLOR_WHITE, .opi_env$KowaAP7000$COLOR_GREEN, .opi_env$KowaAP7000$COLOR_BLUE, or .opi_env$KowaAP7000$COLOR_RED

Value

A list containing

• err String message or NULL for no error.

• seen 1 if seen, 0 otherwise.

• time Reaction time (if seen).

• pupilX

• pupilY

• purkinjeX

• purkinjeY

Description

Implementation of opiPresent for the Kowa AP7000 machine. Version for opiKineticStimulus.

This is for internal use only. Use opiPresent() with stim as class opiStaticStimulus and you will get the Value back.

Description

An implementation of Kontsevich and Tyler (Vis Res 39 (1999) pages 2729–2737 default parameterised for Standard Automated Perimetry. based on A. Turpin, D. Jankovic and A.M. McKendrick, "Identifying Steep Psychometric Function Slope Quickly in Clinical Applications", Vision Research, 50(23). November 2010. Pages 2476-2485

Usage

KTPsi(
  domains = list(slopes = 1:5, thresholds = 20:40, fps = c(0, 0.025, 0.05, 0.1, 0.2), fns
    = c(0, 0.025, 0.05, 0.1, 0.2)),
  priors = list(slopes = rep(1, length(domains$slopes))/length(domains$slopes),
    thresholds = rep(1, length(domains$thresholds))/length(domains$thresholds), fps =
    rep(1, length(domains$fps))/length(domains$fps), fns = rep(1,
    length(domains$fns))/length(domains$fns)),
  stimValues = 17:40,
  stopType = "N",
  stopValue = 140,
  maxPresentations = 200,
  minInterStimInterval = NA,
  maxInterStimInterval = NA,
  verbose = 0,
  makeStim,
  ...
)

KTPsi.start(
  domains = list(slopes = 1:5, thresholds = 20:40, fps = c(0, 0.025, 0.05, 0.1, 0.2), fns
    = c(0, 0.025, 0.05, 0.1, 0.2)),
  priors = list(slopes = rep(1, length(domains$slopes))/length(domains$slopes),
    thresholds = rep(1, length(domains$thresholds))/length(domains$thresholds), fps =
    rep(1, length(domains$fps))/length(domains$fps), fns = rep(1,
    length(domains$fns))/length(domains$fns)),
  stimValues = 17:40,
  stopType = "N",
  stopValue = 140,
  maxPresentations = 200,
  minInterStimInterval = NA,
  maxInterStimInterval = NA,
  verbose = 0,
  makeStim,
  ...
)

KTPsi.step(state, nextStim = NULL, fixedStimValue = NA)

KTPsi.final(state, method = "expectation")

KTPsi.stop(state)

Arguments

Details

The assumed psychometric function is the cumulative Gaussian:

$\mbox{fp}+(1-\mbox{fp}-\mbox{fn})(1-\mbox{pnorm}(x, \mbox{threshold}, \mbox{slope})$

hence domain$slopes are standard deviations and domain$thresholds are the mean.

While it is assumed that domains$thresholds and stimValues are in dB, this need not be the case. As long as the makeStim function converts stimValues into cd/$\mbox{m}^2$ for the opiPresent function, then any units should work.

The checkFixationOK function is called (if present in stim made from makeStim) after each presentation, and if it returns FALSE, the pdf for that state is not changed (ie the presentation is ignored), but the stim, number of presentations etc is recorded in the state.

If more than one KTPsi is to be interleaved (for example, testing multiple locations), then thePsi KTPsi.start, KTPsi.step, KTPsi.stop and KTPsi.final calls can maintain the state of the KTPsi after each presentation, and should be used. If only a single KTPsi is required, then the simpler KTPsi function can be used, which is a wrapper for the four functions that maintain state. See examples below.

Value

Single location

KTPsi returns a list containing

• npres Total number of presentations used.

• respSeq Response sequence stored as a matrix: row 1 is dB values of stimuli, row 2 is 1/0 for seen/not-seen, row 3 is fixated 1/0 (always 1 if checkFixationOK not present in stim objects returned from makeStim).

• pdfs If verbose is bigger than 0, then this is a list of the pdfs used for each presentation, otherwise NULL.

• final The mean/median/mode of the final pdf, depending on stimChoice, which is the determined threshold.

• opiResp A list of responses received from each successful call to opiPresent within KTPsi.

Multilple locations

KTPsi.start returns a list that can be passed to KTPsi.step, KTPsi.stop, and KTPsi.final. It represents the state of a KTPsi at a single location at a point in time and contains the following.

• name KTPsi

• A copy of all of the parameters supplied to KTPsi.start: domains, priors, stimValues, stopType, stopValue, maxPresentations, makeStim and opiParams.

• psi A matrix where psi[domain_index, stim] is the probability of seeing stim assuming the psychometric function for the domain index domain_index.

• labels A text representation of psi[domain_index, ], or the the psychometric function for the domain index domain_index.

• pdf Current pdf: vector of probabilities the same length as product of lengths of domain elements.

• numPresentations The number of times KTPsi.step has been called on this state.

• stimuli A vector containing the stimuli used at each call of KTPsi.step.

• responses A vector containing the responses received at each call of KTPsi.step.

• responseTimes A vector containing the response times received at each call of KTPsi.step.

• fixated A vector containing TRUE/FALSE if fixation was OK according to checkFixationOK for each call of KTPsi.step (defaults to TRUE if checkFixationOK not present).

• opiResp A list of responses received from each call to opiPresent within KTPsi.step.

KTPsi.step returns a list containing

• stat: The new state after presenting a stimuli and getting a response.

• resp The return from the opiPresent call that was made.

KTPsi.stop returns TRUE if the KTPsi has reached its stopping criteria, and FALSE otherwise.

KTPsi.final returns an estimate of threshold based on state based on its parameter.

TRUE if the state has reached its stopping criteria, and FALSE otherwise.

References

Kontsevich and Tyler. Vision Research 39 (1999) pages 2729–2737.

A. Turpin, D. Jankovic and A.M. McKendrick, "Identifying Steep Psychometric Function Slope Quickly in Clinical Applications", Vision Research, 50(23). November 2010. Pages 2476-2485

A. Turpin, P.H. Artes and A.M. McKendrick "The Open Perimetry Interface: An enabling tool for clinical visual psychophysics", Journal of Vision 12(11) 2012.

See Also

dbTocd, opiPresent

Examples

chooseOpi("SimGaussian")
if(!is.null(opiInitialize(sd = 2)$err))
  stop("opiInitialize failed")

     # This section is for single location KTPsi
     # Stimulus is Size III white-on-white as in the HFA
makeStim <- function(db, n) {
  s <- list(x=9, y=9, level=dbTocd(db), size=0.43, color="white",
            duration=200, responseWindow=1500, checkFixationOK=NULL)
  class(s) <- "opiStaticStimulus"
  return(s)
}

KTPsi(makeStim = makeStim, stopType="H", stopValue=  3, tt=30, fpr=0.03)
KTPsi(makeStim = makeStim, stopType="N", stopValue= 27, verbose = 0, tt=30, fpr=0.03)

     # For multiple locations...
## Not run: 
states <- lapply(1:10, function(loc) KTPsi.start(makeStim = makeStim))
unfinished <- 1:10
while (length(unfinished) > 0) {
     loc <- unfinished[[1]]
     states[[loc]] <- KTPsi.step(states[[loc]])$state
     if (KTPsi.stop(states[[loc]]))
         unfinished <- unfinished[-1]
}

## End(Not run)

Description

MOCS performs either a yes/no or n-interval-forced-choice Method of Constant Stimuli test

Usage

MOCS(
  params = NA,
  order = "random",
  responseWindowMeth = "constant",
  responseFloor = 1500,
  responseHistory = 5,
  keyHandler = function(correct, ret) return(list(seen = TRUE, time = 0, err = NULL)),
  interStimMin = 200,
  interStimMax = 500,
  beep_function,
  makeStim,
  stim_print,
  ...
)

Arguments

Details

Whether the test is yes/no or forced-choice is determined by the number of columns in params. The code simply presents all columns from 5 onwards and collects a response at the end. So if there is only 5 columns, it is a yes/no task. If there are 6 columns it is a 2-interval-forced-choice. Generally, an nIFC experiment has 4+n columns in params.

Note that when the order is "random", the number of trials in the test will be the sum of the 3rd column of params. When the order is "fixed", there is only one presentation per row, regardless of the value in the 3rd column of params.

If a response is received before the final trial in a nIFC experiment, it is ignored.

If the checkFixationOK function is present in a stimulus, then it is called after each presentation, and the result is “anded” with each stimulus in a trial to get a TRUE/FALSE for fixating on all stimuli in a trial.

Value

Returns a data.frame with one row per stimulus copied from params with extra columns appended: checkFixation checks, and the return values from opiPresent() (see example). These last values will differ depending on which machine/simulation you are running (as chosen with chooseOpi().

• column 1: x

• column 2: y

• column 3: location number

• column 4: number of times to repeat this stim

• column 5: correct stimulus index

• column 6: TRUE/FALSE was fixating for all presentations in this trial according to checkFixationOK

• column 7...: columns from params

• ...: columns from opiPresent return

References

A. Turpin, P.H. Artes and A.M. McKendrick. "The Open Perimetry Interface: An enabling tool for clinical visual psychophysics", Journal of Vision 12(11) 2012.

See Also

dbTocd, opiPresent

Examples

# For the Octopus 900
# Check if pupil centre is within 10 pixels of (160,140)
checkFixationOK <- function(ret) return(sqrt((ret$pupilX - 160)^2 + (ret$pupilY - 140)^2) < 10)

# Return a list of opi stim objects (list of class opiStaticStimulus) for each level (dB) in
# p[5:length(p)]. Each stim has responseWindow BETWEEN_FLASH_TIME, except the last which has
# rwin. This one assumes p is on old Octopus 900 dB scale (0dB == 4000 cd/m^2).
makeStim <- function(p, rwin) {
  BETWEEN_FLASH_TIME <- 750   # ms
  res <- NULL
  for(i in 5:length(p)) {
    s <- list(x=p[1], y=p[2], level=dbTocd(p[i],4000/pi), size=0.43, duration=200,
              responseWindow=ifelse(i < length(p), BETWEEN_FLASH_TIME, rwin),
              checkFixationOK=NULL)
    class(s) <- "opiStaticStimulus"
    res <- c(res, list(s))
  }
  return(res)
}

################################################################
# Read in a key press 'z' is correct==1, 'm' otherwise
#    correct is either 1 or 2, whichever is the correct interval
#
# Return list(seen={TRUE|FALSE}, time=time, err=NULL))
#        seen is TRUE if correct key pressed
################################################################
## Not run: 
  if (length(dir(".", "getKeyPress.py")) < 1)
    stop('Python script getKeyPress.py missing?')

## End(Not run)

keyHandler <- function(correct, ret) {
  return(list(seen=TRUE, time=0, err=NULL))
  ONE <- "b'z'"
  TWO <- "b'm'"
  time <- Sys.time()
  key <- 'q'
  while (key != ONE && key != TWO) {
    a <- system('python getKeyPress.py', intern=TRUE)
    key <- a # substr(a, nchar(a), nchar(a))
    print(paste('Key pressed: ',key,'from',a))
    if (key == "b'8'")
      stop('Key 8 pressed')
  }
  time <- Sys.time() - time
  if ((key == ONE && correct == 1) || (key == TWO && correct == 2))
    return(list(seen=TRUE, time=time, err=NULL))
  else
    return(list(seen=FALSE, time=time, err=NULL))
}

################################################################
# Read in return value from opipresent with F310 controller.
# First param is correct, next is 1 for left button, 2 for right button
# Left button (LB) is correct for interval 1, RB for interval 2
#    correct is either 1 or 2, whichever is the correct interval
#
# Return list(seen={TRUE|FALSE}, time=time, err=NULL))
#        seen is TRUE if correct key pressed
################################################################
F310Handler <- function(correct, opiResult) {
  z <- opiResult$seen == correct
  opiResult$seen <- z
  return(opiResult)
}

################################################################
# 2 example beep_function
################################################################
## Not run: 
  require(beepr)
  myBeep <- function(type='None') {
    if (type == 'correct') {
      beepr::beep(2)  # coin noise
      Sys.sleep(0.5)
    }
    if (type == 'incorrect') {
      beepr::beep(1) # system("rundll32 user32.dll,MessageBeep -1") # system beep
      #Sys.sleep(0.0)
    }
  }
  require(audio)
  myBeep <- function(type="None") {
    if (type == 'correct') {
      wait(audio::play(sin(1:10000/10)))
    }
    if (type == 'incorrect') {
      wait(audio::play(sin(1:10000/20)))
    }
  }

## End(Not run)

################################################################
# An example stim_print function
################################################################
## Not run: 
  stim_print <- function(s, ret) {
    sprintf("%4.1f %2.0f",cdTodb(s$level,10000/pi), ret$seen)
  }

## End(Not run)

Description

Implementation of opiPresent for Kinetic stimuli on the Octopus090 machine.

This is for internal use only. Use opiPresent() with these Arguments and stim as class opiKineticStimulus and you will get the Value back.

Arguments

Details

stim is a list containing at least the following 3 elements:

• path, A list of $(x,y)$ coordinates in degrees that is usable by xy.coords().

• sizes, A list where sizes[i] is the size of stimulus (diameter in degrees) to use for the section of path specified by path[i]..path[i+1]. Rounded to nearest Goldmann size.

• levels A list where levels[i] is the stimulus luminance in cd/$\mbox{m}^2$ to use for the section of path specified by path[i]..path[i+1].

• speeds A list where speeds[i] is the speed in degrees per second to use for the section of path specified by path[i]..path[i+1].

Value

A list containing

• err String message or NULL for no error.

• seen 1 if seen, 0 otherwise.

• time Reaction time (if seen).

• x Coordinate where button was pressed (degrees - i guess).

• y Coordinate where button was pressed (degrees - i guess).

Description

Implementation of opiPresent for the Octopus090 machine. Version for opiStaticStimulus.

This is for internal use only. Use opiPresent() with stim as class opiStaticStimulus and you will get the Value back.

Arguments

Details

stim is a list containing at least the following 3 elements:

• x, x-coordinate in degrees (floating point) (range $[-30,30]$).

• y, y-coordinate in degrees (floating point) (range $[-30,30]$).

• level is luminance in cd/$\mbox{m}^2$, and is rounded to the nearest whole dB for display (range 0 to 50). 0dB is 10000aps.

It can also contain:

• responseWindow from start of stimulus presentation in milliseconds (default is 1500).

• duration of stimulus on in milliseconds (default 200).

• color one of .opi_env$O900$STIM_WHITE, .opi_env$O900$STIM_BLUE or .opi_env$O900$STIM_RED. It must be same as that initialised by opiSetup() or opiInitialize() (default .opi_env$O900$STIM_WHITE).

• size of stimulus diameter in degrees (default Size III == 0.43). This is rounded to the nearest support Goldmann size.

If responses are taken from the F310 Controller then

• If the L button is pressed, seen is set to 1.

• If the R button is pressed, seen is set to 2.

• If no button is pressed within responseWindow, then seen is set to 0.

If stim is null, always return err = NULL status.

Value

A list containing

• err String message or NULL for no error.

• seen 1 if seen, 0 otherwise. (See details for F310)

• time Reaction time (if seen).

Description

Implementation of opiPresent for Temporal stimuli on the Octopus090 machine.

This is for internal use only. Use opiPresent() with these Arguments (stim as class opiTemporalStimulus) and you will get the Value back.

Arguments

Details

stim is a list containing at least the following 3 elements:

• x, x-coordinate in degrees (floating point) (range $[-30,30]$).

• y, y-coordinate in degrees (floating point) (range $[-30,30]$).

• rate is frequency in Hz.

It can also contain:

• responseWindow from start of stimulus presentation in milliseconds (default is 1500).

• duration of stimulus on in milliseconds (default 200).

• size of stimulus diameter in degrees (default Size III == 0.43). This is rounded to the nearest support Goldmann size.

If stim is null, always return err = NULL status.

Value

A list containing

• err String message or NULL for no error.

• seen 1 if seen, 0 otherwise.

• time Reaction time (if seen).

Description

Internal use only.

Usage

open_socket(ip, port, machineName)

Arguments

Value

Socket or NULL on error

Description

Specific parameters and return values can be seen in the machine specific versions listed below in the ’See Also’.

Usage

opiClose()

Value

Each implementation should(!) return a list with at least the following elements:

• err NULL if no error, otherwise a string describing the error.

See Also

opiClose_for_Compass(), opiClose_for_Octopus900(), opiClose_for_ImoVifa(), opiClose_for_PhoneHMD(), opiClose_for_Display(), opiClose_for_PicoVR(), opiClose_for_SimNo(), opiClose_for_SimYes(), opiClose_for_SimHenson(), opiClose_for_SimGaussian()

Description

This is for internal use only. Use opiClose() with these Arguments and you will get the Value back.

Value

A list with elements

• err, which is an error code, NULL for no error

• fixations, which is a matrix one row per fixation and three columns:

  • time (same as time_hw in opiPresent)

  • x (degrees relative to the centre of the image returned by opiInitialise - not the PRL)

  • y (as for x)

Description

This is for internal use only. Use opiClose() after chooseOPI("Display") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiClose()

Examples

chooseOpi("Display")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiClose()

Description

This is for internal use only. Use opiClose() after chooseOPI("ImoVifa") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiClose()

Examples

chooseOpi("ImoVifa")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiClose()

Description

This is for internal use only. Use opiClose() with the same parameters.

Value

Returns list(err = NULL).

Description

This is for internal use only. Use opiClose() with the same parameters.

Value

Returns list(err = NULL).

Description

This is for internal use only. Use opiClose() with the same parameters.

Value

Returns list(err = NULL).

Description

This is for internal use only. Use opiClose() after chooseOPI("PhoneHMD") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiClose()

Examples

chooseOpi("PhoneHMD")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiClose()

Description

This is for internal use only. Use opiClose() after chooseOPI("PicoVR") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiClose()

Examples

chooseOpi("PicoVR")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiClose()

Description

Does nothing.

Value

A list with elements:

• error Always FALSE.

• msg A string "Close OK".

Description

Does nothing.

Value

A list with elements:

• err Always NULL.

Description

Does nothing.

Value

A list with elements:

• err Always NULL.

Description

Does nothing.

Value

A list with elements:

• error Always FALSE.

• msg A string "Close OK".

Description

Does nothing.

Value

A list with elements:

• err Always NULL.

Description

Establishes connection with the device and a Monitor (aka Server) if appropriate. Sends any startup parameters that might be needed by the machine. Specific parameters and return values can be seen in the machine specific versions listed below in the ’See Also’.

Usage

opiInitialise(...)

opiInitialize(...)

Arguments

Value

A list containing at least the following elements:

• err NULL if no error, otherwise a string describing the error.

See Also

opiInitialise_for_ImoVifa(), opiInitialise_for_PhoneHMD(), opiInitialise_for_Display(), opiInitialise_for_PicoVR(), opiInitialise_for_Octopus900(), opiInitialise_for_Compass(), opiInitialise_for_SimNo(), opiInitialise_for_SimYes(), opiInitialise_for_SimHenson(), opiInitialise_for_SimGaussian()

Description

This is for internal use only. Use opiInitialise() with these Arguments and you will get the Value back.

Arguments

Details

Warning: this returns a list, not a single error code.

Value

A list with elements:

• err NULL if successful, not otherwise.

• prl A pair giving the (x,y) in degrees of the Preferred Retinal Locus detected in the initial alignment.

• onh a pair giving the (x,y) in degrees of the ONH as selected by the user.

• image raw bytes being the JPEG compressed infra-red image acquired during alignment.

Examples

## Not run: 
  # Set up the Compass
  chooseOpi("Compass")
  result <- opiInitialize(ip = "192.168.1.7", port = 44965)
  if (is.null(result$err))
    print(result$prl)

## End(Not run)

Description

This is for internal use only. Use opiInitialise() after chooseOPI("Display") to call this function.

Arguments

Details

port can take on values in the range [0, 65535].

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiInitialise()

Examples

chooseOpi("Display")
result <- opiInitialise(address = list(port = 50001, ip = "localhost"))

Description

This is for internal use only. Use opiInitialise() after chooseOPI("ImoVifa") to call this function.

Arguments

Details

port can take on values in the range [0, 65535].

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiInitialise()

Examples

chooseOpi("ImoVifa")
result <- opiInitialise(address = list(port = 50001, ip = "localhost"))

Description

This is for internal use only. Use opiInitialise() with these Arguments and you will get the Value back.

Arguments

Details

If the chosen OPI implementation is KowaAP7000, then you must specify the IP address and port of the AP-7000 server.

Value

list(err = NULL)

Examples

## Not run: 
  # Set up the Kowa AP-7000
  chooseOpi("KowaAP7000")
  opiInitialize(ip="192.168.1.7", port=44965)

## End(Not run)

Description

This is for internal use only. Use opiInitialise() with these Arguments and you will get the Value back.

Arguments

* 0 is off
* 1 is eye blink
* 2 is eye blink, forehead rest, fixation control
* 3 is eye blink, forehead rest, fixation control, fast eye movements

Details

The default background and stimulus setup is to white-on-white perimetry.

Uses port 50000 on the O600.

Examples

## Not run: 
  # Set up the O600
  chooseOpi("O600")
  opiInitialize(ip="192.168.1.7", eye = "left")

## End(Not run)

Description

This is for internal use only. Use opiInitialise() with these Arguments and you will get the Value back.

Arguments

Details

If the chosen OPI implementation is Octopus900, then you must specify a directory and the eye to be tested.

serverPort is the TCP/IP port on which the server is listening (on localhost).

eyeSuiteSettingsLocation is the folder name containing the EyeSuite setting files, and should include the trailing slash.

eye must be either "left" or "right".

gazeFeed is the name of an existing folder into which the video frames of eye tracker are recorded. Set to NA for no recording.

bigWheel is FALSE for a standard Octopus 900 machine. Some research machines are fitted with an alternate aperture wheel that has 24 sizes, which are accessed with bigWheel is TRUE. The mapping from size to 'hole on wheel' is hard coded; see code for details.

If pres_buzzer is greater than zero, a buzzer will sound with each stimuli presented.

If resp_buzzer is greater than zero, a buzzer will sound with each button press (response). The volume can be one of 0 (no buzzer), 1, 2, or 3 (max volume). If both buzzers are more than zero, the maximum of the two will be used as the volume.

If zero_dB_is_10000_asb is TRUE then 0 dB is taken as 10000 apostilbs, otherwise 0 dB is taken as 4000 apostilbs.

Value

A list containing err which is

• NULL if successful

• 1 if Octopus900 is already initialised by a previous call to opiInitialize

• 2 if some error occurred that prevented initialisation.

Examples

## Not run: 
  chooseOpi("Octopus900")
  res <- opiInitialize(serverPort = 50001,
                       eyeSuiteSettingsLocation = "C:/ProgramData/Haag-Streit/EyeSuite/",
                       eye = "", gazeFeed = "", bigWheel = FALSE,
                       pres_buzzer = 0, resp_buzzer = 0, zero_dB_is_10000_asb = TRUE)
  if (!is.null(res$err))
      stop("opiInitialize failed")

## End(Not run)

Description

This is for internal use only. Use opiInitialise() after chooseOPI("PhoneHMD") to call this function.

Arguments

Details

port can take on values in the range [0, 65535].

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiInitialise()

Examples

chooseOpi("PhoneHMD")
result <- opiInitialise(address = list(port = 50001, ip = "localhost"))

Description

This is for internal use only. Use opiInitialise() after chooseOPI("PicoVR") to call this function.

Arguments

Details

port can take on values in the range [0, 65535].

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiInitialise()

Examples

chooseOpi("PicoVR")
result <- opiInitialise(address = list(port = 50001, ip = "localhost"))

Description

Simulates responses using a Frequency of Seeing (FoS) curve.

The FoS is modelled as a cumulative Gaussian function with standard deviation equal to sd as provided and the mean as the true threshold given as tt opiPresent. All values are in dB relative to maxStim.

This is for internal use only, use opiInitialize().

Arguments

Value

A list with elements:

• err NULL if initialised, a message otherwise

Examples

# Set up a simple simulation for white-on-white perimetry
chooseOpi("SimGaussian")
res <- opiInitialize(sd = 2.5)
if (!is.null(res$err))
  stop(paste("opiInitialize() failed:", res$msg))

Description

Simulates responses using a Frequency of Seeing (FoS) curve.

For internal use only, use opiInitialize().

The FoS is modelled as a cumulative Gaussian function over dB with standard deviation equal to min(cap, exp( A * t + B)), where t is the threshold/mean of the FoS in dB. All values are in dB relative to maxStim.

Arguments

Value

A list with elements:

• err NULL if initialised, msg otherwise

Examples

# Set up a simple simulation for white-on-white perimetry
chooseOpi("SimHenson")
res <- opiInitialize(type = "C", cap = 6)
if (!is.null(res$err))
  stop(paste("opiInitialize() failed:", res$err))

Description

Simulates responses using a Frequency of Seeing (FoS) curve.

For internal use only, use opiInitialize().

The FoS is modelled as a cumulative Gaussian function over dB with standard deviation equal to min(cap, exp( A * t + B)), where t is the threshold/mean of the FoS in dB. All values are in dB relative to maxStim.

Arguments

Details

If the chosen OPI implementation is SimHensonRT, then the first six parameters are as in SimHenson, and rtData is a data frame with at least 2 columns: "Rt", response time; and "Dist", signifying that distance between assumed threshold and stimulus value in your units.

This package contains RtSigmaUnits or RtDbUnits that can be loaded with the commands data(RtSigmaUnits) or data(RtDbUnits), and are suitable to pass as values for rtData.

rtFp gives the vector of values in milliseconds from which a response time for a false positive response is randomly sampled.

Value

A list with elements:

• err NULL if initialised, string msg otherwise

Examples

# Set up a simple simulation for white-on-white perimetry
# and display the stimuli in a plot region and simulate response times
chooseOpi("SimHensonRT")
data(RtSigmaUnits)
oi <- opiInitialize(type="C", cap=6, display=NA, rtData=RtSigmaUnits, rtFP=1:100)
if (!is.null(oi$err))
  stop("opiInitialize failed")

    # Set up a simple simulation for white-on-white perimetry
chooseOpi("SimHenson")
res <- opiInitialize(type = "C", cap = 6)
if (!is.null(res$err))
  stop(paste("opiInitialize() failed:", res$err))

Description

Does nothing.

Arguments

Value

A list with elements:

• err Always NULL.

Description

Does nothing.

Arguments

Value

A list with elements:

• err Always NULL.

Description

For backwards compatibility. Used by Octopus900 and KowaAP7000.

Usage

opiKineticStimulus()

Description

Specific parameters and return values can be seen in the machine specific versions listed below in the ’See Also’.

Usage

opiPresent(stim, ...)

Arguments

Value

Each implementation should(!) return a list with at least the following elements:

• err NULL if no error, otherwise a string describing the error.

• seen TRUE if stimulus seen, FALSE otherwise

• time Response time from onset of stimulus in milliseconds.

See Also

opiPresent_for_Compass(), opiPresent_for_Octopus900(), opiPresent_for_ImoVifa(), opiPresent_for_PhoneHMD(), opiPresent_for_Display(), opiPresent_for_PicoVR(), opiPresent_for_SimNo(), opiPresent_for_SimYes(), opiPresent_for_SimHenson(), opiPresent_for_SimGaussian()

Description

This is for internal use only. Use opiSetup() with these Arguments and you will get the Value back.

Arguments

Details

If the chosen OPI implementation is Compass, then nextStim is ignored. Note that the dB level is rounded to the nearest integer.

If tracking is on, then this will block until the tracking is obtained, and the stimulus presented.

stim is a list containing some or all of the following elements:

• x, x-coordinate in degrees (floating point) (range $[-30,30]$).

• y, y-coordinate in degrees (floating point) (range $[-30,30]$).

• level is luminance in cd/$\mbox{m}^2$, and is rounded to the nearest whole dB for display (range 0 to 50). 0dB is 10000aps.

• responseWindow is in milliseconds (range 0 to 2680).

Stimulus duration is assumed to be 200ms, and size is assumed to be Goldmann III (0.43 degrees diameter), color is assumed to be white. These cannot be changed.

Value

A list containing:

• err 0 all clear, >= 1 some error codes (eg cannot track, etc) (integer)

• seen FALSE for not seen, TRUE for seen (button pressed in response window)

• time Response time in ms (integer) since stimulus onset, -1 for not seen

• time_rec Time since epoch when command was received at Compass (integer ms)

• time_pres Time since epoch that stimulus was presented (integer ms)

• num_track_events Number of tracking events that occurred during presentation (integer)

• num_motor_fails Number of times motor could not follow fixation movement during presentation (integer)

• pupil_diam Pupil diameter in mm (float)

• loc_x Pixels integer, location in image of presentation (integer)

• loc_y Pixels integer, location in image of presentation (integer)

Examples

## Not run: 
  # Set up the Compass
  chooseOpi("Compass")
  result <- opiInitialize(ip = "192.168.1.7", port = 44965)
  if (!is.null(result$err))
    stop("Initialisation failed")

   #' @param x X location of stim in degrees
   #' @param y Y location of stim in degrees
   #' @param size If 3, Goldmann III, else V
   #' @param db Value in dB
   #' @return stim object ready for opiPresent
   makeStim <- function(x, y, size, db) {
       s <- list(x = x, y = y, level = dbTocd(db, 10000 / pi),
           size = ifelse(size == 3, 0.43, 1.77),
           duration = 200, responseWindow = 1500)
       class(s) <- "opiStaticStimulus"
       return(s)
   }

   result <- opiPresent(makeStim(9, 9, 3, 10))

## End(Not run)

Description

This is for internal use only. Use opiPresent() after chooseOPI("Display") to call this function.

Arguments

Details

Elements in lum can take on values in the range [0.0, 1.0E10].

stim.length can take on values in the range [1, 2147483647].

Elements in color1 can take on values in the range [0.0, 1.0].

Elements in sx can take on values in the range [0.0, 180.0].

Elements in sy can take on values in the range [0.0, 180.0].

Elements in eye can take on values in the set {"left", "right", "both", "none"}.

Elements in t can take on values in the range [0.0, 1.0E10].

w can take on values in the range [0.0, 1.0E10].

Elements in x can take on values in the range [-90.0, 90.0].

Elements in y can take on values in the range [-90.0, 90.0].

Elements in envSdx can take on values in the range [-1.0E10, 1.0E10].

Elements in envSdy can take on values in the range [-1.0E10, 1.0E10].

Elements in envRotation can take on values in the range [-1.0E10, 1.0E10].

Elements in type can take on values in the set {"flat", "checkerboard", "sine", "squaresine", "g1", "g2", "g3", "text", "image"}.

Elements in frequency can take on values in the range [0.0, 300.0].

Elements in color2 can take on values in the range [0.0, 1.0].

Elements in fullFoV can take on values in the range [-1.0E10, 1.0E10].

Elements in phase can take on values in the range [0.0, 1.0E10].

Elements in shape can take on values in the set {"triangle", "square", "polygon", "hollow_triangle", "hollow_square", "hollow_polygon", "cross", "maltese", "circle", "annulus", "optotype", "text", "model"}.

Elements in rotation can take on values in the range [0.0, 360.0].

Elements in texRotation can take on values in the range [0.0, 360.0].

Elements in defocus can take on values in the range [0.0, 1.0E10].

Elements in envType can take on values in the set {"none", "square", "circle", "gaussian"}.

Elements in contrast can take on values in the range [0.0, 1.0].

Elements in optotype can take on values in the set {"a", "b", "c", "d", "e", "f", "g", "h", "i", "j", "k", "l", "m", "n", "o", "p", "q", "r", "s", "t", "u", "v", "w", "x", "y", "z"}.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

• time Response time from stimulus onset if button pressed (ms).

• seen '1' if seen, '0' if not.

See Also

opiPresent()

Examples

chooseOpi("Display")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiPresent(stim = list(lum = list(300.0), stim.length = 1, color1 = list(list(1.0,
                  1.0, 1.0)), sx = list(1.72), sy = list(1.72),
                  eye = list("LEFT"), t = list(200.0), w = 1500.0, x = list(0.0), y = list(0.0)))

Description

This is for internal use only. Use opiPresent() after chooseOPI("ImoVifa") to call this function.

Arguments

Details

Elements in lum can take on values in the range [0.0, 1.0E10].

stim.length can take on values in the range [1, 2147483647].

Elements in color1 can take on values in the range [0.0, 1.0].

Elements in sx can take on values in the range [0.0, 180.0].

Elements in sy can take on values in the range [0.0, 180.0].

Elements in eye can take on values in the set {"left", "right", "both", "none"}.

Elements in t can take on values in the range [0.0, 1.0E10].

w can take on values in the range [0.0, 1.0E10].

Elements in x can take on values in the range [-90.0, 90.0].

Elements in y can take on values in the range [-90.0, 90.0].

Elements in envSdx can take on values in the range [-1.0E10, 1.0E10].

Elements in envSdy can take on values in the range [-1.0E10, 1.0E10].

Elements in envRotation can take on values in the range [-1.0E10, 1.0E10].

Elements in type can take on values in the set {"flat", "checkerboard", "sine", "squaresine", "g1", "g2", "g3", "text", "image"}.

Elements in frequency can take on values in the range [0.0, 300.0].

Elements in color2 can take on values in the range [0.0, 1.0].

Elements in fullFoV can take on values in the range [-1.0E10, 1.0E10].

Elements in phase can take on values in the range [0.0, 1.0E10].

Elements in shape can take on values in the set {"triangle", "square", "polygon", "hollow_triangle", "hollow_square", "hollow_polygon", "cross", "maltese", "circle", "annulus", "optotype", "text", "model"}.

Elements in rotation can take on values in the range [0.0, 360.0].

Elements in texRotation can take on values in the range [0.0, 360.0].

Elements in defocus can take on values in the range [0.0, 1.0E10].

Elements in envType can take on values in the set {"none", "square", "circle", "gaussian"}.

Elements in contrast can take on values in the range [0.0, 1.0].

Elements in optotype can take on values in the set {"a", "b", "c", "d", "e", "f", "g", "h", "i", "j", "k", "l", "m", "n", "o", "p", "q", "r", "s", "t", "u", "v", "w", "x", "y", "z"}.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

• eyedStart Diameter of pupil at stimulus onset (pixels).

• eyexEnd x co-ordinate of pupil at button press (pixels from image centre. Image is 640x480, left < 0). Note that for multi-part stimuli (t=0), the eye taken is the last eye in the list of components.

• eyeyEnd y co-ordinate (pixels). See eyexEnd for more details. (up > 0)

• eyedEnd Diameter of pupil at button press or response window expiry (pixels).

• eyexStart x co-ordinates of pupil at stimulus onset (pixels from image centre. Image is 640x480, left < 0). For a multi-part stimulus (t=0), the eye taken is the first eye in the list of components.

• time Response time from stimulus onset if button pressed (ms).

• seen '1' if seen, '0' if not.

• eyeyStart y co-ordinates (pixels from image centre). See eyexStart for more details. (up > 0)

See Also

opiPresent()

Examples

chooseOpi("ImoVifa")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiPresent(stim = list(lum = list(300.0), stim.length = 1, color1 = list(list(1.0,
                  1.0, 1.0)), sx = list(1.72), sy = list(1.72),
                  eye = list("LEFT"), t = list(200.0), w = 1500.0, x = list(0.0), y = list(0.0)))

Description

This is for internal use only. Use opiPresent() with the same arguments and the class of stim as one of opiStaticStimulus, opiTemporalStimulus, or opiKineticStimulus.

Arguments

Value

See kowa.presentStatic, kowa.presentTemporal, or kowa.presentKinetic.

See Also

kowa.presentStatic, kowa.presentTemporal, kowa.presentKinetic

Examples

## Not run: 
  chooseOpi("KowaAP7000")
  if (!is.null(opiInitialize()$err))
      stop("opiInitialize failed")
  s <- list(x=9, y=9, level=dbTocd(db), size=0.43, color="white",
            duration=200, responseWindow=1500, checkFixationOK=NULL)
  class(s) <- "opiStaticStimulus"
  print(opiPresent(s, NULL))

## End(Not run)

Description

This is for internal use only. Use opiPresent() with the same arguments.

Arguments

Value

A list containing

• err String message or NULL for no error.

• seen 1 if seen, 0 otherwise.

• time Reaction time (if seen).

Examples

## Not run: 
  chooseOpi("O600")
  if (!is.null(opiInitialize()$err))
      stop("opiInitialize failed")
  s <- list(x=9, y=9, level=dbTocd(db), size=0.43, color="white",
            duration=200, responseWindow=1500, checkFixationOK=NULL)
  print(opiPresent(s, NULL))

## End(Not run)

Description

This is for internal use only. Use opiPresent() with the same arguments and the class of stim as one of opiStaticStimulus, opiTemporalStimulus, or opiKineticStimulus.

Arguments

See Also

octo900.presentStatic, octo900.presentTemporal, octo900.presentKinetic

Examples

## Not run: 
  chooseOpi("Octopus900")
  if (!is.null(opiInitialize()$err))
      stop("opiInitialize failed")
  s <- list(x=9, y=9, level=dbTocd(db), size=0.43, color="white",
            duration=200, responseWindow=1500, checkFixationOK=NULL)
  class(s) <- "opiStaticStimulus"
  print(opiPresent(s, NULL))

## End(Not run)

Description

This is for internal use only. Use opiPresent() after chooseOPI("PhoneHMD") to call this function.

Arguments

Details

Elements in lum can take on values in the range [0.0, 1.0E10].

stim.length can take on values in the range [1, 2147483647].

Elements in color1 can take on values in the range [0.0, 1.0].

Elements in sx can take on values in the range [0.0, 180.0].

Elements in sy can take on values in the range [0.0, 180.0].

Elements in eye can take on values in the set {"left", "right", "both", "none"}.

Elements in t can take on values in the range [0.0, 1.0E10].

w can take on values in the range [0.0, 1.0E10].

Elements in x can take on values in the range [-90.0, 90.0].

Elements in y can take on values in the range [-90.0, 90.0].

Elements in envSdx can take on values in the range [-1.0E10, 1.0E10].

Elements in envSdy can take on values in the range [-1.0E10, 1.0E10].

Elements in envRotation can take on values in the range [-1.0E10, 1.0E10].

Elements in type can take on values in the set {"flat", "checkerboard", "sine", "squaresine", "g1", "g2", "g3", "text", "image"}.

Elements in frequency can take on values in the range [0.0, 300.0].

Elements in color2 can take on values in the range [0.0, 1.0].

Elements in fullFoV can take on values in the range [-1.0E10, 1.0E10].

Elements in phase can take on values in the range [0.0, 1.0E10].

Elements in shape can take on values in the set {"triangle", "square", "polygon", "hollow_triangle", "hollow_square", "hollow_polygon", "cross", "maltese", "circle", "annulus", "optotype", "text", "model"}.

Elements in rotation can take on values in the range [0.0, 360.0].

Elements in texRotation can take on values in the range [0.0, 360.0].

Elements in defocus can take on values in the range [0.0, 1.0E10].

Elements in envType can take on values in the set {"none", "square", "circle", "gaussian"}.

Elements in contrast can take on values in the range [0.0, 1.0].

Elements in optotype can take on values in the set {"a", "b", "c", "d", "e", "f", "g", "h", "i", "j", "k", "l", "m", "n", "o", "p", "q", "r", "s", "t", "u", "v", "w", "x", "y", "z"}.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

• time Response time from stimulus onset if button pressed (ms).

• seen '1' if seen, '0' if not.

See Also

opiPresent()

Examples

chooseOpi("PhoneHMD")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiPresent(stim = list(lum = list(300.0), stim.length = 1, color1 = list(list(1.0,
                  1.0, 1.0)), sx = list(1.72), sy = list(1.72),
                  eye = list("LEFT"), t = list(200.0), w = 1500.0, x = list(0.0), y = list(0.0)))

Description

This is for internal use only. Use opiPresent() after chooseOPI("PicoVR") to call this function.

Arguments

Details

Elements in lum can take on values in the range [0.0, 1.0E10].

stim.length can take on values in the range [1, 2147483647].

Elements in color1 can take on values in the range [0.0, 1.0].

Elements in sx can take on values in the range [0.0, 180.0].

Elements in sy can take on values in the range [0.0, 180.0].

Elements in eye can take on values in the set {"left", "right", "both", "none"}.

Elements in t can take on values in the range [0.0, 1.0E10].

w can take on values in the range [0.0, 1.0E10].

Elements in x can take on values in the range [-90.0, 90.0].

Elements in y can take on values in the range [-90.0, 90.0].

Elements in envSdx can take on values in the range [-1.0E10, 1.0E10].

Elements in envSdy can take on values in the range [-1.0E10, 1.0E10].

Elements in envRotation can take on values in the range [-1.0E10, 1.0E10].

Elements in type can take on values in the set {"flat", "checkerboard", "sine", "squaresine", "g1", "g2", "g3", "text", "image"}.

Elements in frequency can take on values in the range [0.0, 300.0].

Elements in color2 can take on values in the range [0.0, 1.0].

Elements in fullFoV can take on values in the range [-1.0E10, 1.0E10].

Elements in phase can take on values in the range [0.0, 1.0E10].

Elements in shape can take on values in the set {"triangle", "square", "polygon", "hollow_triangle", "hollow_square", "hollow_polygon", "cross", "maltese", "circle", "annulus", "optotype", "text", "model"}.

Elements in rotation can take on values in the range [0.0, 360.0].

Elements in texRotation can take on values in the range [0.0, 360.0].

Elements in defocus can take on values in the range [0.0, 1.0E10].

Elements in envType can take on values in the set {"none", "square", "circle", "gaussian"}.

Elements in contrast can take on values in the range [0.0, 1.0].

Elements in optotype can take on values in the set {"a", "b", "c", "d", "e", "f", "g", "h", "i", "j", "k", "l", "m", "n", "o", "p", "q", "r", "s", "t", "u", "v", "w", "x", "y", "z"}.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

• time Response time from stimulus onset if button pressed (ms).

• seen '1' if seen, '0' if not.

See Also

opiPresent()

Examples

chooseOpi("PicoVR")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiPresent(stim = list(lum = list(300.0), stim.length = 1, color1 = list(list(1.0,
                  1.0, 1.0)), sx = list(1.72), sy = list(1.72),
                  eye = list("LEFT"), t = list(200.0), w = 1500.0, x = list(0.0), y = list(0.0)))

Description

Determine the response to a stimuli by sampling from a cumulative Gaussian Frequency-of-Seeing (FoS) curve (also known as the psychometric function).

The FoS has formula

$\mbox{fpr}+(1-\mbox{fpr}-\mbox{fnr})(1-\mbox{pnorm}(x, \mbox{tt}, \mbox{sd})$

where $x$ is the stimulus value in dB, and sd is set by opiInitialize and tt, fpr and fnr are parameters.

This is for internal use only, use opiPresent().

Arguments

Value

A list with elements:

• err NULL if no error, a string message otherwise.

• seen TRUE or FALSE.

• time Always NA.

Examples

# Stimulus is Size III white-on-white as in the HFA
chooseOpi("SimGaussian")
res <- opiInitialize(sd = 1.6)
if (!is.null(res$err))
  stop(paste("opiInitialize() failed:", res$err))

result <- opiPresent(stim = list(level = dbTocd(20)), tt = 30, fpr = 0.15, fnr = 0.01)
print(paste("Seen:", result$seen, quote = FALSE))

res <- opiClose()
if (!is.null(res$err))
  warning(paste("opiClose() failed:", res$err))

Description

Determine the response to a stimuli by sampling from a cumulative Gaussian Frequency-of-Seeing (FoS) curve (also known as the psychometric function).

For internal use only, use opiPresent().

Arguments

Details

The FoS formula is

$\mbox{fpr}+(1-\mbox{fpr}-\mbox{fnr})(1-\mbox{pnorm}(x, \mbox{tt}, \mbox{pxVar})$

where x is the stimulus value in dB, and pxVar is

$\min(\mbox{cap}, e^{A\times\mbox{tt}+B}).$

The ceiling cap is set with the call to opiInitialize, and A and B are from Table 1 in Henson et al (2000), also set in the call to opiInitialise using the type parameter.

Value

A list with elements:

• err NULL if no error, a string message otherwise.

• seen TRUE or FALSE.

• time Always NA.

Examples

# Stimulus is Size III white-on-white as in the HFA
chooseOpi("SimHenson")
res <- opiInitialize(type = "C", cap = 6)
if (!is.null(res$err))
  stop(paste("opiInitialize() failed:", res$err))

result <- opiPresent(stim = list(level = dbTocd(20)), tt = 30, fpr = 0.15, fnr = 0.01)
print(paste("Seen:", result$seen, quote = FALSE))

res <- opiClose()
if (!is.null(res$err))
  stop(paste("opiClose() failed:", res$err))

Description

Determine the response to a stimuli by sampling from a cumulative Gaussian Frequency-of-Seeing (FoS) curve (also known as the psychometric function).

For internal use only, use opiPresent().

Arguments

Details

As the response time returned for a false positive is determined separately from a positive response, we first check for a false response. If there is no false response, we use the FoS formula

$1-\mbox{pnorm}(x, \mbox{tt}, \mbox{pxVar})$

where x is the stimulus value in dB, and pxVar is

$\min(\mbox{cap}, e^{A\times\mbox{tt}+B}).$

The ceiling cap is set with the call to opiInitialize, and A and B are from Table 1 in Henson et al (2000), also set in the call to opiInitialise using the type parameter.

Thus, this function is the same as for SimHenson, but reaction times are determined by sampling from rtData as passed to opiInitialize. The dist parameter is the distance of the stimulus level from the true threshold, and should be in the same units as the Dist column of rtData. The default is just the straight difference between the stimulus level and the true threshold, but you might want it scaled somehow to match rtData.

Value

A list with elements:

• err NULL if no error, a string message otherwise.

• seen TRUE or FALSE.

• time The response time.

Examples

# Stimulus is Size III white-on-white as in the HFA
chooseOpi("SimHensonRt")
data(RtSigmaUnits)
res <- opiInitialize(type = "C", cap = 6, rtData = RtSigmaUnits)
if (!is.null(res$err))
  stop(paste("opiInitialize() failed:", res$err))

dist <- (10 - 30) / min(exp(-0.098 * 30 + 3.62), 6)
result <- opiPresent(stim = list(level = dbTocd(20)), tt = 30, fpr = 0.15, fnr = 0.01, dist=dist)
print(result, quote = FALSE)

res <- opiClose()
if (!is.null(res$err))
  stop(paste("opiClose() failed:", res$err))


if (!is.null(opiClose()))
  warning("opiClose() failed")

Description

Always respond 'not seen' to any parameter. No checking is done on the validity of stim.

Arguments

Value

A list with elements:

• err Always NULL.

• time Always NA.

Description

Always respond 'yes' immediately to any parameter.

Arguments

Value

A list with elements:

• err Always FALSE.

• seen Always TRUE.

• time Always NA.

Description

Returns a list that describes the current state of the machine. Specific parameters and return values can be seen in the machine specific versions listed below in the ’See Also’.

Usage

opiQueryDevice()

Value

A list specific to each machine.

See Also

opiQueryDevice_for_ImoVifa(), opiQueryDevice_for_Compass(), opiQueryDevice_for_Octopus900(), opiQueryDevice_for_PhoneHMD(), opiQueryDevice_for_Display(), opiQueryDevice_for_PicoVR(), opiQueryDevice_for_SimNo(), opiQueryDevice_for_SimYes(), opiQueryDevice_for_SimHenson(), opiQueryDevice_for_SimGaussian()

Description

This is for internal use only. Use opiQueryDevice() with these Arguments and you will get the Value back.

Value

A list containing constants and their values used in the OPI Compass module.

Description

This is for internal use only. Use opiQueryDevice() after chooseOPI("Display") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiQueryDevice()

Examples

chooseOpi("Display")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiQueryDevice()

Description

This is for internal use only. Use opiQueryDevice() after chooseOPI("ImoVifa") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

• rightEyex x co-ordinates of right pupil (pixels from image centre, left < 0)

• rightEyey y co-ordinates of right pupil (pixels from image centre, up > 0)

• leftEyex x co-ordinates of left pupil (pixels from image centre, left < 0)

• leftEyey y co-ordinates of left pupil (pixels from image centre, up > 0)

• leftEyed Diameter of left pupil (pixels)

• rightEyed Diameter of right pupil (pixels)

See Also

opiQueryDevice()

Examples

chooseOpi("ImoVifa")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiQueryDevice()

Description

This is for internal use only. Use opiQueryDevice() with these Arguments and you will get the Value back.

Details

Prints defined constants in OPI package pertaining to Kowa AP7000.

Value

List containing

• isSim = FALSE.

• pupilX, the x-coordinate of the pupil position in pixels.

• pupilY, the y-coordinate of the pupil position in pixels.

• purkinjeX, the x-coordinate of the purkinje position in pixels.

• purkinjeY, the y-coordinate of the purkinje position in pixels.

Description

This is for internal use only. Use opiQueryDevice() with these Arguments and you will get the Value back.

Details

Prints defined constants in OPI package pertaining to O600.

Value

Returns a list of 10 items:

* \code{answerButton} 0 = not pressed, 1 = pressed
* \code{headSensor} 0 = no forehead detected, 1 = forehead detected
* \code{eyeLidClosureLeft} 0 = eye is open, 1 = eye is closed
* \code{eyeLidClosureRight} 0 = eye is open, 1 = eye is closed
* \code{fixationLostLeft} 1 = eye pos lost, 0 = eye pos ok
* \code{fixationLostRight} 1 = eye pos lost, 0 = eye pos ok
* \code{pupilPositionXLeft} (in px)
* \code{pupilPositionYLeft} (in px)
* \code{pupilPositionXRight} (in px)
* \code{pupilPositionYRight} (in px)

Description

This is for internal use only. Use opiQueryDevice() with these Arguments and you will get the Value back.

Details

Prints defined constants in OPI package pertaining to Octopus 900.

Value

List containing isSim = FALSE.

Description

This is for internal use only. Use opiQueryDevice() after chooseOPI("PhoneHMD") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiQueryDevice()

Examples

chooseOpi("PhoneHMD")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiQueryDevice()

Description

This is for internal use only. Use opiQueryDevice() after chooseOPI("PicoVR") to call this function.

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiQueryDevice()

Examples

chooseOpi("PicoVR")
opiInitialise(list(port = 50001, ip = "localhost"))
opiSetup(list(eye = "BOTH"))
result <- opiQueryDevice()

Description

Returns a simple list.

Value

A list with elements:

• err Always NULL

• machine that is set to "SimGaussian".

Description

Returns name of the machine.

Value

A list with elements:

• isSim Always TRUE.

• machine that is set to "SimHenson".

Description

Returns name of the machine.

Value

A list with elements:

• isSim Always TRUE.

• machine that is set to "SimHensonRT".

Description

Returns name of the machine.

Value

A list with elements:

• machine that is set to "SimNo".

• isSim that is set to TRUE.

Description

Returns name of the machine.

Value

A list with elements:

• isSim Always TRUE.

• machine that is set to "SimYes".

Deprecated. Use opiSetup().

Description

In older OPIs it set background color and luminance in both eyes. Deprecated for OPI >= v3.0.0 and replaced with opiSetup().

See Also

opiSetup()

Description

Specific parameters and return values can be seen in the machine specific versions listed below in the ’See Also’.

Usage

opiSetup(settings)

Arguments

Value

Each implementation should(!) return a list with at least the following elements:

• err NULL if no error, otherwise a string describing the error.

See Also

opiSetup_for_Compass(), opiSetup_for_Octopus900(), opiSetup_for_ImoVifa(), opiSetup_for_PhoneHMD(), opiSetup_for_Display(), opiSetup_for_PicoVR(), opiSetup_for_SimNo(), opiSetup_for_SimYes(), opiSetup_for_SimHenson(), opiSetup_for_SimGaussian()

Description

This is for internal use only. Use opiSetup() with these Arguments and you will get the Value back.

Arguments

Details

Note: tracking will be relative to the PRL established with the fixation marker used at setup (call to OPI-OPEN), so when tracking is on you should use the same fixation location as in the setup.

Value

A list containing err which is NULL for success, or some string description for fail.

Description

This is for internal use only. Use opiSetup() after chooseOPI("Display") to call this function.

Arguments

Details

eye can take on values in the set {"left", "right", "both", "none"}.

fixShape can take on values in the set {"triangle", "square", "polygon", "hollow_triangle", "hollow_square", "hollow_polygon", "cross", "maltese", "circle", "annulus", "optotype", "text", "model"}.

fixLum can take on values in the range [0.0, 1.0E10].

fixType can take on values in the set {"flat", "checkerboard", "sine", "squaresine", "g1", "g2", "g3", "text", "image"}.

fixCx can take on values in the range [-90.0, 90.0].

fixCy can take on values in the range [-90.0, 90.0].

Elements in fixCol can take on values in the range [0.0, 1.0].

bgLum can take on values in the range [0.0, 1.0E10].

tracking can take on values in the range [0, 1].

Elements in bgCol can take on values in the range [0.0, 1.0].

fixSx can take on values in the range [0.0, 1.0E10].

fixSy can take on values in the range [0.0, 1.0E10].

fixRotation can take on values in the range [0.0, 360.0].

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiSetup()

Examples

chooseOpi("Display")
opiInitialise(list(port = 50001, ip = "localhost"))
result <- opiSetup(settings = list(eye = "BOTH"))

Description

This is for internal use only. Use opiSetup() after chooseOPI("ImoVifa") to call this function.

Arguments

Details

eye can take on values in the set {"left", "right", "both", "none"}.

fixShape can take on values in the set {"triangle", "square", "polygon", "hollow_triangle", "hollow_square", "hollow_polygon", "cross", "maltese", "circle", "annulus", "optotype", "text", "model"}.

fixLum can take on values in the range [0.0, 1.0E10].

fixType can take on values in the set {"flat", "checkerboard", "sine", "squaresine", "g1", "g2", "g3", "text", "image"}.

fixCx can take on values in the range [-90.0, 90.0].

fixCy can take on values in the range [-90.0, 90.0].

Elements in fixCol can take on values in the range [0.0, 1.0].

bgLum can take on values in the range [0.0, 1.0E10].

tracking can take on values in the range [0, 1].

Elements in bgCol can take on values in the range [0.0, 1.0].

fixSx can take on values in the range [0.0, 1.0E10].

fixSy can take on values in the range [0.0, 1.0E10].

fixRotation can take on values in the range [0.0, 360.0].

Value

A list containing:

• err NULL if there was no error, a string message if there is an error.

See Also

opiSetup()

Examples

chooseOpi("ImoVifa")
opiInitialise(list(port = 50001, ip = "localhost"))
result <- opiSetup(settings = list(eye = "BOTH"))

Description

Implementation of opiSetup for the Kowa AP7000 machine.

This is for internal use only. Use opiSetup() with these Arguments.

Usage

opiSetup_for_KowaAP7000(lum = NA, color = NA, fixation = NA)

Arguments

Details

If lum is 10 and color is not set, then .opi_env$KowaAP7000$BACKGROUND_WHITE is assumed.