This problem can be conquered by combining Stars and Bars theory with Inclusion-Exclusion, as is done in this answer. Consider the following enumeration problem:
$~x_1 + x_2 + \cdots + x_{k+1} = (n-k).$
$~x_1, ~x_2, ~\cdots, ~x_{k+1} \in \mathbb{Z_{\geq 0}}.$
$~x_1, ~x_2, ~\cdots, ~x_{k+1} \in \mathbb{Z_{\leq t-1}}.$
To see that this represents the problem, consider the following tableau:
- - - | - - - | - - - ... - - - | - - -
The $~k~$ selected marbles are represented by the bars in the above tableau. These bars create $~(k+1)~$ islands that each represent a group of consecutive marbles that remain. Reading these islands from left to right, and letting $~x_1, ~x_2, ~\cdots, ~x_{k+1}~$ represent the sizes of each of the islands, none of the sizes is allowed to be $~\geq t.$ Since $~k~$ marbles have been selected, the size of the islands must sum to $~(n-k).$
If you ignore the third constraint above (i.e. the upper bound constraint), then by Stars and Bars theory, the enumeration would be
$$\binom{[n-k] + [(k+1)-1]}{[k+1]-1} = \binom{n}{k}, \tag1 $$
as expected.
Let $~r~$ denote $~\displaystyle \left\lfloor ~\frac{n-k}{t} ~\right\rfloor.$ This means that for those solutions where the third constraint is ignored, no more that $~r~$ of the variables can be $~\geq t.~$
For $~j \in \{1,2,\cdots,r\},~$ let $~f(j)~$ denote the number of solutions where the variables $~x_1, x_2, \cdots, x_j~$ are in violation of the upper bound constraint, and where the other $~k+1 - j~$ variables may or may not be in violation of the upper bound constraint.
Then $~f(j)~$ equals the enumeration of the number of solutions to
$~x_1 + x_2 + \cdots + x_{k+1} = (n-k).$
$~x_1, ~x_2, ~\cdots, ~x_{k+1} \in \mathbb{Z_{\geq 0}}.$
$~x_1, ~x_2, ~\cdots, ~x_j \in \mathbb{Z_{\geq t}}.$
$~f(j)~$ is easily computed by setting $~y_i = x_i - t ~: ~i \in \{1,2,\cdots,j\}.~$ Then, the above enumeration is equivalent to the enumeration of
$~y_1 + y_2 + \cdots + y_j + x_{j+1} + x_{j+2} + x_{k+1} = (n-k) - jt.$
$~y_1, ~y_2, ~\cdots, ~y_{j} \in \mathbb{Z_{\geq 0}}.$
$~x_{j+1}, ~x_{j+2}, ~\cdots, ~x_{k+1} \in \mathbb{Z_{\geq 0}}.$
By Stars and Bars theory,
$$f(j) = \binom{[n-k-jt] + [(k+1)-1]}{[k+1]-1} = \binom{n - jt}{k}. \tag2 $$
Now, Inclusion-Exclusion theory makes the final computation routine.
Per (1) above, let
$$~T_0 ~\text{denote} ~\binom{n}{k}. \tag3 $$
Per (2) above, consider that instead of the $~j~$ variables being represented specifically by $~\{x_1,x_2,\cdots,x_j\},~$ there are $~\displaystyle \binom{k+1}{j}~$ ways that the $~j~$ variables could have been designated.
So, for $~j \in \{1,2,\cdots,r\},~$ let
$$T_j ~\text{denote} ~f(j) \times \binom{k+1}{j} = \binom{n - jt}{k} \times \binom{k+1}{j}. \tag4 $$
Then, the final computation is
$$\sum_{j=0}^r (-1)^j T_j. \tag5 $$