r - 特定の要因によるRのt検定

Question

いくつかの変数を含むデータセットがあります:

X は数値変数、Y と Z は 2 つの因子のみを含む因子変数 (Y=1,2 Z=3,4)

             x y z
1  -0.59131983 1 3
2   1.51800178 1 3
3   0.03079412 1 3
4  -0.43881764 1 3
5  -1.44914000 1 3
6  -1.33483914 1 4
7   0.25612595 1 4
8   0.12606742 1 4
9   0.44735965 1 4
10  1.83294817 1 4
11 -0.59131983 2 3
12  1.51800178 2 3
13  0.03079412 2 3
14 -0.43881764 2 3
15 -1.44914000 2 3
16 -1.33483914 2 4
17  0.25612595 2 4
18  0.12606742 2 4
19  0.44735965 2 4
20  1.83294817 2 4

因子変数が Y (t.test(X~Y)) の場合、t 検定は簡単に実行できます。しかし、たとえば、Z (3 と 4) の間で Y==2 の X 値のみを比較する t 検定を行う方法がわかりませんか?

私の言い方が正しいかどうか自信がないので、表で見た方が分かりやすいかもしれません。そこで、因子変数が Z で Y==2 である X のテストを行いたいと思います。どうすればこれを行うことができますか? STATA では簡単です: ttest var1 if var3==3, by(var2)

しかし、私はRでそれを取得しません:(

         x     y  z
11 -0.59131983 2 3
12  1.51800178 2 3
13  0.03079412 2 3
14 -0.43881764 2 3
15 -1.44914000 2 3
16 -1.33483914 2 4
17  0.25612595 2 4
18  0.12606742 2 4
19  0.44735965 2 4
20  1.83294817 2 4