我已经运行了48次t检验(手工编码而不是编写循环),并希望将这些t.tests的某些结果拼接出来,以创建一个我最感兴趣的事物表.
具体来说,我想保留这48个测试中每一个的p值,置信区间,x的平均值和y的平均值,然后构建一个结果表.
除了这里详细介绍的最佳答案之外,还有一种优雅,快速的方法吗?其中我会进行所有48次测试,并使用类似的东西获取所有三个所需的输出ttest$p.value?也许是循环?
下面是一个t检验的编码输入样本,然后是R提供的输出.
# t.test comparing means of Change_Unemp for 2005 government employment (ix)
lowgov6 <- met_res[met_res$Gov_Emp_2005 <= 93310, "Change_Unemp"]
highgov6 <- met_res[met_res$Gov_Emp_2005 > 93310, "Change_Unemp"]
t.test(lowgov6,highgov6,pool.sd=FALSE,na.rm=TRUE)
Welch Two Sample t-test
data: lowgov6 and highgov6
t = 1.5896, df = 78.978, p-value = 0.1159
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-0.1813909 1.6198399
sample estimates:
mean of x mean …Run Code Online (Sandbox Code Playgroud) 我有一个包含 99 行的两列数据框,我已“按升序排序”,但现在我只想根据第二列保留数据框的第一行和最后十行(1:10、90:99) (“改变”)。
例如
ID_NUM Change
1 -55223
42 -2321
6 -201
20 17
99 93
53 1009
Run Code Online (Sandbox Code Playgroud)
...ETC。
如何创建一个新的数据框,从现有数据框中排除中间的 11:89 行?