我有以下data.frame(附加到问题)。我想做一个t.test比较每个变量(100%,75%,66%,50%和25%)的等级= 0和等级= 1的数量SBR的平均值。理想情况下,我想要一个给出mean(SBR,class=0), mean(SBR,class=1), variable, p-value)
的data.frame。目前,我可以获得每个计数密度的p值,但实际上非常繁琐)
例如
t.test(jr_all_melt[jr_all_melt$Class == 0 & jr_all_melt$variable == "100%",5], jr_all_melt[jr_all_melt$Class == 1 & jr_all_melt$variable == "100%",5])
我确信我可以使用plyr
等软件包来完成,但无法找到可行的示例。
structure(list(Age = c(50.2916666666667, 50.2916666666667, 81.7666666666667,
81.7666666666667, 72.4333333333333, 72.4333333333333, 62.1916666666667,
62.1916666666667, 79.3694444444444, 79.3694444444444, 81.0583333333333,
81.0583333333333, 67.4194444444445, 67.4194444444445, 82.6361111111111,
82.6361111111111, 49.2638888888889, 49.2638888888889, 79.5416666666667,
79.5416666666667, 53.8833333333333, 53.8833333333333, 53.7416666666667,
53.7416666666667, 91.1, 91.1, 79.6138888888889, 79.6138888888889,
64.8388888888889, 64.8388888888889, 81.0444444444445, 81.0444444444445,
64.3305555555555, 64.3305555555555, 67.0083333333333, 67.0083333333333,
56.3972222222222, 56.3972222222222, 83.1777777777778, 83.1777777777778,
75.9027777777778, 75.9027777777778, 79.7944444444445, 79.7944444444445,
56.4277777777778, 56.4277777777778, 67.5083333333333, 67.5083333333333,
70.9888888888889, 70.9888888888889, 67.5805555555555, 67.5805555555555,
67.7388888888889, 67.7388888888889, 58.3083333333333, 58.3083333333333,
68.2805555555556, 68.2805555555556, 52.6055555555556, 52.6055555555556,
84.2861111111111, 84.2861111111111, 70.7861111111111, 70.7861111111111,
75.5222222222222, 75.5222222222222, 82.8361111111111, 82.8361111111111,
80.2361111111111, 80.2361111111111, 76.3888888888889, 76.3888888888889,
73.4, 73.4, 59.3055555555556, 59.3055555555556, 74.6972222222222,
74.6972222222222, 76.025, 76.025, 65.0944444444444, 65.0944444444444,
64.1888888888889, 61.9722222222222, 52.0472222222222, 64.8055555555556,
72.5083333333333, 83.0416666666667, 72.9, 68.9666666666667, 60.3694444444444,
80.3027777777778, 70.9333333333333, 68.4555555555555, 70.825,
73.9972222222222, 60.0916666666667, 77.2833333333333, 73.0722222222222,
66.4861111111111, 66.4777777777778, 41.7722222222222, 78.1861111111111,
74.8305555555555, 72.1638888888889, 71.725, 78.7416666666667,
65.0055555555556, 77.7083333333333, 54.1722222222222, 62.8583333333333,
73.0361111111111, 78.4944444444444, 77.4027777777778, 75.2472222222222,
68.475, 64.7222222222222, 81.1972222222222, 51.5861111111111,
60.8833333333333, 73.3972222222222, 76.8694444444444, 53.8805555555556,
63.3083333333333, 78.1083333333333, 78.4666666666667, 50.2916666666667,
50.2916666666667, 81.7666666666667, 81.7666666666667, 72.4333333333333,
72.4333333333333, 62.1916666666667, 62.1916666666667, 79.3694444444444,
79.3694444444444, 81.0583333333333, 81.0583333333333, 67.4194444444445,
67.4194444444445, 82.6361111111111, 82.6361111111111, 49.2638888888889,
49.2638888888889, 79.5416666666667, 79.5416666666667, 53.8833333333333,
53.8833333333333, 53.7416666666667, 53.7416666666667, 91.1, 91.1,
79.6138888888889, 79.6138888888889, 64.8388888888889, 64.8388888888889,
81.0444444444445, 81.0444444444445, 64.3305555555555, 64.3305555555555,
67.0083333333333, 67.0083333333333, 56.3972222222222, 56.3972222222222,
83.1777777777778, 83.1777777777778, 75.9027777777778, 75.9027777777778,
79.7944444444445, 79.7944444444445, 56.4277777777778, 56.4277777777778,
67.5083333333333, 67.5083333333333, 70.9888888888889, 70.9888888888889,
67.5805555555555, 67.5805555555555, 67.7388888888889, 67.7388888888889,
58.3083333333333, 58.3083333333333, 68.2805555555556, 68.2805555555556,
52.6055555555556, 52.6055555555556, 84.2861111111111, 84.2861111111111,
70.7861111111111, 70.7861111111111, 75.5222222222222, 75.5222222222222,
82.8361111111111, 82.8361111111111, 80.2361111111111, 80.2361111111111,
76.3888888888889, 76.3888888888889, 73.4, 73.4, 59.3055555555556,
59.3055555555556, 74.6972222222222, 74.6972222222222, 76.025,
76.025, 65.0944444444444, 65.0944444444444, 64.1888888888889,
61.9722222222222, 52.0472222222222, 64.8055555555556, 72.5083333333333,
83.0416666666667, 72.9, 68.9666666666667, 60.3694444444444, 80.3027777777778,
70.9333333333333, 68.4555555555555, 70.825, 73.9972222222222,
60.0916666666667, 77.2833333333333, 73.0722222222222, 66.4861111111111,
66.4777777777778, 41.7722222222222, 78.1861111111111, 74.8305555555555,
72.1638888888889, 71.725, 78.7416666666667, 65.0055555555556,
77.7083333333333, 54.1722222222222, 62.8583333333333, 73.0361111111111,
78.4944444444444, 77.4027777777778, 75.2472222222222, 68.475,
64.7222222222222, 81.1972222222222, 51.5861111111111, 60.8833333333333,
73.3972222222222, 76.8694444444444, 53.8805555555556, 63.3083333333333,
78.1083333333333, 78.4666666666667, 50.2916666666667, 50.2916666666667,
81.7666666666667, 81.7666666666667, 72.4333333333333, 72.4333333333333,
62.1916666666667, 62.1916666666667, 79.3694444444444, 79.3694444444444,
81.0583333333333, 81.0583333333333, 67.4194444444445, 67.4194444444445,
82.6361111111111, 82.6361111111111, 49.2638888888889, 49.2638888888889,
79.5416666666667, 79.5416666666667, 53.8833333333333, 53.8833333333333,
53.7416666666667, 53.7416666666667, 91.1, 91.1, 79.6138888888889,
79.6138888888889, 64.8388888888889, 64.8388888888889, 81.0444444444445,
81.0444444444445, 64.3305555555555, 64.3305555555555, 67.0083333333333,
67.0083333333333, 56.3972222222222, 56.3972222222222, 83.1777777777778,
83.1777777777778, 75.9027777777778, 75.9027777777778, 79.7944444444445,
79.7944444444445, 56.4277777777778, 56.4277777777778, 67.5083333333333,
67.5083333333333, 70.9888888888889, 70.9888888888889, 67.5805555555555,
67.5805555555555, 67.7388888888889, 67.7388888888889, 58.3083333333333,
58.3083333333333, 68.2805555555556, 68.2805555555556, 52.6055555555556,
52.6055555555556, 84.2861111111111, 84.2861111111111, 70.7861111111111,
70.7861111111111, 75.5222222222222, 75.5222222222222, 82.8361111111111,
82.8361111111111, 80.2361111111111, 80.2361111111111, 76.3888888888889,
76.3888888888889, 73.4, 73.4, 59.3055555555556, 59.3055555555556,
74.6972222222222, 74.6972222222222, 76.025, 76.025, 65.0944444444444,
65.0944444444444, 64.1888888888889, 61.9722222222222, 52.0472222222222,
64.8055555555556, 72.5083333333333, 83.0416666666667, 72.9, 68.9666666666667,
60.3694444444444, 80.3027777777778, 70.9333333333333, 68.4555555555555,
70.825, 73.9972222222222, 60.0916666666667, 77.2833333333333,
73.0722222222222, 66.4861111111111, 66.4777777777778, 41.7722222222222,
78.1861111111111, 74.8305555555555, 72.1638888888889, 71.725,
78.7416666666667, 65.0055555555556, 77.7083333333333, 54.1722222222222,
62.8583333333333, 73.0361111111111, 78.4944444444444, 77.4027777777778,
75.2472222222222, 68.475, 64.7222222222222, 81.1972222222222,
51.5861111111111, 60.8833333333333, 73.3972222222222, 76.8694444444444,
53.8805555555556, 63.3083333333333, 78.1083333333333, 78.4666666666667,
50.2916666666667, 50.2916666666667, 81.7666666666667, 81.7666666666667,
72.4333333333333, 72.4333333333333, 62.1916666666667, 62.1916666666667,
79.3694444444444, 79.3694444444444, 81.0583333333333, 81.0583333333333,
67.4194444444445, 67.4194444444445, 82.6361111111111, 82.6361111111111,
49.2638888888889, 49.2638888888889, 79.5416666666667, 79.5416666666667,
53.8833333333333, 53.8833333333333, 53.7416666666667, 53.7416666666667,
91.1, 91.1, 79.6138888888889, 79.6138888888889, 64.8388888888889,
64.8388888888889, 81.0444444444445, 81.0444444444445, 64.3305555555555,
64.3305555555555, 67.0083333333333, 67.0083333333333, 56.3972222222222,
56.3972222222222, 83.1777777777778, 83.1777777777778, 75.9027777777778,
75.9027777777778, 79.7944444444445, 79.7944444444445, 56.4277777777778,
56.4277777777778, 67.5083333333333, 67.5083333333333, 70.9888888888889,
70.9888888888889, 67.5805555555555, 67.5805555555555, 67.7388888888889,
67.7388888888889, 58.3083333333333, 58.3083333333333, 68.2805555555556,
68.2805555555556, 52.6055555555556, 52.6055555555556, 84.2861111111111,
84.2861111111111, 70.7861111111111, 70.7861111111111, 75.5222222222222,
75.5222222222222, 82.8361111111111, 82.8361111111111, 80.2361111111111,
80.2361111111111, 76.3888888888889, 76.3888888888889, 73.4, 73.4,
59.3055555555556, 59.3055555555556, 74.6972222222222, 74.6972222222222,
76.025, 76.025, 65.0944444444444, 65.0944444444444, 64.1888888888889,
61.9722222222222, 52.0472222222222, 64.8055555555556, 72.5083333333333,
83.0416666666667, 72.9, 68.9666666666667, 60.3694444444444, 80.3027777777778,
70.9333333333333, 68.4555555555555, 70.825, 73.9972222222222,
60.0916666666667, 77.2833333333333, 73.0722222222222, 66.4861111111111,
66.4777777777778, 41.7722222222222, 78.1861111111111, 74.8305555555555,
72.1638888888889, 71.725, 78.7416666666667, 65.0055555555556,
77.7083333333333, 54.1722222222222, 62.8583333333333, 73.0361111111111,
78.4944444444444, 77.4027777777778, 75.2472222222222, 68.475,
64.7222222222222, 81.1972222222222, 51.5861111111111, 60.8833333333333,
73.3972222222222, 76.8694444444444, 53.8805555555556, 63.3083333333333,
78.1083333333333, 78.4666666666667, 50.2916666666667, 50.2916666666667,
81.7666666666667, 81.7666666666667, 72.4333333333333, 72.4333333333333,
62.1916666666667, 62.1916666666667, 79.3694444444444, 79.3694444444444,
81.0583333333333, 81.0583333333333, 67.4194444444445, 67.4194444444445,
82.6361111111111, 82.6361111111111, 49.2638888888889, 49.2638888888889,
79.5416666666667, 79.5416666666667, 53.8833333333333, 53.8833333333333,
53.7416666666667, 53.7416666666667, 91.1, 91.1, 79.6138888888889,
79.6138888888889, 64.8388888888889, 64.8388888888889, 81.0444444444445,
81.0444444444445, 64.3305555555555, 64.3305555555555, 67.0083333333333,
67.0083333333333, 56.3972222222222, 56.3972222222222, 83.1777777777778,
83.1777777777778, 75.9027777777778, 75.9027777777778, 79.7944444444445,
79.7944444444445, 56.4277777777778, 56.4277777777778, 67.5083333333333,
67.5083333333333, 70.9888888888889, 70.9888888888889, 67.5805555555555,
67.5805555555555, 67.7388888888889, 67.7388888888889, 58.3083333333333,
58.3083333333333, 68.2805555555556, 68.2805555555556, 52.6055555555556,
52.6055555555556, 84.2861111111111, 84.2861111111111, 70.7861111111111,
70.7861111111111, 75.5222222222222, 75.5222222222222, 82.8361111111111,
82.8361111111111, 80.2361111111111, 80.2361111111111, 76.3888888888889,
76.3888888888889, 73.4, 73.4, 59.3055555555556, 59.3055555555556,
74.6972222222222, 74.6972222222222, 76.025, 76.025, 65.0944444444444,
65.0944444444444, 64.1888888888889, 61.9722222222222, 52.0472222222222,
64.8055555555556, 72.5083333333333, 83.0416666666667, 72.9, 68.9666666666667,
60.3694444444444, 80.3027777777778, 70.9333333333333, 68.4555555555555,
70.825, 73.9972222222222, 60.0916666666667, 77.2833333333333,
73.0722222222222, 66.4861111111111, 66.4777777777778, 41.7722222222222,
78.1861111111111, 74.8305555555555, 72.1638888888889, 71.725,
78.7416666666667, 65.0055555555556, 77.7083333333333, 54.1722222222222,
62.8583333333333, 73.0361111111111, 78.4944444444444, 77.4027777777778,
75.2472222222222, 68.475, 64.7222222222222, 81.1972222222222,
51.5861111111111, 60.8833333333333, 73.3972222222222, 76.8694444444444,
53.8805555555556, 63.3083333333333, 78.1083333333333, 78.4666666666667
), Sex = c("M", "M", "M", "M", "F", "F", "M", "M", "M", "M",
"M", "M", "M", "M", "F", "F", "F", "F", "M", "M", "F", "F", "M",
"M", "F", "F", "F", "F", "F", "F", "M", "M", "F", "F", "M", "M",
"M", "M", "F", "F", "M", "M", "M", "M", "M", "M", "F", "F", "F",
"F", "F", "F", "F", "F", "M", "M", "M", "M", "F", "F", "F", "F",
"M", "M", "F", "F", "M", "M", "M", "M", "M", "M", "M", "M", "M",
"M", "F", "F", "F", "F", "F", "F", "F", "M", "M", "M", "M", "F",
"M", "F", "M", "M", "M", "M", "M", "M", "F", "M", "M", "M", "F",
"M", "F", "M", "M", "M", "F", "M", "F", "F", "F", "M", "F", "M",
"M", "F", "F", "M", "F", "F", "M", "M", "M", "M", "F", "F", "M",
"M", "M", "M", "F", "F", "M", "M", "M", "M", "M", "M", "M", "M",
"F", "F", "F", "F", "M", "M", "F", "F", "M", "M", "F", "F", "F",
"F", "F", "F", "M", "M", "F", "F", "M", "M", "M", "M", "F", "F",
"M", "M", "M", "M", "M", "M", "F", "F", "F", "F", "F", "F", "F",
"F", "M", "M", "M", "M", "F", "F", "F", "F", "M", "M", "F", "F",
"M", "M", "M", "M", "M", "M", "M", "M", "M", "M", "F", "F", "F",
"F", "F", "F", "F", "M", "M", "M", "M", "F", "M", "F", "M", "M",
"M", "M", "M", "M", "F", "M", "M", "M", "F", "M", "F", "M", "M",
"M", "F", "M", "F", "F", "F", "M", "F", "M", "M", "F", "F", "M",
"F", "F", "M", "M", "M", "M", "F", "F", "M", "M", "M", "M", "F",
"F", "M", "M", "M", "M", "M", "M", "M", "M", "F", "F", "F", "F",
"M", "M", "F", "F", "M", "M", "F", "F", "F", "F", "F", "F", "M",
"M", "F", "F", "M", "M", "M", "M", "F", "F", "M", "M", "M", "M",
"M", "M", "F", "F", "F", "F", "F", "F", "F", "F", "M", "M", "M",
"M", "F", "F", "F", "F", "M", "M", "F", "F", "M", "M", "M", "M",
"M", "M", "M", "M", "M", "M", "F", "F", "F", "F", "F", "F", "F",
"M", "M", "M", "M", "F", "M", "F", "M", "M", "M", "M", "M", "M",
"F", "M", "M", "M", "F", "M", "F", "M", "M", "M", "F", "M", "F",
"F", "F", "M", "F", "M", "M", "F", "F", "M", "F", "F", "M", "M",
"M", "M", "F", "F", "M", "M", "M", "M", "F", "F", "M", "M", "M",
"M", "M", "M", "M", "M", "F", "F", "F", "F", "M", "M", "F", "F",
"M", "M", "F", "F", "F", "F", "F", "F", "M", "M", "F", "F", "M",
"M", "M", "M", "F", "F", "M", "M", "M", "M", "M", "M", "F", "F",
"F", "F", "F", "F", "F", "F", "M", "M", "M", "M", "F", "F", "F",
"F", "M", "M", "F", "F", "M", "M", "M", "M", "M", "M", "M", "M",
"M", "M", "F", "F", "F", "F", "F", "F", "F", "M", "M", "M", "M",
"F", "M", "F", "M", "M", "M", "M", "M", "M", "F", "M", "M", "M",
"F", "M", "F", "M", "M", "M", "F", "M", "F", "F", "F", "M", "F",
"M", "M", "F", "F", "M", "F", "F", "M", "M", "M", "M", "F", "F",
"M", "M", "M", "M", "F", "F", "M", "M", "M", "M", "M", "M", "M",
"M", "F", "F", "F", "F", "M", "M", "F", "F", "M", "M", "F", "F",
"F", "F", "F", "F", "M", "M", "F", "F", "M", "M", "M", "M", "F",
"F", "M", "M", "M", "M", "M", "M", "F", "F", "F", "F", "F", "F",
"F", "F", "M", "M", "M", "M", "F", "F", "F", "F", "M", "M", "F",
"F", "M", "M", "M", "M", "M", "M", "M", "M", "M", "M", "F", "F",
"F", "F", "F", "F", "F", "M", "M", "M", "M", "F", "M", "F", "M",
"M", "M", "M", "M", "M", "F", "M", "M", "M", "F", "M", "F", "M",
"M", "M", "F", "M", "F", "F", "F", "M", "F", "M", "M", "F", "F",
"M", "F", "F", "M", "M", "M", "M", "F", "F"), Class = c(0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), variable = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L), .Label = c("100%", "75%", "66.6%", "50%",
"25%"), class = "factor"), SBR = c(12.681, 13.04, 11.72, 10.124,
10.69, 11.764, 8.561, 6.438, 4.263, 7.568, 9.231, 8.922, 12.651,
12.43, 8.281, 7.087, 12.944, 12.538, 11.443, 10.491, 12.792,
13.091, 7.091, 7.082, 5.406, 5.044, 11.007, 11.679, 9.909, 10.501,
10.268, 9.312, 16.485, 16.838, 11.743, 9.811, 13.447, 13.355,
5.138, 5.475, 8.883, 8.323, 9.054, 9.72, 21.538, 23.007, 8.842,
9.027, 5.757, 4.229, 5.293, 5.56, 9.714, 9.559, 13.698, 13.965,
7.364, 9.345, 9.309, 9.641, 6.893, 5.765, 7.45, 7.414, 3.651,
3.594, 8.243, 9.381, 7.068, 6.376, 7.011, 8.016, 9.391, 8.884,
8.24, 8.142, 9.045, 8.306, 5.427, 6.813, 10.123, 12.302, 3.953,
9.049, 4.292, 4.182, 5.084, 1.839, 2.693, 3.541, 4.079, 5.09,
5.655, 2.735, 8.321, -0.761, 4.331, 3.822, 1.403, 8.115, 7.87,
8.127, 3.931, 7.447, 1.362, 4.431, -1.583, 6.545, 2.851, 4.957,
6.518, 5.894, 3.874, 8.008, 4.72, 7.651, 8.212, 6.915, 9.625,
7.094, 2.291, 8.13, 7.71, 6.319, 2.494, 2.002, 12.478, 13.372,
12.212, 12.551, 10.571, 11.693, 8.674, 6.894, 5.061, 7.714, 8.512,
8.35, 12.251, 12.375, 7.898, 6.654, 12.89, 12.83, 12.24, 10.638,
13.615, 13.877, 6.845, 7.28, 5.224, 4.883, 10.742, 11.542, 10.308,
10.896, 9.356, 11.338, 15.793, 16.179, 12.567, 11.429, 13.431,
13.99, 6.123, 5.167, 8.709, 7.829, 7.599, 8.454, 19.393, 20.918,
8.33, 9.136, 4.824, 3.95, 6.425, 6.67, 8.22, 7.626, 13.909, 13.971,
8.425, 8.845, 8.887, 9.942, 7.078, 5.547, 6.837, 6.094, 4.576,
4.528, 8.619, 10.478, 3.125, 3.382, 7.22, 8.293, 11.238, 10.548,
7.197, 7.097, 9.659, 8.863, 5.577, 6.025, 10.77, 11.954, 3.814,
9.114, 4.435, 3.592, 5.265, 1.093, 3.903, 3.581, 3.602, 4.468,
5.451, 1.995, 7.03, -0.539, 4.875, 3.812, 1.476, 8.107, 8.628,
9.107, 2.54, 6.902, 1.553, 5.64, -2.762, 6.625, 2.715, 5.475,
6.085, 6.385, 2.794, 7.813, 4.801, 5.883, 7.837, 4.8, 9.527,
8.804, 2.988, 5.729, 7.275, 8.683, 2.627, 1.882, 11.884, 12.791,
10.861, 9.382, 9.525, 10.237, 9.299, 7.338, 5.075, 7.583, 10.324,
9.143, 14.003, 12.861, 7.687, 6.719, 13.422, 13.595, 11.828,
10.943, 12.296, 12.836, 7.069, 7.523, 5.691, 5.309, 11.494, 12.253,
8.711, 8.991, 9.009, 8.809, 15.387, 15.631, 11.948, 12.295, 15.034,
13.889, 5.454, 6.518, 9.535, 8.777, 8.436, 8.396, 18.508, 20.444,
9.505, 8.998, 4.705, 3.47, 5.202, 4.86, 7.532, 7.451, 14.499,
14.114, 6.872, 9.877, 9.077, 9.33, 6.009, 4.583, 5.988, 6.252,
4.766, 4.273, 11.44, 10.971, 6.393, 5.551, 7.107, 6.959, 11.844,
10.795, 9.183, 8.909, 9.142, 7.753, 6.055, 7.152, 9.566, 11.637,
3.668, 8.697, 4.701, 3.218, 4.292, 0.273, 2.713, 3.579, 5.496,
6.322, 4.852, 3.553, 6.82, -0.745, 4.581, 4.12, 1.886, 7.633,
8.07, 8.381, 3.627, 7.156, 1.669, 2.754, -2.56, 5.437, 2.587,
4.361, 4.959, 4.653, 3.2, 7.943, 5.922, 6.379, 5.403, 4.698,
9.665, 6.82, 2.432, 8.587, 6.4, 6.454, 2.054, 0.999, 10.129,
11.129, 13.195, 11.133, 10.702, 12.333, 9.4, 7.665, 5.538, 8.751,
8.868, 8.25, 12.695, 12.2, 8.346, 6.802, 12.007, 11.245, 11.456,
9.58, 12.941, 12.826, 6.239, 6.829, 5.887, 4.896, 10.303, 11.478,
8.972, 9.636, 9.538, 8.34, 16.686, 16.453, 10.827, 8.923, 13.867,
13.964, 5.016, 5.51, 7.635, 7.577, 7.725, 9.532, 21.083, 21.567,
9.417, 8.435, 5.635, 4.206, 4.299, 6.018, 8.642, 8.961, 12.613,
12.665, 8.861, 10.567, 9.365, 9.341, 6.435, 5.435, 7.099, 7.285,
5.467, 4.823, 8.959, 10.568, 6.085, 4.575, 7.267, 8.206, 9.589,
9.242, 7.504, 7.43, 9.171, 7.967, 3.426, 5.843, 10.341, 12.598,
4.001, 9.796, 4.188, 5.943, 6.093, 0.063, 2.512, 3.094, 3.106,
5.35, 6.468, 3.677, 7.964, -0.355, 4.088, 3.271, 0.608, 7.864,
7.941, 7.409, 1.796, 6.667, 0.731, 5.275, -1.904, 5.552, 2.487,
4.371, 6.181, 6.772, 2.993, 7.824, 4.659, 5.242, 5.33, 5.572,
9.409, 6.97, 2.791, 7.5, 6.753, 6.453, 1.955, 1.61, 13.193, 13.31,
12.485, 9.241, 11.773, 12.114, 9.747, 6.474, 5.956, 9.922, 11.264,
10.741, 15.046, 14.575, 6.828, 5.505, 13.366, 12.29, 10.432,
10.353, 12.386, 11.557, 6.68, 7.691, 5.158, 5.553, 10.371, 11.818,
9.635, 9.746, 8.955, 7.844, 18.448, 17.837, 11.19, 11.44, 13.722,
13.807, 7.387, 6.504, 8.724, 8.63, 9.441, 11.918, 19.074, 21.111,
8.849, 8.582, 3.558, 3.453, 4.963, 5.903, 9.605, 9.247, 12.929,
14.318, 7.057, 8.258, 9.982, 9.598, 8.709, 6.397, 7.797, 5.743,
5.742, 5.503, 7.563, 8.864, 5.685, 5.297, 7.63, 7.469, 10.426,
9.266, 7.451, 7.535, 9.941, 8.376, 4.793, 5.155, 9.867, 12.821,
3.052, 5.617, 2.7, 4.793, 6.38, 0.828, 2.281, 2.713, 4.508, 5.245,
5.325, 3.849, 6.326, 0.197, 5.909, 3.548, 1.077, 8.596, 10.599,
6.268, 1.86, 7.651, 1.083, 4.869, -2.972, 4.83, 2.146, 5.991,
4.333, 7.388, 2.181, 6.888, 4.85, 5.289, 5.878, 4.706, 4.609,
6.398, 2.106, 10.709, 6.243, 7.267, 1.431, 1.706)), row.names = c(NA,
-630L), .Names = c("Age", "Sex", "Class", "variable", "SBR"), class = "data.frame")
答案 0 :(得分:2)
library(dplyr)
library(purrr)
library(broom)
library(tidyr)
df %>%
group_by(variable) %>% # group by your desired variable
nest() %>% # split data using that variable
mutate(res = map(data, ~ tidy(t.test(SBR~Class, data = .)))) # apply t.test to each subset and keep results in a tidy format / tidy dataset
# # A tibble: 5 x 3
# variable data res
# <fctr> <list> <list>
# 1 100% <tibble [126 x 4]> <data.frame [1 x 10]>
# 2 75% <tibble [126 x 4]> <data.frame [1 x 10]>
# 3 66.6% <tibble [126 x 4]> <data.frame [1 x 10]>
# 4 50% <tibble [126 x 4]> <data.frame [1 x 10]>
# 5 25% <tibble [126 x 4]> <data.frame [1 x 10]>
到目前为止,您拥有单个数据集所需的一切。您有变量来分割数据(variable
),每次使用的数据集(data
)和t.test(res
)的结果。
您可以unnest
结果并保留您喜欢的列:
df %>%
group_by(variable) %>%
nest() %>%
mutate(res = map(data, ~ tidy(t.test(SBR~Class, data = .)))) %>%
unnest(res) %>%
select(variable, estimate1, estimate2, p.value)
# # A tibble: 5 x 4
# variable estimate1 estimate2 p.value
# <fctr> <dbl> <dbl> <dbl>
# 1 100% 9.484195 5.018591 1.146367e-12
# 2 75% 9.429390 4.895682 1.558697e-12
# 3 66.6% 9.421500 4.629750 2.110422e-14
# 4 50% 9.268463 4.683409 3.325438e-13
# 5 25% 9.580171 4.482977 7.025362e-15
如果您希望将p值设置为非科学格式,则可以在运行代码之前使用options(scipen = 999)
:
# # A tibble: 5 x 4
# variable estimate1 estimate2 p.value
# <fctr> <dbl> <dbl> <dbl>
# 1 100% 9.484195 5.018591 0.000000000001146367091
# 2 75% 9.429390 4.895682 0.000000000001558696898
# 3 66.6% 9.421500 4.629750 0.000000000000021104215
# 4 50% 9.268463 4.683409 0.000000000000332543837
# 5 25% 9.580171 4.482977 0.000000000000007025362