  We systematically performed numerical-relativity simulations for black hole
(BH) - neutron star (NS) binary mergers with a variety of the BH spin
orientation and equations of state (EOS) of the NS. The initial misalignment
angles of the BH spin are chosen in the range of i_tilt,0 = 30--90[deg.]. We
employed four models of NS EOS with which the compactness of the NS is in the
range of C = M_NS/R_NS = 0.138--0.180, where M_NS and R_NS are the mass and the
radius of the NS, respectively. The mass ratio of the BH to the NS, Q =
M_BH/M_NS, and the dimensionless spin parameter of the BH, chi, are chosen to
be Q = 5 and chi = 0.75, together with M_NS = 1.35 M_sun. We obtain the
following results: (i) The inclination angle of i_tilt,0 < 70[deg.] and
i_tilt,0 < 50[deg.] are required for the formation of a remnant disk with its
mass larger than 0.1 M_sun for the case C = 0.140 and C = 0.160, respectively,
while the disk mass is always smaller than 0.1M_sun for C = 0.175. The ejecta
with its mass larger than 0.01 M_sun is obtained for i_tilt,0 < 85[deg.] with C
= 0.140, for i_tilt,0 < 65[deg.] with C = 0.160, and for i_tilt,0 < 30[deg.]
with C = 0.175. (ii) The rotational axis of the dense part of the remnant disk
is approximately aligned with the remnant BH spin for i_tilt,0 = 30[deg.]. On
the other hand, the disk axis is misaligned initially with ~ 30[deg.] for
i_tilt,0 = 60[deg.], and the alignment with the remnant BH spin is achieved at
~ 50--60 ms after the onset of merger. The accretion time scale of the remnant
disk is typically ~ 100 ms and depends only weakly on the misalignment angle
and the EOS. (iii) The ejecta velocity is typically ~ 0.2--0.3c and depends
only weakly on i_tilt,0 and the EOS of the NS, while the morphology of the
ejecta depends on its mass. (iv) The gravitational-wave spectra contains the
information of the NS compactness in the cutoff frequency for i_tilt,0 <
60[deg.].
