Introduction {#s1}
============

One of the noted behavioral features of primates is their ability to use their hands to interact with objects in various situations. Motor outputs that control the hands are generated based on environmental information collected through sensory inputs. The central nervous system (CNS) computes a suitable transformation of these sensory inputs to motor outputs that allows motor performance required in daily life. One of these performances is the reach-to-grasp movement, in which the primate extends the arm toward an object placed in front of it, and then grasps the object with its fingers.

Jeannerod ([@B16], [@B17]) has investigated behavioral properties associated with these prehension movements. In his experiments, human participants were instructed to perform natural reach-to-grasp movements toward visually presented target objects of several sizes placed at several distances. He found some basic features: the velocity profile of the hand exhibits one peak in the first half of the movement duration, while the grip aperture demonstrates one peak in the second half. The peak grip aperture (PGA) is scaled to the size of the target object (e.g., Marteniuk et al., [@B22]). Appropriate transformation from the visual perception of the target object into the generation of the PGA is required for successful prehension movement.

How the CNS computes this transformation from sensory input of the visual object size to motor output in the form of a PGA is not yet clear. Previous studies have reported that the PGA is larger when vision is not available during the movement (e.g., Wing et al., [@B40]; Jakobson and Goodale, [@B15]; Fukui and Inui, [@B6]) or when the movement speed is particularly rapid (Wing et al., [@B40]). One interpretation of these observations is that the grip aperture is controlled to prevent inappropriate collisions of the fingers with the target and that the PGA becomes larger when visual uncertainty and/or motor variability are increased because of visual occlusion and/or faster movement (Wing et al., [@B40]). Actually, if the target object is presented in an eccentric view and its actual position is uncertain, the PGA will increase linearly with the eccentricity of the view (Schlicht and Schrater, [@B27]), which indicates that visual uncertainty is influential in the mechanism for grip aperture control. Elucidating the underlying mechanism for this grip aperture control would be possible with the use of appropriate modeling studies. Although previous models have described the reach-to-grasp movement (Hoff and Arbib, [@B11]; Haggard and Wing, [@B9]; Zaal et al., [@B42]; Smeets and Brenner, [@B29]; Ulloa and Bullock, [@B33]; Simmons and Demiris, [@B28]), no model has yet satisfactorily explained the effect of visual uncertainty on the generation of the PGA.

Traditionally, one of the major problems in modeling motor control is determining one trajectory of movement in the presence of motor redundancy. A prevailing idea is that the CNS solves this problem by minimizing costs to generate these movement trajectories, and several costs have been proposed (e.g., Flash and Hogan, [@B5]; Uno et al., [@B34]; Harris and Wolpert, [@B10]; Todorov and Jordan, [@B31]). One index of the cost minimizations that successfully explains human reaching trajectory is movement smoothness (Flash and Hogan, [@B5]; Uno et al., [@B34]). However, why the CNS adopts such costs is not clear. Harris and Wolpert ([@B10]) focused on the importance of the cost of the final variance, because it is clearly related to the achievement of the task. If the final body state is less variable, the movement is more likely to be successful. Motor control also needs to optimize the possibility of task achievement in the presence of motor variability (Miyamoto et al., [@B24]; Todorov, [@B30]; Trommershäuser et al., [@B32]). Thus, recent computational models of motor control have emphasized the importance of both sensory and motor variability (e.g., van Beers et al., [@B35], [@B36]; Guigon et al., [@B8]).

In a reach-to-grasp movement, the compensation for motor variability would result in adjustment of the PGA. A minimal variance model proposed by Simmons and Demiris ([@B28]) indeed explains the emergence of PGA in reach-to-grasp movement. Specifically, their model reported that a higher movement speed in a reach-to-grasp movement creates a larger motor variability, and the grip aperture increases in a compensatory manner to prevent undesirable collisions of the hand with the target. However, since visual occlusion of the target and/or of the hand would affect sensory uncertainty, not motor variability, the manner how sensory feedback contributes to variability-based motor control should be revealed.

To clarify this manner, Todorov and Jordan ([@B31]) proposed a computational model for motor control based on stochastic control theory. In their model, called optimal feedback control, a copy of the motor command and the noisy sensory feedback are used for internal estimation of the current state of the body. The optimal motor command that minimizes the task cost is then computed from this state estimation. The body state is estimated based on both motor and sensory variability, and the motor command is generated by taking into account the uncertainty of the state estimation. The concept underlying this computational model could explain the effect of online vision during reach-to-grasp movement, as it gives sensory feedback the role of compensating for motor variability and improving the estimation of the current body state. Specifically, when visual feedback is absent, the motor variability would directly cause uncertainty in the estimated body state due to the loss of the compensation by sensory feedback. Our hypothesis is that a larger PGA would appear as a result of compensating for the increased uncertainty of state estimation.

In the present paper, we propose a stochastic control model for grasping that can simulate the effects of online vision. Although, Todorov and Jordan ([@B31]) verified their own model by simulating arm reaching (pointing) movement, application of their model directly to reach-to-grasp movement is difficult because the movement is too complicated to be described by the task cost that is available in their model. Instead, we avoid optimal control and consider "compensation for sensory and motor variability" for calculating motor command, thereby preserving the concept of stochastic control.

An argument could be made that the effect of visual occlusion during grasping might be explained without invoking an online control mechanism that includes motor and sensory variability; i.e., some strategic effect against removing vision (e.g., the visual feedback schedule) might be causing a larger PGA (Jakobson and Goodale, [@B15]). However, Fukui and Inui ([@B6]) investigated the temporal and spatial effects of visual occlusion on the grip configuration when several occlusion conditions were randomly presented (see also Whitwell et al., [@B39]) and found that removing the target vision in the early phase of the movement enlarged the PGA, which indicated that the PGA was modified by the availability of online vision, not only by the motor strategy. Here, we verify that the effect of visual occlusion on the PGA is due to an online control mechanism that includes motor and sensory variability by simulating the experiments of Fukui and Inui ([@B6]) using the following model.

Materials and methods {#s2}
=====================

Overview of the aperture control model
--------------------------------------

Reach-to-grasp movement is traditionally thought to have two control components: a transportation (reaching) component and an aperture (manipulation) component (Jeannerod, [@B16]). Although the kinematics of trajectories of each finger are different, which may result from detailed control of each finger (Smeets and Brenner, [@B29]), the grip aperture size that represents the spatial relationship of these fingers shows significant correlation with the object size (Marteniuk et al., [@B22]).

Schlicht and Schrater ([@B27]) suggested, based on the principal components analysis, that only one principal component, which can be represented by the aperture size, would explain the difference in the trajectory of the fingers during reach-to-grasp movements where the target is placed in different visual eccentricities. Therefore, we assume that modeling grip aperture between thumb and index finger, and not the kinematics of each finger, is sufficient for revealing the grip control mechanism influenced by online vision (see also Vilaplana and Coronado, [@B37]).

An outline of the model is shown in Figure [1](#F1){ref-type="fig"}. The output motor command that controls grip aperture is transmitted to the end effector (hand) with delay and noise, and this changes the grip aperture. At the same time, the efference copy of the motor command is transmitted to the State Estimator, with no delay, and the body state (grip aperture) is estimated (or predicted) by the forward model. The estimated aperture is compared with the sensory feedback derived from vision and proprioception of the actual grip aperture, and the estimation error is then used to correct the body state estimation. The next motor command is generated based on the estimated (predicted) body state and the uncertainty of the estimation. The size and uncertainty of the target object that is observed by vision are then used to calculate the motor command. Further details are described below.

![**Outline of the model**. The output motor command that controls grip aperture is transmitted with delay and noise. At the same time, an efference copy of the motor command is transmitted to the State Estimator, and the body state (grip aperture) is estimated (or predicted) by the forward model. The sensory prediction is computed from the estimated aperture, and is compared with the sensory feedback (vision and proprioception) of the actual grip aperture. The Kalman filter corrects the estimated body state based on the sensory prediction error. The next motor command is generated based on the estimated (predicted) body state and the variability (variance) of estimation. The size and variability (variance) of the target object that are observed by vision with noise and delay are then used to calculate the motor command. Although the transport component was assumed to contribute to the aperture control, the result of parameter fitting showed that it was not effective.](fncom-09-00143-g0001){#F1}

Aperture size to compensate for sensorimotor variability
--------------------------------------------------------

In the model setting, the hand travels from the start position to the object position along the y-axis, as indicated in Figure [2A](#F2){ref-type="fig"}. This hand transportation is not modeled here. The transportation component (the profile of hand position along the y-axis) is realized by the data from the experiment by Fukui and Inui ([@B6]). The finger aperture opens along the x-axis because the target object is a simple cylinder and the grip orientation has little effect on the aperture size. When the transportation component properly controls the hand position so that the center point between the thumb and finger matches the center of the target object, the problem becomes how the aperture size is controlled as the hand approaches the object.

![**Hand and target object. (A)** Hand transportation follows the *y*-axis, and grip opening follows the *x*-axis. **(B)** For successful grasping, the grip aperture *a* must be sufficiently larger than object size *s*~*o*~, considering the variability of grip aperture $\sigma_{\textit{a}}^{\text{2}}$, hand position $\sigma_{\textit{h}}^{\text{2}}$, and target object $\sigma_{\textit{o}}^{\text{2}}$.](fncom-09-00143-g0002){#F2}

To avoid an undesirable collision during a reach-to-grasp movement, the grip aperture has to include a margin of safety that accounts for uncertainty of the hand and the object. For example, consider the situation shown in Figure [2B](#F2){ref-type="fig"} to represent the beginning of the finger closure phase of the movement. If the position on the x-axis and the diameter of the object are denoted by *x*~*o*~, *s*~*o*~ respectively, the position of the right edge of the object is *x*~*o*~ + *s*~*o*~∕2. Similarly, the position of the finger at the right side of the object (i.e., the index finger) is *x*~*h*~ + *a*∕2, if the hand position (the center point between the thumb and the finger) and the aperture size are denoted by *x*~*h*~ and *a*, respectively. To avoid collision of the finger with the object, the finger position has to be to the right of the object edge, i.e., *x*~*o*~ + *s*~*o*~∕2 \< *x*~*h*~ + *a*∕2. We assume that *x*~*o*~ − *x*~*h*~ = 0 because the hand position should be controlled to match the object position by the transportation component. However, the motor variability and sensory uncertainty of hand position relative to the object position still exists, and here, the standard deviation (SD) of the variability (uncertainty) is denoted by σ~*h*~. The variability (uncertainty) of the aperture size and the object size (the SDs of them are denoted by σ~*a*~ and σ~*o*~) also take effect. The situation before successful grasping has to satisfy *a* − *s*~*o*~ \>0 under the condition where the variability of *a* − *s*~*o*~ is characterized by the variance $\sigma^{2} = \sigma_{a}^{2} + \sigma_{h}^{2} + \sigma_{o}^{2}$, if Gaussian distribution and the independence of these variabilities is assumed. The same situation takes place at the left side of the object. Here, the probability Φ (*a, s*~*o*~, σ) that the target object is inside the grip aperture is formulated by:

Φ

(

a

,

s

o

,

σ

)

=

1

2

π

σ

2

∫

0

∞

exp

\-

(

a

′

\-

(

a

\-

s

o

)

)

2

2

σ

2

d

a

′

This is the cumulative distribution function that describes the probability that the grip aperture *a* is larger than target object *s*~*o*~ in the presence of Gaussian noise.

In order to execute a grasping movement when uncertainty is involved, the CNS has to take the risk of collision into account. If we consider a particular probability ϕ, which the CNS adopts as a success ratio, the grip aperture to achieve is then calculated by $a^{*} = \Phi^{- 1}\left( {\phi,s_{o},\sigma} \right)$. Note that ϕ is the probability that the CNS predicts before executing the movement with the uncertainty and variability (i.e., σ^2^) that the CNS estimates, and that ϕ does not represent the actual success rate of the grasping. Here, we call the aperture *a*^\*^ the *target aperture*. Once the target aperture is determined, the controller generates a motor command that moves the current grip aperture toward the target aperture.

Among infinite possibilities of motor commands that make the current grip aperture approach the target aperture, one motor command should be determined to generate an actual movement. In human arm control, the movement trajectory is determined by minimizing a certain criterion (Flash and Hogan, [@B5]; Uno et al., [@B34]; Harris and Wolpert, [@B10]). Smoothness of movement is an important factor in determining the trajectory; the trajectory of point-to-point movement starting and ending with zero velocity and zero acceleration is well explained by the minimum jerk model (Flash and Hogan, [@B5]). However, the start point of the trajectory formulation required in the present model is not the movement start point but the *current* point during movement, and the initial velocity and acceleration of the trajectory are not zero. Minimum acceleration, which is another criterion for kinematic smoothness of movement, also well explains the reaching trajectory if boundary conditions (i.e., the initial and final velocity and acceleration) are specified (Ben-Itzhak and Karniel, [@B1]). Here, we adopted the minimum acceleration criterion for generating a motor command (Note that the model only generates the *next* motor command from the *current* grip aperture, not the whole trajectory).

To achieve successful grasping, the thumb and fingers have to attain a suitable contact with the surface of the object. At the same time, they have to avoid undesirable collision with the object. The relative timing of these requirements is different: first, for avoiding collision, and then, for contacting the object. Winges et al. ([@B41]) investigated finger formation for different shaped objects during reach-to-grasp movement under several visual conditions, and they concluded that visual occlusion prolonged the final low-speed phase of reaching when the precise finger formation occurred. Our model is focused on the mechanism for generating peak grip aperture, not for contacting the object; therefore, we did not implement the mechanism for closing fingers upon suitably touching the object surface. Details are mentioned in the section that covers the simulation (Section Simulation).

Formulation of aperture size control
------------------------------------

The movement of grip aperture is modeled by a state-space representation:

x

t

\+

Δ

t

=

F

x

t

\+

D

u

t

\+

G

w

t

F

=

\[

1

Δ

t

0

1

\]

D

=

\[

0

Δ

t

\]

G

=

\[

0

1

\]

Here, $\textit{x}_{t} = \left\lbrack {a_{t},ȧ_{t}} \right\rbrack^{T}$, where *a*~*t*~ and ȧ~*t*~ is grip aperture and its velocity at time *t*, respectively, and \[□\]^*T*^ denotes transposition of a vector or matrix. Motor command *u*~*t*~ is acceleration of grip aperture and *w*~*t*~ is the motor noise generated from *N*(0, *Q*~*t*~), while Δ*t* is a discrete time step in simulation. The controller generates motor commands to achieve the target aperture *a*^\*^ at the movement end time *T*. As a requirement in the experiment, the movement duration is predefined to be 1 second (Fukui and Inui, [@B6]). In order to simulate a smooth movement that achieves target aperture *a*^\*^ and target aperture velocity ȧ^\*^ = 0 at the movement end, we adopted the minimum acceleration criterion. The velocity at time τ (*t* ≤ τ ≤ *T*) in the minimum acceleration profile is described by:

ȧ

τ

=

ȧ

t

\+

b

1

(

τ

\-

t

)

\+

b

2

(

τ

\-

t

)

2

b

1

=

6

(

a

\*

\-

a

t

)

(

T

\-

t

)

2

\-

4

ȧ

t

T

\-

t

b

2

=

\-

6

(

a

\*

\-

a

t

)

(

T

\-

t

)

3

\-

3

ȧ

t

(

T

\-

t

)

2

(See Appendix [A](#A1){ref-type="app"}). What the controller determines here is not the whole trajectory, but the instant acceleration at time *t*. Based on Equation (3), the aperture velocity at time *t* + Δ*t* is:

ȧ

t

\+

Δ

t

=

ȧ

t

\+

b

1

Δ

t

\+

b

2

Δ

t

2

Since ȧ~*t*\ +\ Δ*t*~ = ȧ~*t*~ + *u*~*t*~Δ*t* comes from Equation (2), the motor command at time *t* in the minimum acceleration criterion is:

u

t

=

b

1

\+

b

2

Δ

t

Although, the motor command is calculated from the current and target aperture, the CNS does not know the true state of the current body. The information about body state has to be obtained by the sensory system. Here, grip aperture is observed by vision and proprioception:

y

t

=

H

x

t

\+

v

t

H

=

\[

1

0

1

0

\]

For simplicity, we assumed that only the grip aperture (not the aperture velocity) is used as input information. Here, the first element of ***y***~*t*~ is for vision and the second one is for proprioception, where ***v***~*t*~ is observation noise generated from *N*(0, *R*~*t*~). The Kalman filter (Appendix [B](#A2){ref-type="app"}) estimates the body state from sensory observation. If the estimated body state is ${\overset{\hat{}}{\textit{x}}}_{t} = \left\lbrack {â_{t},{\overset{\hat{}}{ȧ}}_{t}} \right\rbrack^{T}$ then *a*~*t*~ and ȧ~*t*~ in Equations (4) and (5) are replaced by â~*t*~ and ${\overset{\hat{}}{ȧ}}_{t}$, respectively:

b

\^

1

=

6

(

a

\*

\-

â

t

)

(

T

\-

t

)

2

\-

4

ȧ

\^

t

T

\-

t

b

\^

2

=

\-

6

(

a

\*

\-

â

t

)

(

T

\-

t

)

3

\-

3

ȧ

\^

t

(

T

\-

t

)

2

These equations generate the motor command based on the estimated body state.

The actual brain-body system undergoes transmission delays in the motor command and the sensory feedback pathway. When motor delay and sensory delay are *d*~*m*~ and *d*~*s*~, respectively, Equations (2) and (8) are modified to:

x

t

\+

Δ

t

=

F

x

t

\+

D

u

t

\-

d

m

\+

G

w

t

y

t

=

H

x

t

\-

d

s

\+

v

t

The motor command generated at time *t* has to be based on the body state (estimation) at time *t* + *d*~*m*~. Here, Equations (9) and (10) are modified to:

b

\^

1

=

6

(

a

\*

\-

â

t

\+

d

m

)

(

T

\-

(

t

\+

d

m

)

)

2

\-

4

ȧ

\^

t

\+

d

m

T

\-

(

t

\+

d

m

)

b

\^

2

=

\-

6

(

a

\*

\-

â

t

\+

d

m

)

(

T

\-

(

t

\+

d

m

)

)

3

\-

3

ȧ

\^

t

\+

d

m

(

T

\-

(

t

\+

d

m

)

)

2

These equations imply that the "current" body representation in the brain precedes the actual body state. This function is necessary for predictive remapping (Colby and Duhamel, [@B2]; Melcher, [@B23]), where neural representations of planned saccade movements are modified to compensate for future changes of eye position. Similarly, the grip aperture also has to be predicted before the motor command reaches the effector during reach-to-grasp control.

Prediction of uncertainty
-------------------------

As mentioned in Section Aperture Size to Compensate for Sensorimotor Variability, the prediction error variance at the end of the movement is necessary for determining the target aperture. The variance of grip aperture $\sigma_{a}^{2}$ is obtained by prediction of the Kalman filter (Kalman, [@B18]; see Appendix [B](#A2){ref-type="app"}). The Kalman filter predicts the future body state following the state-space representation, Equation (11), as well as the variance of prediction error. Although, in the Kalman filter, the variance is typically used to determine filter gain for correction of an estimated body state based on sensory feedback, the proposed model also uses the calculated estimation variability for aperture control.

The prediction of the SD of hand position σ~*h*~ is assumed to be proportional to the distance of residual hand transportation. That is:

σ

h

=

α

(

h

\*

\-

h

t

)

,

where *h*^\*^ is the distance between the initial hand position and the target object, *h*~*t*~ is the distance between the initial hand position and the hand position at time *t*, and α is the proportionality coefficient. This assumption is based on the idea that increases in distance to the target position give rise to more uncertainty in the prediction of final hand position. This agrees with the fact that travel distance and the target width are proportional in Fitts\' law (Fitts, [@B4]). If target width is interpreted as "tolerable size of variability of the final hand position," then Equation (15) will be derived. As the hand approaches the target position, the prediction of the SD decreases because the residual distance to be traveled decreases.

In a reach-to-grasp movement, the effect of the aperture component on the transportation component is small (Paulignan et al., [@B26]), especially before the movement end point. Consequently, the use of a fixed transportation component is reasonable for the simulation of the generation of peak grip aperture. Here, we used the wrist position profile (Figure [3A](#F3){ref-type="fig"}) that is obtained by temporally normalizing and averaging the experimental data of Fukui and Inui ([@B6]).

![**Profiles of variables in the model**. Variables obtained from a simulation with target object of 4 cm in Full Vision (FV) and 150S (visual occlusion from 150 ms after movement onset) conditions are shown. Black solid and dashed lines denote FV and 150S conditions, respectively. Gray vertical lines denote onset of visual occlusion in 150S condition. **(A)** Hand distance traveled from the start position. **(B)** Predicted SD of grip aperture σ~*a*~. **(C)** Target aperture *a*^\*^. **(D)** Grip aperture. **(E)** SD of object size σ~*o*~.](fncom-09-00143-g0003){#F3}

Finally, because the target object is only observed visually, the variability is determined by visual noise. The variance of uncertainty of the target object size $\sigma_{o}^{2}$ is obtained from the vision element of the covariance matrix of sensory noise *R*~*t*~ (see Section Simulation).

Simulation
----------

Using the proposed model, we simulated a visual occlusion experiment by Fukui and Inui ([@B6]). In the simulation, movement end time *T* and a discrete time step Δ*t* were set to 1000 and 10 ms, respectively. Both visual and motor delays were set to 100 ms (Huettel et al., [@B12]). Because the movement onset in the experiment was defined by wrist release from a start position, the grip aperture and the aperture velocity at the movement onset were not zero. Therefore, in the simulation we set the initial grip aperture to be an average grip aperture at the movement onset in all trials recorded in the experiment. Additionally, the initial motor command (acceleration) in the simulation was set to be a value that achieved the average velocity at the movement onset in the experiment. The initial motor command in the simulation was generated taking into account the *d*~*m*~ (motor delay) period preceding movement onset.

For the simulation of visual occlusion during movement, the variance of visual feedback noise was assumed to be larger than usual when vision was unavailable. On the other hand, proprioceptive noise was assumed to be constant. If the variance of visual noise during occlusion (no vision) and during non-occlusion (full vision) are denoted by $\sigma_{NV}^{2}$ and $\sigma_{FV}^{2}$, respectively, and the variance of proprioceptive noise is denoted by $\sigma_{P}^{2}$, the covariance matrix of sensory noise is:

R

t

=

\[

σ

N

V

2

0

0

σ

P

2

\]

(

during occlusion

)

R

t

=

\[

σ

F

V

2

0

0

σ

P

2

\]

(

during non-occlusion

)

No correlation is assumed between visual noise and proprioceptive noise. Viewed from an engineering perspective, visual occlusion should have been simulated by simply setting the Kalman gain to zero. This implementation means that the system estimates the target state only from Equation (2), but we did not adopt this implementation in the current study. Westwood and Goodale ([@B38]) have demonstrated that the grip aperture during visually occluded grasping was affected by a size-contrast illusion, while the grip aperture in normal vision was not. We assume that this may reflect switching of the visual processing pathways from the dorsal "how" stream to the ventral "what" stream (Goodale and Milner, [@B7]). Therefore, the two kinds of noise assumed in the present model Equation (16) represent these two different types of visual processing.

As described in Equation (2), motor noise is assumed to be a single dimension and affects the system in the same domain as the motor command. This implies that the motor noise originates from motor command. However, the variance of motor noise is assumed to be constant; that is, $Q_{t} = \sigma_{M}^{2}$.

Here, the proposed model has six parameters: the variance of motor noise $\sigma_{M}^{2}$, the variance of visual noise during occlusion and non-occlusion $\sigma_{NV}^{2}$, $\sigma_{FV}^{2}$, the variance of proprioceptive noise $\sigma_{P}^{2}$, the proportionality coefficients α between distance and the variance of hand position, and the probability ϕ that the object has to be inside the aperture. For successful grasping, the probability ϕ should be high. However, if ϕ is too high, the noise variances have to be very small in order to achieve the appropriate grip aperture, and small variances make it difficult to differentiate the visual conditions. Therefore, ϕ was fixed a priori at 0.9 based on the empirical inspection that people could fail prehension performance once out of 10 trials. The other parameters were determined by fitting the model outputs to the experimental data.

A probability ϕ of 0.9 might be rather small for successful grasping. However, the proposed model was not designed to explain precise control of the final phase of the movement, where the successful achievement of grasping is realized by the mechanism for adjusting the aperture size to the object size. Rather, this model aims to demonstrate how the generation of the PGA is influenced by online vision during the movement. We assume that the probability ϕ, which determines the PGA, does not have to be strictly high during the movement. In order to examine the effect of ϕ on the aperture profile, the simulation with various values for ϕ (= \[0.8, 0.85, 0.9, 0.95, 0.99\]) was also conducted.

Parameter fitting
-----------------

Five parameters, σ~*M*~, σ~*P*~, σ~*FV*~, σ~*NV*~ and α, were fixed by fitting the grip aperture profile obtained from the simulation to that in the experiment by Fukui and Inui ([@B6]). In the experiment, a visual occlusion experiment was conducted while the reach-to-grasp movement was in process with two experimental paradigms: a shutting paradigm (SP) and a re-opening paradigm (RP). In the SP, vision was available during target presentation and reaction time phases, and occlusion started immediately after movement onset (NV), or 150 ms (150S), 350 ms (350S), 500 ms (500S), or 700 ms (700S) after movement onset; or vision is never occluded during movement (FV). On the other hand, in RP, occlusion once started at movement onset, and vision was recovered immediately after movement onset (FV), or 150 ms (150R), 350 ms (350R), 500 ms (500R), or 700 ms (700R) after movement onset; or vision was never recovered (NV). Both paradigms had conditions in which vision during the movement was constantly available (FV) or constantly unavailable (NV). Model parameters were fitted against the aperture profile from these constant conditions. For the simulation and comparison of simulated aperture profiles to the ones obtained from the experiment, the movement time of each trial of experimental data was normalized to 1 s, and the grip aperture and wrist position were resampled at 100 Hz by interpolating the experimental data. Note that the participants in the experiment were instructed to keep the movement duration of the trials at approximately 1000 ms. Then, the temporally normalized aperture profiles of all subjects and trials were averaged for each condition.

Fitting was performed by minimizing the squared sum of errors at each time point between the aperture profile of the experiment and the simulations for two constant vision conditions (i.e., FV and NV conditions). Minimization was achieved by the Nelder-Mead simplex method (function *fminsearch* of MATLAB optimization toolbox). The fitting result is shown in Table [1](#T1){ref-type="table"}.

###### 

**Results of parameter fitting**.

  **σ~***M***~ (cm/ms**^2^**)**   **σ~*P*~ (cm)**   **σ~*FV*~ (cm)**   **σ~*NV*~ (cm)**   **α**
  ------------------------------- ----------------- ------------------ ------------------ --------
  0.3824                          0.6043            0.0450             0.3187             0.0000

*Note that the degrees of variability are standard deviation, not variance*.

Results {#s3}
=======

As a result of fitting, parameter α was converged to zero. This means that we do not have to assume the relationship between the transport and aperture component in the variability domain. In other words, the aperture component considers the graspability of the target object only with respect to the variability (uncertainty) of the grip aperture and the object size (i.e., σ~*a*~ and σ~*o*~). This result is in line with the visuomotor channel hypothesis (Jeannerod, [@B16], [@B17], discussed below). The grip aperture profiles of the experiment and the simulation in the FV and NV conditions after parameter fitting are shown in Figure [4](#F4){ref-type="fig"}. The whole aperture profile is well-fitted to the experimental data.

![**Grip aperture for parameter fitting**. Profile of grip aperture from the experiment (Fukui and Inui, [@B6]; re-opening paradigm: blue lines; shutting paradigm: red lines) and the simulation (black) for 4 cm (top panels) and 6 cm (bottom panels) target objects in Full Vision (FV; left panels) and No Vision (NV; right panels) conditions.](fncom-09-00143-g0004){#F4}

The temporal development of the variables in the simulation of grasping toward a 4 cm target object in FV and 150S conditions is shown in Figure [3](#F3){ref-type="fig"}. As the hand approached the target object (Figure [3A](#F3){ref-type="fig"}), the predicted SD of the grip aperture decreased as the movement continued (Figure [3B](#F3){ref-type="fig"}) because the time for future prediction (i.e., the time left until the movement ends) decreased. Consequently, the target aperture decreased (Figure [3C](#F3){ref-type="fig"}) because the uncertainty of the final posture and the "safety margin" decreased. The profile of the grip aperture (Figure [3D](#F3){ref-type="fig"}) was generated by smooth following of the target aperture, exhibiting a peak at around 600 ms after movement onset. After the vision was occluded (Figure [3](#F3){ref-type="fig"}; the gray vertical line denotes 150-ms shutting), the visual noise increased (Figure [3E](#F3){ref-type="fig"}, dashed line) compared to the values obtained in FV conditions after sensory delay (i.e., 100 ms), and the predicted uncertainty of the final grip aperture (Figure [3B](#F3){ref-type="fig"}, dashed line), the target aperture (Figure [3C](#F3){ref-type="fig"}, dashed line) and the grip aperture (Figure [3D](#F3){ref-type="fig"}, dashed line) increased compared to the values obtained in FV condition.

The output aperture profiles of the model\'s simulation in all visual conditions, with the parameters described in Table [1](#T1){ref-type="table"}, are shown in Figure [5](#F5){ref-type="fig"}. For all conditions in both paradigms, the peak grip aperture appeared around 550--650 ms, as it did in the experimental results. Despite the lack of a mechanism for closing fingers upon touching the object, the model first opened the fingers larger than the object size and then closed them. This is because the predicted variability at the final phase of the movement decreased as the movement progressed. This decrease in the predicted variability arose primarily from the decrease in the duration to the end of the movement. Since, the Kalman filter calculates future aperture variability by repeatedly applying the incremental forward model (Appendix [B](#A2){ref-type="app"}), the later prediction points (i.e., prediction in earlier phase of the movement) will be more variable.

![**Grip aperture profiles for simulation of visual occlusion**. **Upper left**: target object size is 4 cm in the shutting paradigm (SP), **Upper right**: 4 cm in the re-opening paradigm (RP), **Lower left**: 6 cm in SP, and **Lower right**: 4 cm in RP. In FV conditions, vision was constantly available during the movement, while in NV conditions, vision was constantly unavailable. In 150S, 350S, 500S, and 700S conditions, visual occlusion started 150, 350, 500, and 700 ms after movement onset, respectively. In 150R, 350R, 500R, and 700R conditions, vision was first unavailable and then recovered 150, 350, 500, and 700 ms after movement onset, respectively. For each panel, aperture profiles within gray rounded rectangles are magnified for precise view.](fncom-09-00143-g0005){#F5}

Peak grip apertures in both the simulation and the experiment are shown in Figure [6](#F6){ref-type="fig"}. In the simulation, the peak grip apertures in 350S, 500S, 700S, and FV (shutting paradigm) conditions were remarkably similar value, as well as in 350R, 500R, 700R, NV (re-opening paradigm). The peak grip apertures observed in the experiment in these conditions did not show statistical difference, indicating that visual availability after 350 ms does not affect the peak grip aperture. The simulation replicated the experimental result that showed that vision in the early phase of the movement is important for grasping. Note that the parameters in the model were determined by fitting to the aperture profile of the NV and FV conditions. Other conditions were not considered in the fitting. Therefore, no temporal information about the change in visual availability was contained in the data used for fitting. Even then, the simulation still reproduced the effect of early phase visual availability on the PGA. This occurred because the predicted variability mostly decreased in 350 ms in the model, and consequently, the visual conditions became insensitive to the aperture control.

![**Peak grip aperture in the simulation and the experiment**. **Upper left**: experimental result in shutting paradigm, **Upper right**: experimental result re-opening paradigm, **Lower left**: simulation result in shutting paradigm, and **Lower right**: simulation result in re-opening paradigm. Experimental data is from Fukui and Inui ([@B6]). Error bars indicate the standard errors of the values between participants.](fncom-09-00143-g0006){#F6}

In order to verify the effect of online vision on the aperture profile, we investigated whether the simulated aperture velocity profiles between 150S (150R) and 350S (350R) conditions diverge as shown in the experiment of Fukui and Inui ([@B6]). In the experiment, the difference between the aperture velocities for the 150S (150R) and 350S (350R) conditions was statistically significant after 400--475 ms following movement onset. The aperture velocity profile of the model\'s output in 150S, 150R, 350S, and 350R conditions for 4 and 6 cm target objects is shown in Figure [7](#F7){ref-type="fig"}. The time points when the aperture velocities diverged between 150S (150R) and 350S (350R) conditions were around 400--450 ms. Specifically, the motor response to the onset (recover) of visual occlusion at 150 ms in 150S (150R) conditions was expected to start at 350 (i.e., 150 + 200) ms after movement onset. Considering that the response becomes noticeable only after the divergence is large enough, the result that the divergence time is similar in the simulation and experiment (i.e., 400--450 ms) indicates that the sensorimotor delay used in the simulation (=200 ms) is reasonable in this prehension task.

![**Simulated velocity profile of grip aperture for conditions of 150S, 350S, 150R, and 350R**. Divergence between 150S (150R) and 350S (350R) appeared at approximately 400--450 ms after movement onset, as found in the experiment.](fncom-09-00143-g0007){#F7}

To verify the robustness of the model to the variability of sensorimotor delay, we performed the simulations with different sensorimotor delays (100 and 300 ms). We found similar effects of (non-)available duration of online vision on the PGA. Specifically, the PGA in 150S and 350S conditions increased with less sensorimotor delay (100 ms), and the PGA in 150R and 350R conditions increased with more sensorimotor delay (300 ms; Figure [8](#F8){ref-type="fig"}). The time of divergence in aperture velocity was dependent on sensorimotor delay (Figure [9](#F9){ref-type="fig"}). Among these three conditions, the simulation in the 200-ms delay condition was best fitted to the experimental result in terms of temporal divergence, and this 200-ms delay would be consistent with neural response latency suggested by Huettel et al. ([@B12]).

![**Peak grip aperture in simulations with various sensorimotor delays**. **Top**: 100 ms sensorimotor delay. **Middle**: 200 ms sensorimotor delay. **Bottom**: 300 ms sensorimotor delay. **Left**: shutting paradigm. **Right**: re-opening paradigm. Note that the figures in the middle panels are identical to those shown in Figure [6](#F6){ref-type="fig"}.](fncom-09-00143-g0008){#F8}

![**Simulated velocity profile of the grip aperture for conditions of 150S, 350S, 150R and 350R, with 100, 200 and 300 ms sensorimotor delays**. The target object is 4 cm. **Top**: 100 ms sensorimotor delay. **Middle**: 200 ms sensorimotor delay. **Bottom**: 300 ms sensorimotor delay. Note that the figure in the middle panel is identical to that shown in Figure [7](#F7){ref-type="fig"}.](fncom-09-00143-g0009){#F9}

These additional simulations indicated that this emerging effect of online vision on the PGA was not due to sensorimotor delay but to the control mechanisms of grip aperture with the predicted variability assumed in the present model.

The grasping probability parameter ϕ had a strong effect on the peak grip aperture. According to the simulation results in which various values were set to ϕ (= \[0.8, 0.85, 0.9, 0.95, 0.99\]) without changing other parameters, the grip aperture was larger when the probability requirement was larger (Figure [10](#F10){ref-type="fig"}). This is consistent with the intuition that a larger safety margin is preserved when humans grasp things more carefully.

![**Aperture profiles for various grasping probability ϕ**. Simulation result with ϕ = \[0.8, 0.85, 0.9, 0.95, 0.99\]. **Left panel**, 4 cm; **Right panel**, 6 cm.](fncom-09-00143-g0010){#F10}

Furthermore, we performed the simulations with minimum jerk instead of minimum acceleration because the minimum jerk was generally used as a criterion of movement smoothness in previous works. We found a slower decrease of the grip aperture in the final phase of the movement and an earlier time to PGA than observed in the experiment.

Discussion {#s4}
==========

Sensory and motor variability in models for reach-to-grasp movement
-------------------------------------------------------------------

We proposed a model for reach-to-grasp movement that controlled grip aperture to compensate for predicted sensory uncertainty and motor variability. This is the first model that includes a control principle to explain how the grip aperture is arranged to grasp an object based on sensory uncertainty and motor variability. The formation of grip aperture was modeled in several previous studies. Hoff and Arbib ([@B11]) divided the aperture formation into preshaping and finger closure phases. In the preshaping phase, aperture is opened toward the peak grip aperture that scaled with object size. They provide no explanation for how the peak grip aperture is actually determined. Smeets and Brenner ([@B29]) modeled grasping as pointing of two fingers to the opposite side of the object, with the trajectory formed by minimum jerk principle. The overshoot of the grip aperture to the object size was generated by the constraint that the fingers had to approach the object orthogonally to the surface of the object. No explanation was given for determination of this "approaching parameter." Ulloa and Bullock ([@B33]) proposed a dynamic system model that included an equation to describe the influence of the transportation component on the aperture component. These models did not capture the function of sensory and motor variability in reach-to-grasp movement.

The fitting parameter α in the present model resulted in zero. This means that the prediction of wrist position variability is not necessary to determine the grip aperture in reach-to-grasp movements. However, it does not mean that there is no connection between the transport and aperture components; instead, a relationship between them in the time domain is implemented in the present model. Specifically, the grip configuration progresses by predicting future state variability at the movement end, which was pre-determined following the experimental instruction of movement time (i.e., 1000 ms). The movement time in the experiment was defined by the duration from the time when the wrist was released from the table to the time when the wrist started to move the target object approximately 30 cm towards the body. In other words, the grip configuration was partially regulated by the movement time, which is defined by the transport component. This type of temporal relationship between the two components has been modeled previously (Hoff and Arbib, [@B11]). Despite the assumption that the connection of these components was taken into account in the variability domain in the present model (i.e., the function of parameter α), the parameter that represents this connection resulted in zero, suggesting the visuomotor channel hypothesis (Jeannerod, [@B16], [@B17]), which assumes an independence of each component in visuomotor transformation. Although the visuomotor hypothesis does not incorporate the framework of variability-based motor control, the present model suggests that this independence of each component could be applicable to the level of variability-based motor control. The assumption of the temporal coordination of these visuomotor channels (Hoff and Arbib, [@B11]) is also implemented in the present model.

One possibility for the result that α became zero after fitting could be the inappropriate assumption of the proportionality between the residual distance and the predicted variability of the hand position (Equation 15). However, the fact that the grip aperture profile was successfully reproduced without assuming the predicted variability of the hand position indicates that the variability-based model of reach-to-grasp movement without assuming a transport component has the capacity to parsimoniously explain the effect of online vision on the PGA shown by Fukui and Inui ([@B6]). Nevertheless, incorporating a predicted variability of the hand position (which might not be describable by a simple function of the residual distance) into the model will contribute to a better understanding of reach-to-grasping movements. Further investigations are required to identify the relationship between the aperture and transport components in the variability domain.

As mentioned in the introduction, Simmons and Demiris ([@B28]) have shown the importance of variability in the control of grip aperture. They modeled aperture control as a movement that goes through via-points for each finger. They found that when motor variability was dependent on the amplitude of a motor command, the via-points that resulted in the smallest variance at the end point of the movement were chosen. However, their model only focused on motor variability. To explain the effect of visual feedback in grasping, the function of sensory uncertainty also has to be modeled. Our model includes the Kalman filter for predicting future variability of the aperture, so that sensory uncertainty also plays a role in the control.

The present model suggests a new computational view expanding the conventional Bayesian sensorimotor framework: the online prediction of motor variability at the task end is used for feedback control. Previous studies have demonstrated that motor variability and sensory uncertainty are integrated on the basis of the Bayesian framework (or the Kalman filter) in order to enhance the task performance in the presence of sensory and motor noise during motor planning (Harris and Wolpert, [@B10]; van Beers et al., [@B36]; Trommershäuser et al., [@B32]), sensorimotor learning (Körding and Wolpert, [@B20]) and online control (Izawa and Shadmehr, [@B14]). In the present model, the target aperture is determined using the predicted variability of grip aperture at the movement end. This means that the online feedback controller changes the "immediate goal (target)" necessary for achieving the task depending on both the Bayesian estimation of the current state of the body and the environment and the Bayesian prediction of a future state. This feature of the present model explains how the effect of online vision during prehension varies according to the phase of the movement. In the early phase of the movement (0--350 ms), where the prediction of a far future (\> 650 ms) is required, the predicted task end state is uncertain and the sensory feedback is important to reduce this uncertainty. We consider that this online variability prediction plays an important role in performing a grasping movement.

In the present model, motor noise is assumed to be a constant during the movement rather than signal-dependent noise (Harris and Wolpert, [@B10]). The controller used the predicted variability to determine a target aperture (*a*^\*^), and it was impossible to predict the future variability of the grip aperture with signal-dependent noise, because the noise affects the target aperture, the motor commands used to achieve this target aperture, and, consequently, the control noise itself. A cost function is necessary to fix the motor command in the presence of signal-dependent noise. The parameter ϕ, which was set to 0.9 in the present study, might be determined with such a cost function and signal-dependent noise in future studies.

Neural computation of variability
---------------------------------

We calculated the future variability of the aperture and updated the variability due to sensory feedback by implementing the Kalman filter in our model. During the estimation of the visual image, forward and inverse calculations similar to the Kalman filter appear to be conducted during the visual process (Kawato et al., [@B19]). In a reaching movement, the states of the body and environment are estimated in a manner similar to the Kalman filter, in off-line adaptation (Izawa et al., [@B13]) as well as in online processing (Izawa and Shadmehr, [@B14]). The Bayesian inference necessary for the Kalman filter can be computed by the neural population that fires in a probabilistic manner (Ma et al., [@B21]), and a neural implementation of the Kalman filter has been suggested (Denève et al., [@B3]). However, the way that the Kalman gain is computed still remains to be resolved.

The neural substrates related to the calculation necessary for the Kalman filter have been speculated (Ogawa et al., [@B25]). During performance of a tracing or tracking movement using a computer mouse with an artificial delay introduced occasionally, the right posterior parietal cortex was activated in relation to the function that detects the error between the target and the mouse cursor. The right temporo-parietal junction was involved in state estimation of self-movement during visually guided movement. Unfortunately, these studies on neural computation related to the Kalman filter only investigated reaching movements, not grasping. However, motor control of the body is clearly related to the sensory and motor variability, and a reach-to-grasp movement also has to be explained in this manner. The present study demonstrated that a computational model based on mechanisms of variability prediction with the Kalman filter explains the online effect of vision during reach-to-grasp movement.

Conclusions {#s5}
===========

We have proposed a computational model for a reach-to-grasp movement, where the state of the hand is estimated and predicted by Kalman filters, and a motor command is generated that establishes a target grip aperture that is sufficiently large, in the stochastic manner, in relation to the target object size. The simulations of the model reproduced the effect of visual occlusion on the PGA during grasping. Online control of the movement would therefore require: (i) internal prediction of future states of the body and its variability, and (ii) motor commands based on the prediction and task constraints. The model was constructed within the framework of optimal feedback control in the sense of predictive stochastic control under sensory and motor variability, although the control has not yet been fully optimized.

Conflict of interest statement
------------------------------

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

This research was supported by Japan Science and Technology Agency, ERATO, Asada Synergistic Intelligence Project and by MEXT Grant-in-Aid for Scientific Research on Innovative Areas "Constructive Developmental Science" (15H01580).

Trajectory in minimum acceleration principle
============================================

Considering the trajectory from a current state *x*~*t*~ at current time *t* to a target state *x*~*T*~ at movement end time *T*, the cost to minimize is described by the following:

C

=

1

2

∫

t

T

(

d

2

x

t

′

d

t

′

2

)

2

d

t

′

In the similar manner as discussed by Flash and Hogan (1985) for the minimum jerk principle, the solution is given by a third order polynomial:

x

τ

=

a

0

\+

a

1

(

τ

\-

t

)

\+

a

2

(

τ

\-

t

)

2

\+

a

3

(

τ

\-

t

)

3

where τ denotes the arbitral time during movement (*t* ≤ τ ≤ *T*).

The solution of (A2) can be obtained by applying boundary conditions, which specify the position *x*~*t*~ and the velocity ẋ~*t*~ at time *t*:

a

0

=

x

t

a

1

=

ẋ

t

a

2

=

3

(

x

T

\-

x

t

)

(

T

\-

t

)

2

\-

ẋ

T

\+

2

ẋ

t

(

T

\-

t

)

a

3

=

\-

2

(

x

T

\-

x

t

)

(

T

\-

t

)

3

\-

ẋ

T

\+

ẋ

t

(

T

\-

t

)

2

Taking the derivative of (A2),

x

τ

=

a

1

\+

2

a

2

(

τ

\-

t

)

\+

3

a

3

(

τ

\-

t

)

2

The velocity profile in the minimum acceleration principle, Equations (4) and (5), is obtained from (A3) and (A4).

The kalman filter with time delay
=================================

The system dynamics are described by Equations (11) and (12). The estimation of the state at time *t*~1~ based on the observation until time *t*~2~ is denoted by $x_{t_{1}|t_{2}}$, and the covariance of the estimation error is denoted by $P_{t_{1}|t_{2}}$. Because of sensory delay, the observation at time *t* is used to correct the estimation of state at time *t* − *d*~*s*~:

K

=

P

t

\-

d

s

\|

t

\-

1

H

T

\[

H

P

t

\-

d

s

\|

t

\-

1

H

T

\+

R

t

\]

\-

1

x

\^

t

\-

d

s

\|

t

=

x

\^

t

\-

d

s

\|

t

\-

1

\+

K

(

y

t

\-

H

x

\^

t

\-

d

s

\|

t

\-

1

)

P

t

\-

d

s

\|

t

=

P

t

\-

d

s

\|

t

\-

1

\-

K

H

P

t

\-

d

s

\|

t

\-

1

Here, \[□\]^*T*^ denotes a transposition of the vector or matrix, while \[□\]^−1^ denotes an inverse matrix. The estimation of the state at *t* is calculated by corrected estimation ${\hat{x}}_{t - d_{s}|t}$ and a motor command that has already been output, by applying the following in order of *i* = −*d*~*s*~, ..., −1:

x

\^

t

\+

i

\+

1

\|

t

=

F

x

\^

t

\+

i

\|

t

\+

D

u

t

\+

i

\-

d

m

P

t

\+

i

\+

1

\|

t

=

F

P

t

\+

i

\|

t

F

T

\+

G

Q

t

\+

i

\-

d

m

G

T

The motor command generated at time *t* − 1 arrives at the end effector at time *t* + *d*~*m*~ − 1 because of the motor delay. Therefore, the state estimation ${\hat{x}}_{t + d_{m}|t}$ at time *t* + *d*~*m*~ can be calculated at time *t* by applying (B4) and (B5) for *i* = 0, ..., *d*~*m*~ − 1. As mentioned in Section Formulation of Aperture Size Control, this estimation is necessary for generation of the motor command *u*~*t*~. Equation (B5) is also used to predict future variability at the end of movement. As variability is predicted farther into the future, the predicted variability increases because repetitive application of Equation (B5) accumulates the motor variability *Q*~*t*~.

[^1]: Edited by: Vincent C. K. Cheung, The Chinese University of Hong Kong, Hong Kong

[^2]: Reviewed by: Robert Van Beers, VU University Amsterdam, Netherlands; Yutaka Sakaguchi, University of Electro-Communications, Japan

[^3]: †Present Address: Naohiro Takemura, Center for Information and Neural Networks, National Institute of Information and Communications Technology, Suita City, Japan;

[^4]: Takao Fukui, Department of Rehabilitation for Brain Functions, Research Institute of National Rehabilitation Center for Persons with Disabilities, Tokorozawa, Japan;

[^5]: Toshio Inui, Department of Psychology, Otemon Gakuin University, Ibaraki, Japan
